ML022030334
| ML022030334 | |
| Person / Time | |
|---|---|
| Site: | Cook  |
| Issue date: | 07/16/2002 |
| From: | Greenlee S Indiana Michigan Power Co |
| To: | Document Control Desk, Office of Nuclear Reactor Regulation |
| References | |
| AEP:NRC:2520, TAC MB3603, TAC MB3604 SE 99-4 | |
| Download: ML022030334 (168) | |
Text
ATTACHMENT 5 TO AEP:NRC:2520 SOLVIA ENGINEERING REPORT SE 99-4, "SOLVIA VERIFICATION MANUAL LINEAR EXAMPLES"
/
SOLVIA Finite Element System Version 99.0 SOLVIA Verification Manual Linear Examples Report SE 99-4 SOLVIA Engineering AB
© SOLVIA Engineering AB, 1987 - 2000. All rights reserved.
SOLVIA Verification Manual, Linear Examples.
This manual, as well as the software described in it, is furnished under license and may be used or copied only in accordance with the terms of such license. The content of this manual is furnished for informational use only, is subject to change without notice, and should not be construed as a commitment by SOLVIA Engineering AB. SOLVIA Engineering AB assumes no responsibility or liability for any errors, inaccuracies or defects that may appear in this manual.
Except as permitted by such license, no part of this manual may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, recording, or otherwise, without the prior written permission of SOLVIA Engineering AB.
SOLVIA is a trademark or registered trademark of SOLVIA Engineering AB. Other tradenames, trademarks and registered trademarks included herein are property of their respective holders.
SOLVIA Engineering AB Tel +46-21-144050 Trefasgatan 3 Fax +46-21-188890 SE-721 30 Vasterts engineering@solvia.se Sweden www.solvia.com
SOLVIA Verification Manual INTRODUCTORY REMARKS The objective with this report is to present example solutions obtained with the SOLVIA-PRE, SOLVIA and SOLVIA-POST computer programs (the SOLVIA System) that verify and demonstrate their usage. Solutions to linear analyses are presented in this report. Nonlinear example solutions are presented in the companion report SE 99-5.
Since the aim is to compare the analysis results with analytical solutions, relatively small problems are solved, that also allow insight into the results. The analyses reported upon can be directly rerun with version 99.0 of the SOLVIA System. Complete input data for SOLVIA-PRE and SOLVIA POST is given for each example. All plot pictures created by the input data have been output to Microsoft Word in the PostScript graphical language.
We intend to update this report with further example solutions as we continue our work on the SOLVIA System. If you have any suggestions regarding the example solutions presented in this manual or suggestions on additional problems, we would be glad to hear from you.
Version 99.0 Linear Examples I
SOLVIA Verification Manual CONTENTS Al Thick cylinder under internal pressure A2 Axisymmetric shell under internal pressure A3 Simply supported circular plate under pressure load A4 Circular cylindrical shell under pressure load A5 Circular cylindrical shell under line load, SHELL A6 Circular cylindrical shell under line load, PLATE A7 Cantilever beam under tip loads, BEAM element A8 Cantilever beam under tip loads, ISOBEAM element A9 Pinched circular ring, ISOBEAM elements A10 Pinched circular ring, BEAM elements All Curved beam under out-of-plane load, ISOBEAM elements A12 Curved beam under out-of-plane load, BEAM elements A13 Cantilever truss structure under concentrated load A14 Cantilever under distributed load using skew systems, PLANE STRESS A15 Cantilever under distributed load using skew systems, PLANE STRAIN A16 Cantilever under distributed load using skew systems, SOLID A17 Cantilever under distributed load using skew systems, BEAM A18 Cantilever under distributed load using skew systems, ISOBEAM A19 Cantilever under distributed load using skew systems, PLATE A20 Cantilever under distributed load using skew systems, SHELL A21 Planar truss A22 Tapered cantilever under tip load A23 Stiffened plate cantilever under tip load A24 Analysis of spherical dome under self weight A25 Clamped square plate under pressure load A26 Material damping in modal superposition A27 Simply supported square plate under pressure load A28 Plate under uniform twisting A29 Edge bending and twisting of a triangular plate on corner supports A30 Scordelis-Lo cylindrical roof, cubic SHELL A31 Scordelis-Lo cylindrical roof, PLATE A32 Pinched cylindrical shell, SHELL elements A33 Pinched cylindrical shell, PLATE elements A34 Analysis of concentric fluid-filled cylinders A35 Z-section cantilever under distributed edge load A36 Cylindrical pressure vessel with hemispherical ends A37 Frequency analysis of cantilever with off-center masses A38 Simply supported skew plate under pressure load A39 Orthotropic plate under pressure load A40 Hemispherical shell under point loads A41 Beam on elastic foundation A42 Perforated tension strip A43 Fracture mechanics analysis of a tensile specimen A44 Fundamental frequency of cantilever, PLANE STRESS A45 Fundamental frequency of cantilever, PLANE STRAIN Version 99.0 Linear Examples 2
SOLVIA Verification Manual A46 Fundamental frequency of cantilever, SOLID A47 Fundamental frequency of cantilever, BEAM A48 Fundamental frequency of cantilever, ISOBEAM A49 Fundamental frequency of cantilever, PLATE A50 Fundamental frequency of cantilever, SHELL A51 Fundamental frequency of a simply supported plate A52 Wave propagation in a rod A53 Wave propagation in a water column A54 Cantilever subjected to ground motion A55 Cylindrical tube under step loading A56 Frequencies of a water-filled acoustic cavity A57 Thick-walled curved beam under end load A58 Beam subjected to a travelling load A59 Frequencies of a clamped thin rhombic plate A60 Transient heat conduction in a semi-infinite solid A61 Steady-state heat conduction in a square column A62 Steady-state heat conduction in hollow cylinder A63 Change in electric potential due to crack growth A64 Thermal eigenvalues and mode shapes A65 Torsional shear stress in a T-section beam A66 Fundamental frequency of a cantilever, 4-node SHELL A67 Scordelis-Lo cylindrical roof, 4-node SHELL A68 Analysis of a flanged elbow A69 Analysis of a rotating tubular shaft A70 Wave propagation in a rod A71 Transient analysis of a comer (Backward-Euler)
A72 Transient analysis of a comer (Trapezoidal rule)
A73 Semi-infinite region subjected to constant heat flux A74 Heat generation in semi-infinite solid (Trapezoidal)
A75 Heat generation in semi-infinite solid (Forward-Euler)
A76 Response spectrum analysis of a simply supported beam A77 Harmonic response of a two-degree-of-freedom system A78 Response spectrum of two load pulses A79 Earthquake excitation of a beam structure, floor response spectra A80 Case-combinations with different boundary conditions A81 Harmonic response of a simply supported beam A82 Stiffened plate cantilever under tip moment A83 Frequency analysis of a torsion spring system A84 Two-dimensional heat transfer with convection A85 Harmonic vibration of a damped plate A86 Harmonic response of a damped two-degree-of-freedom system A87 Modal combination methods in response spectrum analysis A88 Simply supported beam conditions using the end release option A89 Laminated strip A90 Sandwich shell Version 99.0 Linear Examples 3
SOLVIA Verification Manual A91 Laminated square plate under normal pressure A92 Dynamic analysis of a beam A93 Mode superposition analysis of a beam A94 Orthotropic cylinder under internal pressure A95 Material damping in complex-harmonic analysis A96 Frequencies of a beam with U cross-section A97 Segment of a pipe partly filled with water A98 Pipe joint A99 Cantilever beam, standard T cross-section AIOO Frequency analysis of a spherical dome, PLANE STRESS3 AlO Dynamic excitation of a beam structure A102 Conditional case combinations for a multispan beam Version 99.0 Linear Examples 4
SOLVIA Verification Manual EXAMPLE Al THICK CYLINDER UNDER INTERNAL PRESSURE Objective To verify the stress and displacement variation of the PLANE AXISYMMETRIC element and the application of axisymmetric pressure loading.
Physical Problem The figure below shows the cylinder to be analyzed. The cylinder is assumed to be guided so that no axial displacement can occur and it is acted upon by internal pressure.
E=2.0.10" N/rm2 v = 0.3 p = 10.0_ 106 N/rm2 a-=0.5m b= 1.0m Finite Element Model The top figure on page A 1.3 shows the finite element model. It consists of 8 PLANE AXISYMMETRIC elements having parabolic displacement variation. The element pressure is applied along the generated line between nodes 1 and 4.
Solution Results The theoretical stress solution is given in [1] p. 60:
a 2p (1 _ b2) y 27 c*x b2 _a2 az = v(Fx +oy)
The displacement in the radial direction can be calculated from the circumferential strain:
uy =y'x=Y(yx-V(a, + z))
Version 99.0 Linear Examples A1.1
SOLVIA Verification Manual or a2 p(y 2 +b 2 -2y 2 v)(l+v)
Uy
=
Ey(b 2 -a2)
The SOLVIA numerical solution obtained using the input data on pages A1.4 and Al.5 is as follows:
Stresses in 106 N/m 2 at node 1 and node 2:
y-coord.
(in)
SOLVIA Theory SOLVIA Theory SOLVIA Theory 0.5 16.73 16.67
-9.85
-10.00 2.066 2.000 1.0 6.672 6.667 0.012
- 0.
2.005 2.000 Radial displacements:
The plots from SOLVIA-POST on pages Al.3 and A1.4 show the variation of radial displacement, circumferential stress, radial stress and axial stress along a line in the radial direction.
User Hints
"* Symmetry gives that no variation of stresses can occur in the axial (Z) direction. A linear dis placement assumption in the axial direction would therefore be sufficient and the higher order displacement assumption in the axial direction can not improve the results.
"* Note that MYNODES is by default set to 100. The generated node numbers start with 101 as seen in the figure on page A1.3.
Reference
[1]
Timoshenko, S., Goodier, J.N., Theory of Elasticity, Second Edition, McGraw-Hill, 1951.
Version 99.0 y-coord.
SOLVIA Theory (in)
(in)
(in) 0.5 4.7667.10-5 4.7667.10-'
1.0 3.0333.10-s 3.0333.10-5 A1.2 Linear Examples
At THICK CYLINDER UNDER INTERNAL PRESSURE TIM*
i
- Tll, T
If *,,
o I
I
, I*a m
N -
cy)~U
+--
0.0 0.1 0.2 0.3 0.q 0.5 0.0 0
0.2 0.3
- 0.
- 0.
RADIUS NODE 1 -2 CYLINDER SOLVIA-POST 99.0 SOLVIA ENGINEERING AB Version 99.0 A1.3 SOLVIA Verification Manual Linear Examples No 0
x Co o
SOLVIA Verification Manual Linear Examples Lii u) cc a 0
Ai THICK CYLINDER UNDER INTERNAL PRESSURE S0-i c
/
N t
N U)
/
U)
/
F
/
/
N
/
/
00 01 0.2 0.3
- 0. 4
- 0.
0.0 01 0.2 I
TIMt I CYLINDER SOLVIA ENGINEERING AB SOLVIA-POST 99.0 SOLVIA-PRE input HEAD
'Al THICK CYLINDER UNDER INTERNAL PRESSURE' DATABASE CREATE MASTER IDOF=I01111 COORDINATES ENTRIES NODE Y
1 0.5 2
- 1.
3
- 1.
4 0.5 z
0.125 0.125 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
PLANE AXISYMMETRIC RESULTS=NSTRESSES GSURFACE 1 2 3 4 EL1=8 EL2=1 NODES=8 LOADS ELEMENT INPUT=LINE 1 4 10.E6 ic 0
CYLINDER 0
0 0.3
- 0. 4 0.1 SUBFRAME 12 MESH ENUMBER=YES MESH NSYMBOL=YES SOLVIA END VECTOR=LOAD NNUMBER=YES GSCALE=OLD Version 99.0 A1.4
SOLVIA Verification Manual SOLVIA-POST input Al THICK CYLINDER UNDER INTERNAL PRESSURE DATABASE CREATE WRITE FILENAME='al.lis' EPLINE NAME=CYLINDER 1
1 5 2 TO 8
1 5 2 NPLINE NAME=RADIUS 1
101 TO 115 2
AXIS ID=I VMIN=0 VMAX=100.E5 LABEL='STRESS-ZZ' NLINE LINENAME=RADIUS DIRECTION=2 OUTPUT=ALL SUBFRAME=21 ELINE LINENAME=CYLINDER KIND=SXX OUTPUT=ALL ELINE LINENAME=CYLINDER KIND=SYY OUTPUT=ALL SUBFRAME=21 ELINE LINENAME=CYLINDER KIND=SZZ YAXIS=l OUTPUT=ALL END Version 99.0 Linear Examples A1.5
SOLVIA Verification Manual EXAMPLE A2 AXISYMMETRIC SHELL UNDER INTERNAL PRESSURE Objective To verify the PLANE AXISYMMETRIC element under distributed loading when used for axisymmetric shell bending problems.
Physical Problem The figure below shows the cylinder to be analyzed. The radius to thickness ratio is 20 so thin shell theory is applicable. The cylinder is fixed at one end and it is loaded by internal pressure. The cylinder is long so that the effects of bending at one end do not affect the other end.
Ri-h E = 2.0.1011 N/mn 2
v = 0.3 Ri =0.5 m h = 0.025 m p = 1.0_10 6 N/nm2 Portion of interest Finite Element Model The figure on page A2.4 shows the finite element model. The 8-node PLANE AXISYMMETRIC element with 2x2 integration is used to describe the bending behaviour. Since most of the bending is concentrated towards the end a finer mesh is used for that part.
Solution Results The behaviour of the cylindrical shell under edge loading and internal pressure is described for example in [1] p. 140.
Version 99.0 Linear Examples A2.1
SOLVIA Verification Manual Using 4 3(1-v)2 Mo= P-a-C=-2 n
h2 a2 2n 2 )
D Eh' T. =2Pa
]
D 12(l_ V2)
T=P-an
- 2)
C a Pa =P a-h/2
-E2 where E = Young's modulus v =
Poisson's ratio a =
mean radius h =
thickness PA = internal pressure applied at radius a - h/2. Note that the factor C is introduced to calculate equivalent loads acting at the midsurface.
we obtain the radial displacement as w =
Pa C(
-) _-x (cos nx+sin nx) 4n4D V) and the axial bending moment M =
~+M. ) sin nx +M. cos nx) e-n The stresses in the axial (z) and circumferential (x) directions are paa + 6Mr zz 2h
- h-3 wE xx= -
+ V Yxx a
where r is the distance from the mean radius of the cylinder.
Version 99.0 Linear Examples A2.2
SOLVIA Verification Manual Using the input data on pages A2.6 and A2.7 we obtain the following results from SOLVIA:
Radial displacement (10-6 m) z-coord.
SOLVIA Theory (m) 0.025 3.76 2.91 0.050 10.74 9.51 0.075 18.73 17.41 0.100 26.32 25.08 0.200 43.72 43.21 0.300 45.81 45.53 0.400 44.61 44.29 0.500 44.07 43.71 0.900 44.11 43.75 The SOLVIA results are taken from the nodes on the mean surface of the cylinder. The displacements are also shown by the deformed mesh on page A2.5.
Stresses at inner Gauss points (MN/m2) z-coord.
oXYzz SOLVIA Theory SOLVIA Theory 0.0026 7.16 7.99 26.07 26.59 0.0099 7.26 7.36 23.39 23.87 0.0151 7.22 7.06 21.59 22.02 0.0223 6.92 6.84 19.31 19.71 0.0526 8.08 7.73 12.15 12.35 0.1151 13.81 13.32 6.38 6.34 0.2151 19.86 19.52 7.63 7.57 0.9000 20.29 20.00 9.76 9.76 The stress results for the inner Gauss points are also plotted on page A2.5.
User Hints
"* The 8-node PLANE element describes the bending behaviour quite good. The 4-node element, which has a linear assumption on the displacement behaviour in the radial and the axial direction, would not be satisfactory.
"* The portion of most interest in edge bending of an axisymmetric cylindrical shell extends approximately the length 2. --4-h in the axial direction, where a is the mean radius and h is the thickness. This portion should, therefore, be modeled with a finer mesh.
Version 99.0 Linear Examples A2.3
Linear Examples SOLVIA Verification Manual 0 In thin shell theory the pressure loading is applied at the midsurface of the shell. For PLANE elements, however, the pressure loading is applied at the element boundaries. A correction factor C is, therefore, introduced in the formulas obtained using thin shell theory in order to calculate equivalent loads acting at the midsurface.
Reference
[1]
Kraus, H., Thin Elastic Shells, John Wiley & Sons, 1967 A2 AXISYMMETRIC SHELL UNDER INTERNAL PRESSURE ORIGINAL 0.1 TIME 1Z Z
PRESSURE
- 9. 7561E6 ORIGINAL 0.02 ZONE FIXSIDE MAST ER 10011i B 1t0111 C IIIItI SOLVIA ENGINEERING AB SOLVIA-PRE 99.0 Version 99.0 T~-
F IZ A2.4
SOLVIA Verification Manual A2 0.0 0.2 0.4 AXIA AXISYMMETRIC SHELL I
TIMý i
0.6 0.8 1.,
UNDER INTERNAL PRESSURE
ý0
,I
. I,,
CO
- I N
N UU)
Co 04 04 0
to4-TIME I
.0 0.2
- 0. 4 0.6 0.8 i. I AXIAL SOLVIA ENGINEERING AB SOLVIA-POST 99.0 Version 99.0 A2 AXISYMMETRIC SHELL UNDER INTERNAL PRESSURE ORIGINAL H
0.05 MAX DISPL.
i 4.888SE-S TIME I
Z R
SOLVIA-POST 99.0 SOLVIA ENGINEERING AB to 04 004 "N
0D4 to X
X I-u cM 04l
/
/ /
/
i..
I I
Linear Examples A2.5
SOLVIA Verification Manual Linear Examples SOLVIA-PRE input HEADING
'A2 AXISYMMETRIC SHELL UNDER INTERNAL PRESSURE' DATABASE CREATE MASTER IDOF=I00111 COORDINATES ENTRIES NODE Y
Z 1
0.5 TO 3
0.525 4
0.5 0.225 5
0.525 0.225 6
0.5 0.9 7
0.525 0.9 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
PLANE AXISYMMETRIC INT=2 GSURFACE 6 4 5 7 EL1=27 EL2=1 NODES=8 GSURAFCE 4 1 3 5 EL1=18 EL2=1 NODES=8 LOADS ELEMENT INPUT=LINE 1 4 1.E6 4 6 1.E6 6 7 -9.7561E6 3
FIXBOUNDARIES 3
/
1 3 FIXBOUNDARIES 23
/
2 VIEW ID=1 XVIEW=1 ROTATION=-90 SET NSYMBOLS=MYNODES VIEW=1 MESH NNUMBERS=MYNODES VECTOR=LOAD SUBFRAME=12 ZONE NAME=FIXSIDE INPUT=GLOBAL-LIMITS ZMAX=0.230 MESH ZONENAME=FIXSIDE ENUMBER=YES BCODE=ALL SOLVIA END Version 99.0 A2.6
SOLVIA Verification Manual SOLVIA-POST input A2 AXISYMMETRIC SHELL UNDER INTERNAL PRESSURE DATABASE CREATE WRITE FILENAME='a2.1is' VIEW ID=1 XVIEW=1 ROTATION=-90 SET ORIGINAL=DASHED DEFORMED=YES VIEW=1 MESH OUTLINE=YES DMAX=1 EPLINE NAME=AXIAL 45 2 4 TO 1
2 4 ELINE LINENAME=AXIAL KIND=SXX OUTPUT=ALL SUBFRAME=21 ELINE LINENAME=AXIAL KIND=SZZ OUTPUT=ALL NLIST KIND=DISPLACEMENT NLIST KIND=REACTION END Version 99.0 Linear Examples A2.7
SOLVIA Verification Manual EXAMPLE A3 SIMPLY SUPPORTED CIRCULAR PLATE UNDER PRESSURE LOAD Objective To verify the PLANE AXISYMMETRIC element under distributed loading when used for plate bending problems.
Physical Problem The figure below shows the circular plate to be analyzed. It is a thin plate since the diameter to thickness ratio is 20. The boundary of the plate is simply supported.
a P
I I1L4I 9
I-------
5; t
E =2.0. 10" N/rm2 v =0.3 a=0.5m h =0.05 m p =1.0.10 6 N/rm2 Finite Element Model The top figure on page A3.3 shows the finite element model. Ten 8-node PLANE AXISYMMETRIC elements are used.
Solution Results The theoretical solution is given for example in [1] p. 56 as follows:
Center deflection:
Wmx(5 v)pa 4 =0.1739.10-T m 64(l+V)D D=
Eh 3 12(1-v
- 2)
Moment:
My =
(3+v)(a2 -y2)
Mx = P(a' (3+v)- y2 (1+ 3v))
16 Version 99.0 I I I I I Linear Examples I
A3.1
SOLVIA Verification Manual Stresses at x=0 and y--0:
6My y=-
= 123.75.106 Nm2
= 6Mx = 123.75.106 Nm2 Using the input data on page A3.5 the following results are obtained:
Center deflection SOLVIA Theory
-0.1752.102
-0.1739.10-2 Stresses of element 10 at node 1 (at y = z = 0):
Oxx (N / m2) ar (N/m2)
SOLVIA Theory SOLVIA Theory 124.02 123.75 124.02 123.75 The distribution of the von Mises effective stress and the stresses 0. and Yyy along the lower surface of the plate as obtained from SOLVIA are shown in the bottom figure on page A3.3 and the top figure on page A3.4.
User Hints
"* The 8-node PLANE element should be used in order to model the bending behaviour.
"* Note that the PLANE AXISYMMETRIC element extends 1 radian in the circumferential direction.
Reference
[1]
Timoshenko, S., Woinowsky-Krieger, S., Theory of Plates and Shells, Second Edition, McGraw-Hill, 1959.
Version 99.0 Linear Examples A3.2
SOLVIA Verification Manual Linear Examples A3 SIMPLY SUPPORTED CIRCULAR PLATE UNDER PRESSURE LOAD ORIGINAL 0.05 Z
LY R MAS TER 100/
B 10 1 C i10'"11 z
LY R
ORIGINAL i 0.05 TIME I
PRESSURE 1E6 SOLVIA ENGINEERING AB SOLVIA-PRE 99.0 Version 99.0 AS SIMPLY SUPPORTED CIRCULAR PLATE UNDER PRESSURE LOAD ORIGINAL ý-
0.05 Z
MAX DISPL.
1.752E-3 Y
TIME 1
R J
LOAD 15833 MAX I..2400E8
- 1. 1628E8 1 008SE8
- 8. S421E 7
- 6. 9992E7 S 4562E7
- 3. 9132E7
- 2. 3703E7
- 8. 2728E6 MIN S.S800ES SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A3.3 I.
t 6J 1
SOLVIA Verification Manual Linear Examples Version 99.0 A3 SIMPLY SUPPORTED CIRCULAR PLATE UNDER PRESSURE LOAD TI I.
ITIM TIIMý I
xx I-F
\\N o
c
)....
.o i..
.i.
0.0 1
0!2 0!3 0.4 0.0 0.0'
. 0!2 0!3 0.4 O.5 RADIAL RADIAL SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A3.4
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A3 SIMPLY SUPPORTED CIRCULAR PLATE UNDER PRESSURE LOAD?
DATABASE CREATE MASTER IDOF=100111 COORDINATES ENTRIES NODE Y
Z 1
TO 3
0.0 0.05 4
0.5 TO 6
0.5 0.05 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
PLANE AXISYMMETRIC RESULTS=NSTRESSES GSURFACE 6 3 1 4 EL1=10 EL2=1 NODES=8 LOADS ELEMENT INPUT=LINE 3
6 1.E6 FIXBOUNDARIES 2
/
1 2 3 FIXBOUNDARIES 3
/
5 SET NSYMBOLS=MYNODES SUBFRAME 12 MESH ENUMBER=YES BCODE=ALL MESH NNUMBER=MYNODES VECTOR=LOAD SOLVIA END SOLVIA-POST input A3 SIMPLY SUPPORTED CIRCULAR PLATE UNDER PRESSURE LOAD DATABASE CREATE WRITE FILENAME='a3.1is' EPLINE NAME=RADIAL 10 3 7 4 TO 1 3 7 4 SET NSYMBOLS=MYNODES MESH ORIGINAL=YES CONTOUR=MISES VECTOR=LOAD SUBFRAME 21 ELINE LINENAME=RADIAL KIND=SXX OUTPUT=ALL ELINE LINENAME=RADIAL KIND=SYY OUTPUT=ALL NMAX KIND=DISPLACEMENT NUMBER=3 SUMMATION KIND=REACTION DETAILS=YES END Version 99.0 Linear Examples A3.5
SOLVIA Verification Manual EXAMPLE A4 CIRCULAR CYLINDRICAL SHELL UNDER PRESSURE LOAD Objective To verify the membrane behaviour of a curved thin SHELL element under distributed loading.
Physical Problem The pressure vessel shown in the figure below is loaded by internal pressure. A section of the shell subjected only to membrane action is to be analyzed.
E = 6.625-10'0 N/rm2 v =0.33 a= 2.0000 m h = 0.0010 m p = 0.1 -10 N/rm 2
Finite Element Model The figure below shows the finite element model. We use one 9-node SHELL element extending 15 degrees and skew systems oriented in the radial and circumferential directions for ease of applying the boundary conditions. The axial loading is simulated by SHELL edge loading. The numbering of nodes and the boundary conditions are shown in the top figure on page A4.3.
Skew systems defined at nodes Version 99.0 Linear Examples A4.1
SOLVIA Verification Manual Solution Results The radial deformation is w Eh(
V and the stresses are:
_ p a 2h a~hoop = pa h
The SOLVIA numerical results obtained using the input data shown on page A4.4 are as follows:
Displacement Stresses w(i0-,M)
Ora,(16 NI/M2) 7hOp (106 N /rn)
SOLVIA Theory SOLVIA Theory SOLVIA Theory 0.5041 0.5041 10.0 10.0 20.0 20.0 The deformed mesh is shown in the bottom figure on page A4.3.
User Hints
"* The quadratic SHELL element is chosen in order to approximate the circular shell. This approxi mation leads to some very small variation in the displacement and stress results since the curvature of the element varies slightly over the element.
"* Note that SOLVIA-PRE automatically assigns GLOBAL rotations for a SHELL midsurface node when the node has a boundary condition for rotation specified [1].
"* If the 9-node SHELL element shall model bending then three or more elements are recommended for the 15 degree portion, thus at least 5 degrees per element.
Reference
[1]
SOLVIA-PRE 99.0 Users Manual, Stress Analysis, Report SE 99-1, p. 7.23 - 7.34.
Version 99.0 Linear Examples A4.2
SOLVIA Verification Manual Linear Examples Version 99.0 A4 CIRCULAR CYLINDRICAL SHELL UNDER PRESSURE LOAD ORIGINAL 0.2 Z
ORIGINAL 0.1 z
TIME I
EFORCE 10000 PRESSURE 10000 C
MASTER 000000 a
b 010a 1 NAXES=SKEW C
110111 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A4 CIRCULAR CYLINDRICAL SHELL UNDER PRESSURE LOAD ORIGINAL
- O.S z
MAX DISPL.
-H S.0482E-4 TIME I
X y
ý tt REACTION N
/6666.S t
r s
EAXES STRESS-RST SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A4.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A4 CIRCULAR CYLINDRICAL SHELL UNDER PRESSURE LOAD' DATABASE CREATE SYSTEM 1 CYLINDRICAL COORDINATES ENTRIES NODE R
THETA XL 1
- 2.
- 90.
- 0.
TO 3
- 2.
- 75.
- 0.
4
- 2.
- 90.
0.5 TO 6
- 2.
- 75.
0.5 SKEWSYSTEM EULERANGLES 1
-7.5 2 -15.
MATERIAL 1
ELASTIC E=6.625E10 NU=0.33 EGROUP 1
SHELL STRESSREFERENCE=ELEMENT THICKNESS 1 1.E-3 GSURFACE 3 1 4 6 EL1=1 EL2=1 NODES=9 SYSTEM=1 LOADS ELEMENT TYPE=PRESSURE INPUT=ELEMENTS 1 -T I.E4 LOADS ELEMENT TYPE=FORCE INPUT=LINE 4
6 OUT 1.E4 NSKEWS INPUT=NODES
/
2 1
/
5 1 NSKEWS INPUT=LINES
/
3 6 2
FIXBOUNDARIES 12456 INPUT=LINE 1 3 FIXBOUNDARIES 2456 INPUT=LINE
/
1 4 4 6 3 6 SET SMOOTHNESS=YES MESH NSYMBOLS=YES NAXES=SKEW SUBFRAME=21 MESH VECTOR=LOAD BCODE=ALL SOLVIA END SOLVIA-POST input A4 CIRCULAR CYLINDRICAL SHELL UNDER PRESSURE LOAD DATABASE CREATE WRITE FILENAME='a4.1is' SET SMOOTHNESS=YES MESH VECTOR=REACTION ORIGINAL=DASHED EAXES=STRESS-RST NLIST DIRECTION=123 EMAX END Version 99.0 Linear Examples A4.4
SOLVIA Verification Manual EXAMPLE A5 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, SHELL Objective To verify the bending behaviour of the curved thin SHELL element.
Physical Problem The figure below shows the circular cylindrical shell to be analyzed. It is acted upon by equal and opposite line loads. The solution is sought for the part of the shell which is not disturbed by end effects.
z I TL I ThL r =2.00 mn h =0.02 mF b=0.50 m
= 400 N/rm 2
= 2.0.10" N/im2 i = 0.3 Finite Element Model The figure on page A5.3 shows the finite element model used. Symmetry gives that only a 900 portion of the cylinder need be modelled. All 4 boundary faces have symmetry boundary conditions. In addition, the displacements in the X-direction are zero. The 16-node SHELL element with Gauss integration orders 4x4x2 as well as 3x3x2 have been used for comparison.
Solution Results The theoretical solution is given in [I] p. 381 for a thin ring under equal and opposite forces. A correction to account for the plane strain condition in the axial direction gives then 6, =-(1-c-2 bPr 3 2k1v 4 7t)
El g 1/2(-2 (2-i brcLT
- 2)
El where 5z is the radial displacement at the application of the load and By is the radial displacement 90' from the load (along the Y-axis).
Version 99.0 Linear Examples A5.1
SOLVIA Verification Manual The following displacement results are obtained using the input data shown on page A5.7:
Theory SOLVIA SOLVIA SOLVIA 6 elements 12 elements 6 elements 4x4x2 int.
4x4x2 int.
3x3x2 int.
8 (mm)
-1.625
-1.607
-1.625
-1.625 5Y (mm) 1.492 1.480 1.492 1.492 Compare also the bending stress at the outer face for ý=O and 04=90'.
The theoretical bending moment is ([1] p. 380):
M= bPr Cos
-2) 2
(
Ir)
M((=0)= 72.676[Nm]
(145.352Nm/m)
M(p = 90°) =-127.324[Nm]
(-254.648Nm/m)
The following results are obtained Bending stress [MPa]
Theory SOLVIA SOLVIA SOLVIA 6 elements 12 elements 6 elements 4x4x2 int.
4x4x2 int.
3x3x2 int.
2.180 2.249 2.182 2.264 S-90°
-3.820
-3.915
-3.871
-4.383 The variations of bending stress (STRESS-RR) and axial stress (STRESS-SS) along the shell are shown in figures on pages A5.4 and A5.5 for 4x4x2 integration and in figures on page A5.6 for 3x3x2 integration.
User Hints
"* Note that the stress distribution resulting from the analysis using 4x4x2 integration is more accurate than the 3x3x2 stresses.
"* Note that SOLVIA-PRE automatically assigns GLOBAL rotations for a SHELL midsurface node when the node has a boundary condition for rotation specified [2].
References
[1]
Timoshenko, S., Strength of Materials, Part I. Elementary Theory and Problems, Third Edition, D.Van Nostrand, 1955.
[2]
SOLVIA-PRE 99.0 Users Manual, Stress Analysis, Report SE 99-1, p. 7.23 - 7.34.
Version 99.0 Linear Examples A5.2
SOLVIA Verification Manual Linear Examples Version 99.0 AS CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, SHELL ORIGINAL 0-i O.S z
TIME I
x EFORCE 200 r2 EAXES=
STRESS-RST SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A5.3
SOLVIA Verification Manual Linear Examples Version 99.0 AS CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, SHELL ORIGINAL H--
0.5 z
MAX DISPL.
1.6066E-3 TIME I
X Y
LOAD 37.5 "MRR-SECTION MAX 117.2S I123. 14 26.694
-21. 528
-69.749 I*,*q-117.97
-166.19
- t.
-214.41 MIN-238.52 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A5.4
SOLVIA Verification Manual Linear Examples ASA CIRCULAR CYLINDRICAL SHELL i
TIME I I
I I
/
'7'
/
/~
/
12, 0[
2 3
UNDER LINE LOAD, 12 CUBIC SHELL TI CO o I..
.1T I 0
LO V)
/
en/ V)
F--c FI SHELLSURFACE TOP SOLVIA-POST 99.0 FI SHELLSURFACE TOP SOLVIA ENGINEERING AB Version 99.0 ASA CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, 12 CUBIC SHELL ORIGINAL F-
-H D.5 Z
MAX DISPL.
1.624SE-3 TIME I
X Y
LDAD 37.5 6
MRR-SECTION MAX 146.28 121.23 71.1I26 21.026
-29.073
- -79.173 r
-129.27 I
-179.37 i
-229.47 MIN-254.52 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB C
C en If LJ A5.5
SOLVIA Verification Manual Linear Examples ASB CIRCULAR CYLINDRICAL SHELL TTMr I
IC I
0
-I Cu /
0 4
FI SHELLSURFACE TOP SOLVIA-POST 99.0 UNDER LI C>
U)
U)
U)
(y I
U)
U) 0 C
U)
NE LOAD,
- SHELL, int. 3*3*2 1
TItE 1 1
2 3
FI SHELLSURFACE TOP SOLVIA ENGINEERING AB Version 99.0 ASB CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD,
- SHELL, int.
3*3*2 ORIGINAL H--
-- 1 0.5 Z
MAX DISPL.
1.6249E-3 TIME I
X Y
LOAD 37.S MRR-SECTION MAX 145.76 120.69 70.543 20.400
-29.744
-79.888
-130.03
-180.17 b
j-230.32 MIN-255.39 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB 0
C my I--
//
Ai
-I-i i
i I
1 1
A5.6
SOLVIA Verification Manual SOLVIA-PRE input (6 elements, 4x4x2 int.)
HEAD
'AS CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, SHELL' DATABASE CREATE SYSTEM 1
CYLINDRICAL COORDINATES ENTRIES NODE R
THETA XL 1
- 2.
2
- 2.
3
- 2.
4 2.
- 0.
0.5
- 0.
- 0.
- 90.
0.5
- 90.
0.
MATERIAL 1
ELASTIC E=2.E11 NU=.3 EGROUP 1 SHELL RINT=4 RESULTS=NSTRESSES STRESSREFERENCE=ELEMENT THICKNESS 1
.02 GSURFACE 3 1 2 4 ELT=6 EL2=1 NODES=16 SYSTEM=1 LOADS ELEMENT TYPE=FORCE INPUT=LINE 3 4 THICKNESS -200 FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES 345 INPUT=LINES 246 INPUT=LINES 156 INPUT=LINES
/
1 2
/
3 4
/
1 3
/
2 4 SET NSYMBOLS=MYNODES NNUMBERS=MYNODES MESH VECTOR=LOAD EAXES=STRESS-RST SOLVIA END SOLVIA-POST input (6 elements, 4x4x2 int.)
A5 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, SHELL DATABASE CREATE WRITE FILENAME='a5.1is' VIEW ID=1 XVIEW=1.
YVIEW=0.75 ZVIEW=0.5 SET VIEW=1 ORIGINAL=DASHED MESH CONTOUR=MRR-SECTION NSYMBOLS=MYNODES VECTOR=LOAD EPLINE NAME=FI 1
1 5 9 2 TO 6
ELINE LINENAME=FI ELINE LINENAME=FI ZONE NLIST NLIST END 1592 KIND=SRR OUTPUT=ALL SUBFRAME=21 KIND=SSS OUTPUT=ALL NAME=EDGES INPUT=NODES
/
1 TO 4 ZONENAME=EDGES DIRECTION=123 KIND=REACTION DIRECTION=34 Version 99.0 Linear Examples A5.7
SOLVIA Verification Manual EXAMPLE A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE Objective To verify the bending behaviour of the PLATE element when applied to a curved shell structure.
Physical Problem The cylindrical shell to be analyzed is the same as in example A.5. It is shown again in the figure below.
Ph
_ _F-y p
r =2.00 m h =0.02 m b =0.50 m P = 400 N/m E= 2.0.l10" N/m 2 v =0.3 Finite Element Model As in the previous example only a 900 portion of the cylinder is modeled due to symmetry reasons, see figures on page A6.3. Note that the boundary conditions are such that the displacements in the X-direction are zero. In addition, all 4 boundary faces have symmetry boundary conditions.
Solution Results The theoretical solution is the same as for example A.5. Using the input data on pages A6.5 and A6.6 the following results are obtained:
Radial displacements:
8, (mm)4 = 90 degrees 58y (mm)q = 0 degrees Theory SOLVIA Theory SOLVIA
-1.625
-1.627 1.492 1.492 The distribution of radial displacements along a line from node 3 to node 5 is shown in figure on page A6.4.
The theoretical bending moment per unit length is ([1] p. 380):
Pr (
- 2) 2 (
IrJ Version 99.0 b
Linear Examples I
A6.1
SOLVIA Verification Manual The bending moment at the point of application of the line load (=--90 degrees) can then be compared with the results obtained using the input shown on pages A6.5 and A6.6.
Bending moment (Nm/im)
Theory SOLVIA el.45 pt 2 el.46 pt 4
-255
-276
-251 The distribution of the bending moment per unit length along two lines is shown on page A6.4. The line "MID" is located in the symmetry plane of the model (X=0.25) while the line "SIDE" is located along one side of the model (X=O).
The distribution of the membrane force per unit length in the circumferential direction is also shown on page A6.4.
A Local Cylindrical System is used in displaying the bending moment (Ml 1) and the membrane force (F1l1) per unit length. The x,-direction of the Local System is in the circumferential direction.
User Hints
"* The PLATE element does not include shear deformation effects; hence even when the shell modelled is thick, shear deformation effects are not included in the model the element can reliably be employed for thin and even very thin plates/shells.
" The membrane action in the PLATE element corresponds to a constant strain assumption. Hence, the membrane forces are constant over each element, see for example the variation of the Fl11 force per unit length on page A6.4.
" Note that significantly larger moment jumps between neighboring elements occur for the line "SIDE" than for the line "MID".
Reference
[1]
Timoshenko, S., Strength of Materials, Part I. Elementary Theory and Problems, Third Edition, D Van Nostrand, 1955.
Version 99.0 Linear Examples A6.2
SOLVIA Verification Manual Linear Examples A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE ORIGINAL O.S Z
X 121 120 18 t17 16 1 is t13 MASTER B
10011 B
C 1iO11t D I 9 0 I SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE ORIGINAL F
0.5 TIME I
SOLVIA-PRE 99.0 SOLVIA E z
X Y
FORCE so NGINEERING AS Version 99.0 A6.3
SOLVIA Verification Manual Version 99.0 Linear Examples A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE TIME I
TI1E I o
I
/
°.
,/
/
U)
/.
0
/
//
U ° U )
.0
/
03 1
2 3I I /
0 a
/
LU)
S9 0
E/
/ii
(
i
/
L._
o o/
L-0 1
2 3
01 2
3 LINE3S NODE 3 S
SIDE SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE 0.TIrE I
I TIPE I
//V 7
z
/
ft o
/f SIDE MID SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A6.4
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE' DATABASE CREATE MASTER IDOF=100011 COORDINATES 1
.5
/
2
/
3
.5 0 2
/
4 0 0 2
/
5
.5 2
/
6 0 2 LINE ARC N1=5 N2=3 NCENTER=*
EL=I2 LINE ARC N1=6 N2=4 NCENTER=2 EL=12 MATERIAL 1
ELASTIC E=2.E1J NU=.3 EGROUP 1
PLATE RESULTS=TABLES STRESSTABLE 1
1 2 3 4 5 6 GSURFACE 6 4 3 5 EL1=12 EL2=1 EDATA
/
1 0.02 LOADS CONCENTRATED 3 3
-50 4 3
-50 3
FIXBOUNDARIES 34
/
5 6 FIXBOUNDARIES 24
/
3 4 FIXBOUNDARIES 234
/
1 2 MESH NSYMBOL=YES NNUMBER=YES OUTLINE=YES BCODE=ALL MESH ENUMBER=YES VECTOR=LOAD SOLVIA END Version 99.0 Linear Examples A6.5
SOLVIA Verification Manual SOLVIA-POST input A6 CIRCULAR CYLINDRICAL SHELL UNDER LINE LOAD, PLATE DATABASE CREATE SYSTEM 1 CYLINDRICAL WRITE FILENAME='a6.1is' NPLINE NAME=LINE35 3 111 STEP -1 TO 101 5 EPLINE NAME=SIDE 45 2 4 1 STEP -4 TO 1
2 4 1 SUBFRAME 21 NLINE LINENAME=LINE35 DIRECTION=3 OUTPUT=ALL SYSTEM=1 ELINE LINENAME=SIDE KIND=F11 OUTPUT=ALL SYSTEM=1 SUBFRAME 21 ELINE LINENAME=SIDE KIND=M11 OUTPUT=ALL SYSTEM=1 EPLINE NAME=MID 46 4 3 /48 3 4 /42 4 3 /44 3 4 38 4 3 /40 3 4 /34 4 3 /36 3 4 30 4 3 /
32 3 4 / 26 4 3 / 28 3 4 22 4 3 /24 3 4 /18 4 3 /20 3 4 14 4 3 / 16 3 4 /
10 4 3 / 12 3 4 643/
834/
243/
434 ELINE LINENAME=MID KIND=M11 OUTPUT=ALL SYSTEM=1 NLIST KIND=REACTION DIRECTION=34 END Version 99.0 Linear Examples A6.6
SOLVIA Verification Manual EXAMPLE A7 CANTILEVER BEAM UNDER TIP LOADS, BEAM ELEMENT Objective To verify the three-dimensional action of the straight BEAM element under end loads.
Physical Problem A cantilever beam as shown in figure below under transverse end load and torsion.
L z
Pz X
P, MY 0 Y z
x lb a
L=1 m a=0.01 m b = 0.02 m E=2.0-10 1 N/im 2 v = 0.3 P*I N
P =10 N My=1 Nm Finite Element Model The behaviour of the cantilever is described using one BEAM element only.
Solution Results The theoretical solution is for example given in [1] for the end displacement 8 and for the end rotations Ob due to bending and 0, due to torsion.
PL' PL 3EI AsG pL2 Pb 2EI (bending) t MGJ (torsion)
Version 99.0 Linear Examples A7.1
SOLVIA Verification Manual The following solution is obtained using the input data on page A7.4:
Displacements (mm)
Linear Examples 3.7500 Reference
[1]
Roark, R.J.,
Theory
- 4) 2.8424 Rotations (radians. 10-3)
-7.5000 3.7500 Formulas for Stress and Strain, Fourth Edition, McGraw Hill, 1965.
Version 99.0 SOLVIA 2 y 2.8424
-7.5000 A7.2
SOLVIA Verification Manual Linear Examples Version 99.0 A7 CANTILEVER BEAM UNDER TIP LOADS, BEAM ELEMENT ORIGINAL 0.2 z
ORIGINAL 0.1 Z
TIME I
MOMENT t
1 MASTER r
S 000000 EAXES=RST B 111111 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A7 CANTILEVER BEAM UNDER TIP LOADS, BEAM ELEMENT ORIGINAL 0.2 Z
ORIGINAL 0.2 Z
TIME I Ly TIME I
Ly MOMENT-S MOMENT-T MAX 10.000 MAX 5.0000 S9.3750 4.687S 8.1250 4.062S 6.87S0 3.4375
-~5.6250 2.8125 S4.37S0
- 2. 187S 3.1250 1.5625 1.87S0 0.93750 0.62S00 0.31250 MIN 0 MIN 0
SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A7.3
SOLVIA Verification Manual Linear Examples SOLVIA-PRE input HEAD
'A7 CANTILEVER BEAM UNDER TIP LOADS, BEAM ELEMENT' DATABASE CREATE COORDINATES 1
2
- 0.
- 1.
MATERIAL 1
ELASTIC E=2.E11 NU=.3 EGROUP 1
BEAM RESULT=FORCES SECTION 1
RECTANGULAR WTOP=0.01 D=0.02 BEAMVECTOR 1
- 1.
ENODES 1
-1 1 2 LOADS CONCENTRATED 2 1
- 5.
2 3 10.
2 5
- 1.
FIXBOUNDARIES
/
1 SET NSYMBOLS=YES NNUMBERS=YES MESH EAXES=RST SUBFRAME=21 MESH BCODE=ALL VECTOR=MOMENT SOLVIA END SOLVIA-POST input A7 CANTILEVER BEAM UNDER TIP LOADS, BEAM ELEMENT DATABASE CREATE WRITE FILENAME='a7.1is' SET VIEW=X ORIGINAL-YES DEFORMED=NO MESH CONTOUR=MS SUBFRAME=21 MESH CONTOUR=MT NLIST KIND=DISPLACEMENT NLIST KIND=REACTION ELIST END Version 99.0 A7.4
SOLVIA Verification Manual EXAMPLE A8 CANTILEVER BEAM UNDER TIP LOADS, ISOBEAM ELEMENT Objective To verify the three-dimensional action of the straight ISOBEAM element under end loads.
Physical Problem Same as in Example A7.
Finite Element Model One 4-node ISOBEAM element as shown in the top figure on page A8.3 is chosen to model the cantilever.
Solution Results The theoretical solution is the same as given in Example A7. The following SOLVIA results are obtained using the input data on page A8.4.
Displacements (mm):
Theory SOLVIA 5.0004 2.5008 5.0003 2.5007 Rotations (radians.l0-3)"
User Hints
- The 4-node ISOBEAM element has cubic assumptions for all displacements along the beam and can in this case describe the bending behaviour exactly. Special assumptions are used for the warping displacements due to torsion.
Version 99.0 Theory SOLVIA 3.7500 2.8424
-7.5000 3.7500 2.8188
-7.5000 Linear Examples A8.1
SOLVIA Verification Manual
" The ISOBEAM element is more expensive to use than the corresponding linear BEAM element.
The reason is that the ISOBEAM stiffness matrix and stresses are evaluated by numerical integration while the linear BEAM is formulated in closed form.
" Comparing the performance of the 4-node ISOBEAM element versus the BEAM element, we may note for linear analysis of straight beam sections, always use the BEAM element; for linear analysis of curved beams the ISOBEAM element can be effective (because transverse and longitudinal displacements are interpolated to the same order), but an alternative is to model a curved beam as an assemblage of straight BEAM elements (see Examples A9 and A10);
for nonlinear analysis the ISOBEAM element can be more effective, because in large displacement analysis geometry changes are modelled more accurately, (and transverse and longitudinal displacements are interpolated to the same order);
the ISOBEAM element is compatible with the isoparametric SHELL element and can be employed to model stiffeners.
Note that the nodal point forces and moments output by SOLVIA-POST are virtual work equivalent nodal point forces/moments, F=JBTt'dV v
where B is the strain-displacement matrix and t is the vector of element stresses. This vector F is equal to the externally applied forces/moments. Hence the forces/moments output by SOLVIA POST corresponding to the element internal nodes are zero, because no loads are applied (in this problem) at these nodes.
Version 99.0 Linear Examples A8.2
SOLVIA Verification Manual Linear Examples A8 CANTILEVER BEAM UNDER TIP LOADS, ISOBEAM ELEMENT ORIGINAL 0.2 Z
ORIGINAL 0.2 TIME IIy TIME I z
MOMENT-S MAX 10.000 E 9.3750 8.1250
- .8750 5.6250
- 4. 3750
- 3. 1250 1.8750 0.62S00 MIN 8.2072E-12 SOLVIA-POST 99.0 MOMENT-T MAX S.0000 E4.6875 4.062S 3.4375 2.8125 S2.187S I.S62S 0.93750 0.312S0 MIN 9.2786E-12 SOLVIA ENGINEERING AS Version 99.0 A8.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A8 CANTILEVER BEAM UNDER TIP LOADS, ISOBEAM ELEMENT' DATABASE CREATE COORDINATES 1
TO 4
- 0.
- 1. /
5 0.1 MATERIAL 1
ELASTIC E=2.E11 NU=.3 EGROUP 1
ISOBEAM RESULTS=FORCES SECTION 1
SDIM=0.01 TDIM=0.02 ENODES 1
51423 LOADS CONCENTRATED 4 1
- 5.
4 3 10.
4 5
- 1.
FIXBOUNDARIES
/
1 5 SET NNUMBERS=YES NSYMBOLS=YES MESH EAXES=RST SUBFRAME=21 MESH BCODE=ALL VECTOR=MOMENT SOLVIA END SOLVIA-POST input A8 CANTILEVER BEAM UNDER TIP LOADS, ISOBEAM ELEMENT DATABASE CREATE WRITE FILENAME='a8.1is' SET VIEW=X ORIGINAL=YES DEFORMED=NO MESH CONTOUR=MS SUBFRAME=21 MESH CONTOUR=MT NLIST DIRECTION=13456 NLIST DIRECTION=13456 KIND=REACTION ELIST END Version 99.0 Linear Examples A8.4
SOLVIA Verification Manual EXAMPLE A9 PINCHED CIRCULAR RING, ISOBEAM ELEMENTS Objective To verify the curved ISOBEAM element under bending.
Physical Problem A circular ring of square cross section is subjected to equal and opposite concentrated forces, see figure below. The thickness of the ring and the mean radius are the same as for the shell in Examples A5 and A6.
z r = 2.000 m a=0.020 m
E =2.0.0lolN/m 2
v =0.3 P =8.000 N
t A
Y Finite Element Model Symmetry gives that only a quarter of the ring needs to be modeled, see top figures on page A9.3 Four cubic ISOBEAM elements are used.
Solution Results The theoretical solution is given in [1] p. 381 as follows:
_=-1 (7c_ 2)Pr 3 2z 2(4 7c EI
=1(2 _1 Pr 3 2Trc
- 2) El where 8, is the radial displacement at the application of the load and 5Y is the radial displacement 900 from the load (along the y-axis).
Version 99.0 Linear Examples A9.1
SOLVIA Verification Manual The following displacement results are obtained:
8, (mm) 8y (mm)
Theory SOLVIA Theory SOLVIA
-1.785
-1.786 1.639 1.640 The SOLVIA results are obtained using the input data on pages A9.4 and A9.5.
The variation of radial displacements along a line from node 2 to node 1 is shown in the top figure on page A9.4.
The theoretical bending moment at the application of the load is ([1] p. 380):
Pr Comparison with the SOLVIA results yields:
Bending moment (Nm)
Theory SOLVIA
-5.09
-5.09 The variation of the bending moment is shown in the top figure on page A9.4.
User Hints
"* The SOLVIA results are obtained using three point integration along the centroidal axis of each element. The four point integration results give about 4% less displacements.
"* The present example structure is in a state of plane stress while the previous Examples A5 and A6 were plane strain solutions. As for these previous examples the corresponding 8-node PLANE element could also effectively be employed.
"* Radial displacements are calculated in SOLVIA-POST using the x3 displacements of a Local Cylindrical System.
Reference
[1]
Timoshenko, S., Strength of Materials, Part I, Elementary Theory and Problems, Third Edition, D. Van Nostrand, 1955.
Version 99.0 Linear Examples A9.2
SOLVIA Verification Manual Linear Examples A9 PINCHED CIRCULAR
- RING, ORIGINAL F-0.2 TIME I A9 PINCHED CIRCULAR RING, ISOBEAM ELEMENTS ORIGINAL 0.2 Z
L C
2y 10091 18 10111 107 106 De31 6
C11011 10111111 ALC A1 B
D.3tll SOLVA-PR 990 SOVIA NGIEAERING A A9 PINCHED CIRCULAR RING, ISOBEAM ELEMENTS ORIGINAL I-H 0.2 MAX DISPL.
I.78S6E-3 TIME I
SOLVIA-POST 99.0 SOLV ISOBEAM ELEMENTS z Ly FORCE A
SOLVIA ENGINEERING AB z
LOAD 4
IA ENGINEERING AB Version 99.0 SOLVIA-PRE 99.0 A9.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A9 PINCHED CIRCULAR RING, ISOBEAM ELEMENTS' DATABASE CREATE MASTER IDOF=100011 COORDINATES ENTRIES NODE Y
Z 1
2 2
0 2
3 LINE ARC N1=1 N2=2 NCENTER=3 EL=4 MIDNODES=2 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
SECTION 1
GLINE 1 2 ISOBEAM RESULTS=FORCES SDIM=0.02 TDIM=0.02 3
EL=4 NODES=4 LOADS CONCENTRATED 2 3 -4 FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES 34
/
24
/ /
Version 99.0 1
2 3
Linear Examples A9.4
SOLVIA Verification Manual SOLVIA-PRE input (cont.)
SET VIEW=X SMOOTHNESS=YES NSYMBOLS=YES PLOTORIENTATION=PORTRAIT MESH NNUMBERS=YES EAXES=RST BCODE=ALL MESH ENUMBER=YES VECTOR=LOAD SOLVIA END SOLVIA-POST input A9 PINCHED CIRCULAR RING, ISOBEAM ELEMENTS DATABASE CREATE SYSTEM 1 CYLINDRICAL WRITE FILENAME='a9.1is' MESH VIEW=X SMOOTHNESS=YES ORIGINAL=DASHED VECTOR=LOAD NPLINE NAME=LINE
/
2 111 STEP -1 TO 101 1
SUBFRAME 21 NLINE LINENAME=LINE DIRECTION=3 SYMBOL=1 OUTPUT=ALL SYSTEM=1 EPLINE NAME=RING 42/41/32/31 22/21/12/11 ELINE LINENAME=RING KIND=MT OUTPUT=ALL END Version 99.0 Linear Examples A9.5
SOLVIA Verification Manual EXAMPLE A10 PINCHED CIRCULAR RING, BEAM ELEMENTS Objective To verify the BEAM element when applied to a curved structure.
Physical Problem A circular ring of square cross-section subject to equal and opposite concentrated forces, see figure below, is considered. The problem is the same as in Example A9.
z I
A Y
r = 2.000 m a= 0.020 m E = 2.0- 10" N/m 2 v = 0.3 P = 8.000 N Finite Element Model Symmetry gives that only a quarter of the ring need be modelled, see the bottom figures on page A 10.2. Twelve BEAM elements are used.
Solution Results The theoretical solution is the same as for Example A9. The following SOLVIA results are obtained using the input data shown on page A 10.4.
Displacements:
8, = radial displacement at the point of load application 8y = radial displacement 900 from the point of load application The radial displacements of the ring are shown in the left bottom figure on page A10.3. The x3 direction in the Local Cylindrical System is in the radial direction.
Version 99.0 Linear Examples A10.1
SOLVIA Verification Manual Linear Examples Forces and moments:
A contour plot of the bending moment is shown in the top figure on page A10.3. The variation of the axial force is shown in the bottom figure on page A10.3.
User Hints
" The BEAM element is straight and each element has a cubic displacement assumption for bending and a linear displacement assumption for axial displacements. Hence, several elements are re quired to model a curved structure. The linear stiffness matrix is evaluated in closed form which is effective, and even when many elements are used it results in a relatively low cost of analysis.
- The axial force shown in the top figure on page A10.3 refers to the element coordinate system.
Discontinuities, therefore, occur between elements due to the different orientation of the elements.
AIO PINCHED CIRCULAR RIN ORIGINAL 0.2 TIME I SOLVIA-PRE 99.0 G, BEAM ELEMENTS z Ly 6
S 4
3 2
FORCE 4
SOLVIA ENGINEERING AB Version 99.0 AIO PINCHED CIRCULAR RING, BEAM ELEMENTS ORIGINAL ý--
0.2 Z
L C
2*
0B 10111 C2 C11011 S10301 102 101 EAXES=RST MASTER 100011 B 1011111 C 1 10 11 D
.111l SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A10.2
SOLVIA Verification Manual Version 99.0 Linear Examples A10.3 AID PINCHED CIRCULAR RING, BEAM ELEMENTS o
TIME I
TIE El U)j 2
3 o
U)o N
I-I
/~
~
I-"-
fl
?-
/
L LU/
SOVI0 PS 990IL NGNEIGA
/7 P
-]
I 0i 2
3 0
1 2
3 LINE NODE 2 -
I RING SOLVIA-POST 99.0 SOLVIA ENGINEERING AS
SOLVIA Verification Manual SOLVIA-PRE input:
HEAD
'A10 PINCHED CIRCULAR RING, BEAM ELEMENTS' DATABASE CREATE MASTER IDOF=100011 COORDINATES 1
- 0.
- 2.
/
2
- 0.
- 0.
- 2.
/
3
- 0.
LINE ARC N1=1 N2=2 NCENTER=3 EL=12 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
BE2 SECTION 1
RI GLINE 1 2 3 FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES AM RESULT=FORCES ECTANGULAR WTOP=0.02 D=0.02 EL=12 34 24 LOADS CONCENTRATED 2 3 -4 SET MESH MESH
/ /
/
1 2
3 VIEW=X PLOTORIENTATION=PORTRAIT NNUMBER=YES NSYMBOL=YES EAXES=RST BCODE=ALL ENUMBER=YES VECTOR=LOAD SOLVIA END SOLVIA-POST input:
A10 PINCHED CIRCULAR RING, BEAM ELEMENTS DATABASE CREATE SYSTEM 1 CYLINDRICAL WRITE FILENAME='al0.1is' MESH VIEW=X CONTOUR=MT ORIGINAL=YES DEFORMED=NO NPLINE NAME=LINE EPLINE NAME=RING
/ /
SUBFRAME 21 NLINE LINENAME=LINE ELINE LINENAME=RING ELINE LINENAME=RING END 2
111 STEP -1 TO 101 12 2 1 TO 1
2 1 1
DIRECTION=3 SYMBOL=1 OUTPUT=ALL SYSTEM=1 KIND=MT OUTPUT=LIST KIND=FR OUTPUT=ALL Version 99.0 Linear Examples A10.4
SOLVIA Verification Manual EXAMPLE All CURVED BEAM UNDER OUT-OF-PLANE LOAD, ISOBEAM ELEMENTS Objective To verify the three-dimensional action of curved ISOBEAM elements including torsion.
Physical Problem The figure below shows the curved beam to be analyzed. The beam is relatively slender and fixed at one end. The other end is loaded by a concentrated force in the out-of-plane direction.
E=2.0.10" N/m 2 v=0.3 R=l.0m a = 0.02 m b=0.02m P =100 N Y
Finite Element Model The figure on page Al 1.3 shows the finite element model. Six 3-node ISOBEAM elements are used to approximate the 90-degree circular bend.
Solution Results The theoretical solution is given in [1] p. 412.
Tip displacement in the direction of force:
Q PR3 (T El,(3iTY ox __W
-+cs 2 ))I El, ~4 C
where C =0.141-a 4 -G
([1] p. 290)
Torsional moment Mr and bending moment M.:
IM I=PR(1-sin4)
[M3[ = P R'cos ý Version 99.0 z
P Linear Examples All.1
SOLVIA Verification Manual A comparison with results from SOLVIA using the input data on page Al 1.5 yields:
M, (Nm)
M, (Nm)
(degree)
Theory SOLVIA Theory SOLVIA 0
100.0 99.9 100.0 100.1 15 74.1 74.1 96.6 96.6 30 50.0 50.0 86.6 86.6 45 29.3 29.3 70.7 70.7 60 13.4 13.4 50.0 50.0 75 3.4 3.4 25.9 25.9 90 0
0 0
0 The deformed shape, contour plot of torsional moment and the variation of the moments Mr and M, along the curved beam are shown in figures on page A 11.4.
User Hints
"* Torsional effects are very important in this example. The ISOBEAM element contains special displacement assumptions for the torsional behaviour.
" The nodal forces and moments at an internal node of the ISOBEAM element are zero unless the node is acted upon by an externally applied force or moment. The forces and moments of the element end nodes are, however, those which balance the forces and moments of the adjoining elements and the applied loads, see Example A8.
"* The constant curvature of the circular bend is only approximated by the 3-node (parabolic)
ISOBEAM elements. Since six elements are used for the bend the solution still agrees very well with the exact solution.
Reference
[1]
Timoshenko, S., Strength of Materials, Part I. Elementary Theory and Problems, Third Edition, D. Van Nostrand, 1955.
Version 99.0 Linear Examples Al11.2
SOLVIA Verification Manual All CURVED BEAM UNDER ORIGINAL
-i 0.2 OUT-OF-PLANE LOAD, ISOBEAM ELEMENTS z
ORIGINAL I 0.2 3
x y
4 B
FORCE t00 MASTER 000000 B ll11 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB Version 99.0 z
x 'l y
Linear Examples All1.3
SOLVIA Verification Manual Linear Examples All CURVED BEAM UNDER OUT-OF-PLANE
- LOAD, TIME I
In 0
z In ISOBEAM ELEMENTS T IPE I CA co 0
C10.
No 2.0 1.0 0.5 1.0 1.5 BEAM BEAM SOLVIA ENGINEERING AB SOLVIA-POST 99.0 0.0 0.5 1.0 1.5 2.0 Version 99.0 All CURVED BEAM UNDER OUT-OF-PLANE LOAD, ISOBEAM ELEMENTS ORIGINAL I-
-- 0*1 Z
MAX DISPL.
0.049976 TIME I
X y
LOAD LOAO 9O0 ROLVIA-POST 99.I SOLVIA ENGINEERING AR All CURVED BEAM UNDER OUT-OF-PLANE LOAD. ISOBEAM ELEMENTS ORIGINAL *
.- 0.
Z TIME I
Ly MOMENT-R MAX 99 944 93.698 81.205 68.712 43.726 31.233 18.740 6.2465 MIN-I.0301E-It SOLVIA-POST 99.0 SOLVIA ENGINEERING AB
.I/
S.
..../
/ /
/
/
/ /
/
/
/ /
/
z 0
Y:
0
(NJ o
o_
All1.4
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'All CURVED BEAM UNDER OUT-OF-PLANE LOAD, ISOBEAM ELEMENTS' DATABASE CREATE SYSTEM 1
CYLINDRICAL COORDINATES
/
ENTRIES NODE R
THETA 1
/
2
- 1.
/
3
- 1.
- 90.
MATERIAL 1
ELASTIC E=2.Ell NU=0.3 EGROUP 1
ISOBEAM RESULTS=FORCES SECTION 1
SDIM=0.02 TDIM=0.02 GLINE N1=2 N2=3 AUX=l EL=6 NODES=3 SYSTEM=l FIXBOUNDARIES
/
1 2 LOADS CONCENTRATED 3 1 100.
SET SMOOTHNESS=YES NSYMBOLS=YES SUBFRAME 21 MESH VECTOR=LOAD BCODE=ALL MESH ENUMBER=YES NNUMBERS=MYNODES SOLVIA END SOLVIA-POST input All CURVED BEAM UNDER OUT-OF-PLANE LOAD, ISOBEAM ELEMENTS DATABASE CREATE WRITE FILENAME='all.Iis' SET PLOTORIENTATION=PORTRAIT MESH ORIGINAL=DASHED VECTOR=LOAD SMOOTHNESS=YES MESH CONTOUR=MR ORIGINAL=YES DEFORMED=NO VIEW=X EPLINE NAME=BEAM 6
2 1 TO 1
2 1 SET PLOTORIENTATION=LANDSCAPE ELINE LINENAME=BEAM KIND=MR OUTPUT=ALL SUBFRAME=21 ELINE LINENAME=BEAM KIND=MS OUTPUT=ALL NLIST ZONENAME=N3 DIRECTION=I56 NLIST DIRECTION=-56 KIND=REACTION END Version 99.0 Linear Examples All1.5
SOLVIA Verification Manual EXAMPLE A12 CURVED BEAM UNDER OUT-OF-PLANE LOAD, BEAM ELEMENTS Objective To verify the three-dimensional action of the BEAM element, including torsion, when applied to a curved structure.
Physical Problem The curved beam extends 900 and is loaded by a concentrated force at one end and is fixed at the other end. The problem is the same as in Example A11, see figure below.
z
\\b E =2.0 -10"1 N/rn2 v=0.3 R=l.0m a
aa=0.02m b=0.02m P=100N
/Y Finite Element Model The figure on page A12.2 shows the finite element model. Twelve BEAM elements are used.
Solution Results The theoretical solution is given in Example A11. The following SOLVIA results are obtained using the input data on page A 12.4.
Tip displacement (in) in the direction of the force:
Version 99.0 Linear Examples A 12.1
SOLVIA Verification Manual Moments (Nm) at the built-in end transformed to the Z-axis (torsion) and the Y-axis (bending):
A contour plot of the torsional moment is shown in the top figure on page A 12.3. The variation of the torsion and the bending moments are shown in the bottom figure on page A 12.3.
User Hints
" The stiffness matrix of the linear BEAM element is evaluated in closed form, which is effective computationally. The element forces/moments or stresses are referred to the element coordinate system. Since the element is straight there are in general discontinuities in the variation of these quantities along BEAM elements modeling a curved structure.
"* Note that the torsional moment in each BEAM element is constant. The bending moment distribution is then forced to be of saw tooth form to maintain equilibrium with the moment formed by the applied force.
A12 CURVED BEAM UNDER OUT-OF-PLANE
- LOAD, BEAM ELEMENTS ORIGINAL 0.2 TIME I z
x y
FORCE 100 SOLVIA-PRE 99.0 ORIGINAL 0.2
- -,12 z
X y
\\S
\\3
- I MASTER 000000 B 1It11 SOLVIA ENGINEERING AB Version 99.0 Torsional Moment Bending Moment Theory SOLVIA Theory SOLVIA 100.0 100.0 100.0 100.0 Linear Examples A12.2
SOLVIA Verification Manual Version 99.0 Linear Examples A12 CURVED BEAM UNDER OUT-OF-PLANE LOAD, BEAM ELEMENTS ORIGINAL 0.1 Z
TIME I
L1 MOMENT-R "NO AVERAGING MAX 93.246 87.418 75.762 64.106 5.451 40.795 29.139 17.481 5.8278 MIN-2.8422E-13 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A12 CURVED BEAM UNDER OUT-OF-PLANE LOAD, BEAM ELEMENTS TI ME 1
TIME 1 I
0 A
C 0
0 It
4 0.0 O.S 1.0 I.S 2.0 0.0
- 0. S 1.0 1.5 2.0 BEAM BEAM SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A12.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A12 CURVED BEAM UNDER OUT-OF-PLANE LOAD, BEAM ELEMENTS' DATABASE CREATE SYSTEM 1 CYLINDRICAL COORDINATES
/
ENTRIES NODE R
THETA 1
/
2
- 1.
/
3
- 1. 90.
MATERIAL 1
ELASTIC E=2.E1T NU=0.3 EGROUP 1
BEAM RESULTS=FORCES SECTION 1
RECTANGULAR WTOP-0.02 D=0.02 GLINE 2 3 1 EL=12 SYSTEM=1 FIXBOUNDARIES
/
1 2 LOADS CONCENTRATED 3 1 100.
SET NSYMBOLS=YES SUBFRAME 21 MESH NNUMBERS=MYNODES VECTOR=LOAD MESH ENUMBER=YES BCODE=ALL SOLVIA END SOLVIA-POST input A12 CURVED BEAM UNDER OUT-OF-PLANE LOAD, BEAM ELEMENTS DATABASE CREATE WRITE FILENAME='a12.1is' CONTOUR AVERAGE=NO MESH CONTOUR=MR ORIGINAL=YES DEFORMED=NO VIEW=X EPLINE NAME=BEAM 12 2 1 TO 1
2 1 ELINE LINENAME=BEAM KIND=MR OUTPUT=ALL SUBFRAME=21 ELINE LINENAME=BEAM KIND=MS OUTPUT=ALL NLIST ZONENAME=N3 DIRECTION=156 NLIST DIRECTION=156 KIND=REACTION END Version 99.0 Linear Examples A12.4
SOLVIA Verification Manual EXAMPLE A13 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LOAD Objective To verify the TRUSS element when used with SKEW degree-of-freedom directions.
Physical Problem A cantilever truss structure under concentrated end load is considered as shown in the figure below.
z L = 5.000 m h = 0.500 m I = 0.500 m A = 0.0001 m2 P = 1000 N E =2.0.10" N/m 2 d =
h'-- + 2 Y
Finite Element Model The truss structure is modelled with 2-node TRUSS elements as shown in the figures on page A 13.3.
The structure is inclined in the global coordinate system so it is convenient to employ SKEW degree of-freedom directions.
Solution Results The structure is statically determinate and the internal forces are as shown in the figure below.
-o-p
-0L
-2P, hP Version 99.0 Linear Examples A13.1
SOLVIA Verification Manual The displacement of the point of load application (node 22) in the direction of the cantilever is:
9ipl_
= 1.125.10-3 (im) i=1 hAE The SOLVIA numerical solution obtained using input data on pages A13.5 and A13.6 is as follows:
Internal element forces (N):
Element Theory SOLVIA 1
-10000.00
-10000.00 11 9000.00 9000.00 21 1414.21 1414.21 Displacement 8 (mm):
The deformed mesh is shown in the top figure on page A 13.4. A contour plot of the internal element forces, R-force, is shown in the bottom figure on page A 13.4.
User Hints
- The SKEW degree-of-freedom System used for each node is oriented in the principal directions of the cantilever. The input concentrated force and all nodal results are then referred to this coordinate system.
Version 99.0 Linear Examples A13.2
SOLVIA Verification Manual Linear Examples AD AI3 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LO ORIGINAL i
0.5S Z
0 40 0
9 10I' MASTER 100111 B l1lltl AREA S E IE-4 SOLVIA ENGINEERING AB SOLVIA-PRE 99.0 Version 99.0 ELEMENTS A13 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LOAD ORIGINAL F--
0.5 Z
TIME I
L S~NODES FORCE 1000 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A13.3
SOLVIA Verification Manual Linear Examples A13 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LOAD z
ORIGINAL H-
--- O.S MAX DISPL..
0.017743 TIME I
LOAD 1000 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB Version 99.0 A13 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LOAD MAX DISPL.
0.017743 Z
TIME I
L FORCE-R NO AVERAGING MAX 9000.0
- 7812.S 5437.5 3062.5
'0,00 687.50
- lip,
-1687.5
, -4062.S
-88 12.5 MIN-1OO0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A13.4
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A13 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LOAD' DATABASE CREATE MASTER IDOF=-00111 SKEWSYSTEMS EULERANGLES 1 30 SYSTEM 1
CARTESIAN PHI=30 COORDINATES ENTRIES NODE YL ZL 1
- 0.
- 0.
TO 11
- 5.
- 0.
12
- 0.
0.5 TO 22
- 5.
0.5 NSKEWS 1 1 TO 22 1 MATERIAL 1
ELASTIC E=2.E11 NU=0.
EGROUP 1 TRUSS ENODES ENTRIES EL NI N2 1
1 2
TO 10 10 11 11 12 13 TO 20 21 22 21 12 2
TO 30 21 11 31 13 2
TO 40 22 11 EDATA ENTRIES EL AREA 1
.0001 LOADS CONCENTRATED 22 3 -1000 FIXBOUNDARIES
/
1 12 SET VIEW=X HEIGHT=0.25 MESH NNUMBERS=YES NSYMBOLS=YES VECTOR=LOAD TEXT STRING='NODES' XPT=10.
YPT=2.
MESH ENUMBERS=YES BCODE=ALL CONTOUR=AREA TEXT STRING='ELEMENTS' XPT=10.
YPT=2.
SOLVIA END Version 99.0 Linear Examples A13.5
SOLVIA Verification Manual SOLVIA-POST input A13 CANTILEVER TRUSS STRUCTURE UNDER CONCENTRATED LOAD DATABASE CREATE WRITE FILENAME='a13.1is' SET VIEW=X MESH ORIGINAL=DASHED VECTOR=LOAD CONTOUR AVERAGE=NO MESH CONTOUR=FR NLIST NLIST ELIST END ZONENAME=N22 KIND=REACTIONS Version 99.0 Linear Examples A13.6
SOLVIA Verification Manual EXAMPLE A14 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLANE STRESS Objective To verify the PLANE STRESS element subjected to distributed loading and when employing SKEW degree of freedom systems.
Physical Problem:
A cantilever beam of rectangular cross-section is considered when subjected to a triangular distributed load, see figure below.
L = 1.000 m p = 0.200.10 6 N/m2 h = 0.100mE=2.0-10.N/m 2
b = 0.050 my = 0.3 Finite Element Model The figures on page A 14.3 shows the finite element model. The upper face of the elements is acted upon by a linearly varying pressure load. The model is inclined in the Global System.
Solution Results Using beam theory the theoretical solution for the end displacement is:
11 PL3 5 PL 60 El 6 AG Version 99.0 Linear Examples A14.1
SOLVIA Verification Manual Linear Examples The input data, as given on pages A14.5 and A14.6 yields the following solution:
End displacement (mm):
Theory SOLVIA
-1.113
-1.110 The axial stress (N / m2, in the direction of the beam) for element 10 and node 3 is:
Theory SOLVIA 40.00.106 39.77.106 User Hints
"* The exact displacement solution according to beam theory for a linearly varying distributed load contains the 5-th power of the coordinate. Several 8-node (parabolic) PLANE elements are, therefore, used to obtain a good approximation to the analytical solution.
" For the PLANE STRESS 2 element the stress results are calculated in the global Y and Z directions. By specifying STRESSREFERENCE = ELEMENT in the DATABASE command in SOLVIA-POST the program transforms element stress results to the Element Stress System for PLANE and SOLID elements [1]. Alternatively, the stresses could be referenced to a Local System defined by the SYSTEM command, see example A15.
"* Note that the summation of reaction forces and applied loads in SOLVIA-POST is made in the global coordinate system, i.e., nodal forces referring to skew systems are transformed to the global coordinate system before summation.
Reference
[1]
SOLVIA-POST 99.0 Users Manual, Report SE 99-1, pp. 5.1-5.6.
Version 99.0 A14.2
SOLVIA Verification Manual A14 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS.
PLANE STRESS z
ORIGINAL ý-
0.1 TIME I
S 4
I tA PRESSURE 200000 z
ORIGINAL
ý-H 0.1 L
y B,
MASTER 100t11 SOLVIA-PRE 99.0 B 110111 C 111111 SOLVIA ENGINEERING AB Version 99.0 Linear Examples LY A 14.3
SOLVIA Verification Manual Version 99.0 Linear Examples A14.4
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A14 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLANE STRESS' DATABASE CREATE MASTER IDOF=100111 SET MYNODES=10 SKEWSYSTEMS EULERANGLES 1 30 SYSTEM 1 CARTESIAN PHI=30 COORDINATES ENTRIES NODE YL ZL 1
0 0
2 0
0.05 3
0 0.1 4
1 0
5 1
0.1 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
PLANE STRESS2 RESULTS=NSTRESSES GSURFACE 5 3 1 4 EL1=10 EL2=1 NODES=8 EDATA ENTRIES EL THICK 1
0.05 LOADS ELEMENT INPUT=LINE 3
5
- 0.
2.E5 FIXBOUNDARIES 23 /
2 FIXBOUNDARIES 2
/
1 3 NSKEWS INPUT=SURFACE 5314 1
SET PLOTORIENTATION=PORTRAIT NSYMBOLS=MYNODES SUBFRAME 12 MESH NNUMBERS=MYNODES VECTOR=LOAD MESH ENUMBERS=YES BCODE=ALL SOLVIA END Version 99.0 Linear Examples A14.5
SOLVIA Verification Manual SOLVIA-POST input A14 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLANE STRESS DATABASE CREATE STRESSREFERENCE=ELEMENT WRITE FILENAME='a14.1is' SET PLOTORIENTATION=PORTRAIT VIEW ID=1 XVIEW=I ROTATION=-30 SET VIEW=1 SUBFRAME 13 MESH VECTOR=LOAD MESH ORIGINAL=DASHED GSCALE=OLD EPLINE NAME=TOP 10 2 5 1 TO 1
2 5 1 ELINE LINENAME=TOP KIND=SRR OUTPUT=ALL SUMMATION KIND=LOAD SUMMATION KIND=REACTION DETAILS=YES NLIST ZONENAME=EL1 EMAX SELECT=S-EFFECTIVE END Version 99.0 Linear Examples A14.6
SOLVIA Verification Manual Linear Examples EXAMPLE A15 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLANE STRAIN Objective To verify the PLANE STRAIN element subjected to distributed loading and when employing a SKEW degree-of-freedom System.
Physical Problem Same as in Example A14, see figure on page A14-1.
Finite Element Model Same as in Example A 14 except that PLANE STRAIN elements are used, see figure on page A 15.2.
Solution Results To obtain the same theoretical solution as for the plane stress case of Example A 14 the following material data is used:
E*= 1+2v E=1.89349"10 11N/m 2
(1+ v V*=
= 0.230769 1+v where E and v are the Young's modulus and Poisson's ratio used in Example A14.
The input data on page A 15.4 gives the following solutions:
End displacement (mm):
The axial stress (N/m 2 in the direction of the beam) for element 10 at node 3 is:
Version 99.0 A15.1
SOLVIA Verification Manual Since the PLANE STRAIN element in SOLVIA gives stress results in the Global System, a Local System with the x,-direction oriented in the axial direction was defined using the SYSTEM command.
Alternatively, SOLVIA-POST could be employed to transform the global stresses to the Element Stress System.
User Hints
- Note that the thickness of the PLANE STRAIN element is always equal to unity. The pressure is, therefore, here acting on a 20 times wider cantilever than in Example A14, giving a factor of 20 times larger reaction forces than for Example A 14.
AIS CANTILEVER UNDER DISTRIBUTED LOAD ORIGINAL 0
1 TIME I
SOLVIA-PRE 99.0 USING SKEW SYSTEMS, PLANE STRAIN z
L-y PRESSURE 200000 A
r S
EAXES=
STRESS-RST SOLVIA ENGINEERING AB Version 99.0 Linear Examples A15.2
SOLVIA Verification Manual Version 99.0 Linear Examples A15.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A15 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLANE STRAIN' DATABASE CREATE MASTER IDOF=100111 SKEWSYSTEMS EULERANGLES 1 30 SYSTEM 1 CARTESIAN PHI=30 COORDINATES 1
TO 3
- 0.
- 0.
0.1 4
- 0.
- 1. /
5
- 0.
- 1. 0.1 MATERIAL 1
ELASTIC E=1.89349E11 NU=0.230769 EGROUP 1
PLANE STRAIN RESULTS=NSTRESSES GSURFACE 5 3 1 4 EL1=10 EL2=1 NODES=8 LOADS ELEMENT INPUT=LINE 3 5
- 0.
2.E5 FIXBOUNDARIES 23 2
FIXBOUNDARIES 2
/
3 NSKEWS INPUT=SURFACE 5314 1
SET NSYMBOLS=MYNODES MESH NNUMBERS=MYNODES EAXES=STRESS-RST VECTOR=LOAD SOLVIA END Version 99.0 Linear Examples A15.4
SOLVIA Verification Manual SOLVIA-POST input A15 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLANE STRAIN DATABASE CREATE SYSTEM 1 CARTESIAN NA=4 NB=3 WRITE FILENAME='alS.1is' MESH ORIGINAL=DASHED VECTOR=LOAD CONTOUR=S11 SYSTEM=1 EPLINE NAME=TOP 10 2 5 1 TO 1
2 5 1 ELINE LINENAME=TOP KIND=S11 OUTPUT=ALL SYSTEM=1 SUBFRAME=21 ELINE LINENAME=TOP KIND=S33 SYSTEM=1 SUMMATION KIND=LOAD SUMMATION KIND=REACTION DETAILS=YES NLIST ZONENAME=EL1 EMAX SELECT=S-EFFECTIVE END Version 99.0 Linear Examples A 15.5
SOLVIA Verification Manual EXAMPLE A16 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SOLID Objective To verify the SOLID element subjected to pressure loading and when employing SKEW degree-of freedom Systems.
Physical Problem Same problem as shown in Example A14, see figure on page A14-1.
Finite Element Model The model is shown in the figure on page A 16.2 and the top figure on page A 16.3. The elements are 20-node SOLID elements which have a quadratic assumption for the displacements, also in the transverse direction of the cantilever.
Solution Results The theoretical solution for the end displacement is the same as given in Example A 14.
The input data shown on pages A16.4 and A16.5 gives the following results:
End displacement (mm):
The axial stress (N/m2, in the direction of the beam) for element 10 at node 2 is:
Since SOLVIA gives the stresses referred to the Global System, SOLVIA-POST is employed to transform the global stress components to the Element Stress System by specifying STRESSREFERENCE = ELEMENT in the DATABASE command. Alternatively, the stresses could be referenced to a local system defined by the SYSTEM command, see example A15.
User Hints
- The midside nodes in the transverse direction of the beam are included in order to model the anticlastic curvature. For a description of the phenomenon of anticlastic bending, see for example
[1] p. 175.
Version 99.0 Linear Examples A16.1
SOLVIA Verification Manual Linear Examples If 2x2x2 integration is employed spurious (zero energy) modes are excited in the transverse direction. These spurious modes could be avoided by introducing additional fixed boundary conditions for the nodes at the built-in end of the beam. However, this change would also introduce additional stresses in the beam.
Reference
[1]
Oden, J.T., Ripperger, E.A., Mechanics of Elastic Structures, Second Edition, McGraw-Hill, 1981.
A16 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SOLID ORIGINAL I O.OS Z
x y
a b
NAXES=SKE\\W SOLVIA ENGINEERING AB SOLVIA-PRE 99.0 Version 99.0 A16.2
SOLVIA Verification Manual Version 99.0 A16 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS.
SOLID MAX DISPL.
1.l13iE-3 S
TIME I MISES MAX 4.042SE7
'3.7904E7 2.7820E7 2.2862E7 2.2777E7 S1.773SE7
- i1.2693E7 S7.6506E6 2.6084E6 MIN 87291 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB I
Linear Examples A16.3
SOLVIA Verification Manual Linear Examples SOLVIA-PRE input HEAD
'A16 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SOLID' DATABASE CREATE MASTER IDOF=000111 SYSTEM 1
CARTESIAN PHI=30 COORDINATES 1
-0.025 1. 0.1 2
-0.025 0. 0.1 3
0.025 0. 0.1 4
0.025 1. 0.1 5
-0.025
- 1. 0.
6
-0.025
- 0. 0.
7 0.025 0. 0.
8 0.025 1. 0.
SKEWSYSTEMS 1 30 EULERANGLES MATERIAL 1
ELASTIC E=2.E11 NU=0.3 Version 99.0 A16 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SOLID 10TIMý 1I TINý
\\\\
\\
\\,
I
\\
ofo 0.0 0.2 0.4 0!6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 10 EDGE EDGE SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A16.4
SOLVIA Verification Manual SOLVIA-PRE input (cont.)
EGROUP 1
SOLID RESULTS=NSTRESSES GVOLUME 1 2 3 4 5 6 7 8
EL1=10 EL2=1 EL3=1 NODES=20 LOADS ELEMENT INPUT=SURFACE 2
3 4 1
- 0.
- 0.
2.E5 2.E5 NSKEWS INPUT=SURFACE 2367 1
1458 1
FIXBOUNDARIES 123 INPUT=SURFACE
/
2 3 6 7 FREEBOUNDARIES 3
INPUT=LINES
/
2 3
/
6 7 FREEBOUNDARIES 1
INPUT=LINES
/
2 6
/
3 7 SET NSYMBOLS=MYNODES MESH NNUMBERS=MYNODES NAXES=SKEW MESH CONTOUR=PRESSURE ENUMBERS=YES BCODE=ALL SOLVIA END SOLVIA-POST input AI6 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SOLID DATABASE CREATE STRESSREFERENCE=ELEMENT WRITE FILENAME='a16.1is' MESH CONTOUR=MISES NSYMBOLS=MYNODES OUTLINE=YES EPLINE NAME=EDGE 10 2 9 1 TO 1
2 9 1 SUBFRAME 21 ELINE LINENAME=EDGE KIND=SRR OUTPUT=ALL ELINE LINENAME=EDGE KIND=ERR ZONE NAME=TIPNODES INPUT=NODES 1 4 5 8 NLIST ZONENAME=TIPNODES NLIST KIND=REACTIONS EMAX NUMBER=3 END Version 99.0 Linear Examples A16.5
SOLVIA Verification Manual EXAMPLE A17 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, BEAM Objective To verify the BEAM element subjected to distributed loading and when employing skew coordinate systems.
Physical Problem Same as in Example A14, see figure on page A14-1.
Finite Element Model The model is shown in the top figure on page A17.2. Ten BEAM elements are used to model the cantilever.
Solution Results The theoretical solution for the end displacement is the same as given in Example Al14.
The input data on page A17.3 gives the following results:
End displacement (mm):
The axial stress ( N/m 2, in the direction of the beam) on the top surface, closest to the fixed end is:
User Hints The BEAM elements do not predict the exact displacement in this example since its displacement assumption along the beam is a cubic polynomial while the exact displacement varies with the 5th power of the axial coordinate. However, 10 BEAM elements give a very good solution to this example problem.
- Note that the distributed load for the BEAM element is force per unit length.
Version 99.0 Linear Examples 17.1
SOLVIA Verification Manual Version 99.0 A17 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, BEAM ORIGINAL 0
0.1 Z
TIME i L,
8 6
7 EFORCE 10000 t
r S
EAXES=RST InI \\VTA-PRF 9Q 0 SOLVIA ENGINEERING AB A17 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, BEAM ORIGINAL
-- 0.1 z
MAX DISPL.
-i 1.1104E-3 TIME I
y LOAD 900
.TIM I
TIM I
0.0 0.2 0.4 0,6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 TOP MID SOLVIA-POST 99.0 SOLVIA ENGINEERING AB Linear Examples 17.2
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A17 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, BEAM' DATABASE CREATE SKEWSYSTEMS EULERANGLES 1 30 SYSTEM 1
CARTESIAN PHI=30 COORDINATES 1
/
2
- 0. 1.
MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
BEAM RESULT=STRESSES SECTION 1
RECTANGULAR WTOP=0.05 D=0.1 BEAMVECTOR 1
-1.
GLINE N1=1 N2=2 AUX=-1 EL=10 LOADS ELEMENT TYPE=FORCE INPUT=LINE 1 2
-T
- 0.
1.E4 NSKEWS INPUT=LINE 12 1
FIXBOUNDARIES
/
1 SET VIEW=X NSYMBOLS=MYNODES MESH ENUMBERS=YES NNUMBERS=MYNODES EAXES=RST VECTOR=LOAD SOLVIA END SOLVIA-POST input A17 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, BEAM DATABASE CREATE WRITE FILENAME='a17.1is' VIEW ID=1 XVIEW=1 ROTATION=-30 SET VIEW=1 MESH VECTOR=LOAD ORIGINAL=DASHED SUBFRAME=12 EPLINE NAME=TOP
/
1 1 9 TO 10 1 9 EPLINE NAME=MID
/
1 6 14 TO 10 6 14 ELINE LINENAME=TOP KIND=SRR OUTPUT=ALL SUBFRAME=2121 ELINE LINENAME=MID KIND=SRT SUBFRAME=2221 NLIST DIRECTION=34 NLIST KIND=REACTIONS DIRECTION=34 SUMMATION KIND=LOAD END Version 99.0 Linear Examples 17.3
SOLVIA Verification Manual EXAMPLE A18 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, ISOBEAM Objective To verify the ISOBEAM element subjected to distributed loading and when employing SKEW degree-of-freedom Systems.
Physical Problem Same as in Example A14, see figure on page A14-1.
Finite Element Model The model is shown in the figures on page A18.2. Ten parabolic ISOBEAM elements are used to model the cantilever.
Solution Results The theoretical solution for the end displacement is the same as given in Example A14.
The input data shown on pages A 18.4 and A 18.5 gives the following results:
End displacement (mm):
The axial stress (N/m 2, in the direction of the beam) on the top surface, closest to the fixed end is:
SOLVIA-POST is employed to calculate bending and shear stresses from element bending moments and shear forces.
The deformed mesh and a contour plot of the bending moment and the variation of axial stress and shear stress along the top surface of the beam are shown on page A 18.3.
User Hints
" Note that the SOLVIA predicted displacements agree well with the analytical solution even though 3-node ISOBEAM elements were used (with a parabolic displacement assumption in each beam element).
"* The forces and moments calculated by SOLVIA for the internal nodes of the elements only balance the external loads applied to these nodes. They are, therefore, not equal to the stress resultants obtained by integrating the stresses through the beam thickness (see Example A8).
Version 99.0 Linear Examples A18.1
SOLVIA Verification Manual Linear Examples The forces and moments calculated for element end nodes are, however, in equilibrium with the externally applied forces/moments plus those supplied by adjoining elements. Note that the exact forces and moments for the node at the fixed end are calculated.
A18 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, ISOBEAM ORIGINAL i
-i 0.1 Z Ly t
r S
EAXES=RST SOLVIA ENGINEERING AB AI8 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, ISOBEAM ORIGINAL I
0.1 TIME I SOLVIA-PRE 99.0 y
KY "to 0ý7 6
S-FORCE MAX 10000 S9375.0 8125.0 687S.0 S.. 62S.0 S4375.i50 1875 0 625.00 MIN 0 SOLVIA ENGINEERING AB
.3 SOLVIA-PRE 99.0 Version 99.0 A 18.2 03 3
Version 99.0 A18.2
SOLVIA Verification Manual Linear Examples Version 99.0 AI8 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, ISOBEAM ORIGINAL
-- 1 0.1 z
MAX DISPL.
H 1.1087E-3 TIME I
. Y LOAD 633.33 ORIGINAL
- 0. 1 MAX DISPL.
H 1.i087E-3 MOMENT-T TIME I
MAX 3333.3 S3125.0 2708.3 2291.7 1875.0 1458.3 t1041.7 625.00 208.33 MIN-1.1416E-I1 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A18.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A18 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, ISOBEAM' DATABASE CREATE SKEWSYSTEMS EULERANGLES 1 30 SYSTEM 1
CARTESIAN PHI=30 COORDINATES 1
/
2
- 0.
- 1.
/
3
- 0.
- 0.
-0.1 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
ISOBEAM RESULTS=FORCES SECTION 1
SDIM=0.1 TDIM=0.05 GLINE 1 2 3 EL=10 NODES=3 LOADS ELEMENT TYPE=FORCE INPUT=LINE 1 2 S
- 0.
1.E4 NSKEWS INPUT=LINE 12 1
FIXBOUNDARIES
/
1 3 SET VIEW=X NSYMBOLS=YES NNUMBERS=MYNODES MESH EAXES=RST MESH ENUMBERS=YES CONTOUR=SFORCE SOLVIA END Version 99.0 Linear Examples A18.4
SOLVIA Verification Manual SOLVIA-POST input Ai8 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, ISOBEAM DATABASE CREATE WRITE FILENAME='aI8.1is' VIEW ID=1 XVIEW=1 ROTATION=-30 SET VIEW=1 ORIGINAL=DASHED MESH ORIGINAL=DASHED VECTOR=LOAD SUBFRAME=12 MESH CONTOUR=MT AXES=NO EPLINE NAME=TOP 1
1 2 TO 10 1 2 EVARIABLE NAME=TMOMENT TYPE=ISOBEAM KIND=MT EVARIABLE NAME=SFORCE TYPE=ISOBEAM KIND=FS CONSTANT NAME=W VALUE=83.3333E-6 CONSTANT NAME=A VALUE=5.OOOOE-3 RESULTANT NAME=STRESSBB STRING='
ABS(
TMOMENT /
W )'
RESULTANT NAME=STRESSBC STRING='-ABS(
SFORCE
/
A )'
AXIS ID=1 LABEL='LENGTH COORDINATE' AXIS ID=2 LABEL='STRESS-BB' AXIS ID=3 LABEL='STRESS-BC' SUBFRAME 21 RLINE LINENAME=TOP RESULTANTNAME=STRESSBB XAXIS=1 YAXIS=2, OUTPUT=ALL RLINE LINENAME=TOP RESULTANTNAME=STRESSBC XAXIS=1 YAXIS=3, OUTPUT=ALL NLIST ZONENAME=ELT0 DIRECTION=34 NLIST KIND=REACTIONS DIRECTION=34 END Version 99.0 Linear Examples A18.5
SOLVIA Verification Manual EXAMPLE A19 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLATE Objective To verify the PLATE element subjected to distributed loading and when employing SKEW degree-of freedom Systems.
Physical Problem Same problem as in Example A14, see figure A14-1.
Finite Element Model The model is shown on page A19.2. Twenty PLATE elements are used to model the cantilever.
Solution Results The theoretical solution for the end displacement is the same as given in Example A14.
The input data on pages A 19.4 and A 19.5 give the following results:
The moments at the fixed end are:
End displacement (mm):
The deformed mesh and the variation of the bending moment along the beam (at stress points at the middle and at the edge) are shown in the figures on page A19.3.
User Hints
"* Since the PLATE element model does not include the strain energy due to transverse shear strains in the formulation and the Kirchhoff hypothesis of zero transverse shear strains is imposed at discrete points, the element is applicable to thin (including very thin) plates and shells.
"* Note that smaller stress jumps are obtained for the line along the edge than for the line along the middle of the cantilever, see figures on page A 19.3. A finer mesh is necessary to reduce the stress jumps.
"* Section forces/moments per unit length are obtained as the basic element result for the PLATE element. Using the RESULTANT command a transformation to the axial stress values is achieved.
Version 99.0 Theory SOLVIA (Nm)
(Nm/m)
(Nm) 3333 67528 3376 (element 4, at node 1) 3333 64442 3222 (element 1, at node 1) 3333 64442 3222 (element 1, at node 2)
Linear Examples A19.1
SOLVIA Verification Manual Linear Examples A19 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLATE ORIGINAL 0.05 TIME I
.4 z
X ',Y PRESSURE 200000 SOLVIA ENGINEERING AB SOLVIA-PRE 99.0 Version 99.0 A19 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLATE ORIGINAL 0.02 SOLVIA-PRE 99.0 SOLVIA ENGINEE x
Y MASTER 000000 B 011i0i C 1111 1 RING AB A19.2
A19 CANTILEVER UNDER DISTRIBUTED I...
TIME I
0.2 0.4 EDGE 0.6 0.8 1.(
LOAD U)
IL of N
0 SOLVIA-POST 99.0 USING SKEW SYSTEMS, PLATE TIME 1 I
I I
.0 0.2 0.4 0.6 0.8 I.C MIDDLE SOLVIA ENGINEERING AB Version 99.0 SOLVIA Verification Manual Linear Examples AI9 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLATE ORIGINAL h-
-- 0.05 Z
MAX DISPL.
1.1012E-3 TIME I X'
Y LOAD 600 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB CD u)
C)0 C')
o LU V-)
I C,)
N 041 0.0
ý -I--
0 I
A19.3
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A19 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLATE' DATABASE CREATE COORDINATES 1
- 0.
2 0.05 3
0.05 0.8660254 0.5 4
- 0.
0.8660254 0.5 SKEWSYSTEM EULERANGLES 1 30 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
PLATE GSURFACE 1 2 3 4 EL1=1 EL2=5 STRESSTABLE 1
1 2 3 4 5 6 7 EDATA
/
1 0.1 LOADS ELEMENT INPUT=SURFACE 1 2 3 4 T
- 0.
- 0.
2.E5 2.E5 NSKEWS INPUT=SURFACE 1234 1
FIXBOUNDARIES 12346 1
FIXBOUNDARIES 2346
/
2 SET NSYMBOLS=MYNODES MESH VIEW=I OUTLINE=YES VECTOR=LOAD NNUMBERS=MYNODES VIEW ID=1 XVIEW=1 YVIEW=3 ZVIEW=1 MESH VIEW=1 ENUMBER=YES BCODE=ALL SOLVIA END Version 99.0 Linear Examples A19.4
SOLVIA Verification Manual SOLVIA-POST input A19 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, PLATE DATABASE CREATE SYSTEM 1 CARTESIAN NA=2 NB=4 WRITE FILENAME='a19.1is' SET NSYMBOLS=MYNODES MESH ORIGINAL=DASHED VECTOR=LOAD EPLINE NAME=EDGE 4
2 4 1 STEP 4 TO 20 2 4 1 EPLINE NAME=MIDDLE 1
473
/
3 374 5
473
/
7 374 9
4 7 3
/
11 3 7 4 13 4 7 3
/
15 3 7 4 17 4 7 3
/
19 3 7 4 F
EVARIABLE FORCE PLATE F22 EVARIABLE MOMENT PLATE M22 CONSTANT THICK 0.1 CONSTANT SIX 6.0 RESULTANT STRESS22
'FORCE/THICK+SIX*MOMENT/(THICK*THICK)'
RLINE LINENAME=EDGE RESULTANT=STRESS22 OUTPUT=ALL SYSTEM=1 SUBFRAME=21 RLINE LINENAME=MIDDLE RESULTANT=STRESS22 OUTPUT=ALL SYSTEM=1 ELIST ZONENAME=ELI SYSTEM=1 NLIST DIRECTION=34 NLIST KIND=REACTION DIRECTION=34 END Version 99.0 Linear Examples A19.5
SOLVIA Verification Manual EXAMPLE A20 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SHELL Objective To verify the SHELL element when subjected to distributed loading and when employing SKEW degree-of-freedom Systems.
Physical Problem Same problem as in Example A14, see figure A14-1.
Finite Element Model The model is shown in the figure on page A20.2. It consists of five 9-node SHELL elements.
Solution Results The theoretical solution for the end displacement is the same as given in A14.
The input data on pages A20.4 and A20.5 gives the following results:
Displacement (mm):
The axial stress (N/m 2 ), in the direction of the beam) at node 3 is:
The deformed mesh and the variation of the axial stress and strain along the beam is shown on page A20.3.
User Hints
- All shell director vectors are in this example generated automatically by SOLVIA-PRE and the shell rotations of nodes with rotation boundary conditions specified are automatically referenced to the Global (or SKEW) directions [1].
- In this example the Element Stress System is oriented in the principal directions of the beam, since all shell director vectors are orthogonal to the mid-surface. The stresses are, therefore, conveniently requested to be given in the Element Stress System.
Version 99.0 Linear Examples A20.1
SOLVIA Verification Manual
- If 2x2x2 integration is employed, spurious modes are present in the model, compare Example A16. The default integration orders must, therefore, be used in this example and are also recommended for general use.
Reference
[1]
SOLVIA-PRE 99.0 Users Manual, Stress Analysis, Report SE 99-1, p. 7.25.
Version 99.0 A20 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SHELL ORIGINAL C.
Z TIME I
x Y
PRESSURE 200000 MASTER 000000 B li110t C 111ilt SOLVIA-PRE 99-0 SOLVIA ENGINEERING AB Linear Examples A20.2
SOLVIA Verification Manual Version 99.0 Linear Examples A20.3 A20 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SHELL
.1.
- TIMý I
I TIMý I 0
C!UC)
-ý C!)
Z H -
i.
. i
.I
0-In 0.0 0.2 0.1
. 0.6 0.8 1.0 0 0 0.2 0,'
0.6 0.8 1.0 TOP SHELLSURFACE TOP TOP SHELLSURFACE TOP SOLVIA-POST 99.0 SOLVIA ENGINEERING AB
SOLVIA Verification Manual SOLVIA-PRE input HEAD
'A20 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SHELL' DATABASE CREATE SYSTEM 1 CYLINDRICAL COORDINATES ENTRIES NODE R
THE¶ 1
1 2
1 3
0 TA XL 30 0.05 30 0
30 0
4 0
30 0.05 5
0 30 0.025 SKEWSYSTEMS EULER-ANGLES 1
30 MATERIAL 1
ELASTIC E=2.E11 NU=0.3 EGROUP 1
SHELL STRESSREFERENCE=ELEMENT RESULTS=NSTRESSES THICKNESS 1
0.1 GSURFACE 3 4 1 2 EL1=1 EL2=5 NODES=9 LOADS ELEMENT TYPE=PRESSURE INPUT=SURFACE 3 4 1 2 T
- 0.
- 0.
2.E5 2.E5 NSKEWS INPUT=SURFACE 1234 1
FIXBOUNDARIES 12346 FIXBOUNDARIES SET NSYMBOLS=MYNODES MESH NNUMBER=MYNODES SOLVIA END
/
3 4
/5 ENUMBERS=YES BCODE=ALL VECTOR=LOAD Version 99.0 Linear Examples A20.4
SOLVIA Verification Manual SOLVIA-POST input A20 CANTILEVER UNDER DISTRIBUTED LOAD USING SKEW SYSTEMS, SHELL DATABASE CREATE WRITE FILENAME='a20.1is' MESH ORIGINAL=DASHED CONTOUR=MISES EPLINE NAME=TOP 1
1 8 4 TO 5
1 8 4 ELINE LINENAME=TOP KIND=SSS OUTPUT=ALL SUBFRAME=21 ELINE LINENAME=TOP KIND=ESS NLIST DIRECTION=34 NLIST KIND=REACTION DIRECTION=34 EMAX SELECT=S-EFFECTIVE END Version 99.0 Linear Examples A20.5
SOLVIA Verification Manual EXAMPLE A21 PLANAR TRUSS Objective To verify the TRUSS element when applied to a 2-dimensional structure.
Physical Problem The planar truss structure shown in figure below is considered.
Pz
-py I =20 ft h= 15 ft E-A = 30-10 4 kips-ft2/ft2 P, = 10 kips Y
PY = 20 kips Finite Element Model The model is shown in the top figure on page A21.2. The 2-node TRUSS element is used to model the elements in this problem.
Solution Results The theoretical displacement solution is given in [1]I p. 257. Using the input data on page A21.3 the following results are obtained:
Nodal displacements (10-3 feet):
Node Theory SOLVIA y-dir.
z-dir.
y-dir.
z-dir.
1 0
0 0
0 2
1.48
-4.38 1.48
-4.38 3
2.96
-2.50 2.96
-2.50 4
3.70 0
3.70 0
5 3.00
-4.38 3.00
-4.38 6
3.59
-2.08 3.59
-2.08 A contour plot of the element forces is shown in the bottom figure on page A21.2.
Reference
[1]
Tuma, J.J., Munshi, R.K., Theory and Problems of Advanced Structural Analysis, Schaum's Outline Series, McGraw-Hill.
Version 99.0 T
Ih 74,I 7ý Linear Examples A21.1 Z
z,
SOLVIA Verification Manual Linear Examples Version 99.0 A21 PLANAR TRUSS ORIGINAL F--
S.
Z 8
1 2
3 MASTER 100111 B 101111 C 11ll1 AREA I
SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A21 PLANAR TRUSS z L_¥ MAX DISPL. '
5.3085E-3 TIME I
5.
FORCE-R NO AVERAGING MAX 22.222 1i9.S83 13.929 8.2738 2.6190 MI-3.03S7
- -8.6905
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A21 PLANAR TRUSS' DATABASE CREATE MASTER IDOF=100111 COORDINATES 1
TO 4
- 0.
- 60.
/
5
- 0.
- 20.
15.
/
6
- 0.
- 40.
15.
MATERIAL 1
ELASTIC E=30.E4 EGROUP 1
TRUSS ENODES 1
1 2 TO 3
3 4 4
5 1 TO 6
5 3 7
63
/
8 65 EDATA
/
1
- 1.
FIXBOUNDARIES FIXBOUNDARIES LOADS 53 62 23 3
/
9 6 4
/1
/
4 CONCENTRATED 10 20 MESH VIEW=X NSYMBOLS=YES NNUMBERS=YES ENUMBERS=YES, BCODE=ALL CONTOUR=AREA SOLVIA END SOLVIA-POST input A21 PLANAR TRUSS DATABASE CREATE WRITE FILENAME='a21.lis' NLIST NLIST KIND=REACTION SET VIEW=X CONTOUR AVERAGE=NO VMIN=-20 MESH CONTOUR=FR END Version 99.0 Linear Examples A21.3
SOLVIA Verification Manual EXAMPLE A22 TAPERED CANTILEVER UNDER TIP LOAD Objective To verify the SHELL element with variable thickness and when used as a transition element connected to a SOLID element.
Physical Problem A cantilever of variable thickness and subjected to a moment load at the end, see figure below, is considered.
z L
L= 100 in.
b=10 in.
M 12 h = 10/(l+9.y/L) in.
=3.0.10 7 psi h
V=0 m = 1000 in.Ibf Finite Element Model The model is shown in the figures on page A22.3. It consists of seven SHELL elements and two SOLID elements. The SHELL elements use a quadratic displacement assumption. A linear thickness variation of each SHELL element is assumed. Reduced integration in the rs-plane and the closed Newton-Cotes method in the thickness direction are used.
Solution Results Beam theory gives the following displacement and top surface bending stress solution along the cantilever:
+( 9y1-05 45y I
5.81.104 100) 100 6.M y (y) --b.h 2 The SOLVIA numerical solution obtained using the input data shown on pages A22.6 and A22.7 is as follows:
End displacement (inch):
Version 99.0 Linear Examples A22.1
SOLVIA Verification Manual A vector plot of principal stresses for the top layer of integration points and a contour plot of the von Mises effective stress are also shown on page A22.4. The top figure on page A22.5 shows the devia tion of effective stress from the nodal mean value.
Bending and shear stresses along the shell elements using reduced and full integration are shown in the bottom figures on page A22.5. The bending stress is shown together with the theoretical solution.
User Hints
"* Poisson's ratio is set to zero in order to simulate the plane stress condition associated with beam theory.
" The material law for the SHELL element is such that the stress in the thickness direction is zero.
The stiffness in the thickness direction of a transition SHELL element is, therefore, small at the transition nodes. However, an adjoining SOLID element provides stiffness in the thickness direction.
" SOLVIA-PRE automatically assigns GLOBAL rotations for SHELL midsurface nodes with applied moments specified. The applied moments at the end nodes are then ensured to be acting about the global X-direction [1].
" Since the SHELL element midsurface nodes have no rotational stiffness about the thickness direction, the rotational stiffness about the Z-direction at the end nodes must be fixed [l].
" The plotting of contours for SHELL elements is always performed on a surface formed by the nodes defining the shell elements. When only midsurface nodes are used, contour plots are always located on the SHELL midsurface although the contour results are valid for the TOP, MID or BOTTOM shell surface. For a SHELL transition element, thus an element employing top and bottom nodes at some locations as in this example, then the contours are plotted on the surface formed by the most distant of the nodes defining the SHELL element. For further information, see the commands CONTOUR and SHELLSURFACE in the SOLVIA-POST Users Manual.
" The calculated stresses are less accurate than the calculated displacements. The exact displace ments vary as a 5th order polynomial, while the SOLID and SHELL elements are only capable to describe a parabolic variation.
The approximation in stresses is illustrated by the level of stress jumps (stress deviations) between elements, see the figure on page A22.5. Note, however, that the calculated stresses at the tip of the cantilever are also approximated although no stress jump can be shown there since it is part of the boundary of the structure.
" Linear variation of the SHELL thickness is used in the model to simplify the input data. The use of a complete SHELL thickness table will increase the accuracy of the solution.
" In this example the reduced integration of the SHELL elements improves the stress distribution.
Reference
[1]
SOLVIA-PRE 99.0 Users Manual, Stress Analysis, Report SE 99-1, p. 7.25.
Version 99.0 Linear Examples A22.2
SOLVIA Verification Manual Linear Examples A22A TAPERED CANTILEVER UNDER TIP LOAD, INTEGRATION 3*3*3 z
x y
ORIGINAL
- 5.
TIME 1
,01 109 111 I5 EMOMENT 100 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB Version 99.0
,j 3 A22.3
SOLVIA Verification Manual Linear Examples Version 99.0 A22 TAPERED CANTILEVER UNDER TIP LOAD ORIGINAL S
- 5.
Z TIME I X Y SPRINCIPAL SHELL TOP S75.68
-575.68 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A22 TAPERED CANTILEVER UNDER TIP LOAD ORIGINAL S.
Z TIME I X
Y MISES SHELL TOP MAX 597.4S I560.19 485.66 411.13 336.60 262.07 187.54 113.02 38.488 MIN 1.2233 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A22.4
SOLVIA Verification Manual Linear Examples 1'0 20 30 40 SO SHELL DISTANCE 6G 7q TIME I 0
150 20 30 40 slo 0
70 ESHELL SOLVIA-POST 99.0 SOLVIA ENGINEERING AS full integration A22A TAPERED CANTILEVER UNDER TIP LOAD, o4 1
0 2o0 30 40 so 60 7
SHELL DISTANCE ESHELL SOLVIA-POST 99.0 SOLVIA ENGINEERING AB Version 99.0 A22 TAPERED CANTILEVER UNDER TIP LOAD ORIGINAL 5.z TIME I
DEVIATION OF MISES SHELL TOP MAX 16.521 I
15.488
- 13. 423 11. 358 9.2931 7.2280 5.1628 S3.0977
- 1. 0326 MIN 0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB reduced integration 0 I+
t t
r.
r
- T A22 TAPERED CANTILEVER UNDER TIP LOAD TIME I INTEGRATION 3-3-3 TIME I
7
0
1--
-C o 8.
g-o d A22.5
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A22 TAPERED CANTILEVER UNDER TIP LOAD' DATABASE CREATE SET MYNODES=200 COORDINATES ENTRIES NODE Y
Z 1
- 0.
- 5.
/
2
- 0.
- 0.
4 7.5 2.98507
/
5 7.5
-2.98507 6
- 15.
2.12766
/
7
- 15.
- 0.
9 22.5 1.65289
/
10 22.5 -1.65289 11
- 30.
1.35135
/ 12
- 30.
-1.35135 13
- 35.
- 0.
TO 14
- 40.
- 0.
15 100.
NGENERATION TIMES=2 NSTEP=50 XSTEP=-5 1 3 6 8 11 TO 14 NGENERATION NSTEP=i00 XSTEP=-10 2
4 5 7 9 10 15 MATERIAL 1
ELASTIC E=3.E7 EGROUP 1
SOLID ENODES ENTRIES EL N1 N13 1
106 105 2
111 110 EGROUP 2
ENODES ENTRIES GSURFACE THICKNESS THICKNESS THICKNESS THICKNESS THICKNESS THICKNESS THICKNESS N2 N14 101 53 106 58 N3 Ni5 1
5 6
10 N4 Ni6 6
58 11 62 N5 Ni7 108 107 112 0
N6 NI8 103 102 108 107
/
3
- 0.
-5.
/
8
- 15.
-2.12766 N7 Ni9 3
2 8
7 N8 N9 N20 8
104 7
12 109 N10 Nil N12,
51 4
56 9
SHELL RINT=2 SINT=2 TINT=-3 EL 1
14 1
2 3
4 5
6 7
N1 N2 N3 N4 N6 N8 N17 N2 ii 14 114 i11 64 61 12 6
15 115 114 EL1=6 NODES=9 T1=2. 703 T1=2. 174 T1=I. 818 T1=I 563 TI=1. 370 Ti=I 220 TI=I.099 EDATA
/
ENTRIES EL NTH 1 1 TO 7 7 T2=2. 174 T2=I. 818 T2=I. 563 T2=1. 370 T2=I.220 T2=1 099 T2=I 000 T3=2. 174 T3=I. 818 T3=I. 563 T3=1. 370 T3=1. 220 T3=I. 099 T3=I. 000 24 N20 N5 N13 N7 i2 112 13 63 113 T4=2.703 T4=2.174 T4=1.818 T4=1.563 T4=I 370 T4=1. 220 T4=I 099 FIXBOUNDARIES FIXBOUNDARIES
/
1 TO 4
6 INPUT=LINE /
LOADS ELEMENT TYPE=MOMENT INPUT=LINE 15 115 edge -100 SET MESH MESH MIDSURFACE=NO SMOOTHNESS=YES NSYMBOL=MYNODES NNUMBERS=MYNODES VECTOR=LOAD CONTOUR=THICKNESS BCODE=ALL ENUMBERS=YES SOLVIA END Version 99.0 56
/
61 15 115 Linear Examples A22.6
SOLVIA Verification Manual SOLVIA-POST input A22 TAPERED CANTILEVER UNDER TIP LOAD DATABASE CREATE WRITE FILENAME='a22.1is' SET MIDSURFACE=NO ORIGINAL=YES DEFORMED=NO OUTLINE=YES EGROUP 2 EPLINE ESHELL /
1 12 6 TO 7 12 6 MESH VECTOR=SPRINCIPAL MESH CONTOUR=MISES MESH CONTOUR=SDEVIATION PLINES=ESHELL AXIS 1 VMIN=O VMAX=70 LABEL='SHELL DISTANCE' AXIS 2 VMIN=O.
VMAX=600 LABEL='STRESS-YY' USERCURVE 1 /
READ A22A.DAT SET PLOTORIENTATION=PORTRAIT ELINE ESHELL KIND=SYY SYMBOL=1 XAXIS=1 YAXIS=2 SUBFRAME=12 PLOT USERCURVE 1 XAXIS=-1 YAXIS=-2 SUNFRAME=OLD ELINE ESHELL KIND=SYZ SYMBOL=1 NLIST MYNODES DIRECTION=34 NLIST KIND=REACTION DIRECTION=2 EMAX SELECT=S-EFFECTIVE END Version 99.0 Linear Examples A22.7
SOLVIA Verification Manual EXAMPLE A23 STIFFENED PLATE CANTILEVER UNDER TIP LOAD Objective To verify the use of the ISOBEAM element as a stiffener for the shell element using rigid links.
Physical Problem A cantilever of channel cross-section is loaded by concentrated end loads as shown in the figure below.
STIFFENED P LATE FINITE
/FIXED EDG P/2>
ELEMENT 'Z
/
MODEL p = 200 lbf E = 3.0-.10 7 psi F e m0 Finite Element Model The model is shown in the figure on page A23.2. The nodes 1, 5, 9 and 13 are fixed. Each ISOBEAM node is coupled to the corresponding SHELL node using a rigid link. Using this model with a linear displacement assumption in the X-direction for the SHELL element, no shear-lag effects are modeled, which allows an easy comparison with theoretical results from beam theory. Poisson's ratio is for the same reason set to zero.
Solution Results The theoretical displacement is PL3 3EI Version 99.0 Linear Examples A23.1
Linear Examples SOLVIA Verification Manual The input data on page A23.4 gives the following result:
End displacement 6 (inch):
Theory SOLVIA 0.02520 0.02523 Contour plots of the effective stress and the bending section moment in the SHELL portion are shown in the figures on page A23.3.
User Hints
- Note that SOLVIA-PRE automatically assigns GLOBAL rotations for the shell nodes connected to rigid links and for nodes with a specified boundary condition for rotation. The procedure is further described in the command SHELLNODES of the SOLVIA-PRE 99.0 Users Manual, Stress Analysis, p. 7.23.
Version 99.0 A23 STIFFENED PLATE CANTILEVER UNDER TIP LOAD ORIGINAL 1
i0.
Z TIME I N
Y N
FORCE 100 t
r EAXES=RST SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A23.2
SOLVIA Verification Manual Linear Examples Version 99.0 A23 STIFFENED PLATE CANTILEVER UNDER TIP LOAD MAX DISPL. i- 0.025245 Z
TIME 1
X Y
REACTION 2227 MRR-SECTION MAX 13.122
- 12.302 10.662 9.0217 7.3814 5.7411 2.4605 0.82015 MIN-6.3860E-13 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A23.3
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A23 STIFFENED PLATE CANTILEVER UNDER TIP LOAD' DATABASE CREATE COORDINATES 1
- 0.
- 0.
- 3.
TO 4
- 0.
120.
- 3.
5
- 0.
- 0.
6.5 TO 8
- 0.
120.
6.5 9
- 40.
- 0.
- 3.
TO 12
- 40.
120.
- 3.
13
- 40.
- 0.
6.5 TO 16
- 40.
120.
6.5 MATERIAL 1
ELASTIC E=3.E7 EGROUP 1
SHELL STRESSREFERENCE=GLOBAL THICKNESS 1
1.0 ENODES ENTRIES EL N1 N2 N3 N4 N5 N7 N9 Nil 1
8 5 13 16 7 14 6
15 EGROUP 2
ISOBEAM SECTION 1
SDIM=6.
TDIM=i.
STRESSTABLE 1
1ii 121 112 122, 211 221 212 222, 311 321 312 322 ENODES ENTRIES EL AUX N1 N2 N3 N4 1
5 1
4 2
3 2 13 9 12 10 11 LOADS CONCENTRATED 8 3 -100
/
16 3 -100 FIXBOUNDARIES
/
1 5 13 9 RIGIDLINK 2 6
/
3 7
/
4 8
/
10 14
/
ii 15
/
12 16 MESH EAXES=RST NSYMBOLS=YES NNUMBERS=YES VECTOR=LOAD SOLVIA END SOLVIA-POST input A23 STIFFENED PLATE CANTILEVER UNDER TIP LOAD DATABASE CREATE WRITE FILENAME='a23.1is' MESH CONTOUR=MISES VECTOR=LOAD MESH CONTOUR=MRR VECTOR=REACTION NLIST DIRECTION=234 NLIST KIND=REACTION DIRECTION=234 EMAX SELECT=S-EFFECTIVE END Version 99.0 Linear Examples A23.4
SOLVIA Verification Manual EXAMPLE A24 ANALYSIS OF SPHERICAL DOME UNDER SELF WEIGHT Objective To verify the PLANE STRESS3 (membrane) element to model a spherical surface when subjected to gravity loading.
Physical Problem A spherical roof (hemisphere) subjected to gravity loading is considered as shown in the figure below.
R = 4.5 m h=0.08m E = 2.0.1010 N/m2 v=0.3 p = 3000 kg/m3 g = 9.81 m/s 2 Finite Element Model A sector of 30 degrees of the spherical surface is modeled using 80 parabolic PLANE STRESS3 elements. SKEW degree of freedom Systems are used to define the boundary conditions in the non global circumferential direction as seen in the figures on pages A24.2 and A24.3. A consistent mass matrix assumption is used in the analysis.
Solution Results The analytical solution to this problem can be found in many text books, e.g., [ 1 ]. For the membrane solution of this problem it is important that the support reactions are acting in the tangent direction to the meridians.
A Local Spherical System is defined with the x, and x2 axis in the meridional and the tangential directions, respectively. The input data of pages A24.4 and A24.5 gives the following results:
At the top of the spherical dome (Z = 4.5)
Version 99.0 Stress-11 Stress-22 [kPa]
Analytical
-66.2
-66.2 SOLVIA
-66.4
-66.4 Linear Examples A24.1
SOLVIA Verification Manual Linear Examples At the support of the spherical dome (Z = 0)
Stress-11 Stress-22 [kPa]
Analytical
-132.4 132.4 SOLVIA
-131.4 132.9 Contour plots and stress distributions along the spherical surface can be seen in the figures on pages A24.3 and A24.4.
User Hints
"* A relative large number of elements must be used in the meridional direction to describe the stress variation in this example. In the circumferential direction, however, only one element can be used.
" Note that a flat PLANE STRESS3 element has no stiffness in the direction orthogonal to its plane in a small displacement analysis.
Reference
[1]
Timoshenko, S.P. and Woinowsky-Krieger, S., Theory of Plates and Shells, Second Edition, McGraw-Hill, 1959, p. 436.
Version 99.0 A24 ANALYSIS CF SPHERICAL DOME UNDER SELF WEIGHT ORIGINAL 0.S Z
CE 00Y O*
MASTER BO001111 C 010til D Sill H E 110111 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A21 ANALYSIS OF SPHERICAL DOME UNDER SELF WEIGHT ORIGINAL 0.5 Z
ZONE LINEI-2 I
X,,
t EAXES=
STRESS-RST SOLVIA-PRE 99.0 SOLVIA ENGINEERING AS A24.2
SOLVIA Verification Manual Linear Examples A24 ANALYSIS OF SPHERICAL DOME UNDER SELF WEIGHT ORIGINAL O.S Y
X NAXES=SKEW SOLVIA ENGINEERING AB SOLVIA-PRE 99.0 A24 ANALYSIS OF SPHERICAL DOME UNDER SELF WEIGHT M AXDISPL
-S.
6694E-S z
TIME I x
T STRESS-11 SYSTEM I
MAX-66238.2
-70310.0
-78453.7
-86S97.4
-94741.1
-102885 1-1028 19172
,.127316 MIN-131388 SOLVIA ENGINEERING AB A24 ANALYSIS OF SPHERICAL DONE UNDER SELF WEIGHT MAX DISPL.
5.8694E-S Z
TIME I REACTION 3995.8 STRESS-22 SYSTEM I
MAX 1.3291E5 1.2049ES 70661 45747 20833 4080. 4 28994
-53908 IN-66364 STI VIA-PST 99.0 SOLVIA ENGINEEZTNG AS Version 99.0 SOLVIA-POST 99.0 A24.3
SOLVIA Verification Manual Linear Examples A24 ANALYSIS OF SPHERICAL TIME 1
0 2
4 EDGE SOLVIA-POST 99.0 U)
U)
N N U)
U)
Lu Cr F
U) 6 DOME UNDER SELF WEIGHT I
I T I. E I
0 C,
0
/
/
/
EDGE SOLVIA ENGINEERING AB SOLVIA-PRE input HEADING
'A24 ANALYSIS OF SPHERICAL DOME UNDER SELF WEIGHT' DATABASE CREATE MASTER IDOF=000111 ANALYSIS TYPE=STATIC MASSMATRIX=CONSISTENT SYSTEM 1 SPHERICAL COORD I NATES ENTRIES NODE R
THETA PHI 1
4.5
- 0.
- 0.
2 4.5
- 0.
- 90.
3 4.5
- 30.
- 90.
SKEWSYSTEMS VECTORS 1
0.8660254 0.5 0.
-0.5 0.8660254 0.
MATERIAL 1
ELASTIC E=2.E10 NU=0.3 DENSITY=3000.
EGROUP 1
PLANE STRESS3 RESULTS=NSTRESSES EDATA
/
1 0.08 GSURFACE 1 2 3 1 EL1=20 EL2=4 NODES=8 SYSTEM=1 BLENDING=ANGLES NSKEWS INPUT=LINE
/
1 3 1
Version 99.0 0
0 00 U) 0")
U)
Lu U)
N 0 -t or A24.4
SOLVIA Verification Manual SOLVIA-PRE input (cont.)
FIXBOUNDARIES 2
INPUT=LINES
/
1 2
/
1 3 FIXBOUNDARIES 3
INPUT=LINES
/
2 3 FIXBOUNDARIES 1
INPUT=NODES
/
1 LOADS MASSPROPORTIONAL ZFACTOR=1.
ACCGRA=-9.81 SET NSYMBOLS=MYNODES PLOTORIENTATION=PORTRAIT MESH BCODE=ALL NNUMBERS=MYNODES ZONE NAME=LINE1-2 INPUT=ELEMENTS
/
1 TO 20 MESH ZONENAME=LINE1-2 EAXES=STRESS-RST SET PLOTORIENTATION=LANDSCAPE MESH VIEW=Z NAXES=SKEW SOLVIA END SOLVIA-POST input A24 ANALYSIS OF SPHERICAL DOME UNDER SELF WEIGHT DATABASE CREATE SYSTEM 1 SPHERICAL WRITE FILENAME='a24.1is' EPLINE NAME=EDGE 1
1 5 2 TO 20 1 5 2 SET PLOTORIENTATION=PORTRAIT MESH CONTOUR=S11 SYSTEM=1 MESH CONTOUR=S22 VECTOR=REACTION SYSTEM=1 SET PLOTORIENTATION=LANDSCAPE ELINE LINENAME=EDGE KIND=S11 OUTPUT=ALL SYSTEM=1 SUBFRAME=21 ELINE LINENAME=EDGE KIND=S22 OUTPUT=ALL SYSTEM=4 NLIST ZONENAME=N1 SUMMATION KIND=LOAD MASS-PROPERITES END Version 99.0 Linear Examples A24.5
SOLVIA Verification Manual EXAMPLE A25 CLAMPED SQUARE PLATE UNDER PRESSURE LOAD Objective To verify the behaviour of the SHELL element subjected to pressure load.
Physical Problem A thin plate with clamped edges under pressure load as shown in the figure below is considered. The load is acting in the negative Z-direction.
Y A
B 7
7//,// /7/I/
/,
/
/
/
/
/0
/
/
/
/
E =2.1. 10"N / m 2 v =0.3 a=l.8 m h = 0.01 m a
p = 1.0_10 3 N / M2 (pressure)
/
Finite Element Model Because of symmetry conditions only a quarter of the plate is considered, A-B-C-D in the figure above. Nine cubic SHELL elements are used for the finite element model shown in figures on page A25.3. The model allows no deformation in the X-and Y-directions and no rotation about the Z-axis.
Solution Results The theoretical solution of this problem is given in [1] p. 197 and the expression for the central deflection is given as w = 0.00126-p. a4 /D where p - pressure load a = side length D=
Eh3 12(1-v2 )
Version 99.0 I,
I' I
'I I
I I
1 I
I I
I I
I I
z
/Z (thickness)
A25.1 Linear Examples
////t///
SOLVIA Verification Manual The expression for the top surface stress at points C and D are:
6..paa (Yxxc =yyc =-0.0231.
h---2 C*xxD = 0.0513" 2
The following numerical solutions have been obtained using the input data shown on page A25.5.
Case Displ. (m)
Stress
,,x (N/ m2) point C point C point D Theory
-6.88-10-4
-4.49.106 9.97.106 16-node SHELL, 3x3
-6.91.10-4
-4.61.106 9.23.106 In addition the following solutions have been obtained for comparison:
16-node SHELL, 6x6
-6.91.10-4
-4.48.106 9.79.106 4-node SHELL, 3x3
-6.77-10-4
-4.60.106 5.13.106 4-node SHELL, 6x6
-6.88.10-4
-4.49.106 7.08.106 4-node SHELL, 12x12
-6.90.10-4
-4.46.106 8.38.106 PLATE, 3x3 (quad.)
-7.06.10-4
-4.61.106 10.6.106 PLATE, 6x6 (quad.)
-6.95.10-4
-4.49.106 10.6.106 User Hints
"* Note that the 16-node SHELL elements provide a reliable and accurate solution.
"* The solutions obtained using the 4-node SHELL element are good regarding displacements but more approximate regarding stresses. Considerable stress jumps occur between neighbouring elements since the stresses are rather constant over each element.
"* The solutions obtained using PLATE elements are good but not as accurate as the solutions obtained by the 16-node SHELL elements. The PLATE element behaves well in bending. In membrane action the PLATE is, however, a constant strain element.
"* Note that the 16-node SHELL solutions require much more CPU time than the 4-node SHELL solutions and the PLATE solutions.
Reference
[1]
Timoshenko, S.P., Woinowsky-Krieger, S., Theory of Plates and Shells, Second Edition, McGraw-Hill, 1959.
Version 99.0 Linear Examples A25.2
SOLVIA Verification Manual Linear Examples A25 CLAMPED SQUARE PLATE UNDER PRESSURE LOAD ORIGINAL I
TIME I z
X PRESSURE 1000 SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A2S CLAMPED SQUARE PLATE UNDER PRESSURE LOAD ORIGINAL 0.1 z
X Y
MASTER 1o1o01 B 110011 C I110D1 D 110,,1 E
II1 11 SOLVTA ENGINEERING AB SOLVIA-PRE 99.0 Version 99.0 Y
A25.3
SOLVIA Verification Manual Linear Examples A2S CLAMPED SQUARE TIME I I
I..
I
.1 Co 4
N 0
N B LINE-DC SHELLSURFACE TOP SOLVIA-POST 99.0 0.2 0.4 0.6 0.8 PL 1.0 ATE U 10 C,
Cl)
NDER PRESSURE LOAD TTMF i
0.0 0.2 0.,
0.6 0.8 1.0 LINE-DC SHELLSURFACE FOP SOLVIA ENGINEERING AB Version 99.0 A2S CLAMPED SQUARE PLATE UNDER PRESSURE LOAD ORIGINAL
- 0.
OZ MAX DISPL.
i 6.91E-4 TIME I
X y
MISES
"';*i SHELL TOP MAX 8.2055E6
- 7.6927E6 6.6670E6 S4..
..5.6413E6 4.61S6E6 3.5899E6
- 2.5642E6 i. 5385E6 S. 1285E5 MIN 0
SOLVIA-POST 99.0 SOLVIA ENGINEERING AB w
X X
Li C-1 U)
I-
].0 A25.4
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A25 CLAMPED SQUARE PLATE UNDER PRESSURE LOAD' DATABASE CREATE PARAMETER
$ELX=3
$ELY=3 $SHELL=16 MASTER IDOF=410001 COORDINATES 1
0.9 0.9
/
2
- 0. 0.9
/
3
/
4 0.9 MATERIAL 1
ELASTIC E=2.1E11 NU=0.3 EGROUP 1
SHELL STRESSREFERENCE=GLOBAL RESULTS=NSTRESSES THICKNESS 1
0.01 GSURFACE 1 2 3 4 EL1=$ELX EL2=$ELY NODES=$SHELL FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES 4
5 INPUT=LINES INPUT=LINES INPUT=LINES LOADS ELEMENT TYPE=PRESSURE INPUT=SURFACE 1 2 3 4 T
1000 SET MESH MESH NSYMBOLS=MYNODES VECTOR=LOAD NNUMBERS=MYNODES ENUMBERS=YES BCODE=ALL SOLVIA END SOLVIA-POST input A25 CLAMPED SQUARE PLATE UNDER PRESSURE LOAD DATABASE CREATE WRITE FILENAME='a25.1is' MESH ORIGINAL=DASHED OUTLINE=YES CONTOUR=MISES EPLINE NAME=LINE-DC 7
4 11 7 3
TO 9
4 11 7 3 ELINE LINENAME=LINE-DC KIND=SXX ELINE LINENAME=LINE-DC KIND=SYY NMAX NUMBER=5 END OUTPUT=ALL SUBFRAME=21 OUTPUT=ALL Version 99.0
/
1 2
/
3 4
/
2 3
/
1 4 Linear Examples A25.5
SOLVIA Verification Manual Linear Examples A2SA CLAMPED SQUARE PLATE UNDER PRESSURE
- LOAD, 36 CUBIC SHELL ORIGINAL
-1 0.1 MAX DISPL.
6.9112E-4 TIME I
z X
Y MISES SHELL TOP MAX 8.7019E6 I
8.158iE6 7.0703E6 S5.9826E6 4.8948E6 3.8071E6 2.7194E6 1.6316E6 5.4387E5 MIN 0 SOLVIA ENGINEERING AB SOLVIA-POST 99.0 Version 99.0 A2SA CLAMPED SQUARE PLATE UNDER PRESSURE LOAD, 36 CUBIC SHELL T IMý I
.I TIMý I
\\
\\,
,C\\
I N\\
0.4 0.6 08 L.
0.0 0!2 0.4 0.6 0.81.0 LINE-DC SHELLSURFACL TOP LINE-DC SHELLSURFACE TOP SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A25.6
SOLVIA Verification Manual Linear Examples A2SB CLAMPED SQUARE PLATE UNDER PRESSURE LOAD, 9 4-NODE SHELL x X 0.0 0.
(II 4J 0t 0.0 0.2 0.4 0.6 0.8
[.0 LINE-DC SHELLSURFACE TOP SOLVIA-POST 99.0 0.2 0.4 0.6 0.8 1.0 LINE-DC SHELLSURFACE TOP SOLVIA ENGINEERING AB Version 99.0 A2SB CLAMPED SQUARE PLATE UNDER PRESSURE LOAD, 9 4-NODE SHELL ORIGINAL ý-- -1 0.I z
MAX DISPL.
-6.7703E-4 TIME I
MISES SHELL TOP MAX 4.6023E6 P
4.3146E6 3.7393E6 S3.1641E6 2.5888E6 2.0135E6 1.4382E6 8.6293E5 2.8764E5 MIN 0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A25.7
SOLVIA Verification Manual Version 99.0 Linear Examples A25.8
SOLVIA Verification Manual Version 99.0 Linear Examples A2SD CLAMPED SQUARE PLATE UNDER PRESSURE LOAD, 144 4-NODE SHELL ORIGINAL i-- --1 0.1 z
MAX DISPL.
6.9024E-4 TIME I
X Y
MISES 4ý'*
SHELL TOP MAX 7.44SOE6 S6. 9797E6
'6.0490E6 S.
1.84E6 4.1878E6 3.2572E6 2.3266E6 S1. 3959E6 4.6531E5 MIN 0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A2SD CLAMPED SQUARE PLATE UNDER PRESSURE
- LOAD, 144 4-NODE SHELL ko TIME I
M TIr' O*
0 C'!
l
\\j 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.1 0.6 0
LINE-DC SHELLSURFACE TOP LINE-DC SHELLSURFACE SOLVIA-POST 99.0 SOLVIA ENGINEE 3.8 I.0 TOP
[RING AB A25.9
SOLVIA Verification Manual Version 99.0 Linear Examples A25. 10
SOLVIA Verification Manual Linear Examples Version 99.0 A2SF CLAMPED SQUARE PLATE UNDER PRESSURE LOAD, 6x6 PLATE I---6.9497E-4 MAX DISPL.
TIME I
REACTION 109.58 DISPLACEMENT MAX 6. 9497E-4
- 6.
5153E-4 5.6466E-4 4.7779E-4 3.9092E-4
- 3. 040SE-4 2.1718E-4 1.3031E-4 4.3436E-S MIN 0 NEERING AB SOLVIA ENGI SOLVIA-POST 99.0 A2SF CLAMPED SQUARE PLATE UNDER PRESSURE LOAD, 6x6 PLATE TIMý I D
1D I
I TIMý I I*
CD 0.0 0.2 0.1 0.6 0.8 1.0 0.0 0.2 0.4
- 0.
6
- 0. 8 1.0 LINE-DC LINE-DC SOLVIA-POST 99.0 SOLVIA ENGINEERING AS A25.11
SOLVIA Verification Manual EXAMPLE A26 MATERIAL DAMPING IN MODAL SUPERPOSITION Objective To demonstrate the calculation of weighted modal damping factors from modal damping factors of materials in a structure.
Physical Problem A beam structure with a light top part and a heavy bottom part as shown in the figure below is considered in one plane. The natural frequencies and the weighted modal damping factors are calculated.
Top part of the structure L2 =lm Rectangular section 5 x 50 mm E2 = 2.1010 N/m 2 V2 =0 P2 = 1000 kg/m 3 Modal damping ý2 = 5%
Bottom part of the structure L, =5m Rectangular section 50 x 50 mm E, =2-10" N/m 2 V, =0 Pi = 7800 kg/mi3 Modal damping I,
= 1%
Finite Element Model The model consists of two element groups of ISOBEAM elements as shown on page A26.2.
Solution Results The input data on page A26.3 is used to calculate the natural frequencies and stiffness weighted modal damping factors. The mass weighted modal damping factors are calculated in a separate run.
Version 99.0 L2 2
Mode no.
Frequency Modal damping Modal damping K-weighted M-weighted Hz 1
1.62 1.01 1.07 2
3.63 4.97 4.94 3
10.23 1.02 1.03 4
22.65 4.93 4.94 Linear Examples L
\\\\" *\\x\\
A26.1
SOLVIA Verification Manual Linear Examples The mode shapes are shown below. The bottom part dominates modes I and 3 which consequently have a weighted modal damping factor near 1 %. Similarly, the top part dominates modes 2 and 4 which then have damping factors near 5 %.
A26 MATERIAL DAMPING IN MODAL SUPERPOSITION, K-WEIGHTED ORIGINAL i I1.
B SOLVIA-PRE 99.0 z Iy 2
'V MASTER OOODI I B I11111 SOLVIA ENGINEERING AB A26 MATERIAL DAMPING IN REFERENCE -
-t0O5 MAX DISPL. -
0.35087 MODE I FRED 1.6232 DISPLACEMENT 0.35087 SOLVIA-POST 99.0 MODAL SUPERPOSITION. K-WEIGHTED Z
REFERENCE -
0O.R L
r MAX DrSPL. ý-3.9825 MODE 2 FRED 3.6261 F-K 7 DISPLACEMENT 3.9825 SOLVIA ENGINEERING AR Version 99.0 A26.2
SOLVIA Verification Manual SOLVIA-PRE input DATABASE CREATE HEADING
'A26 MATERIAL DAMPING IN MODAL SUPERPOSITION, K-WEIGHTED' MASTER IDOF=-00011 FREQUENCIES SUBSPACE-ITERATION NEIG=4 MODALDAMPING INPUT=K-WEIGHTED ANALYSIS TYPE=DYNAMIC IMODS=1 NMODES=4 COORDINATES ENTRIES NODE Y Z
1 0
0 2
0 5
3 0
6 4
.1 0 MATERIAL 1
ELASTIC E=2.EII NU=0.
DENSITY=7800.
MODALDAMPING=*.
MATERIAL 2
ELASTIC E=2.El0 NU=0.
DENSITY=7000.
MODALDAMPING=5.
EGROUP 1
ISOBEAM MATERIAL=I SECTION 1
SDIM=0.05 TDIM=0.05 GLINE 1 2 AUX=4 EL=A 0 NODES=4 EGROUP 2
ISOBEAM MATERIAL=2 SECTION 1
SDIM=0.005 TDIM=0.05 GLINE 2 3 AUX=4 EL=2 NODES=4 FIXBOUNDARIES
/
1 SET VIEW=X NSYMBOLS=YES NNUMBER=MYNODES BCODE=ALL MESH SOLVIA END SOLVIA-POST input A26 MATERIAL DAMPING IN MODAL SUPERPOSITION, K-WEIGHTED DATABASE CREATE WRITE 'a26.1is' FREQUENCIES SET PLOTORIENTATION=PORTRAIT SET RESPONSETYPE=VIBRATION ORIGINAL=DASHED SMOOTHNESS=YES SET VECTOR=DISPLACEMENT VIEW=X HEIGHT=0.25 MESH TIME=i SUBFRAME=21 MESH TIME=2 MESH TIME=3 SUBFRAME=21 MESH TIME=4 END Version 99.0 Linear Examples A26.3
SOLVIA Verification Manual EXAMPLE A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD Objective To verify the bending and twisting behaviour of the 16-node SHELL element and the PLATE element when modelling a simply supported plate subjected to uniform pressure loading.
Physical Problem A square plate with simply supported edges under uniform pressure loading as shown in the figure below is considered.
Y T
I E = 2.1.10" N/m 2 v=0.3 a=2.0 in h = 0.05 m (thickness) p=l.0-105 N/m 2 (pressure)
Finite Element Model Because of the symmetry conditions only one quarter of the plate is considered, 1-2-3-4 in the figure above. Two finite element models are used, one with 10x10 16-node SHELL elements, see page A27.5, and one with lOx 10x4 PLATE elements, see page A27.1 1. The SHELL element model uses a finer mesh towards the boundaries. Two sets of boundary conditions are used with both models:
"* Soft boundary conditions, examples A27 and A27B: Both bending and twisting degrees of freedom are free along the simply supported boundary.
"* Hard boundary conditions, examples A27A and A27C: The bending degree of freedom is free but the twisting degree of freedom is fixed along the simply supported boundary.
Symmetrical boundary conditions are employed along the symmetry boundaries. All nodes of the models have fixed X-and Y-translations and Z-rotation.
Solution Results The theoretical solution for a simply supported Kirchhoff plate under the considered loading is given in [1], article no. 30. A Kirchhoff plate has no transverse shear deformation but the effect of transverse shear deformation is discussed at the end of article no. 39.
Version 99.0 Linear Examples i t A27.1
SOLVIA Verification Manual The following theoretical solution applies to a square Kirchhoff plate:
Node 3 deflection Node 3 bending moment Node 4 reactive shear force Node 2 reactive shear force Node 1 concentrated reaction w = 0.00406.pa 4 / D M
= Myy = 0.0479. pa 2 F,, = 0.420- pa Zrz = 0.420. pa R = 0.065-pa 2
E h3 The flexural rigidity of the plate is D = 12(1-v2)
Using the input data shown on pages A27.13 to A27.16 the following results are obtained, which are compared with the corresponding theoretical values:
16-node SHELL PLATE Kirchhoff Soft BC Hard BC Soft BC Hard BC theory A27 A27A A27B A27C w at node 3 (mm)
-2.70
-2.7707
-2.7127
-2.7001
-2.7001 M,, at node 3, (Nm/rnm)
-19200
-19473
-19163
-19136
-19136 Fx, at node 4 (N/m) 84000 84181 67607 80741 **
80632 **
R comer reaction
-26000
-22558
- 1114 *
-25379
-24926 max MY (Nm/rm) 12814 12994 12933 12968 max 0 rmises (MPa) 53.287 54.015
- )
Summation of reactions at the nodes that show downward reaction in case A27.
- )
Calculated from reactions. With the used mesh about 4 percent of the pressure goes directly to the degrees of freedom that are fixed in the Z-direction and is not included in the reactions. A finer mesh increases the support reactions.
The distribution of displacements, twisting moment, bending moments, von Mises effective stress and stress deviation, transverse shear forces and force and moment reactions are shown for the 16-node SHELL models on pages A27.5 to A27. 10. The r-s-t axes of the Element Stress System are orientated in the same directions as the Global X-Y-Z System, see figure on page A27.5.
Result plots for the PLATE element model with soft boundary conditions are shown on pages A27. 11 and A27.12.
The following observations can be made regarding the results:
The deflection is slightly larger for the models with 16-node SHELL elements since transverse shear deformation is included.
Version 99.0 (m)
(Nm/m)
(N/m)
(N/m)
(N)
Linear Examples A27.2
SOLVIA Verification Manual The bending moment values are almost the same for the studied models.
The transverse shear force per unit length at node 4 is almost the same for the 16-node SHELL model with soft boundary conditions and the PLATE models when the PLATE values are calculated from reactions. Some small differences occur because of the mesh size of the PLATE models adjacent to the simply supported boundaries. The portion of the pressure load which acts on the boundary nodes is not included in the reactions. Only pressure load in the directions of active degrees of freedom is considered in the load and reaction calculations.
The transverse (reactive) shear force per unit length calculated for the 16-node SHELL model with hard boundary conditions is significantly less than the Kirchhoff theoretical value. The reason is that the portion associated with the twisting moment goes directly to ground without being transformed to statically equivalent transverse shear forces.
For a Kirchhoff plate the undeformed normals to the plate midsurface remain normals during defor mation. The Kirchhoff plate can be regarded to have an orthotropic material with infinite shear moduli Gx, and Gyz. A consequence is that the effect on the plate of an applied distributed twisting moment is the same as the statically equivalent distributed force couples, see figure below with a portion of the 4 - 1 boundary.
dy dy dy Z I aa aM S o.
X +--ady F MY
,M~x/+day x
Y Y
M\\,
,(+aMg 2Mxya
+ aM y dy Ma M
-y dy Twisting moment per unit length Statically equivalent force couples acting and the corresponding shear stress on two dy portions. The net transverse DMx shear force on one dy portion is am 1y dy ay We note that the net transverse shear force per unit length, which is statically equivalent to the twisting moment, is aMy / ay. Following [1 ] we obtain at node 4 the theoretical value DM'y = 0.082 pa = 16400 N/m ay If this portion is transferred directly to the fixed twisting support then the transverse shear force reaction remaining at node 4 is Fx2 = 84000 -16400 = 67600 N/m which is in excellent agreement with the value 67607 N/m computed for the 16-node SHELL model with hard boundary conditions.
Version 99.0 Linear Examples A27.3
SOLVIA Verification Manual The statically equivalent transverse shear force corresponding to the twisting moment is for a Kirchhoff plate balanced at the comer node 1 by a concentrated reaction. In this example half of the reaction is due to the twisting moment along the 4 - 1 boundary and is There is a good agreement between the concentrated reaction calculated by the PLATE element models and the theoretical value.
The 16-node SHELL model with hard boundary conditions gives no concentrated comer reaction, since the twisting moment is not transformed to statically equivalent transverse shear forces. Using soft boundary conditions there is a distributed shear force close to the comer node acting in the downward direction, but no concentrated comer reaction. The sum of the downward shear forces close to the comer is in reasonable agreement with the concentrated value according to the Kirchhoff theory considering that the shear deformation is a relieving factor.
The maximum value of the twisting moment Mxy for the 16-node SHELL and the PLATE models are in close agreement. However, when soft boundary conditions are used for the 16-node SHELL model there is a transition region along the simply supported boundaries extending into the plate of the order one plate thickness. In this transition region the shear stresses cyu corresponding to the twisting moment are transformed to shear stresses in the thickness direction.
The maximum value of the von Mises effective stress is about the same for the two SHELL models with soft and hard boundary conditions. The stress deviation values indicate that the used mesh is good.
User Hints
"* Soft boundary conditions are in general recommended for a SHELL model with simply supported edges.
" Hard boundary conditions for a SHELL model with simply supported edges allow the twisting moment to be directly resisted by a twisting moment reaction instead of being transformed to statically equivalent transverse shear reactions and concentrated comer reactions. This results in lower values for the vertical reaction forces.
"* If only bending moment and deflection results are of interest then both hard and soft boundary conditions can be used along the simply supported boundary.
"* For a PLATE model of a simply supported plate then soft and hard boundary conditions give about the same results if a reasonably fine mesh is used.
"* The 4-node and the 9-node SHELL elements can also be used successfully for a simply supported plate but the 16-node SHELL element used in this example gives a smoother stress distribution.
Reference
[1]
Timoshenko, S.P. and Woinowsky-Krieger, S., Theory of Plates and Shells, Second Edition, McGraw-Hill, 1959.
Version 99.0 Linear Examples A27.4
SOLVIA Verification Manual Linear Examples ORTED SQUARE PLATE UNDER PRESSURE LOAD, 16 NODE SHELL EL Y
Lx B
C D
E F
G ORIGINAL 0.2 MASTER I1000i 1t0o01 110t01 S101iI 1i1001 111011
!11101 SOLVIA-PRE 99.0 t
r S
EAXES=
STRESS-RST SOLVIA ENGINEERING AB A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE
- LOAD, 16 NODE SHELL EL.
MAX DISPL.
-H 2.7707E-3 TIME 1 Y
MAX DISPL LX TIME I
DISPLACEMENT MAX 2.7707E-3 2 $S976E-3 S2.251t2E-3 1 9049E-3 1.5585E-3 1.2122E-3 8.6586E-4 S.951E-4 1.7317E-4 MIN 0 SOLVIA-POST 99.0 SMRS-SECTION MAX 12814 11949 C10218 8487. S 6756.8 S026.0 32952 i-166.28 MIN-1031.7 SOLVIA ENGINEERING AB Version 99.0 A27 SIMPLY SUPPI ORIGINAL 0.2 2.7707E-3 Y
Lx Y X_
A27.5
SOLVIA Verification Manual Linear Examples A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, 16 NODE SHELL EL.
MAX DISPL.
m-i 2.7707E-3 TIME I
Y L-x Y
MAX DISPL.
i-]
2.7707E-3 L
X TIME I
MRR-SECTION MAX 4415.9
- 2922.9
-63.198
-3049.3
-6035.3
-9021.4
-12007
-14993
- -17980 MIN-19473 SOLVIA-POST 99.0 MSS-SECTION MAX 4415.9 S~
2922.9
-3049.3
-6035.3
-9021.4
- - 12007
-N,
_1 4 9 9 3
ý-17980 MIN-19473 SOLVIA ENGINEERING AB A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, 16 NODE SHELL EL.
MAX DISPL. ý--
2.7707E-3 TIME I
Y Lx Y
MAX DISPL. ý-- 2.7707E-3 L
xTIME I MISES SHELL TOP MAX 5.3287E7 IS.0166E7
,,:,4.392IE7 3.7682E7
- 3.1440E7 2.5198E7 1.8956E7 1.2714E7 6.4725E6 MIN 3.3516E6 SOLVIA-POST 99.0 DEVIATION OF MISES SHELL TOP MAX 2.3933E6
!2.2437E6 i.9446E6 1.64S4E6 1.3462E6
- 1.0471E6 7.4791E5 4.4875E5 1.49S8ES MIN 0 SOLVIA ENGINEERING AB Version 99.0 A27.6
SOLVIA Verification Manual Linear Examples Version 99.0 A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, 16 NODE SHELL EL.
MAX DISPL.
m 2.7707E-3 Y
MAX DISPL.
i 2.7707E-3 Y
TIME I LX TIME I LX FRT-SECTIONF CIO MAX 841[81[
MAX 811[81 70000 70000 I-1 S000o 500000 30000
-30000
-1O000 1-S 000
~ -10000
-00 r~r 30000 S
0000 00
-70000
-70000 MIN-8.2302ES MIN-8.2302ES SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, 16 NODE SHELL EL.
MAX DISPL.
t 2.7707E-3 z
MAX DISPL.
H 2.7707E-3 z
TIME I
Y 4,,,X TIME I
YLX REACTION MREACTION 1167.8 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A27.7
SOLVIA Verification Manual MAX D!SP[
TIME I
2.7127E-3 Y
MAX DISPL. i-H 2.7127E L
xTIME I
Linear Examples 16 NODE SHELL EL.,
HARD BC A27A SIMPLY SUPPORTED SQUARE PLATE, Y L-x MRS-SECTION MAX 12994 S12 1 8 2 1 0SS8 8933.3 7309.0 568
ý 7 1060.4 2436.2 S811.85 MIN-0,24981 SOLVIA ENGINEERING AB SOLVTA-POST 99.0 Version 99.0 DISPLACEMENT MAX 2.7127E-3 S2.543iE-3 2.2040E-3 1.86S0E-3 1.5259E-3 1 1868E-3
- 8. 4771E-4 5.0862E-1
- 1.
I 6954E-4 MIN A27A SIMPLY SUPPORTED SCUARE PLATE, 16 NODE SHELL EL.,
HARD BC ORIGINAL 0.2 Y
ORIGINAL 0.2 Y
Lx Lx E E E E F E E E E E
E EEE EEEEE FI N N >
SNN' N:7",N B !i00[1i C 110 10 1 E 111011 r
S F
l!LOI0 EAXES=
G 11,1111I STRESS-RST SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB A27.8
SOLVIA Verification Manual Linear Examples A27A SIMPLY SUPPORTED SOUARE PLATE, 16 NODE SHELL EL.,
HARD BC MAX DISPL. F-2.7127E-3 TIME I Y Ix Y
MAX DISPL. F-2.7127E-3 IX TIME I
MRR-SECTION MAX 0
-1197.7 3593.0
-*-988.4
-8383.8
-10779
-13174
-ISS70 J-17965 MIN-19163 SOLVIA-POST 99.0 MSS-SECTION
-1197. 7
-5988.4
-8383.8
-;*- 1 0779
-13174
-15570
ýj-17965 MIN-19163 SOLVIA ENGINEERING AB A27A SIMPLY SUPPORTED SQUARE PLATE, 16 NODE SHELL EL.,
HARD BC MAX DISPL. ý-i 2.7127E-3 TIME 1 Y
MAX DISPL. ý--
2.7127E-3 LX TIME I MISES SHELL TOP MAX 5.401SE7 I
S.078SE7 S4.4326E7 3.7867E7 3.1408E7 2.4949E7 1.8490E7 1.2031E7 S.571SE6 MIN 2.3420E6 SOLVIA-POST 99.0 DEVIATION OF MISES SHELL TOP MAX 4691.3 4'398.1 3811.
- ! 322S.2 N"2052.4 11466.0 879.61 293.20 MIN 0 SOLVIA ENGINEERING AB Version 99.0 A27.9
SOLVIA Verification Manual Linear Examples A27A SIMPLY SUPPORTED SQUARE PLATE, 16 NODE SHELL EL.,
HARD BC MAX DISPL. v-2.7127E-3 TIME I
z MAX DISPL.
.-- 2.7127E-3 Yd, /X TIME I
REACTION 4021.9 SOLVIA-POST 99.0 MREACTION 1149.2 SOLVIA ENGINEERING AB Version 99.0 z
Y.,IX A27. 10
SOLVIA Verification Manual Linear Examples A27B SIMPLY SUPPORTED SQUARE PL TIME 1
2>
0 V>
C3
/
("4
/
oD r
0.2 0A' 06 08
[0 LINE-32 SOLVIA-POST 99.0 ATE UNDER PRESSURE LOAD, PLATE EL.
I.
i TIM 1
c2 z
f C')
-I
/>
/
U-I
/J I
/..
37"
.0 0.2 0,4 0-6 0.8 1,
L1NE-3z0 SOLVIA ENGINEERING AB Version 99.0 I
I A27B SIMPLY SU'PORTED SOUARE PLATE UNDER PRESSURE LOAD. PLATE EL.
ORIGINAL 1
ORIGINAL 1
L ZONE OUTER
-*4 llO01 4"*
C\\> NJz I0 li,00 132 0
SOLVIA-PRE 99.0
$OLVIA ENGINEERING AB 0
A27.11I
SOLVIA Verification Manual Linear Examples A27B SIMPLY SUPPORTED SQUARE PLATE UNDE TIME1~
{~~
.0 0.2
- 0. 4 0.6 0.8 1.
C-
-J 0x.
03 Ni
/
C, R PRESSURE
- LOAD, PLATE EL.
TIME I 0.0 LINE-21 0.2 0.4 0.6 0.8
- 1.
LINE-34 NODE 3 -
4 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB Version 99.0
(/)
(/)
z H
C)
U)
(I)
C',
C
('1 0
0 0
0
0
0 0
C',
/<
//
/
/
/7
//
/
o' S....
I....
I
/
A27.12
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, 16 NODE SHELL EL.'
DATABASE CREATE MASTER IDOF=110001 COORDINATES 1
- 1. 1. /
2
- 0. 1.
/
3
- 0.
- 0.
/
4
- 1. 0.
MATERIAL 1
ELASTIC E=2.1E11 NU=0.3 EGROUP 1
SHELL RESULT=NSTRESS SET MIDNODES=2 LINE STRAIGHT 1 2 EL=-0 RATIO=5 LINE STRAIGHT 2 3 EL=10 RATIO=5 GSURFACE 1 2 3 4 THICK 1 0.05 LOADS ELEMENT INPUT=SURFACE 1 2 3 4 T
1.0E5 FIXBOUNDARIES 3
INPUT=LINE 1 2 4 1 FIXBOUNDARIES 4
INPUT=LINE
/
1 4 FIXBOUNDARIES 5
INPUT=LINE
/2 3
- HARD BC USED IN A27A:
- FIXBOUNDARIES 34 INPUT=LINE
/
4 1
- FIXBOUNDARIES 35 INPUT=LINE
/
1 2
- FIXBOUNDARIES 4
INPUT=LINE
/
3 4
- FIXBOUNDARIES 5
INPUT=LINE
/
2 3 SET VIEW=Z NNUMBER=MY NSYMBOL=MY SUBFRAME 21 MESH BCODE=ALL MESH EAXES=STRESS SOLVIA END Version 99.0 Linear Examples A27.13
SOLVIA Verification Manual SOLVIA-POST input A27 SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD 16 NODE SHELL EL.
DATABASE CREATE SET VIEW=Z SUBFRAME 21 MESH CONTOUR=DISP SET OUTLINE=YES MESH CONTOUR=MRS SUBFRAME 21 MESH CONTOUR=MRR MESH CONTOUR=MSS SUBFRAME 21 MESH CONTOUR=MISES MESH CONTOUR=SDEVIATION SET OUTLINE=NO CONTOUR VMAX=70000 VMIN=-70000 SUBFRAME 21 MESH CONTOUR=FRT MESH CONTOUR=FST VIEW 1 1 1 SUBFRAME 21 MESH VECTOR=REACTION MESH VECTOR=MREACTION WRITE a27.1is NLIST KIND=REACTIONS DIR=345 END Version 99.0 Linear Examples A27.14
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A27B SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, PLATE EL.'
DATABASE CREATE MASTER IDOF=110001 COORDINATES
- 11.
1./20.1.
/30.0.
/41.0.
MATERIAL 1
ELASTIC E=2.1E11 NU=0.3 EGROUP 1
PLATE STRESSTABLE 1
1 2 3 4 5 6 7 GSURFACE 1 2 3 4 EL1=10 EL2=10 EDATA
/
1 0.05 LOADS ELEMENT INPUT=SURFACE 1 2 3 4 T
1.E5 FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES HARD BC USED
- FIXBOUNDARIES
- FIXBOUNDARIES
- FIXBOUNDARIES
- FIXBOUNDARIES 3
4 5
INPUT=LINE INPUT=LINE INPUT=LINE IN A27C:
34 INPUT=LINE 35 INPUT=LINE 4
INPUT=LINE 5
INPUT=LINE SET VIEW=Z HEIGHT=0.20 SUBFRAME 21 MESH NNUMBER=YES NSYMBOL=YES OUTI ZONE OUTER GLOBAL ZONE OUTER GLOBAL DELETE XMIN=0.1 MESH MESH OUTER ENUMBER=YES GSCALE=OLD SOLVIA END LINE=YES BCODE=ALL XMAX=0.9 YMIN=0.1 YMAX=0.9 SUBFRAME=OLD Version 99.0
/ /
/
/
/
/
/
1 2 3 4 2 3 4
1 3
2
/41 1
2 4
3 Linear Examples A27.15
SOLVIA Verification Manual SOLVIA-POST input A27B SIMPLY SUPPORTED SQUARE PLATE UNDER PRESSURE LOAD, PLATE EL.
DATABASE CREATE SYSTEM 1 CARTESIAN WRITE FILENAME='a27b.lis' EPLINE LINE-34 399 1 4 2 STEP -4 TO 363 1 4 2 EPLINE LINE-32 398 2 4 1 STEP -40 TO 38 2 4 1 EPLINE LINE-21 37 2 4 1 STEP -4 TO 1 2 4 1 NPLINE LINE-34 3 127 STEP -1 TO 119 4 ELINE LINE-32 KIND=M22 OUTPUT=ALL ELINE LINE-34 KIND=M11 OUTPUT=ALL ELINE LINE-21 KIND=M12 OUTPUT=ALL NLINE LINE-34 DIR=3 OUTPUT=ALL SYSTEM=1 SUBFRAME=21 SYSTEM=1 SYSTEM=I SUBFRAME=21 VIEW 1 1 1 SET OUTLINE=YES ORIGINAL=DASHED MESH VECTOR=REACTIONS SUBFRAME=21 MESH VECTOR=MREACTIONS NLIST KIND=REACTIONS DIRECTION=345 END Version 99.0 Linear Examples A27.16
SOLVIA Verification Manual EXAMPLE A28 PLATE UNDER UNIFORM TWISTING Objective To verify the twisting behaviour of the SHELL and PLATE elements.
Physical Problem A square plate of side lengths a, which is supported at three of its corners as shown in the figure below is considered. Consistent twisting moments along the plate sides are used in the SHELL element models. A concentrated load acting in the negative Z-direction is applied at the free comer on the PLATE element model.
I "
tLI~IImL WILL NIrr E=4.10'0 N/m 2 v =0.2 a
I 1 m t =0.1 m (thickness)
Mt = 5000 Nmrn/m i'
a
. [
Finite Element Model The figure on page A28.3 shows the irregular mesh used with the 16-node SHELL element model.
The irregular mesh of the PLATE element model is shown in the right bottom figure of page A28.4.
Solution Results The theoretical solution for this problem is a constant twist and the central deflection is Mt a 2 E
u =
where G =
and 4GK 2(1l+v at3 3
The input data on pages A28.5 and A28.6 gives the following results:
Central deflection, u (mm):
Version 99.0 r
Linear Examples A28.1
SOLVIA Verification Manual Bending and twisting moments (Nm/m) in the whole model referenced to the global X-Y-Z system:
Theory PLATE element M
=Myy =0 M
=Myy = 0 Mxy = -5000 Mxy = -5000 The bending stresses are zero for the SHELL element models.
The torsion stresses ay in MPa at the TOP surface for the SHELL element models are:
Theory 16-node 9-node 4-node SHELL SHELL SHELL
-3.00
-3.00
-3.00
-3.00 User Hints
" This example may be regarded as a patch test problem since the theoretical solution is a constant twist. The SOLVIA solution obtained for the irregular mesh is exactly equal to the theoretical solution and the SHELL and PLATE elements have, therefore, the ability to represent stresses due to a constant twist.
" For a plate without transverse shear deformation (Kirchhoff plate) there is a correspondence between the transverse shear force and the twisting moment so that a constant twisting moment along an edge can be replaced by comer forces, see figure.
Fx=
x The SHELL element is capable of transverse shear deformation. The above correspondence between transverse shear force and twisting moment is then not valid exactly. If the twisting moment is applied as transverse shear force there will be a transition region along each boundary where transverse shear stresses consistent with the transverse shear force are transformed to in plane shear stresses consistent with a twisting moment.
Version 99.0 Linear Examples A28.2
SOLVIA Verification Manual Linear Examples A28 PLATE UNDER UNIFORM TWISTING, 16-NODE SHELL ORIGINAL 02 Y
ORIGNAL v-ý 0.2 Y
L~x TIME I
Lx MASTER 000000 S 000001 C 011001 D t01001 E
I1001I SOLVIA-PRE 99-0 EMOMENT SO00 SOLVIA ENGINEERING AB Version 99.0 A28.3
SOLVIA Verification Manual A28 PLATE UNDER UNIFORM TWISTING.
16-NODE SHELL ORIGINAL
-- 0.2 MAX DISPL. --
I.8E-3 TIME I
z S
DISPLACEMENT t.8E-3 DISPLACEMENT MAX 1.8000E-3 1.6875E-3 Is 14620E-3 "1-2375E-3 1.012SE-3 S7.8750E-4 5.6250E-4 3.37S5E-4
'1.1250E-4 MIN 0 ORIGINAL
-- 0.2 MAX DISPL.
i-1.8E-3 TIME I
SOLVIA-POST 99.0 2
MLOAD 1312.5 STRESS-Xy 54ELL TOP MAO-3010000.S
\\.x
-2999317 3
.. -29999S5.4
-2999983.5
'-3000016.5
-3000049.6
'-3000082.7
-3000111.8 MIN-3000100.0 SOLVIA ENGINEERING AB Linear Examples Version 99.0 A28A PLATE UNDER UNIFORM TWISTING.
9-NODE SHELL ORIGINAL
ý-
--1 0.2 Z
MAX DISPL.
1i SE-3 Y
TIME I
X DISPLACEMENT 7DISPLACEMENT ASI. 811E-3 MAX 1.8000E-3
.6875E-3 I. 462SE-3 1~*[.2375E-3 S1.Ot25E-3 S7,8750E-4 5.6250E-4 3.3750E-4 1 12SOE-4 MIN 0 ORIGINAL
-- 02 Z
MAX DISPL.
3 ISE3Y TIME I
x MLIAD 2333.3 7
STRESS-XY SHELL TOP MAX-300000.0 2999804.2 2999917.3
.A -2999950.4
-2999983. S
-3000016. S
- -3000049.6 S-3000082.7
-3000115.8 MIN-3000000
.0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A280 PLATE UNDER UNIFORM TWISTING.
4-NODE SHELL ORIGINAL F-0.2 Z
MAX DISPL.
-- 1.8E-3 Y
TIME I x
DISPLACEMENT 7
1.8E-3
~.LDISPLACEMENT MAX 1.8000E-3 1237SE-3 1.0125E-3 7.8710E-4 5.6250E-4 3.3750E-4 M.12SOE-4 MIN 0 ORIGINAL 0-2 Z
MAX DISPL.4-I 1.8E-3 Y
TIME I X
MLOAD 2500 SHELL TOP MAX-IGIOOII.0 2
2999884.2
,-2999917.3
"-2999950.4
-2999883.5
ý;*3000O16.S
-3000049.6
-3000082.7 L-300011S.8 MIN-3000000.0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AR A28C PLATE UNDER UNIFORM TWISTING, PLATE ELEMENT ORIGINAL F-
-0.2 Z
MAX DISPL.
I.8E-3 Y
TIME I x
LOAD 10000 MASTER 000000 B 011000 C 101000 D
111000 ORIGINAL 0-
-- 1.2 Z
MAX DISPL.--4 1.8E-3 Y
TIME I x
DISPLACEMENT 1.8E-3 DISPLACEMENT MAX 5.8000E-3 1.687SE-3 I.462SE-3 1.2375E-3 I10125E-3 7.8750E-4 5.6250E-4 3.3750E-4 I 1250E-4 MIN 0 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB A28.4
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A28 PLATE UNDER UNIFORM TWISTING, DATABASE CREATE 16-NODE SHELL' COORDINATES 1
- 0.
- 0.
6 0.4 1.
/ /
2 7
.3
.0
/
.0
- 1.
/
3
- 1.
.0
/
8
- 0.
.5
/
4 9
- 1.
0.3
.5
.5
/
5
- 1. 1.
MATERIAL 1 ELASTIC E=4.E10 NU=0.2 SET NODES=16
"* 9-node SHELL:
SET NODES=9
"* 4-node SHELL:
SET NODES=4 EGROUP 1 SHELL RESULT=NSTRESS GSURFACE 1 2 9 8 / GSURFACE 2 3 4 9 GSURFACE 4 5 6 9 / GSURFACE 8 9 6 7 THICKNESS 1 0.1 LOAD ELEMENT TYPE=MOMENT INPUT=LINES 1 3 out -5.E3 /
3 5 out 5.E3 /
5 7 FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES FIXBOUNDARIES out -5.E3 /
7 1 out 5.E3 3
/
137 1
/
17 2
/
13 6 LINES /
1 3 /
3 5 /
5 7 /
7 1 SET VIEW=Z NSYMBOL=YES MESH NNUMBER=MYNODES ENUMBER=YES BCODE=ALL MESH VECTOR=LOAD X=0.5 SOLVIA END SOLVIA-POST input A28 PLATE UNDER TWISTING, 16-NODE SHELL DATABASE CREATE WRITE FILENAME='a28.1is' VIEW ID=1 XVIEW=1.
YVIEW=-0.5 ZVIEW=0.5 SET NSYMBOL=MYNODES NNUMBERS=MYNODES ORIGINAL=DASHED SET PLOTORIENTATION=PORTRAIT VIEW=1 SUBFRAME 12 MESH CONTOUR=DISPLACEMENT VECTOR=DISPLACEMENT MESH CONTOUR=SXY VECTOR=MLOAD SHELLSURFACE LIST=TOP NLIST MYNODES ELIST EMAX END Version 99.0 SUBFRAME=21 Linear Examples A28.5
SOLVIA Verification Manual SOLVIA-PRE input HEADING
'A28C PLATE UNDER UNIFORM TWISTING, PLATE ELEMENT' DATABASE CREATE COORDINATES 1
- 0.
- 0.
/
2
.3
.0
/
3
- 1..0
/
4
- 1.
0.3
/
5 1.1.
6 0.4 1.
7
.0
- 1.
/8 0..5
/
9
.5
.5 MATERIAL 1 ELASTIC E=4.E10 NU=0.2 EGROUP 1 PLATE ENODES 1129/2239/3349/4459 5569/6679/8789/9819 EDATA /
1 0.1 LOAD CONCENTRATED 5
3 1E4 FIXBOUNDARIES 3
1 3
7 FIXBOUNDARIES 1
/
1 7 FIXBOUNDARIES 2
/
1 3 SET VIEW=Z NSYMBOL=YES PLOTORIENTATION=PORTRAIT MESH NNUMBER=MYNODES ENUMBER=YES BCODE=ALL SUBFRAME=12 MESH VECTOR=LOAD GSCALE=8 SOLVIA END SOLVIA-POST input A28 PLATE UNDER TWISTING, 16-NODE SHELL DATABASE CREATE WRITE FILENAME='a28.1is' VIEW ID=1 XVIEW=1.
YVIEW=-0.5 ZVIEW=0.5 SET NSYMBOL=MYNODES NNUMBERS=MYNODES ORIGINAL=DASHED SET PLOTORIENTATION=PORTRAIT VIEW=1 SUBFRAME 12 MESH CONTOUR=DISPLACEMENT VECTOR=DISPLACEMENT MESH CONTOUR=SXY VECTOR=MLOAD SHELLSURFACE LIST=TOP NLIST MYNODES ELIST EMAX END Version 99.0 Linear Examples A28.6
SOLVIA Verification Manual EXAMPLE A29 EDGE BENDING AND TWISTING OF A TRIANGULAR PLATE ON CORNER SUPPORTS Objective To verify the behaviour of the SHELL element subjected to edge moment loading.
Physical Problem A triangular plate on three comer point supports as shown in the figure below is considered. The displacement in the Z-direction is constrained to zero at each comer. The plate is loaded along its boundaries with a constant bending moment and a linearly varying twisting moment.
Y 2c 3
6 M,
b/3 E = 207.10 9 N/m 2 v = 0.25 Mb = 300 Nm/m Mt = 194.85 Nm/m t=2.54.10-3 m b=0.24 m c = 0.138564 m Finite Element Model Three 16-node SHELL elements are used to model the plate. The element mesh with the boundary condition and the edge moments is shown in the figure on page A29.3.
Solution Results Theoretical solutions of this problem are presented in references [1] and [2]. Different FEM-solutions can be found in references [3]. The SOLVIA solution presented here has also been presented in [4].
Version 99.0 Linear Examples 1ý M,
A29.1
Linear Examples SOLVIA Verification Manual Centroidal displacement (node 7)
The internal element moment can be formed from the TOP surface stresses:
t 2 Mx= -a-x t2 6y
= -
'y Calculated moments from SOLVIA analysis (Nm/m)
Node 6 Node 7 Node 2 Node I Node 3 Mxx 300.000 187.503
-37.4903 300.000 300.000 Myy 75.0066 187.503 412.500 75.0066 75.0066 Mxy_
FI__ 1.10.
3.10-'
9.10-11
-194.851 194.851 Theoretical moments from reference [1] (NmI/m)
Node 6 Node 7 Node 2 Node 1 Node 3 Mxx 300.000 187.500
-37.5000 300.000 300.000 Myy 75.0000 187.500 412.500 75.0000 75.0000 Mxy _
F 0
0 0
-194.856 194.856 An excellent agreement can be observed.
The stress distributions are shown on page A29.4.
User Hints The loading consists of two parts, namely a constant bending moment and a linearly varying twisting moment. The constant bending moment load produces constant curvatures in the plate and hence a deformed shape that is spherical. The stress field consists only of bending stresses and the bending stress components are equal and constant over the faces of the plate. It is a basic require ment for any plate element that constant bending stresses can be modeled exactly.
The difficulty in the triangular plate problem lies instead in obtaining a good solution to the load ing by the linearly varying twisting moment. The magnitude of this applied twisting moment is such that the edges of the plate, which are deformed by the constant bending moment to lie on a Version 99.0 M "y =-* t xy A29.2
SOLVIA Verification Manual spherical surface, are deformed back to their original straight line shape. This magnitude of the twisting moment makes it easy to check the resulting edge displacements since these should be zero.
The two bending stress components and the in-plane shear stress component are varying linearly due to the applied twisting moment. The shear stress components in the thickness direction are zero. The case of the linearly varying twisting moment can be regarded as a patch test, which illuminates the ability of a plate element to model a linearly varying state of stress.
References
[1]
Robinson, J., "A FEN Project - Triangular Plate-Bending Continuum on Three Point Supports",
Finite Element News, 1992 Issue no. 1.
[2]
Robinson, J., "An Elasticity Solution for a Triangular Plate-Bending Continuum on Three Point Supports", Finite Element News, 1995 Issue no. 2.
[3]
Finite Element News, 1992 Issues no. 2, 4, 5, 6 and 1993 Issue no. I.
[4]
Larsson, G., "A Note on the FEN Project Triangular Plate-Bending on Continuum Supports",
Finite Element News, 1994 Issue No. 3.
Version 99.0 A29 EDGE BENDING AND TW4ISTING OF A TRIANGULAR PLATE ON CORNER SUPPORTS ORIGINAL 0 O.OG Y
ORIGINAL i
O.OS Y
TIME 1
Lx L
x C,3 42 E;
1ASTER B 000001 EMOMENT C C01001 D ilOIC0 3I E EE ilNOO SOLVIA-PRE 99.0 SOLVIA ENGINEERING AB Linear Examples A29.3
SOLVIA Verification Manual Linear Examples A29 EDGE BENDING AND TWISTING OF A TRIANGULAR PLATE ON CORNER SUPPORTS MAX DISPL.
-i 2.1228E-3 TIME i
MAX DISPL.
-H 2-1228E-3 TIME I STRESS-XX SHELL TOP MAX 2.7900E8 M 2.5938E8 i*2.201SEB 1.8092E8 1#i
.4168E8 S.0245E8 S.3217E7 2.3984E7
-. 249E7 MIN-3.4866E7 STRESS-XY SHELL TOP MAX 1.8121E8 S1.58S6E8 1.1326E8 i~*6.79S4E7
- 2.26SIE7 i-2.26SIE7
-6'.7954E7 N-1.I326E8
ýj-i.5856E8 MIN-t.812IE8 SOLVIA-POST 99.0 MAX DISPL.
H 2.1228E-3 TIME 1
STRESS-YY SHELL TOP MAX 3.8362E8 S3.6401E8 3.2477E8 2.8SS4E8
- 2. 4631E8 2.0707E8 1.6784E8 1,286IE8 j 8.9373E7 MIN 6.9756E7 MAX DISPL.
1 2.1228E-3 TIME I
Z-DIR DISPL.
MAX 2.1228E-3 I
1.9901E-3 lA 1.7248E-3
- 1. 4594E-3 K::,,I1.1941E-3
- 9. 2873E-4 S6.6338E-4 3.9803E-4 1.3268E-4 MIN 0 SOLVIA ENGINEERING AB A29 EDGE BENDING AND TWISTING OF A TRIANGULAR PLATE ON CORNER SUPPORTS MAX DISPL. i-2.1228E-3 TIME I MAX DISPL.
i-i 2.1228E-3 PRINCIPAL TIME I
STRESS MAX.
SHELL TOP MAX 3.8362E8 P3.7054E8 3.4439E8 3.1823E8 2.9208E8 2.6592E8 N2.3977E8 2.1361E8 1.8746E8 MIN 1.7438E8 PRINCIPAL STRESS MID.
SHELL TOP MAX I.7438E8
!I. 6348E8 I.4168EB 1.1989E8
- 9. 8088E7 7.6291E7 S.4493E7
- 3. 2696E7 I. 0899E7 MIN 0 MAX DISPL.
.-. 2.1228E-3 TIME I PRINCIPAL ORIGINAL 0.05 MAX DISPL.
-i 2.1228E-3 SHELL TOP TM MAX 3.5622E-7
-2.179IE6
-6.5374E6 1.0896E7
-. 5254E7
-1.9612E7
-2.3970E7 I
28329E7
.2687E7 MIN-3.4866E7 SOLVIA-POST 99.0 SOLVIA ENGINEERING AB Version 99.0 z
X 1,Y REACTION I.5643E-9 A29.4
SOLVIA Verification Manual SOLVIA-PRE input DATABASE CREATE HEADING 'A29 EDGE BENDING AND TWISTING OF A TRIANGULAR PLATE ON CORNER SUPPORTS'.
COORDINATES 1
- 0.
-0.138564 2
0.24 0
3 0
0.138564 4
0.12 -0.069282 5
0.12 0.069282 6
0 0
7 0.08 0
MATERIAL 1 ELASTIC E=207E9 NU=0.25 EGROUP 1 SHELL RESULT=NSTRESS THICKNESS 1 0.00254 GSURFACE 1 4 7 6 GSURFACE 2 5 7 4 GSURFACE 3 6 7 5 LOADS ELEMENT MOMENTINTENSITY INPUT=LINES 1 2 edge -300
-300 2 3 edge -300
-300 3
1 edge -300
-300 1 2 out
-194.85 194.85 2 3 out
-194.85 194.85 3 1 out
-194.85 194.85 FIXBOUNDARIES 123 1
FIXBOUNDARIES 13
/ 2 FIXBOUNDARIES 3
/ 3 FIXBOUNDARIES 6 INPUT=LINE /
2 3 /
3 1 /
1 2 SET NSYMBOL=YES VIEW=Z MESH VECTOR=LOAD SUBFRAME=21 MESH BCODE=ALL NNUMBER=MYNODES ENUMBERS=YES SOLVIA END Version 99.0 Linear Examples A29.5
SOLVIA Verification Manual SOLVIA-POST input A29 EDGE BENDING AND TWISTING OF A TRIANGULAR PLATE ON CORNER SUPPORTS DATABASE CREATE WRITE 'a29.1is' NLIST MYNODES ELIST ELIST SELECT=SEFFECTIVE SET VIEW=Z AXES=NO HEIGHT=0.28 SUBFRAME 22 MESH CONTOUR=SXX MESH CONTOUR=SYY MESH CONTOUR=SXY MESH CONTOUR=DZISPLACEMENT SUBFRAME 22 MESH CONTOUR=SPMAX MESH CONTOUR=SPMID MESH CONTOUR=SPMIN COLOR LINE SET VIEW=I ORIGINAL=YES AXES=YES MESH VECTOR=REACTION END Version 99.0 Linear Examples A29.6