The two most used 2D element types in GSA are the plane stress element and the shell element.
Plane stress element
No variation and no stress through the thickness.
The state of stress at any point represented by the simplified stress tensor:
Shell elements
Assumes no through thickness direct stress.
Includes through thickness shear stress.
The stress tensor reduces to:
The stress tensor varies through the thickness of the element. For a linear model the stress at the mid-plane can be thought of as the in-plane stress and the difference between the top and bottom stresses as the bending stress. For a shell element, it is easier to convert these to ‘forces’ and ‘moments’.
These can be defined either by reference to an axis set (global or user defined axes), using the 2D Element Property ‘Axis’, or topologically (local axes).
2D element axes are defined by projecting the 2D Element Property ‘Axis’ on to the element, or, if set to ‘Local’, set the first edge defining the element x and the first and last edges defining the element xy plane.
Defining 2D element local axes by reference to an axis set results in more consistent local axes in the mesh. When the 2D element property axis is set to other than ‘local’, then the specified axis system is projected on to the element.
Look out for: Local axis definition - adjacent elements may have completely different axis directions making it difficult to interpret results.
Tip: If you aren’t sure about the orientation of the local axes, turn on ‘element axes’ in the Labels and Display Methods dialog. The axes are coloured Red / Green / Blue for x / y / z.
Note: Forces and moments are given per unit length in the relevant direction. These are essentially the appropriate stresses multiplied by the element thickness. For a 0.5 × 0.5m square element, therefore, you need to multiply the value by 0.5 to get the total force in the element.
These correspond to a 2D tensor based on the bending stresses, aggregated over the thickness of the element.
As in-plane forces and moments are tensor quantities we can define principal forces and moment, similar to principal stresses and maximum shear force and twisting moment.
These are the principal demands in whichever direction they occur and can be ‘projected’ onto an axis of interest.
When you bend a 2D element you get moments about the two in-plane axes ( and ) and a twisting moment (). These form the basis of projected moments.
Note: 'Projected’ refers to the way the output axis set is projected on to the surface of the element.
Maximum and minimum moments are the principal moments, (i.e., projected moments rotated so that the twisting moment is 0).
These are ‘derived’; calculated from the basic (projected) moment results and are independent of the axis set.
The same principles apply to the in-plane forces (, & and & ).
Note: The projected moments and are based on the stress in the x and y direction respectively (Timoshenko convention) and are not about the x and y axes.
Projected moments are given relative to the specified axes.
Example: gives moments in the direction of (not about) the x axis. Whether this is the local axis for each element, or the global axis for the whole model is very important. Always check which you are using by clicking the ‘Axes’ button in the contour plot dialogue box, via the contour settings icon to the right of the graphic interface.
If you’re unsure about the orientation of the local axes, turn on ‘element axes’ in the labels dialogue box. Both local and global axes are coloured in the same way.
The local x axis is defined in the direction of the first edge you draw in the 2D element.
Both are signed; the latter refers to principal compression. Look at both results to understand the principal compression/tension regime in each element.
To account for in-plane twisting, use the and results.
This will indicate the maximum force / unit length (GSA doesn’t define an equivalent stress) at all locations. There are likely to be large inaccuracies, especially adjacent to point supports.
Check the results first without averaging stresses and forces at nodes, to ensure the error is small.
If using envelopes, check the ‘envelope method’ used, located near the bottom of the contours dialogue box.
The fundamental result for 2D element is the stress state. This is represented by a tensor which can be thought of as a matrix, where each term corresponds to a force per unit area:
Where the shear terms may be denoted by instead of
We can use Mohr’s circle to understand what is happening in the element, and to the principal stresses.
The three circles correspond to the three pairs of directions x-y, y-z and z-x.
Taking the x-y directions the centre of the circle is:
and the radius:
giving principle stresses:
The same applies to the other two direction pairs leading to , and : maximum, intermediate and minimum principal stresses.
The maximum shear stress is simply the largest value of . Plotting the principal stresses give a good indication of the ‘flow’ of the stresses in the model, but there are still multiple values to consider.
Taking an element with an isotropic material under uniaxial stress we can measure how close to yielding we are by comparing the stress in the element with the yield stress in the material. To compare the general stress state to yield we normally use the von Mises stress.
This is defined as:
and in the uniaxial stress state reduces to:
Note:GSA RC design will check if your walls can be reinforced. See Projected Forces and Moments to do independent checks by hand.
To check shear stresses use the projected stresses. Normally, GSA sets the element axes with y vertical and x horizontal and z out-of-plane. So, the shear stresses that you are interested in are Projected xy. A quick hand calculation to check approximate stresses at, e.g., ground level, should confirm that you are looking at the right thing.