Users:General FEM Analysis/Elements Reference/Beam1

From Carat++ Public Wiki
(Difference between revisions)
Jump to: navigation, search
Line 61: Line 61:
 
</pre>
 
</pre>
  
 +
=== 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.
  
 
== Element Loading ==
 
== Element Loading ==
Line 81: Line 83:
 
Schwarz <ref name="Schwarz1991">Schwarz, H.: Methode der Finiten Elemente, Teubner, 1991</ref> and
 
Schwarz <ref name="Schwarz1991">Schwarz, H.: Methode der Finiten Elemente, Teubner, 1991</ref> and
 
Wunderlich <ref name="Wunderlich">Wunderlich, W.: Statik der Stabtragwerke, Teubner, 2004</ref>.
 
Wunderlich <ref name="Wunderlich">Wunderlich, W.: Statik der Stabtragwerke, Teubner, 2004</ref>.
 
For the correct use of the membrane element and the interpretation of related results, the following aspects should be considered:
 
 
 
=== Material parameter ===
 
With the parameter MAT the material, that would be used in the calculation, for the membrane element is defined. A very popular mistake is that a material number would be used, that isn't defined already. Please ensure that the material that you would like to use at this point are already exits. The following materials are testet for the membrane element:
 
* linear elastic isotropic
 
* linear elastic orthotropic (Münsch-Reinhardt)
 
* multilinear elastic isotropic
 
* elastoplastic isotropic
 
* material on the basis of response functions
 
 
=== Thickness parameter ===
 
The parameter THICKNESS defines the thickness of the membrane. The thickness should be constant over the element.
 
 
=== Prestress directions on the surface ===
 
The parameters ''A_X, A_Y, A_Z and B_X, B_Y, B_Z'' are used to define the prestress directions on the surface. In this approach the principle directions of the prestress are defined in a plane area (see figure below). The definition of the area is given by the two vectors '''f'''<sub>1</sub> and '''f'''<sub>2</sub>. The normal vector of the area can be calculated with the cross product of the in plane vectors '''f'''<sub>3</sub>='''f'''<sub>1</sub> x '''f'''<sub>2</sub>. Afterwards the line of intersection '''T'''<sub>1</sub> of the area which is given by '''f'''<sub>1</sub> and '''f'''<sub>3</sub> and the curved surface can be calculated. In this approach '''T'''<sub>1</sub> is interpreted as the first principle direction of the prestress on the curved surface. With the assumption that '''T'''<sub>3</sub> is equal to the surface normal vector '''G'''<sub>3</sub> (not uniformed), the second direction of the prestress is calculated as '''T'''<sub>2</sub>='''T'''<sub>1</sub> x '''T'''<sub>3</sub>. W.r.t. the parameters for the input file, only the plane area with the vectors '''f'''<sub>1</sub> and '''f'''<sub>2</sub> have to be defined. Referring to the depict approach the vector '''A''' defines the vector '''f'''<sub>1</sub> and the vector '''B''' defines the vector '''f'''<sub>2</sub>.
 
 
[[File:Prestress_Directions_for_Membranes.jpg|600px|up|Definition of the prestress direction for Membrane]]
 
 
=== Prestress state '''&sigma''' ===
 
SIG11, SIG22, SIG12 describes the prestress of the membrane element. The element is based on the plane stress assuption. Due to that only normal and shear stresses in the midplane have to be defined. SIGG11 is the stress acting in '''T'''<sub>1</sub>, SIG22 acting in '''T'''<sub>2</sub> and SIG12 is the in plane shear, whereas SIG12=SIG21 (see figure below).
 
 
[[File:Stress_state_for_membranes.jpg|600px|up|Stress state for Membranes]]
 
 
=== Lagrange type ===
 
With the LARANGE parameter it is possible to switch between form finding and statical/dynamical analysis. For the value UPDATED the element is for form finding and for the value TOTAL the element is for statical/dynamical analysis. It is important that the LANGRANGE parameter match to the type of analysis.
 
 
 
  
 
== References ==
 
== References ==
  
 
<references/>
 
<references/>

Revision as of 09:07, 31 October 2011

Coded, still under development and testing. Missing:

  • documentation in Carat++-Wiki
  • benchmark-examples

Coming soon

Contents

General Description

Element Type

  • This beam element is a linear 3D-beam taking into account shear deformation (Timoshenko-beam element).
  • This beam element has 6 DOFs per node (three translations and three rotations)

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.


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
KY, KZ 0 Shear correction factor; if 0 or no value specified, shear is not taken into account
Optional Parameters

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

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.

Element Loading

Pressure

  • not defined yet

Dead Load

  • not defined yet

Snow Load

  • not defined yet


Theory

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

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