2017 Tutorial InPlane Principal Stress Output
2017 Tutorial InPlane Principal Stress Output
Name: InPlanePrincipals/rosette
Path: /Examples/DianaIE/InPlanePrincipals/rosette
Keywords: analys: linear static.
class: large.
constr: suppor.
elemen: hx24l solid tp18l.
load: edge elemen force node.
materi: elasti isotro.
option: direct.
result: cauchy displa princi stress total.
Outline
Stress-rosettes in finite element analysis are used to represent the stress paths in the models. A rosette is represented in the middle of each element.
For 2D models these representations are easy to read, because there are only two principal directions presented in planar surfaces. This is not the
case for 3D models as the rosettes have 3 orthogonal lines and are located inside the model making cloudy pictures that are very difficult to interpret
(example in [Fig. 1]).
For this reason, DianaIE has an alternative approach: principal stresses can be projected to a surface (an outer surface of a cross-section located
inside the model). 2-dimensional rosettes with two orthogonal lines are presented in this surface resulting into readable pictures that provide the
information about stress orientations in a 3D model. These are the so-called in-plane principal stresses represented with 2-dimensional rosettes and
normal stress component displayed as contour plots.
It is a sectional representation of principal stress components that allows to see stress directions at outer and inside surfaces in the 3D model.
This tutorial explains the meaning of in-plane principal stress output in DianaIE followed by an example.
In a general 3D problem, the principal stresses are the eigenvalues of the stress tensor. The principal directions are the direction of the eigenvalues
(eigenvectors). Shear stresses are null in the principal directions.
zz
3
zx zy 2
xz yz
xx xy yx yy
1
xx
yy 1
zz = 2
= =
xy yx Eigenvalues = principal stresses 3
Eigenvectors = principal directions
yz = zy
zx = xz
1. Average the 3D stress tensor in the element from the stress tensor in the nodes;
3. Rotate the averaged 3D stress tensor such that the z-axis is perpendicular to the surface [Fig. 3];
zz
zx zy y'
xz yz z' x'
z xx xy yx yy
y
x
xx 'xx
yy 'yy
zz Referential transformation 'zz
= ' = Q Q T ' = '
xy xy
yz Q is the transformation matrix
with the director cosines 'yz
zx 'zx
'yy '1
'yx
'xy '2
Y' 'xx
Z ' X'
'xx '1
'IP = 'yy Eigenvalues = principal stresses
'IP = '
2
'xy Eigenvectors = principal directions
'zz
Figure 4: Determination of in-plane principals and normal stress component in the output
The model is a box with dimensions of height=2 m, width=1 m and length=10 m, supported at 2 end bottom edges (at begin length and at end
length) kj 0. The material is linear elastic with Youngs modulus of 30E+7 N/m2 and null Poisson ratio.
In the load case 1 (2D problem) in-plane principal stress components will directly correspond to the solver principal stresses as there is no shear
components.
In the load case 2 (3D problem) the in-plane principal stress components will differ from the solver principal stresses as the shear components are
not zero.
Figure 5: Box (2x1x10m) under edge load case Figure 6: Box (2x1x10m) under point load case
DianaIE
DianaIE
Figure 10: Geometry - Add polyline Figure 11: Geometry - Extrude polyline Figure 12: Geometry - Extrude polyline
Figure 14: Geometry - Add vertex Figure 15: Geometry - Imprint vertex Figure 17: Imprint vertex - Tool
Figure 18: Add edge support Figure 19: View of the model - Edge support
Figure 20: Add edge support Figure 21: View of the model - Edge support
Figure 22: Add point support Figure 23: View of the model - Point support
Figure 24: Assign properties Figure 25: Add new material - Linear elastic Figure 26: Material properties
Figure 27: Edge load Figure 29: Model view - Edge load
Figure 30: Point load Figure 32: Model view - Point load
Figure 35: Analysis window Figure 36: Add command Figure 37: Analysis tree
DianaIE
Analysis browser Analysis1 Structural nonlinear Output Edit properties [Fig. 38]
Properties - OUTPUT Result User selection Modify [Fig. 39]
Results Selection STRESS TOTAL CAUCHY GLOBAL Properties Location Element center [Fig. 40]
<do the same for the STRESS TOTAL CAUCHY PRINCI >
Figure 38: Edit output properties Figure 39: Selection of results Figure 40: Result properties
Y=1
Y=0
X=5
X=6
Results browser Result view setting Properties Contour plot settings Clip settings Add Plane 4x [Fig. 44] [Fig. 45]
<rename the 4 planes as Y=0, Y=1, X=5, X=6 (click in the names and rename)> [Fig. 46]
Figure 44: Add clipping plane Figure 45: Clipping planes Figure 46: Rename clipping planes
DianaIE
Contour plot settings Clip settings Y=0 Location [ 0.0 0.0001 0.0 ] , Normal [ 0 1 0 ] [Fig. 47]
Contour plot settings Clip settings Y=1 Location [ 0.0 0.9999 0.0 ] , Normal [ 0 1 0 ] [Fig. 48]
Figure 47: Clip plane settings Y=0 Figure 48: Clip plane settings Y=1
Contour plot settings Clip settings X=5 Location [ 5.0001 0.0 0.0 ] , Normal [ 1 0 0 ] [Fig. 49]
Contour plot settings Clip settings X=6 Location [ 6.0000 0.0 0.0 ] , Normal [ 1 0 0 ] [Fig. 50]
Figure 49: Clip plane settings X=5 Figure 50: Clip plane settings X=6
Note: For the outer surfaces we dont exactly set Y=0 and Y=1 but an approximate number to place the surfaces inside the model (Y=0.0001 and
Y=0.9999, respectively). For the X=5 surface we also give an approximate value (X=5.0001) in order move away from the symmetry point. For the
outer surfaces we could see the results without the clipping planes. However, the clipping panes we can do the contour plots only for these surfaces
and not for the entire mesh which make it easier to analyze the results.
Figure 51: Clip settings Figure 52: Result tree view - Total stresses Figure 53: SXX in the clipping planes
Contour plot setting Y=0 Enable plane <deactivate (take the tick) from the other planes> [Fig. 54]
Contour plot setting Y=1 Enable plane <deactivate (take the tick) from the other planes> [Fig. 55]
Contour plot setting X=5 Enable plane <deactivate (take the tick) from the other planes> [Fig. 56]
Contour plot setting X=6 Enable plane <deactivate (take the tick) from the other planes> [Fig. 57]
Figure 54: SXX in plane Y=0 Figure 56: SXX in plane X=5
Figure 55: SXX in plane Y=1 Figure 57: SXX in plane X=6
Figure 58: Edit clipping planes Figure 59: Edit clipping plane in the Graphics Window
Results browser Result view setting Properties Tensor plot settings Cutting plane settings Enable [Fig. 60]
Cutting plane settings Location [0.0 0.0001 0.0 ] , Normal [ 0 1 0 ] (plane Y=0) [Fig. 61]
Cutting plane settings Location [0.0 0.9999 0.0 ] , Normal [ 0 1 0 ] (plane Y=1) [Fig. 62]
Cutting plane settings Location [5.0001 0.0 0.0 ] , Normal [ 1 0 0 ] (plane X=5) [Fig. 63]
Cutting plane settings Location [6.0000 0.0 0.0 ] , Normal [ 1 0 0 ] (plane X=6) [Fig. 64]
Figure 61: Cutting plane settings Y=0 Figure 63: Cutting plane settings X=5
Figure 60: Cutting plane settings Figure 62: Cutting plane settings Y=1 Figure 64: Cutting plane settings X=6
Figure 66: In-plane stresses in plane Y=0 Figure 68: In-plane stresses in plane X=5
[Fig. 70] presents the total stresses SXX in the different surfaces (see [Fig. 52]). In the following we will see the global and principal stresses and
rosettes for each surface in detail. For that we activate only one clipping plane at a time (see [Fig. 61]-[Fig. 64] and use left view of the model.
Figure 71: SXX in clipping plane Y=0 Figure 73: S1 in clipping plane Y=0
Figure 72: SZZ in clipping plane Y=0 Figure 74: S3 in clipping plane Y=0
Note: As the Edge Load is symmetric, plane Y=1 (Back face) presents the same results as plane Y=0 (Front face) and is not presented for this case.
As this a load concentration area the in-plane principal stress components are not zero [Fig. 81]. The in-principal stresses displayed as rosettes are
approximately null in the bottom (where we have mainly SXX stresses [Fig. 77]) and approximately equal to SZZ [Fig. 78] in the top of the cross
section.
In the center length we have high SZZ stresses, so we have not null in plane stress components. The normal component show the compressive and
tension zones [Fig. 82].
Here the in-plane principals [Fig. 87] are approximately zero because the stresses appear predominantly in XX direction. So, the normal component
[Fig. 88] is similar to the global stress SXX [Fig. 83], showing the compressive and tension zones.
[Fig. 89] presents the Cauchy Total Stresses SXX in the different surfaces. In the following we will see the global and principal stresses and rosettes
for each surface in detail. For that, we define one clipping plane at a time (see [Fig. 61]- [Fig. 64]) and use left view of the model.
Figure 90: SXX in clipping plane Y=0 Figure 93: S1 in clipping plane Y=0
Figure 91: SYY in clipping plane Y=0 Figure 94: S2 in clipping plane Y=0
Figure 92: SZZ in clipping plane Y=0 Figure 95: S3 in clipping plane Y=0
Figure 98: SXX in clipping plane Y=1 Figure 101: S1 in clipping plane Y=1
Figure 99: SYY in clipping plane Y=1 Figure 102: S2 in clipping plane Y=1
Figure 100: SZZ in clipping plane Y=1 Figure 103: S3 in clipping plane Y=1
Figure 106: SXX in clipping plane X=6 Figure 109: S1 in clipping plan X=6
Figure 107: SYY in clipping plane X=6 Figure 110: S2 in clipping plan X=6