Users:General FEM Analysis/Elements Reference/Beam1

From Carat++ Public Wiki
(Difference between revisions)
Jump to: navigation, search
Line 25: Line 25:
 
=== Orientation of the resultant forces ===
 
=== Orientation of the resultant forces ===
  
The Beam1 element has 6 resultant forces in accordance with the degrees of freedom. The resultant forces, ''[ N V1 V2 | MT M1 M2 ]'', are oriented along the local axes on the positive section of the beam (and opposite on the negative one, obviously). See also the following picture:
+
The Beam1 element has 6 resultant forces in accordance with the degrees of freedom. The resultant forces, ''[ N V1 V2 | MT M1 M2 ]'', are oriented along the local axes on the positive section of the beam. Note that the resultantforces are evaluated in the center of the beam element, so it's the mean value of the two ends. For the orientation, cf. also the orientation sketch:
 
[[File:Beam1-Resultant_forces_orientation.png|200px|frame|right|Resultant forces orientation for element Beam1]]
 
[[File:Beam1-Resultant_forces_orientation.png|200px|frame|right|Resultant forces orientation for element Beam1]]
  

Revision as of 12:03, 1 June 2012


Contents

General Description

Element Type

  • This beam element is a 2 node linear 3D-beam taking into account shear deformation (Timoshenko-beam element).
  • This beam element has 6 DOFs per node (three translations and three rotations)
  • The stiffness matrix is hard-coded, thus not needing any integration (and obviously limited to linear applications).

Degrees of Freedom

For the Beam1 element use the 3 translatoric degrees of freedom DISP_X, DISP_Y, DISP_Z and the 3 rotatoric degrees of freedom ROT_X, ROT_Y, ROT_Z.

Orientation of the local coordinate system

The Beam1 element uses the following definition for the determination of the local coordinate system (needed for the orientation of IYY and IZZ,...):

  • the local x-axis is oriented from node 1 to node 2 of the beam
  • the local y-axis lies in the global XY-plane, such that the local z-axis points in the same half-space as the global Z-axis (mathematically spoken: the local z-axis and the global Z-axis result in a positive dot-product)
  • the local z-axis is perpendicular to the other two local axis, following the right-hand-rule for x-y-z
  • exception: If the local x-axis (i.e. the beam axis) points in the direction of global Z, the local y-axis points in the direction of the global Z-axis. The local z-axis once again follows the right-hand-rule for x-y-z.

In case that a rotation of the local coordinate system is needed (rotated elements, inverse definition of IYY and IZZ,...) an angle THETA has to be specified. This angle rotates the whole coordinate system around the local x-axis, following the right-thumb rule (i.e. the thumb of the right hand points in the direction of the local x-axis).

Orientation of the resultant forces

The Beam1 element has 6 resultant forces in accordance with the degrees of freedom. The resultant forces, [ N V1 V2 | MT M1 M2 ], are oriented along the local axes on the positive section of the beam. Note that the resultantforces are evaluated in the center of the beam element, so it's the mean value of the two ends. For the orientation, cf. also the orientation sketch:

Resultant forces orientation for element Beam1


Input Parameters

Parameter Description

Compulsory Parameters
Parameter Values, Default(*) Description
MAT EL-MAT int Number for the used Material

e.g. MAT=EL-MAT 1

AREA Definition of the cross-sectional area of the beam
IYY, IZZ Definition of the moments of inertia
Optional Parameters
KY, KZ 0 Shear correction factor; if 0 or no value specified, shear is not taken into account (Bernouolli-beam theory)
THETA 0 Angle of rotation of the local coordinate system around the beam axis in degrees

Example of a Complete Input Block

EL-PROP 1 : BEAM1
MAT= EL-MAT 1
AREA=5.555  IYY=1.222  IZZ=1.333
KY=1.4   KZ=1.5
THETA=90

Use of the shear correction factors KY and KZ

The shear correction factors KY and KZ depend on the cross section of the beam. For a rectangle, KY resp. KZ are equal to 1.2, which gives a factor of α=1/Ki=5/6.

Use of the rotation parameter THETA

The rotation parameter THETA and its use is explained in the section concerning the coordinate system above.

The torsional resistance IT

The torsional resistance IT is interpreted as the polar moment Ipp, i.e.: IT = Ipp = IYY + IZZ.

Element Loading

For the moment, only nodal forces in the three global directions can be applied (i.e. Fx, Fy, Fz).

Pressure

  • not defined yet

Dead Load

  • not defined yet

Snow Load

  • not defined yet


Tests and Benchmarks

Static linear analysis

Static test for element Beam1

For the moment, the element Beam1 has successfully been tested in 3D in all its linear static features, including

  • bending, axial deformation, torsion,
  • shear deformation (separately definable for both local axes),
  • rotation around the local axis.

As an example, the structure on the right was part of the final tests.


Theory

The element implementation mainly follows the implementation of a linear 3D-beam element in FELyX [1], Schwarz [2] and Wunderlich [3].

As it is based on the Hermite-form functions, a Bernoulli beam can be modelled by one single element.

References

  1. https://www.rdb.ethz.ch/projects/project.php?proj_id=8314
  2. Schwarz, H.: Methode der Finiten Elemente, Teubner, 1991
  3. Wunderlich, W.: Statik der Stabtragwerke, Teubner, 2004




Whos here now:   Members 0   Guests 0   Bots & Crawlers 1
 
Personal tools
Content for Developers