MSC Nastran 2021.4 Verification Guide
MSC Nastran 2021.4 Verification Guide
tion Guide
Japan Asia-Pacific
KANDA SQUARE 16F 100 Beach Road
2-2-1 Kanda Nishikicho, Chiyoda-ku #16-05 Shaw Tower
1-Chome, Shinjuku-Ku Singapore 189702
Tokyo 101-0054, Japan Telephone: 65-6272-0082
Telephone: (81)(3) 6275 0870 Email: APAC.Contact@mscsoftware.com
Email: MSCJ.Market@mscsoftware.com
Worldwide Web
www.mscsoftware.com, www.hexagon.com
Support
https://simcompanion.hexagon.com/customers/s/article/support-contact-information-kb8019304
Disclaimer
Hexagon reserves the right to make changes in specifications and other information contained in this document without prior notice.
The concepts, methods, and examples presented in this text are for illustrative and educational purposes only, and are not intended
to be exhaustive or to apply to any particular engineering problem or design. Hexagon assumes no liability or responsibility to any
person or company for direct or indirect damages resulting from the use of any information contained herein.
User Documentation: Copyright 2021 Hexagon AB and/or its subsidiaries. Printed in U.S.A. All Rights Reserved.
This notice shall be marked on any reproduction of this documentation, in whole or in part. Any reproduction or distribution of this
document, in whole or in part, without the prior written consent of Hexagon is prohibited.
This software may contain certain third-party software that is protected by copyright and licensed from Hexagon suppliers. Additional
terms and conditions and/or notices may apply for certain third party software. Such additional third party software terms and
conditions and/or notices may be set forth in documentation and/or at http://www.mscsoftware.com/thirdpartysoftware (or successor
website designated by Hexagon from time to time).
PCGLSS 8.0, Copyright © 1992-2016, Computational Applications and System Integration Inc. All rights reserved. PCGLSS 8.0 is
licensed from Computational Applications and System Integration Inc.
MSC, Dytran, Marc, MSC Nastran, Patran, Hexagon logo, MSC Software logo, e-Xstream, Digimat, and Simulating Reality are
trademarks or registered trademarks of Hexagon AB and/or its subsidiaries in the United States and/or other countries.
NASTRAN is a registered trademark of NASA. FLEXlm and FlexNet Publisher are trademarks or registered trademarks of Flexera
Software. All other trademarks are the property of their respective owners.
Use, duplicate, or disclosure by the U.S. Government is subjected to restrictions as set forth in FAR 12.212 (Commercial Computer
Software) and DFARS 227.7202 (Commercial Computer Software and Commercial Computer Software Documentation), as
applicable.
U.S. Patent 9,361,413
December 3, 2021
NA:V2021.4:Z:Z:Z:DC-VER-PDF
Documentation Feedback
At MSC Software, we strive to produce the highest quality documentation and welcome your feedback.
If you have comments or suggestions about our documentation, write to us at: documentation-
feedback@mscsoftware.com.
Please include the following information with your feedback:
Document name
Release/Version number
Chapter/Section name
Topic title (for Online Help)
Brief description of the content (for example, incomplete/incorrect information, grammatical
errors, information that requires clarification or more details and so on).
Your suggestions for correcting/improving documentation
You may also provide your feedback about MSC Software documentation by taking a short 5-minute
survey at: http://msc-documentation.questionpro.com.
Note: The above mentioned e-mail address is only for providing documentation specific
feedback. If you have any technical problems, issues, or queries, please contact Technical
Support.
Contents
MSC Nastran Verification Guide
Contents
Preface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
Main Index
2 MSC Nastran Verification Guide
Main Index
Contents 3
Main Index
4 MSC Nastran Verification Guide
Heat Transfer
Main Index
Contents 5
6.41. Test 21R: Simply Supported Thick Square Plate: Random Forced Vibration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .358
Nonlinear
6.54. Test 1A: 2-D Plane Stress - Uniaxial Load, Secondary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .401
6.55. Test 1B: 2-D Plane Stress - Uniaxial Displacement, Secondary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .405
6.56. Test 2A: 2-D Plane Stress - Biaxial Load, Secondary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .409
6.57. Test 2B: 2-D Plane Stress - Biaxial Displacement, Secondary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .412
6.58. Test 3B: 2-D Plane Stress - Biaxial (-ve) Displacement, Secondary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .416
6.59. Test 8A: 2D Plane Stress - Uniaxial Load, Primary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .420
6.60. Test 8B: 2-D Plane Stress - Uniaxial Displacement, Primary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .423
6.61. Test 8C: 2D Plane Stress - Stepped Load Primary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .426
6.62. Test 9A: 2-D Plane Stress - Biaxial Load, Primary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .429
6.63. Test 9B: 2-D Plane Stress - Biaxial Displacement, Primary Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .432
6.64. Test 9C: 2-D Plane Stress - Biaxial Stepped Load, Primary Creep. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .435
Composite
Main Index
6 MSC Nastran Verification Guide
Contact
6.73. 3DNLG1: Elastic Large Deflection Response of A Z-Shaped Cantilever under End Load . . . . . . . . . . . . . . . . . . . . . 465
6.74. 3DNLG2: Elastic Large Deflection Response of a Pear-Shaped Cylinder under End Shortening. . . . . . . . . . . . . . . . 468
6.75. 3DNLG6: Buckling of a Flat Plate When Subjected To In-Plane Shear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 471
6.76. 3DNLG7: Elastic Large Deflection Response of a Hinged Spherical Shell Under Pressure Loading . . . . . . . . . . . . . 474
6.77. 3DNLG9: Large Elastic Deflection of a Pinched Hemispherical Shell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 477
6.78. 3DNLG10: Elastic-Plastic Behavior of a Stiffened Cylindrical Panel under Compressive End Load. . . . . . . . . . . . . . 480
Main Index
Preface
About this Guide 14
List of MSC Nastran Guides 15
Technical Support 16
Training and Internet Resources 16
Main Index
14 MSC Nastran Verification Guide
Main Index
MSC Nastran Verification Guide 15
Preface
Main Index
16 MSC Nastran Verification Guide
You may find any of these documents from MSC Software at:
https://simcompanion.hexagon.com
Technical Support
For technical support phone numbers and contact information, please visit:
https://simcompanion.hexagon.com/customers/s/article/support-contact-information-kb8019304
Support Center (http://simcompanion.hexagon.com)
The SimCompanion link above gives you access to the wealth of resources for MSC Software products. Here
you will find product and support contact information, product documentations, knowledge base articles,
product error list, knowledge base articles and SimAcademy Webinars. It is a searchable database which allows
you to find articles relevant to your inquiry. Valid MSC customer entitlement and login is required to access
the database and documents. It is a single sign-on that gives you access to product documentation for
complete list of products from MSC Software, allows you to manage your support cases, and participate in
our discussion forums.
MSC Software corporate site with information on the latest events, products and services for the
CAD/CAE/CAM marketplace.
http://simcompanion.hexagon.com
The SimCompanion link above gives you access to the wealth of resources for MSC Software products. Here
you will find product and support contact information, product documentations, knowledge base articles,
product error list, knowledge base articles and SimAcademy Webinars. It is a searchable database which allows
you to find articles relevant to your inquiry. Valid MSC customer entitlement and login is required to access
the database and documents. It is a single sign-on that gives you access to product documentation for
complete list of products from MSC Software, allows you to manage your support cases, and participate in
our discussion forums.
http://www.mscsoftware.com/msc-training
The MSC Training link above will point you to schedule and description of MSC Seminars. Following
courses are recommended for beginning MSC Nastran users.
NAS120 - Linear Static Analysis using MSC Nastran and Patran
This seminar introduces basic finite element analysis techniques for linear static, normal modes, and buckling
analysis of structures using MSC Nastran and Patran. MSC Nastran data structure, the element library,
modeling practices, model validation, and guidelines for efficient solutions are discussed and illustrated with
examples and workshops. Patran will be an integral part of the examples and workshops and will be used to
generate and verify illustrative MSC Nastran models, manage analysis submission requests, and visualize
Main Index
MSC Nastran Verification Guide 17
Preface
results. This seminar provides the foundation required for intermediate and advanced MSC Nastran
applications.
Main Index
18 MSC Nastran Verification Guide
Main Index
Section 1: MacNeal-Harder Problems
Main Index
Main Index
1.1 Straight Cantilever Beam
Problem Description
This problem demonstrates the use of different elements to model a cantilevered beam with a rectangular
cross section subjected to unit forces at free end. The solution is compared with the analytical solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is linear elastic with a Young's modulus, E = 1.0 × 107 and Poisson's ratio, = 0.30.
Geometrical properties
Length = 6.0
Width = 0.2
Depth = 0.1
Loading
A force of one unit applied at the free end.
Main Index
10 MSC Nastran Verification Guide
Boundary Condition
The displacements ux = uy= uz = 0.0 at one end of the beam as shown in Figure 1.1-1.
Figure 1.1-1 Straight cantilever beam. (a) Regular shape elements (b) Trapezoidal shape elements
(c) Parallelogram shape elements
Main Index
CHAPTER 1.1 11
Straight Cantilever Beam
Results
Rectangular Elements
Trapezoidal Elements
Main Index
12 MSC Nastran Verification Guide
Parallelogram Elements
Conclusion
From the above observations, the rectangular element results from MSC Nastran are in better agreement with
analytical results than trapezoidal and parallelogram elements.
Main Index
CHAPTER 1.2 13
Curved Beam
Problem Description
This problem demonstrates the use of different elements to model a 90-degree section of a circular beam.
The beam is analyzed with two types of loading conditions:
a. loaded in its plane in a radial direction on its free end, and
b. loaded in out of plane direction on its free end.
The solution is compared to the analytical solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is linear elastic with a Young's modulus, E = 1.0 × 107 and Poisson's ratio = 0.25.
Geometrical properties
Inner radius = 4.12
Outer radius=4.32
Thickness=0.1
Main Index
14 MSC Nastran Verification Guide
Loading
A force of one unit applied at the tip shown in Figure 1.2-1.
Boundary Condition
The displacements ux = uy= uz = 0.0 at one end of the beam as shown in Figure 1.2-1.
Main Index
CHAPTER 1.2 15
Curved Beam
Results
Conclusion
The results from MSC Nastran are in a good agreement with analytical results.
Main Index
16 MSC Nastran Verification Guide
Problem Description
This problem contains basic patch test cases for Nastran elements. The solution is compared to the analytical
solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is elastic with a Young's modulus, E = 1.0 × 107 and Poisson's ratio = 0.25.
Geometrical properties
Membrane Plate
thickness=0.001
length=0.24
width=0.12
Solid Cube
A cube of unit length is used as shown in Figure 1.3-1.
Boundary Condition
Displacement boundary conditions at all exterior nodes for constant stress loading are given below:
Main Index
CHAPTER 1.3 17
Patch Test
Membrane Plate
u = 10-3 (x+y/2)
v = 10-3 (y+x/2)
Solid Cube
u = 10-3 (2x+y+z)/2
v = 10-3 (x+2y+z)/2
w = 10-3 (x+y+2z)/2
a. plates b. solids
Main Index
18 MSC Nastran Verification Guide
Results
Theoretical Results
Conclusion
All elements show exact solutions except for the QUAD8 which differ from the analytical solution by about
15%.
Main Index
CHAPTER 1.4 19
Rectangular Plate
Problem Description
This problem demonstrates the use of different elements to model a rectangular plate. The plate is analyzed
with two types of loading conditions:
a. uniform pressure, and
b. concentrated central load.
The solution is compared to the analytical solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is elastic with a Young's modulus, E = 1.7472 × 107 and Poisson's ratio = 0.3.
Geometrical properties
Two types of plate geometry are used for the analysis.
1. Square plate: a=2.0
2. Rectangular plate: a=2.0, b=10.0
3. Thickness= 0.0001(Shell element) and 0.01 (Solid element).
Main Index
20 MSC Nastran Verification Guide
Loading
Two types of loadings are used for the analysis
1. Uniform pressure= 10-4
2. Central load = 4.0 × 10-4
Boundary Condition
The model is analyzed using two types of boundary conditions
1. Simply supported
2. Clamped
Main Index
CHAPTER 1.4 21
Rectangular Plate
(b)
Figure 1.4-1 (a) Plate with Aspect Ratio 1.0 (b) Plate with Aspect Ratio 5.0
Main Index
22 MSC Nastran Verification Guide
Main Index
CHAPTER 1.4 23
Rectangular Plate
Main Index
24 MSC Nastran Verification Guide
Results
1. Simply Supported and Uniform Load
Main Index
CHAPTER 1.4 25
Rectangular Plate
Main Index
26 MSC Nastran Verification Guide
Conclusion
The results from MSC Nastran are in a good agreement with analytical results.
Main Index
CHAPTER 1.5 27
Scordelis-Lo-Roof
1.5 Scordelis-Lo-Roof
Problem Description
This problem demonstrates the deformation of thin curved membrane loaded under its own weight. The
computed maximum Z displacement is compared to the theoretical solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is linear elastic with a Young's modulus, E = 1.0 × 107 and Poisson's ratio, = 0.25.
Geometrical properties
Radius = 25.0
length=50.
Thickness=0.25
Angle = 40°
Loading
A force of 90.0 per unit area in the - Z direction is applied on the surface.
Main Index
28 MSC Nastran Verification Guide
Boundary Condition
The displacements boundary conditions are applied as shown in Figure 1.5-1.
Main Index
CHAPTER 1.5 29
Scordelis-Lo-Roof
Main Index
30 MSC Nastran Verification Guide
Results
Conclusion
Except for PENTA elements, all other element results from MSC Nastran are in good agreement with the
theoretical solution.
Main Index
CHAPTER 1.6 31
Spherical Shell
Problem Description
This problem demonstrates the use of different elements to model a thin hemispherical shell subjected to
concentrated loads. The solution is compared to the theoretical solution [1].
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is linear elastic with a Young's modulus, E = 6.825 × 107 and a Poisson's ratio, = 0.3.
Geometrical properties
Radius = 10.0
Thickness = 0.04
Loading
1.0-unit load is applied at quadrant.
Boundary Condition
The displacements boundary conditions are applied as shown in Figure 1.6-1.
Main Index
32 MSC Nastran Verification Guide
Main Index
CHAPTER 1.6 33
Spherical Shell
Main Index
34 MSC Nastran Verification Guide
Results
Conclusion
Except for PENTA and CTRIAR elements, the results from MSC Nastran are in a good agreement with
theoretical results.
Main Index
CHAPTER 1.7 35
Thick Wall Cylinder
Problem Description
This problem demonstrates the use of different elements to model a thick-walled cylinder subjected to
internal pressure. Nearly incompressible material is used. The solution is compared to the analytical solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is linear elastic with a Young's modulus, E = 1.0 × 107 and Poisson's ratio = 0.49, 0.499 and
0.4999.
Geometrical properties
Inner radius = 3.0
Outer radius=9.0
Thickness=1
Loading
A unit pressure is applied at inner radius as shown in Figure 1.7-1.
Main Index
36 MSC Nastran Verification Guide
Boundary Condition
uy= uz=Rx=Ry= Rz= 0.0 for all the nodes.
Results
Main Index
CHAPTER 1.7 37
Thick Wall Cylinder
Conclusion
Except CTRIAR elements, the results from MSC Nastran are in a good agreement with analytical results.
Main Index
38 MSC Nastran Verification Guide
Problem Description
This problem demonstrates the use of different elements to model a twisted beam with a rectangular cross
section subjected to unit forces at free end. The solution is compared with the analytical solution.
Solution Number
SOL 101
Features Used
Linear Statics.
Reference
1. Macneal., R.H., Harder, R.L., 1985. A proposed standard set of problems to test finite element
accuracy. Finite Elem. Anal. Des. 1, 3-20.
Material Properties
The material is linear elastic with a Young's modulus, E = 29.0 × 106 and Poisson's ratio = 0.3.
Geometrical properties
Length = 12.0
Width = 1.1
Depth = 0.32
Twist= 90o (root to tip)
Loading
A force of one unit applied at the free end.
Main Index
CHAPTER 1.8 39
Twisted Beam
Boundary Condition
The displacements ux= uy=uz = 0.0 at one end of the beam as shown in Figure 1.8-1.
Results
Conclusion
The results from MSC Nastran are in a good agreement with analytical results.
Main Index
40 MSC Nastran Verification Guide
Main Index
Section 2: Original Verification
Manual Problems
Main Index
Main Index
2.1 Static Analysis of a Pin-
Jointed Truss
Problem Description
A truss consisting of two pin-jointed axial elements is loaded vertically. Find the vertical deflection of the end
and the internal forces and stresses in the elements.
Solution Number
SOL 101
Features Used
Static Analysis
Reference
1. Egor P. Popov, Introduction to Mechanics of Solids (Englewood Cliffs: Prentice-Hall, Inc. 1968), P.
499
Main Index
42 MSC Nastran Verification Guide
Reference Deck
Model Name from Old Manual Model Description Location
v10101.dat v2401.dat Static analysis of pin jointed truss. tpl\verifman\oldverif\
Results
Theoretical Solution
At point B: U y = – 0.0444 in
Force in Element 1: P 1 = +2500. lb (tension)
Force in Element 2: P 2 = – 2500. lb (compression)
Stress in Element 1: 1 = P 1 A 1 = +16667. lb/in2
Stress in Element 2: 2 = P 2 A 2 = – 10000. lb/in2
Main Index
CHAPTER 2.1 43
Static Analysis of a Pin-Jointed Truss
At point B: U y = – 0.0444 in
Force in Rod 1: +2500. lb (tension)
Force in Element 2: -2500. lb (compression)
Stress in Element 1: +16667. lb/in2
Stress in Element 2: -10000. lb/in2
Comparison of Results
Solution Summary
Element Stress (lb/in2)
Displacement (in) Element Force (lb) Rod 1 Rod 2
Theory -0.0444 2500. +16667. -10000.
MSC Nastran -0.0444 2500. +16667. -10000.
% Difference 0% 0% 0% 0%
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
44 MSC Nastran Verification Guide
Problem Description
A fin of circular cross section is maintained at 250o at one end. The fin is extended into the air. The air
temperature is 70o. Assuming a material conductivity of 132 Btu/hr-ft-oF. Determine the steady-state
temperature distribution along the rod.
Solution Number
SOL 400
Features Used
Steady heat transfer
Reference
1. Alan J. Chapman, Heat transfer, 3rd ed. (New York: MacMillan Publishing Co., Inc., 1974) P. 76,
Ex. 3.8
Main Index
CHAPTER 2.2 45
Heat Transfer from Cooling Fin
Reference Deck
Model Name from Old Manual Model Description Location
v40001.dat v2402.dat Heat transfer from cooling fin tpl\verifman\oldverif\
Results
Theoretical Solution
cosh mL – x
T = T B – T F ---------------------------------- + T F
cosh mL
where m = 4h
------ (for cylindrical rod)
kd
m = 4 1.6 -
--------------------------------
1
0.63 ---------------
2 12
1
m = 1.079 ----
ft
X/L T (oF)
0.25 217.0243
0.50 196.0559
0.75 183.7156
1.00 179.7002
Main Index
46 MSC Nastran Verification Guide
Result comparison of theoretical solution and MSC Nastran results are as follows:
T (oF)
X/L 0.250 0.500 0.750 1.000
Theory 217.62 196.06 183.72 179.7
MSC Nastran 217.8 196.338 184.054 180.056
Difference 0.083% 0.142% 0.182% 0.198%
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
CHAPTER 2.3 47
Truss Analysis using Thermal Loads
Problem Description
A pin joined truss is loaded with a force at one end of the component elements is heated. Considering thermal
effects, find the displacements of the truss joints and the forces and stresses in the axial elements.
Solution Number
SOL 101
Features Used
Static analysis with thermal loads
Reference
1. J. S. Przemieniecki, Theory of Matrix Structural Analysis (New York: McGraw- Hill Book Company,
1968), P. 155
Main Index
48 MSC Nastran Verification Guide
Reference Deck
Model Name from Old Manual Model Description Location
v10103.dat v2403.dat Truss analysis with thermal loads tpl\verifman\oldverif\
Results
Displacement (in)
u1 v1 u2 v2
Theory -1.272727E-03 6.363636E-03 7.272727E-04 3.636364E-03
MSC Nastran -1.272727E-03 6.363636E-03 7.272727E-04 3.636364E-03
Difference 0.00% 0.00% 0.00% 0.00%
Main Index
CHAPTER 2.3 49
Truss Analysis using Thermal Loads
Forces (lb)
p1 p2 p3 p4 p5 p6
Theory -636.36 899.95 363.64 -514.26 363.64 0.00
MSC Nastran -636.36 899.95 363.64 -514.26 363.64 0.00
Difference 0.00% 0.00% 0.00% 0.00% 0.00% 0.00%
Stresses (lb/in2)
1 2 3 4 5 6
Theory -636.36 1272.73 363.64 -727.27 363.64 0.00
MSC Nastran -636.36 1272.73 363.64 -727.27 363.64 0.00
Difference 0.00% 0.00% 0.00% 0.00% 0.00% 0.00%
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
50 MSC Nastran Verification Guide
Problem Description
Assuming an aircraft is making a level landing with the nose gear slightly off the ground, find the resulting
inertia loads on the pilot due to the main gear landing loads.
Solution Number
SOL 101
Features Used
Inertia Relief Analysis
Reference
1. E. F. Bruhn, Analysis and Design of flight Vehicle Structures (Cincinnati: Tri- state Offset Company,
1965) P. A4.11
Main Index
CHAPTER 2.4 51
Inertia Load on Pilot During Landing
Reference Deck
from Old
Model Name Manual Model Description Location
v10104.dat v2404.dat Inertia load on Pilot during tpl\verifman\oldverif\
Landing
Results
Theoretical Solution
Loads on Pilot: Fx = 162.7 lb Fz = 378.84 lb
Aircraft
Ma = (258.8 lb-sec2/in)
Pilot
Mp = 0.4658 (lb-sec2/in)
Also the aircraft must be adjusted for the pilot.
Main Index
52 MSC Nastran Verification Guide
Ma = 258.8 – 0.4658
Ma = 258.3342 (lb-sec2/in)
The MSC Nastran results are:
Fx = 162.62 lb
Fz = 378.46 lb
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
CHAPTER 2.5 53
Cantilever Beam, Statics, Checkpoint
Problem Description
A simple cantilever beam is loaded with a vertical shear load at the free end. The solution results are saved
(checkpoint) for a restart in a modal solution (v10305r).
Solution Number
SOL 101
Features Used
Static Analysis
Reference
1. R.J. Roark and W. C. Young, Formulas for stress and strain (New York: McGraw-Hill Book Co.
1975) P. 96
Main Index
54 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10105d.dat v2405n.dat Cantilever beam, statics, tpl\verifman\oldverif\
v10305r.dat v0305r.dat checkpoint tpl\verifman\oldverif\
Results
Theoretical Solution
1. Displacement at free end:
3
Fl
= ---------
3EI
3
100 30
= --------------------------------------------
7
-
3 1.0 10 0.0390
= 2.3077 in
2. Moment at fixed end:
M = FL
M = 100(30)
M = 3000 lb-in
Main Index
CHAPTER 2.5 55
Cantilever Beam, Statics, Checkpoint
3000 0.50
= -------------------------------
0.0390
= 38462 lb/in2
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
56 MSC Nastran Verification Guide
Problem Description
A coordinate system exists that is rotated 45o with respect to a basic frame of reference in the x-y plane. A bar
located in this system is simply supported by elastic supports. A load is applied to the bar at an intermediate
point. Calculate the reaction forces in the supports.
Solution Number
SOL 101
Features Used
Static Analysis
Reference
1. Joseph Shigley, Mechanical Engineering Design (New York: McGraw-Hill Inc., 1977) P. 642
Main Index
CHAPTER 2.6 57
Spring Supported Bar in a Skewed Coordinate System
Reference Deck
from Old
Model Name Manual Model Description Location
v10106.dat V2406.dat Spring supported bar in tpl\verifman\oldverif\
skewed system
Results
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
58 MSC Nastran Verification Guide
Problem Description
A simple structurally intermediate problem (see schematic below) is analyzed using MSC NASTRAN with
the resulting displacement and stress output punched.
Solution Number
SOL 101
Features Used
Static Analysis
Main Index
CHAPTER 2.7 59
Punched Output of a Static Solution
Reference Deck
Model Name from Old Manual Model Description Location
v10107.dat V2407.dat Punched output of static solution tpl\verifman\oldverif\
Main Index
60 MSC Nastran Verification Guide
Results
Punched output of static solution by MSC NASTRAN is as follows.
Main Index
CHAPTER 2.8 61
Cantilever Beam Constructed of Plates, Stress Contours
Cantilever Beam
2.8 Constructed of Plates,
Stress Contours
Problem Description
A Simple cantilever beam is constructed using plate element and is loaded with vertical shear force at one
end. Find the stress at the outside fibers.
Solution Number
SOL 101
Features Used
Static analysis
Reference
1. Ferdinand L. Singer, Strength of Materials (Harper & Row, 1962) Art. 52, pp. 133.
Main Index
62 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10108.dat V2408a.dat Cantilever beam tpl\verifman\oldverif
constructed of plates, stress
contours
Results
Theoretical Solution
Moment of Inertia:
3 3
th 0.10 0.6
I = ------- = -------------------------- = 1.8E-3 in4
12 12
Main Index
CHAPTER 2.8 63
Cantilever Beam Constructed of Plates, Stress Contours
Stresses (psi)
Fibre Distance C = 0.0 C = 0.10 C = 0.20 C = 0.30
Theory 0 20000 40000 60000
MSC Nastran 0 19400 39350 54444
%Difference 0.00% -3.00% -1.63% -9.26%
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
64 MSC Nastran Verification Guide
Problem Description
A coordinate system exists that is rotated 45o with respect to a basic frame of reference in the x-y plane. A bar
located in this system is simply supported by elastic supports. The grid points are located along the x1 axis
and the bar at an intermediate point. Calculate the reaction forces in the supports.
Solution Number
SOL 101
Features Used
Static Analysis
Reference
1. Joseph Shigley, Mechanical Engineering Design (New York: McGraw-Hill Inc., 1977) P. 642
Main Index
CHAPTER 2.9 65
Bar with Offsets in a Skewed Coordinate System
Reference Deck
from Old
Model Name Manual Model Description Location
v10109.dat V2409.dat Bar with offsets in skewed tpl\verifman\oldverif\
coordinate system
Results
Theoretical Solution
Fb 600 2
R 1 = ------- = ----------------- = 400.0 lb
l 3
Fa 600 1
R 2 = ------- = ----------------- = 200.0 lb
l 3
Result comparison of theoretical solution and MSC Nastran results are as follows:
Main Index
66 MSC Nastran Verification Guide
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
CHAPTER 2.10 67
Thick Walled Cylinder, Internal Pressure
Problem Description
Find the stress distribution for a thick-walled cylinder subjected to internal pressure.
Solution Number
SOL 101
Features Used
Static Analysis
Reference
1. S. H. Crandall, N. C. Dahl, and T. J. Lardner, An Introduction to the Mechanics of Solids, 2nd ed.
(New York: McGraw-Hill Book Company, 1972) P. 293-297
Main Index
68 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10110.dat V2410.dat Thick walled cylinder, tpl\verifman\oldverif\
internal pressure
Results
Theoretical Solution
2
–Pi ro r – 1
r = -------------------------------------------
2
-
ro ri – 1
2
Pi ro r + 1
= ----------------------------------------
2
-
ro ri – 1
Result comparison of theoretical solution and MSC Nastran results are as follows:
r lb/in2 lb/in2
r r/ri Theory MSC Nastran % Error Theory MSC Nastran % Error
6.0 1 -10.00 16.67
6.50 1.083 -8.030 -8.095 0.811% 14.69 14.852 1.104%
7.0 1.167 -6.460 13.13
7.50 1.25 -5.200 -5.238 0.731% 11.87 11.956 0.722%
Main Index
CHAPTER 2.10 69
Thick Walled Cylinder, Internal Pressure
r lb/in2 lb/in2
r r/ri Theory MSC Nastran % Error Theory MSC Nastran % Error
8.0 1.333 -4.170 10.83
8.50 1.417 -3.310 -3.333 0.704% 9.98 10.031 0.509%
9.0 1.5 -2.590 9.26
9.50 1.583 -1.990 -2.000 0.502% 8.65 8.686 0.421%
10.0 1.667 -1.470 8.13
10.50 1.75 -1.020 -1.030 1.009% 7.69 7.710 0.262%
11.0 1.833 -0.630 7.3
11.50 1.917 -0.300 -0.303 1.009% 6.96 6.979 0.271%
12.0 2 0.00 6.67
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
70 MSC Nastran Verification Guide
Problem Description
Find the radial deflection and meridional moments in a cantilevered cylindrical shell.
Solution Number
SOL 101
Features Used
Static Analysis
Reference
1. "Analysis of Axisymmetrical shells by the Direct Stiffness Method" P. E. Grafton and D. R. Strome,
AIAA Journal, 1(10): 2342-2347
2. "Evaluation of NASTRAN" J. W. Jones and H. H. Fong, Structural Mechanics Software Series, Vol.
IV (N. Perrone and W. Pilkey, eds.), 1982
3 3
Bh 1 0.01 –8
Moment of inertia per unit width: I = --------- = --------------------- = 8.33 10
12 12
Main Index
CHAPTER 2.11 71
Cantilevered Cylindrical Shell, Axisymmetric Force
Reference Deck
from Old
Model Name Manual Model Description Location
v10111.dat V2411.dat Cantilevered Cylinder shell tpl\verifman\oldverif\
Results
Theoretical Solution
Q - – X
Radial deflection: U r = -------------
3
e cos X
2 D
Q – X
Meridional moment: M = ---- e sin X
where,
= 1.285 Rt = 5.746695 in-1
Result comparison of theoretical solution and MSC Nastran results are as follows:
Main Index
72 MSC Nastran Verification Guide
Main Index
CHAPTER 2.11 73
Cantilevered Cylindrical Shell, Axisymmetric Force
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
74 MSC Nastran Verification Guide
Problem Description
A truss consisting of two pin-jointed axial elements is loaded vertically. Find the vertical deflection of the end
and the internal forces.
Solution Number
SOL 101
Features Used
Static Analysis using Superelements
Reference
1. Egor P. Popov, Introduction to Mechanics of Solids (Englewood Cliffs: Prentice-Hall, Inc.. 1968) P.
499.
Main Index
CHAPTER 2.12 75
Static Analysis of a Pin-Joined Truss Using Superelements
Reference Deck
from Old
Model Name Manual Model Description Location
v10112s.dat V6101s.dat Static analysis of pin joined truss tpl\verifman\oldverif
Results
Theoretical Solution
At point B: U y = – 0.0444 in
Force in Element 1: +2500. lb (tension)
Force in Element 2: -2500. lb (compression)
At point B: U y = – 0.0444 in
Force in Rod 1: +2500. lb
Force in Rod 1: -2500. lb
Main Index
76 MSC Nastran Verification Guide
Comparison of Results
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.13 77
Beam with Point Masses, Enforced Acceleration
Problem Description
A beam, with equal point masses on each end, is given a unit acceleration on its center point normal to its
axis. Find the reaction force, Q, required to produce this acceleration, and the resulting deflections relative
to the point of application of the force.
Solution Number
SOL 101
Features Used
Superelement Statics with Inertial relief
Reference
1. Raymond J. Roark, Formulas for Stress and Strain, 2nd ed. (New York: mcGraaw-Hill Book
Company, 1965), P. 104.
Main Index
78 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10113s.dat V9101s.dat Beam with point masses, tpl\verifman\oldverif
enforced acceleration
Results
Theoretical Solution
From symmetry arguments, the loading condition is equivalent to a point load on a cantilevered beam of half
the total length of this beam. The reference gives the load-deflection equation:
3
p3 l
u 3 = ----------
3EI
Where P3 is the applied load equal to the unit tip mass times the unit acceleration, or unity.
Therefore,
3
10 - –5
u 3 = ------------------------
7
= 1.111... 10
3 10 3
The reaction force is twice the applied load for one side, therefore,
Q = – 2.0
Main Index
CHAPTER 2.13 79
Beam with Point Masses, Enforced Acceleration
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
80 MSC Nastran Verification Guide
Problem Description
A simple frictionless two degree of freedom system is constructed using springs and concentrated masses. Find
all natural frequencies and mode shapes.
Solution Number
SOL 103
Features Used
Normal Modes Analysis
Inverse Power Method and Lanczos Method
Reference
1. Donald. T. Greenwood, Principles of Dynamic (Englewood Cliffs: Prentice-Hall, Inc., 1965) p.459,
ex. 9-1.
Main Index
CHAPTER 2.14 81
Two Degree of Freedom Modal Analysis
Figure 2.14-1 Two degree of freedom system using spring and concentrated masses
Results
MSC Nastran by using rods (CONROD) and Concentrated mass (CONM2) the inverse power method
gives following results and are compared with theoretical results:
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
82 MSC Nastran Verification Guide
Problem Description
Find all the natural frequencies of the clamped-clamped massless beam with the concentrated masses, shown
below. (Note that the model produces uncoupled modes with integer eigenvalues.)
Solution Number
SOL 103
Features Used
Modal Analysis.
Component Modal Synthesis
Givens Method
Reference
1. J. P. Den Hartog, Mechanical Vibrations 4th ed., (New York: McGraw-Hill Book Company, 1956),
pp. 429-430.
Main Index
CHAPTER 2.15 83
Normal Modes of a Beam Structure using Component Modal Synthesis
Main Index
84 MSC Nastran Verification Guide
Results
Theoretical Solution
For single degree of freedom systems (see (2.15-18), for reference)
K
= ---- (2.15-1)
N
u
f = ------ , cycles/ unit time (2.15-3)
2
Axial Extension Mode (T1)
(Axial springs in parallel)
EA EA 4EA
K 1 = --------- + --------- = ----------- = 4.0 (2.15-4)
l2 l2 l
m 2 = 64 (2.15-7)
192EI 1
K 3 = ------------------ = 384 (2.15-9)
3
l
m 3 = 64 (2.15-10)
l 1 = 1.0 (2.15-13)
Main Index
CHAPTER 2.15 85
Normal Modes of a Beam Structure using Component Modal Synthesis
I 33 = 16 (2.15-16)
16EI 1
K 2 = --------------- = 32 (2.15-18)
3
l
I 22 = 15. (2.15-19)
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
86 MSC Nastran Verification Guide
Problem Description
Find the torsional modes of a shaft fixed at one end.
Solution Number
SOL 103
Features Used
Modal Analysis.
Reference
1. Walter C. Hurty and Moshe F. Rubinstein, Dynamics of Structures (Englewood Cliffs: Prentice-
Hall, Inc.. 1964).
Main Index
CHAPTER 2.16 87
Torsional Vibration of a Shaft with Three Disks
Results
MSC Nastran
Frequency (rad/sec) Theory MGIV LAN
1 89.00 89.008 89.008
2 249.40 249.396 249.396
3 360.40 360.388 360.388
Conclusion
The results from MSC Nastran are comparable with theoretical results.
Main Index
88 MSC Nastran Verification Guide
Problem Description
Find the two lowest natural frequencies of a cantilever beam, Restart the run from a previous statics solution,
v10105d.
Solution Number
SOL 103
Features Used
Eigenvalue solution using Modified Givens method. (Restart)
Reference
1. R.J. Roark and W. C. Young, Formulas for stress and strain (New York : McGraw-Hill Book Co..
1975) P. 576
Main Index
CHAPTER 2.17 89
Cantilever Beam, Normal Modes, Restart
Results
Theoretical Solution
1.732
f 1 = ------------- EIg
---------
2 wl 3
For I2 = 0.0241
1.732
f 1 = ------------- 1.0 + 7 0.0241
---------------------------------------
2 1.0 30
3
f1= 0.823 Hz
For I1 = 0.039
Main Index
90 MSC Nastran Verification Guide
1.732
f 2 = ------------- 1.0 + 7 0.039
------------------------------------
2 1.0 30
3
f2= 1.048 Hz
Frequency (Hz)
Theory MSC Nastran
Mode 1 0.823 0.823583
Mode 2 1.048 1.047686
Conclusion
The results from MSC Nastran are comparable with theoretical results. Note: Modified Givens method is
good for use on small problems, but for large problems, Lanczos method is preferred.
Main Index
CHAPTER 2.18 91
Normal Modes Analysis with Reduction Techniques
Problem Description
Find the two lowest modes of a simply supported shaft. Compare the results of the inverse power iteration
method without reduction, the Givens method with Guyan reduction and the Lanczos method.
Solution Number
SOL 103
Features Used
Modal Analysis
Reference
1. J. P. Den Hartog, Mechanical Vibration 4th ed.. (New York: McGraw- Hill, 1956), P. 432
Main Index
92 MSC Nastran Verification Guide
E = 30E+7 lb/in2
= 0.3
= 7.764E-4 lb sec2 / in4 (Mass density)
Main Index
CHAPTER 2.18 93
Normal Modes Analysis with Reduction Techniques
Results
Theoretical Solution
The reference gives the natural frequency equation as
a n EI
f n = ------ -----------
- (2.18-1)
2 L 4
1
2
where 1 is the mass per unit length, n is the mode index and a n = n . Approximate equations for thin
walled circular cylinder section properties are:
I = r2t = 0.15708 (2.18-2)
A = 2rt =0.31416 (2.18-3)
r E 2
- n = 21.8335n2 cycles/second
f n = ------ ----- (2.18-5)
e 2 2
Conclusion
The model was run using three methods. The inverse power method and Lanczos method involve no
approximation and are used to give a measure of the discretization error caused by using 20 elements with
lumped masses. The error in natural frequency is less than one tenth of one percent, even for worst case.
For Guyan reduction, the error for the first and second mode is well below one percent, the criteria used for
picking the number of A-set points. However, the error rapidly climbs to 60 percent for the sixth mode.
Main Index
94 MSC Nastran Verification Guide
Problem Description
Find the natural modes of fixed circular plate.
Solution Number
SOL 103
Features Used
Modal Analysis.
Inverse Power method with enhancements and Lanczos method
Reference
1. S. Timoshenko, D. H. Young and W. Weaver. Jr., Vibration Problems in Engineering, 4th ed. (New
York: John Wiley and Sons, 1947) p. 501
Main Index
CHAPTER 2.19 95
Modes of Fixed Circular Plate with Trapezoidal Ring Element
Results
Theoretical Solution
D
= ----- -----
- (2.19-1)
2
a h
3
Eh
D = --------------------------- (2.19-2)
2
12 1 –
3
30.E + 6 0.125
D = ------------------------------------------------ (2.19-3)
2
12 1 – 0.3
D = 5365.72 (2.19-4)
5365.72 -
= --------------- ---------------------------------------- (2.19-5)
2 7.764E-4 0.125
3.5
Main Index
96 MSC Nastran Verification Guide
F = (606.988) (2.19-6)
Frequencies for N = 0
S (rad/sec)
0 10.21 6197.35
1 39.78 24145.98
2 88.90 53961.23
Comparison of Results
Conclusion
For MSC Nastran results using Inverse power method with enhancement and Lanzos, error is less than one
percent for first two modes and it rises to approximately two percent for third mode.
Main Index
CHAPTER 2.20 97
Two Degree of Freedom Modal Analysis Using Superelements
Problem Description
A simple frictionless two degree of freedom system is constructed using spring and concentrated masses. Find
all natural frequencies and mode shapes.
Solution Number
SOL 103
Features Used
Modal Analysis
Superelements Method
Reference
1. Donald. T. Greenwood, Principles of Dynamic (Englewood Cliffs: Prentice-Hall, Inc., 1965)
P.459, Ex. 9-1.
Main Index
98 MSC Nastran Verification Guide
Figure 2.20-1 Two degree of freedom system using spring and concentrated masses
Reference Deck
from Old
Model Name Manual Model Description Location
v10308s.dat V6301s.dat Two D.O.F. system modal analysis tpl\verifman\oldverif
using superelement
(Inverse Power Method).
Results
MSC Nastran by using rods (CONROD) and Concentrated mass (CONM2) the inverse power method
gives following results and are compared with theoretical results.
Eigenvectors
Frequency
(rad/sec) Theoretical MSC Nastran (INV)
Mode 1 10.83 1.000 1.000
0.531 0.531
Mode 2 46.18 -0.133 -0.133
1.000 1.000
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.21 99
Modes of Fixed Circular Plate with Superelement Method
Problem Description
Find the natural modes of fixed circular plate.
Solution Number
SOL 103
Features Used
Modal Analysis
Superelements Method
Reference
1. S. Timoshenko, D. H. Young and W. Weaver. Jr., Vibration Problems in Engineering, 4th ed. (New
York: John Wiley and Sons, 1947) P. 501
Main Index
100 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10309s.dat V6302s.dat Modes of fixed circular plate with tpl\verifman\oldverif
trapezoidal ring element using inverse
power method with enhancement
Results
Theoretical Solution
D
= ----- -----
- (2.21-1)
2
a h
3
Eh
D = --------------------------- (2.21-2)
2
12 1 –
3
30.E + 6 0.125
D = ------------------------------------------------ (2.21-3)
2
12 1 – 0.3
D = 5365.72 (2.21-4)
5365.72 -
= --------------- ---------------------------------------- (2.21-5)
2 7.764E-4 0.125
3.5
F = (606.988) (2.21-6)
Main Index
CHAPTER 2.21 101
Modes of Fixed Circular Plate with Superelement Method
Frequencies for N = 0
S (rad/sec)
0 10.21 6197.35
1 39.78 24145.98
2 88.90 53961.23
Comparison of Results
Frequency (rad/sec)
Mode No. Theory MSC Nastran %Difference
1 6197.35 6107.253 -1.4538%
2 24145.98 23568.68 -2.3909%
3 53961.23 52342.36 -3.0001%
Conclusion
MSC Nastran results are comparable with theoretical results using superelement method.
Main Index
102 MSC Nastran Verification Guide
Natural Frequencies of a
2.22 Circular Plate Immersed in
Fluid
Problem Description
Find the natural modes of fixed circular plate immersed in fluid.
Solution Number
SOL 103
Features Used
Modal Analysis
Reference
1. S. Timoshenko, D. H. Young and W. Weaver. Jr., Vibration Problems in Engineering, 4th ed. (New
York: John Wiley and Sons, 1947) P. 501
Fluid:
f = 0.0385 lb/in3 (1/g = 0.002588 sec2/in)
Main Index
CHAPTER 2.22 103
Natural Frequencies of a Circular Plate Immersed in Fluid
Reference Deck
from Old
Model Name Manual Model Description Location
v10310s.dat V6303s.dat Natural frequency of a circular plate tpl\verifman\oldverif
immersed in fluid
Results
Theoretical Solution
First mode:
2
1 10.21 1 Et
f = ------ ------------------ ----- ------------------------------------------------------------------
- (2.22-1)
2 1 + r 2 1 – 2 0.002588 12
p
f r
= 0.6689 ----- - (2.22-2)
p t
0.0385 3.5
= 0.6689 ---------------- ------------- = 2.4036 (2.22-3)
0.3 0.125
Main Index
104 MSC Nastran Verification Guide
Comparison of Results
Frequency (rad/sec)
Mode No. Theory MSC Nastran %Difference
1 534.6 568.16 6.27%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.23 105
DDAM Method of Response Spectra
Problem Description
The purpose of this Verification Problem is to show how to calculate peak structural response estimates in
compliance with the U.S. Navy shock design modal summation convention (knowns as the NRL method).
The NRL convention calls for a choice among modal contributions to peak response
u ai = ai q i
It is noted that the maximum contributing mode for a particular response is not necessarily the same for all
physical responses.
This problem demonstrates input of a shock spectrum to a frame structure, to calculate internal loads and
accelerations at points on the structure.
Solution Number
SOL 103
Features Used
Modal Analysis
Reference
1. "Shock Design Criteria for Surface Ships", NAVSEA 0908-LP-000-3010, May 1976, P. 21.
Main Index
106 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10311d.dat V6304d.dat DDAM method of response spectra tpl\verifman\oldverif
v10311r.dat V6304r.dat tpl\verifman\oldverif
Results
Theoretical Solution
Peak absolute modal responses were calculated for the eight flexible body modes (4-11). They are printed in
the run with the title “MATRIX UHVR”.
Main Index
CHAPTER 2.23 107
DDAM Method of Response Spectra
Estimates for peak acceleration of grid point 6 in the x-direction and bending moment in element 104
(end A) are shown below on a mode-by-mode basis as well as the NRL sum. Corresponding results for the
standard SRSS and ABS conventions are also given.
The 6x i column is extracted from the “REAL EIGENVECTOR” output of the first run.
The u··6x i column can be calculated from the product of the u··h column times the 6x i data.
The M104A column is extracted from the “FORCES IN BAR ELEMENT” output.
The M104A column can be calculated from the product of the uh column times the M104A data.
The enclosed modal quantities represent the largest absolute contributor among selected nodes.
Note that the second mode provides the highest acceleration component.
Main Index
108 MSC Nastran Verification Guide
Comparison of Results
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.24 109
Lateral Buckling of Cantilever Beam
Problem Description
A Cantilever beam of constant cross section is loaded at the free end with a vertical shear load. The critical
buckling condition occurs when the deflected shape in the vertical plane is unstable and lateral buckling
occurs. Determine the critical buckling eigenvalues and loads for the structure.
Solution Number
SOL 105
Features Used
Buckling Analysis
Reference
1. Stephen P. Timoshenko and James M. Gere. Theory of Elastic Stability, 2nd ed. (New York:
Mc Graw-Hill Book Co. 1961), P. 257
Main Index
110 MSC Nastran Verification Guide
Results
Theoretical Solution
P CR = EI GJ - N
----------------------- (2.24-1)
2 CR
L
Lateral Buckling of a Cantilever Beam
NCR = ±2.0063, ±5.1231, ±8.2580 (Derived in the reference) (2.24-2)
PCR = ±11.27, ±28.77, ±46.37 (2.24-3)
MSC Nastran Solution
P CR
Eigenvalue = ---------- , solved from K – K d u = 0
P
Letting P = 50, then
Main Index
CHAPTER 2.24 111
Lateral Buckling of Cantilever Beam
Comparison of Results
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
112 MSC Nastran Verification Guide
Problem Description
A simple frame with lateral motion constrained, is loaded vertically. Using Static analysis, find the vertical
column forces and stresses. Then release the constraint to allow lateral motion and find the lowest critical
buckling load.
Solution Number
SOL 105
Features Used
Static and Buckling Analysis
Reference
1. Stephen P. Timoshenko and James M. Gere. Theory of Elastic Stability, 2nd ed. (New York:
Mc Graw-Hill Book Co. 1961), P. 149
Main Index
CHAPTER 2.25 113
Simple Frame Analysis with Buckling
Results
Theoretical Solution
Static Analysis of Symmetric Deformation Mode
Axial force in vertical column: -500.0 lb (comp.)
Axial stress of vertical column: -500.0 lb/in2
Buckling Analysis of Antisymmetric Deformation Mode
2 PC R
k = ----------
EI
where kL = -6 tan (kL)
lowest k: k1 = 2.71645974 x 10-2 p = k2EI = 737.9 lb
Main Index
114 MSC Nastran Verification Guide
Comparison of Results
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.26 115
Beam with Gap/Lift-Off
Problem Description
A simply supported beam is hinged at one end and supported by lifting rollers at two other locations.
Allowing lift-off to occur at the rollers, determine the vertical deflections under the load points.
Solution Number
SOL 106
Features Used
Nonlinear statics
Reference
1. Jack C. McCormac, Structural Analysis, 3rd ed. (New York: Intext Educational Publishers, 1957) P.
323, ex 17.12.
Main Index
116 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10601s.dat V6601s.dat Beam with gap/lift-off tpl\verifman\oldverif
Results
Theoretical Solution
(Lift-off occurs at point D only)
Point B: Uy = -1.01 in
Point D: Uy = +0.546 in
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
CHAPTER 2.27 117
Cantilever Beam with Plastic Hinge
Problem Description
Given a cantilever beam with a concentrated end loading, find the collapse load using plastic hinge analysis.
Solution Number
SOL 106
Features Used
Nonlinear statics
Reference
1. A. Mendelson, PLASTICITY: Theory and Application (New York: The Macmillan Company)P.
308
Fi = 16.0
F = 0.5
Main Index
118 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10603s.dat V6603s.dat Cantilever beam with tpl\verifman\oldverif
plastic hinge
Results
Theoretical Solution
With a concentrated end load, the maximum moment occurs at the base. Collapse occurs when the moment
at the base causes the root cross section to become fully plastic.
2
Collapse Load: F = 2bh yld L = 25.0
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
CHAPTER 2.28 119
Complex Eigenvalues of Three-by-Three Matrix
Problem Description
Given the real, unsymmetric matrix [A], find its characteristic values (“eigenvalues”) and modal matrix
(“matrix of eigenvectors").
Solution Number
SOL 107
Features Used
Direct Complex Eigenvalues
Reference
1. Paul M. DeRusso, Rob J. Roy, and Charles M. Close, State Variables for Engineers (New York: John
Wiley & Sons, 1967), ex. 4.7-2
Engineering Data
2 –2 3
A= 1 1 1
1 3 –1
Main Index
120 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10701.dat V2801.dat Complex eigenvalues of three-by- tpl\verifman\oldverif
three matrix
Results
Theoretical Solution
Evaluate A – I = 0 to produce the characteristic equation 3 – 2 2 – 5 + 6 = 0 , whose roots are
tabulated below. Their validity can be determined by substitute into the characteristic equation. The reference
calculates the eigenvector matrix [M], tabulated below, by use of adjoint matrix techniques. Its validity can
be checked by back-substitution into the equation A – i I M = 0 , where {Mi} is the ith column of M.
Eigenvectors are arbitrary to a scaling constant. The MSC Nastran results were renormalized by hand to
obtain the results tabulated above.
Conclusion
The closed form and tabulated solutions agree exactly to all printed digits. (MSC NASTRAN outputs 6 or 7
digits depending on machine type.)
Main Index
CHAPTER 2.29 121
Complex Roots as a Function of Gain
Problem Description
Find the root-locus of the control system shown in the figure.
Solution Number
SOL 107
Features Used
Direct Complex Eigenvalues
Reference
1. C. J. Savant, Jr., Basic Feedback Control System Design (New York: McGraw-Hill Book Company,
Inc., 1958)
Main Index
122 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10702.dat V2802.dat Complex roots as a function of gain tpl\verifman\oldverif
Results
Theoretical Solution
The denominator of the transfer function is,
2 1 AK m K s
p + --- p + ------------------ = 0
Solutions of this quadratic equation are tabulated for three values of AKs as given on page 82 of the reference.
They are:
Root Location
AKs P1 P2
1.0 -94.7 -5.3
5.0 -50.0 -50.0
10.0 -50 + j50 -50 - j50
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.30 123
Frequency Response of a Single Degree of Freedom System
Frequency Response of a
2.30 Single Degree of Freedom
System
Problem Description
A 3000-pound weight is attached to a platform by a spring and a damper. If the platform is moved up and
down at the system's resonant frequency with an amplitude of one inch, find the maximum amplitude of the
3000-pound weight.
Solution Number
SOL 108
Features Used
Direct Frequency Response
Reference
1. J.P. Den Hartog, Mechanical vibrations (New York: McGraw-Hill Book Company,1956), pp. 49-55
3000 lb
M 1 = -----------------------------------------
32.2 ft- 12
---------------- inch
------------------
sec
2 1 ft
M 1 = 7.764 lb-sec2/in
Main Index
124 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10801.dat V2601.dat Frequency response of tpl\verifman\oldverif\
single degree of system
Results
Theoretical Solution
The motion of M1 is given by equation 2.28a of the reference (p. 49)
P0 k
u = ----------------------------------------------------------------------------------------------
2 2 2 2
1 – n + 2 c cc n
where,
2 2
P0 = ka 0 + ca 0
c c = 2 km
Main Index
CHAPTER 2.30 125
Frequency Response of a Single Degree of Freedom System
Amplitude (inches)
Theory 2.0731
MSC Nastran 2.073576
% Error 0.023%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
126 MSC Nastran Verification Guide
Transient Response of a
2.31 Rocket - Direct Time
Integration
Problem Description
A rocket is subjected to axial loading for a finite duration of time. Using axial elements in a free-free structure,
find the displacement histories of the midpoint and two ends, considering both the rigid body motion and
the structural vibration.
Solution Number
SOL 109
Features Used
Direct Transient Response
Reference
1. J.S. Przemieniecki, Theory of Matrix structural Analysis (New York: McGraw-Hill Book Co., 1968)
P. 340
Main Index
CHAPTER 2.31 127
Transient Response of a Rocket - Direct Time Integration
Reference Deck
from Old
Model Name Manual Model Description Location
v10901.dat V2701.dat Transient response of a rocket - direct tpl\verifman\oldverif
time integration
Results
Theoretical Solution
Using the matrix forms for the equations of motion, with zero initial conditions, the base response
for t t 0 can be written as:
P 0 t 0 2t – t 0 P 0 3 2 1
- ------ cos 2 t – t 0 – cos 2 t + ------ cos 3 t – t 0 – cos 3 t
u 1 t = -------------------------------- + ------- -----
2AL 4A E 2 3
where,
2 = 12E
---------- ; 3 = 48E
----------
2 2
L L
Main Index
128 MSC Nastran Verification Guide
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.32 129
Response Spectra of a Single Oscillator
Problem Description
A large mass, M1, is subjected to a force doublet. A very small mass, M2, is attached to it with a spring and
damper, tuned to an undamped natural frequency of one Hz. A response spectrum is computed for the large
mass at this same frequency. A variable, v3, is connected by a multipoint constraint and is a model of the
relative displacement between the large and small mass. As the small mass-spring combination has the same
physical model as the response spectrum calculation, it should produce the same results as the inertial motion
component (variable 2) and the relative motion component (variable 3).
Solution Number
SOL 109
Features Used
Direct Transient Response
Reference
1. Harris, Cyril M. and Crede, Charles E., Shock and Vibration Handbook, (New York: McGraw-Hill
Book Company, Inc., 1961)
= 0.125664
0 = 1.0E+6
Main Index
130 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10902.dat V2702.dat Response spectra of a single oscillator tpl\verifman\oldverif
Results
Magnitude of Peak Response, Inertial Coordinates
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.33 131
Transient Response with Combined Loading Functions
Problem Description
Three identical structures are given different time-varying loads. The loads on the third structure are a linear
combination of the loads on the other two structures. It is shown that the response of the third structure is
also the same linear combination of the response of the first two structures.
Solution Number
SOL 109
Features Used
Direct Transient Response
Main Index
132 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10903.dat V2703.dat Transient response with tpl\verifman\oldverif
combined loading function
Results
All three systems are sampled at t = 1.5, with their displacements tabulated below seen that the response of
the third structure is the sum of the other two.
Response at t=1.5
Mass 1 1.576766E-05
Mass 2 -1.361804E-05
sum 2.149620E-06
Mass 3 2.149617E-06
% Difference 0.000%
Conclusion
It is proved by MSC Nastran that response of the third structure is also same linear combination of the
response of the first two structures.
Main Index
CHAPTER 2.34 133
Direct Transient Response to an Internal Pressure
Problem Description
Find the peak displacements in a thick-walled cylinder resulting from a suddenly applied internal pressure
load.
Solution Number
SOL 109
Features Used
Direct Transient Response
Reference
1. S. H. Crandall, N. C. Dahl, and T.J. Lardner, An Introduction to the Mechanics of Solids, 2nd ed.
(New York: McGraw-Hill book Company, 1972) pp. 293-297
Main Index
134 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v10904.dat V2704.dat Direct transient response to an tpl\verifman\oldverif
internal pressure
Results
Theoretical Solution
Radial displacements from a static analysis (SOL 101) were: Grid 1, 3.326E-6 in.; Grid 7, 2.773E-6 in.; and
Grid 11, 2.673E-6 in.
MSC Nastran Solution
The force was applied in an instantaneous manner, as shown below. Several time steps of zero load were input
to ensure simulation of a step function in the loading. The integration time step was chosen to be 1.0
microsecond. Approximately three cycles of response were simulated, and damping was neglected. Note that
Main Index
CHAPTER 2.34 135
Direct Transient Response to an Internal Pressure
for a single degree-of-freedom system, an instantaneously applied load will give a maximum displacement of
twice the static displacement.
Comparison of Results
Conclusion
For a Single degree of freedom system, an instantaneously applied load will give a maximum displacement
of twice the static displacement.
Main Index
136 MSC Nastran Verification Guide
Complex Eigenvalues of
2.35 Three-by-Three Matrix,
Modal Method
Problem Description
Given the same unsymmetric matrix [A] used in problem v10701, find its eigen solution in two steps.
1. Add enough terms to [A] to produce a symmetric matrix. Find the real eigensolution of this matrix
to use as modal coordinates.
2. Add in the terms to reconstruct [A], and compute its eigensolution.
Use the so-called 'particular' matrix [AB] of reference 2 as the symmetric matrix, as it has a closed-form
solution for its eigenvalues. Use the negative of this matrix, as this is more typical of the values encountered
on structures.
Solution Number
SOL 110
Features Used
Modal Complex Eigenvalues
Reference
1. Verification Problem No. v10701
2. Marvin Marcus, 'Basic Theorems in Matrix Theory, "National Bureau of Standards Applied
Mathematics Series, 57 (1964)"
Main Index
CHAPTER 2.35 137
Complex Eigenvalues of Three-by-Three Matrix, Modal Method
[A] is tabulated in v10701. [AB] has -2 on the diagonal, +1 above and below the diagonal, and zeros
elsewhere. The difference matrix [D] is then
[D] = [A] - [-AB] = [A] + [AB]
0 –1 3 2 –2 3 –2 1 0
2 –1 2 = 1 1 1 + 1 –2 1
1 4 –3 1 3 –1 0 1 –2
Reference Deck
from Old
Model Name Manual Model Description Location
v11001.dat V2901.dat Complex eigenvalues of three-by- tpl\verifman\oldverif
three matrix, modal method
Results
Theoretical Solution
The theory of the complex eigensolution is described in v10701. A closed form solution for the real
eigensolution to [AB] is given in Reference 2 as
2
i = – 4 sin i 2 N + 1
2
1 = 4 sin --- = 0.585786
8
2
2 = 4 sin --- = 2.0
4
3 2
3 = 4 sin ------ = 3.414214
8
Main Index
138 MSC Nastran Verification Guide
Comparison of Results
Eigenvalues
Theory MSC Nastran % Difference
Mode 1 0.5858 0.5858 0.00%
Mode 2 2.0000 2.0000 0.00%
Mode 3 3.4142 3.4142 0.00%
Both the intermediate solution to the real eigenvalue problem, tabulated above, and the final complex eigen
solution, tabulated in v10701, agree to all places printed.
Conclusion
The MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.36 139
Modal Frequency Response of a Beam Structure
Problem Description
Find the response at response of the model described in problem v10302 for six values of damping.
Solution Number
SOL 111
Features Used
Modal Frequency Response
Reference
1. J. P. Den Hartog, Mechanical Vibrations, 4th ed. (New York: McGraw- Hill Book Company, 1956)
pp. 47-49
Main Index
140 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v11101.dat V3001.dat Modal frequency response of a beam tpl\verifman\oldverif
structure
Results
Theoretical Solution
The response of a single degree of freedom oscillator at resonance is given in the reference by equation 2.28a
(page 49). Simplifying that equation for our case gives,
P0 k
x = --------------- (2.36-1)
= n 2c c c
Where,
x = Magnitude of displacement
P0 = Magnitude of applied load
c = viscous damping coefficient
cc = viscous damping coefficient required for critical damping
Main Index
CHAPTER 2.36 141
Modal Frequency Response of a Beam Structure
Phase angle data is plotted on page 51 of the reference. The graph shows a value of 900 for all values of
damping.
As shown in problem v0302, each degree of freedom shows different stiffness. Here we apply a load
proportional to each stiffness at its corresponding degree of freedom ( P 0 = k ).
With this assumption, equation reduces to
1 -
x - --------------
------------- = (2.36-2)
x static 2c c c
P0
Where, x static = ------
k
Comparison of Results
Conclusion
The MSC Nastran response solutions agree with the above closed-form solutions to all places printed. The
MSC Nastran phase angles are also in agreement with the theoretical results, after accounting for the MSC
Nastran sign convention.
Main Index
142 MSC Nastran Verification Guide
Problem Description
A simply supported beam is mounted on elastic springs. Given a time-dependent symmetric distributed
loading function, find the transient response of the midspan deflections. Use the Modal method, considering
symmetric modes only.
Solution Number
SOL 112
Features Used
Modal Transient Response
Reference
1. Walter C. Hurty and Moshe F. Rubinstein, Dynamics of Structure, (Englewood Cliffs: Prentice-Hall,
Inc., 1964), P. 294.
t0 = 1.0
P0 = 100.0
Main Index
CHAPTER 2.37 143
Transient Response of a Beam Modal Method
Reference Deck
from Old
Model Name Manual Model Description Location
v11201.dat V3101.dat Modal frequency response tpl\verifman\oldverif
of a beam structure
Results
Theoretical Solution
Using a Rayleigh-Ritz approximation of the first three symmetric modes, an equation for the center
deflection from the modal method is
3
–P0 l
= 0 t = -------------
4
- 2.603D 1 t + 0.050D 2 t (2.37-1)
EI
where,
1
D 1 t = ----------- sin 1 t – t 0 – sin 1 t + cos 1 t – t 0 (2.37-2)
1 t0
1
D 2 t = ----------- sin 2 t – t 0 – sin 2 t + cos 2 t – t 0 (2.37-3)
2 t0
2
EI
1 = 0.725 -----
2
------ = 7.155 (2.37-4)
l m
2
EI
2 = 2.860 -----
2
------ = 28.3 (2.37-5)
l m
Main Index
144 MSC Nastran Verification Guide
Comparison of Results
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
CHAPTER 2.38 145
Deflection of a Circular Plate, Cyclic Symmetry
Problem Description
Find the maximum displacement of a circular plate with a uniform pressure load and fixed edges.
Solution Number
SOL 114
Features Used
Cyclic Static
Reference
1. Fred G. Seely and James O. Smith, Advanced Mechanics of Materials, 2nd ed. (New York: John
Wiley & Sons, 1952) P. 227
Main Index
146 MSC Nastran Verification Guide
Results
Theoretical Solution
For a circular plate with the edges fixed and a uniform load
4
3 2 PR
max = ------ 1 – ---------3- (2.38-1)
16 Et
4
3 2 1.0 11.0
max = ------ 1 – 0.3 ------------------------------------3- (2.38-2)
16 10E+6 0.1
Main Index
CHAPTER 2.38 147
Deflection of a Circular Plate, Cyclic Symmetry
Comparison of Results
z (in)
Theory 0.2500
MSC Nastran 0.2526
Difference 1.03%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
148 MSC Nastran Verification Guide
Problem Description
Find the natural frequencies of a free circular plate.
Solution Number
SOL 115
Features Used
Cyclic Modes
Reference
1. S. Timoshenko, D. H. Young and W. Weaver, Jr., Vibration Problems in Engineering, 4th ed. (New
York: John Wiley & Sons, 1947) P. 501
Main Index
CHAPTER 2.39 149
Normal Modes of a Circular Plate, Cyclic Symmetry
Reference Deck
from Old
Model Name Manual Model Description Location
v11501.dat V4801.dat Normal modes of a circular tpl\verifman\oldverif
plate- cyclic symmetry
Results
Theoretical Solution
Frequency of vibration in rad/sec
D
= ----2- ------ (2.39-1)
a h
3
Eh
D = --------------------------
3
- (2.39-2)
12 1 –
D = 915.75 (2.39-3)
Main Index
150 MSC Nastran Verification Guide
-
---
(258.14) (601.23)
1 9.076 20.52 35.24 52.91
(1893.68) (2937.86)
Comparison of Results
Following table shows frequency (rad/sec) comparison of MSC Nastran and theoretical results.
Main Index
CHAPTER 2.39 151
Normal Modes of a Circular Plate, Cyclic Symmetry
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
152 MSC Nastran Verification Guide
Problem Description
A simply supported beam of constant cross sections is loaded at the ends with an axial force. Determine the
critical buckling load.
Solution Number
SOL 116
Features Used
Cyclic Symmetry Buckling
Reference
1. Stephen P. Timoshenko and James M. Gere, Theory of Elastic Stability, 2nd ed. (New York:
McGrawHill Book Co., 1961) P.49.
Main Index
CHAPTER 2.40 153
Euler Buckling of a Simply Supported Beam, Cyclic Symmetry
Reference Deck
from Old
Model Name Manual Model Description Location
v11601.dat V7701s.dat Euler Buckling of a simply tpl\verifman\oldverif
supported beam
Results
Theoretical Solution
2
EI-
P CR = -----------
2
l
2 6
10.0 10 0.08333 -
P CR = ----------------------------------------------------------
2
100
P CR = 822.467 lb
P applied = 822.467
= 0.99982
P CR = 0.99982 822.467
P CR = 822.319 lb
Main Index
154 MSC Nastran Verification Guide
Comparison of Results
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
CHAPTER 2.41 155
Frequency Response of a Beam-Mass Structure using Dihedral Cyclic Symmetry
Frequency Response of a
2.41 Beam-Mass Structure using
Dihedral Cyclic Symmetry
Problem Description
Find the response at resonance for the model described in v10302, modeling one-half of the model, using
cyclic symmetry analysis.
Solution Number
SOL 118
Features Used
Cyclic Frequency Response
Reference
1. J. P. Den Hartog, Mechanical Vibrations, 4th ed. (New York: McGraw-Hill Book Company, Inc.,
1956) pp. 429-431.
Main Index
156 MSC Nastran Verification Guide
6. The default format for frequency response analysis output is the Cartesian form of the complex
number (“RE_IMAG”). By requesting (“PHASE”) on the DISP case control card a polar form
(“MAG_PHASE”) is output instead. The SORT2 option on the DISP card produces a more compact
output.
The model details are given below:
Engineering Data
Beam properties:
A = 1.0 (Extensional Area)
I1 = 2. 0 (Bending moment of inertia in 1-3 plane)
I2 = 1.0 (Bending Moment of inertia in 1-2 plane)
J = 1.25 (Moment of inertia for Torsion)
E = 1.0 (Young's Modulus)
G = 1.0 (Shear Modulus)
L = 1.0 (Length)
Mass properties
M1 = 1.0 (Axial)
M2, M3 = 64.0
I11 = 1.0
I22, I33 = 16.0
Main Index
CHAPTER 2.41 157
Frequency Response of a Beam-Mass Structure using Dihedral Cyclic Symmetry
Reference Deck
from Old
Model Name Manual Model Description Location
v11801.dat V7801s.dat Frequency response of a beam-mass tpl\verifman\oldverif
structure using Dihedral cyclic symmetry
Results
Theoretical Solution
The response of a single degree of freedom oscillator at resonance is given in equation 2.28a (page 49 of the
reference) as
P0 k
X 0 = n = --------------- (2.41-1)
2c c c
The phase angle between applied force and displacement is plotted on page 51 and shows a 900 value for all
values of damping. The symbols are defined as:
X0 = Magnitude of displacement
P0 = Magnitude of applied load
k = Stiffness of oscillator
c = Viscous damping coefficient
Main Index
158 MSC Nastran Verification Guide
The expected solutions should all have the magnitude of one-half and a phase of 2700.
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
CHAPTER 2.42 159
Stress Wave Propagation
Problem Description
A rod of uniform cross section is fixed at one end and a constant force is suddenly applied to its free end. A
stress wave propagates along the length of the rod. Analyze this wave propagation phenomenon by SOL 129
for transient response analysis.
Solution Number
SOL 129
Features Used
Transient Response
Reference
1. Stephen P. Timoshenko and J. N. Goodier, Theory of Elastic Stability, 2nd ed. (New York:
McGraw-Hill Book Co. 1951), pp. 492-496.
2. R. C. Juvinall, Engineering Consideration of Stress, Strain and Strength, McGraw-Hill, New York,
1967, pp. 185-188.
Main Index
160 MSC Nastran Verification Guide
Results
Theoretical Solution
Transient responses at two points of interest, the free end and the fixed end points are analyzed.
Due to constant force F, applied at the free end, the stress at the free end = – F A is constant at all times.
Before the wave front first reaches the fixed end, the stress at the fixed end is zero. As soon as the wave reaches
the fixed end, it is reflected and produces a stress of 2 at the fixed end. The stress remains unchanged until
the wave front returns to this point. The compressive wave leaves the fixed end and upon reaching the free
end it is reflected as a tensile stress wave. When this tensile wave arrives at the fixed end and is reflected there,
the stress at this point becomes zero. The zero stress is again unchanged for a complete cycle of the stress wave
motion.
The displacement-time history at the free end can be expressed, in general, as:
L
x t
ut = ----------------
E
dx
0
FL
u t = u 2 + ------- t – 2 for 2 t 4
AE
etc.
Main Index
CHAPTER 2.42 161
Stress Wave Propagation
This shows that the displacement for the free end is a linear function of time in each time interval.
Where, n max is the highest natural frequency in cycles/sec of all the interested modes of the corresponding
model. Introduction of small damping in the rod elements is recommended to reduce the high frequency
oscillations in the time history plots.
Comparison of Results
The stress time history plots given in Figure 2.42-1 and Figure 2.42-2 for the free end and fixed end, respectively
show that the MSC Nastran results are consistent with the theoretical ones in their global trends and
characters. The oscillation about the exact solution may be attributed to higher natural frequency modes.
Main Index
162 MSC Nastran Verification Guide
Conclusion
The displacement-time plot for the free end (Figure 2.42-3) Shows a close agreement with the theoretical
solution.
Main Index
CHAPTER 2.43 163
Impact Analysis
Problem Description
A rod with fixed end is struck by a moving mass at the other end. Let v0 be the initial velocity of the mass.
Consider the mass body as absolutely rigid. The velocity of particles at the end of the rod at the instant of
impact (t = 0) is v0. Determine the time history of stress and displacement at various cross sections of the rod
during the impact. Find the duration of the impact and the maximum stress in the rod for various values of
mass ratios.
Solution Number
SOL 129
Features Used
Transient Response
Reference
1. Stephen P. Timoshenko and J. N. Goodier, Theory of Elastic Stability, 2nd ed. (New York:
McGraw-Hill Book Co. 1951), pp. 497-504.
2. April 1986 Application Note, Application Manual Section.
Main Index
164 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
v12902.dat V9902s.dat Impact analysis tpl\verifman\oldverif
Results
MSC Nastran Solution
The rod is modeled with twenty CONRODs and a GAP element. One end of the rod is constrained
completely and the other is left free axially. One CONM2 is used to represent the mass body, and the mass
value is varied for various ratios of M to m. A gap element connects the CONM2 and the free end of the rod
with a zero initial opening. An initial velocity of v0 = -0.1 m/sec is specified for both the CONM2 and the
free end grid point using TIC cards. A time step is chosen so that twenty vibration modes of the rod model
can be adequately represented. A small damping (0.4% at resonance) is introduced in the rod elements to
smooth out the high frequency oscillations. The consistent mass matrix option (PARAM, COUPMASS) is
recommended for more stable responses.
Comparison of Results
The instant when the stress at the free end (or force in the gap element) vanishes signifies the end of the
impact. The durations of the impact for four distinct values of M/m are determined from MSC Nastran
analyses and compared with the calculations of Saint-Venant in the following table.
M/m 1 2 4 6
Theory 3.068 4.708 5.900 7.419
MSC Nastran 3.065 4.708 6.700* 7.420
% Difference 0.1% 0.00% 13.56% 0.01%
*The discrepancy is due to a second contact on the rebound which is not accounted for in the theory as
explained in the references.
Main Index
CHAPTER 2.43 165
Impact Analysis
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
166 MSC Nastran Verification Guide
Problem Description
A two-dimensional airfoil is mounted on bending and torsion springs that are aft of the aerodynamic center
as shown in the figure below. A critical velocity will cause either a static instability (torsional divergence) or
an oscillatory instability (flutter). Find both the divergence and flutter speeds for the system.
Solution Number
SOL 145
Features Used
Flutter
Reference
1. Raymond L. Bisplinghoff, Holt Ashley, and Robert L. Halfman, Aeroelasticity (Reading: Addison-
Wesley Publishing Co., Inc., 1955) pp. 538-541
Main Index
CHAPTER 2.44 167
Flutter of a Two-Dimensional Airfoil
Reference Deck
Model Name from Old Manual Model Description Location
v14501q.dat V7501q.dat Flutter of two dimensional Airfoil tpl\verifman\oldverif
Results
Theoretical Solution
The data correspond to Case (r) on page 538 of the reference. The divergence speed is given there as
U D b = 2.88 (2.44-1)
The flutter speed is found from the graph (r) on page 541 of the reference. For X = 0.1 and frequency ratio
h = 0.4 , we find
U F b = 2.25 (2.44-3)
from which
U F = 169 fps (2.44-4)
Main Index
168 MSC Nastran Verification Guide
The solution is carried out by the k-method using the exact Theodorsen function. The plotted output is
shown in the following figure. Interpolation of the printed output for the flutter speed at zero damping gives
VF = 166 fps; extrapolation of the printed output for the divergence speed at zero damping gives VD = 216
fps.
Comparison of Results
VF VD
Theory 169 216
MSC Nastran 166 216
% Difference 1.8% 0.0%
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
CHAPTER 2.44 169
Flutter of a Two-Dimensional Airfoil
Main Index
170 MSC Nastran Verification Guide
Problem Description
Two Parallel 'grey' planes have emissivities of 0.8 and 0.7 and are maintained at 500oF and 1000oF
respectively. Assuming all energy radiated from each plate is irradiated onto the other (View factors = 1.0)
Calculate the net radiation exchange.
Solution Number
SOL 153
Features Used
Nonlinear steady state Heat transfer
Reference
1. Alan J., Chapman, Heat Transfer, 3rd ed. (New York: Macmillan Publishing Co., Inc., 1974) P. 447,
ex. 11.4
Main Index
CHAPTER 2.45 171
Radiation Exchange Between Two Panels
Reference Deck
from Old
Model Name Manual Model Description Location
v15301s.dat V7401s.dat Radiation exchange between tpl\verifman\oldverif
two panels
Results
Theoretical Solution
The net radiation exchange between two “gray” surfaces is given by,
4 4
T1 – T2
q = ---------------------------------------------------- (2.45-1)
1 1 + 1 2 – 1
Note: This is the exact answer obtained when above values are used in the above equation. Text example
contains some round-off error.
Conclusion
The results from MSC Nastran are comparable to the theoretical results.
Main Index
172 MSC Nastran Verification Guide
Problem Description
A 10-inch-thick slab, initially at a uniform temperature of 180o F, suddenly has its end surfaces reduced to
60o F. Find the temperature at a plane four inches from the end surface after an elapsed time of two hours.
Solution Number
SOL 159
Features Used
Transient Heat Transfer
Reference
1. Alan J., Chapman, Heat Transfer, 3rd ed. (New York: Macmillan Publishing Co., Inc., 1974) P. 126,
ex. 4.2.
Tinitial = 1800F
T ends time = 0
= 600F
Main Index
CHAPTER 2.46 173
Transient Temperature Distribution in a Slab
Results
Theoretical Solution
For the given conditions, from a graph in the reference it can be determined that
T–T
----------------s = 0.685 (2.46-1)
Ti – Ts
Main Index
174 MSC Nastran Verification Guide
Comparison of Results
Temperature
Theory 142.2
MSC Nastran 142.05
% Difference -0.10549%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.47 175
Transient Temperature Distribution in a Slab with Convective Effects
Transient Temperature
2.47 Distribution in a Slab with
Convective Effects
Problem Description
A same slab as in problem v15901s (10 inch thick, initially at 180oF) has a fluid at 60oF placed in contact
with its end surfaces at time equal 0.0. Assuming a convective film coefficient of 4.0 Btu/hr-ft-F, find the
temperature at a depth of one inch after 10 hours have passed.
Solution Number
SOL 159
Features Used
Transient Heat Transfer
Reference
1. Alan J., Chapman, Heat Transfer, 3rd ed. (New York: Macmillan Publishing Co., Inc., 1974) P. 126,
ex. 4.3.
Main Index
176 MSC Nastran Verification Guide
Tinitial = 1800F
T ends time = 0
= 600F
Reference Deck
Model Name from Old Manual Model Description Location
v15902s.dat V8902s.dat Transient temperature distribution in tpl\verifman\oldverif
a slab with convective effects
Results
Theoretical Solution
For the given conditions, it can be graphically determined that
T–T
----------------f = 0.13 (2.47-1)
Ti – Tf
Main Index
CHAPTER 2.47 177
Transient Temperature Distribution in a Slab with Convective Effects
Comparison of Results
Temperature
Theory 75.6
MSC Nastran 71.89
% Difference -4.90741%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
178 MSC Nastran Verification Guide
Problem Description
A square plate with clamped edges is loaded uniformly such that the center deflection is more than the plate
thickness. Considering large deflection (small strain theory) find the deflection of the center of the plate.
Solution Number
SOL 400
Features Used
Geometric Nonlinearity
Reference
1. S. Timoshenko and S. Woinowsky - Krieger, Theory of plates and Shells (New York: McGraw-Hill
Book Company, Inc., 1959), P. 422.
Main Index
CHAPTER 2.48 179
Large Deflection of Uniformly Loaded Plate
Reference Deck
Model Name from Old Manual Model Description Location
V40002.dat V6401.dat Large deflection of uniformly loaded tpl\verifman\oldverif
plate
Results
Theoretical Solution
4
For qb Dh = 109.3 , W max h = 1.20
4
For qb Dh = 218.6 , W max h = 1.66
Main Index
180 MSC Nastran Verification Guide
Comparison of Results
Wmax/h
4
qb Dh 109.3 218.6
Theory 1.2 1.66
MSC Nastran 1.26 1.75
% Difference 5.00% 5.42%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 2.49 181
Large Deflection of a Beam
Problem Description
Find the displaced position of a post-buckled column, fixed at one end with an applied axial load at other
end.
Solution Number
SOL 400
Features Used
Geometric Nonlinearity
Reference
1. Stephen P. Timoshenko and James M. Gere, "Theory of Elastic Stability" (New York: McGraw-Hill
Book Company, Inc., 1961), P. 48.
Main Index
182 MSC Nastran Verification Guide
Reference Deck
from Old
Model Name Manual Model Description Location
V40003.dat V6402s.dat Large deflection of a beam tpl\verifman\oldverif
Results
Theoretical Solution
The critical Euler buckling load is found using equation 2-4 on page 48 of the reference.
2
EI-
P CR = -----------
2
(2.49-1)
4L
2 7 –4
2.0 10 2.8125 10 -
P CR = -------------------------------------------------------------------------
2
(2.49-2)
4 12
P CR = 96.38286 (2.49-3)
For a load greater than P CR , the following tip rotations can be determined:
Main Index
CHAPTER 2.49 183
Large Deflection of a Beam
R3
Subcase Load ID (rad) (deg)
14 200 1.04385 59.808
26 300 1.3919 79.750
Comparison of Results
(deg)
P/PCR 1.152 1.293
Theory 600 800
MSC Nastran 59.808 79.750
% Difference 0.32% 0.31%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
184 MSC Nastran Verification Guide
Main Index
Section 3: Rotordynamics Problems
Main Index
Main Index
3.1 Flexible Rotor Critical Speed
and Mode Shape
Problem Description
Complex eigenvalue analysis for Rotors.
Solution Number
SOL 107
Features Used
The analysis is performed in fixed coordinate systems. For the purpose of comparison, same rotor is
modeled using 1-D elements (ROTORG), axisymmetric harmonic elements (ROTORAX), and
solid elements (ROTOR).
Reference
1. Geradin M. and Kill N., (1984), “A new Approach to Finite Element Modeling of Flexible Rotors.”
Engineering Computations, Vol. 1, 52-64.
Main Index
188 MSC Nastran Verification Guide
Results
Eigenvalue comparison to the paper [1]:
Main Index
CHAPTER 3.1 189
Flexible Rotor Critical Speed and Mode Shape
Main Index
190 MSC Nastran Verification Guide
Main Index
CHAPTER 3.1 191
Flexible Rotor Critical Speed and Mode Shape
Conclusion
The results obtained here compare well with those obtained by Geradin et al. (Geradin and Kill, 1984).
Main Index
192 MSC Nastran Verification Guide
Problem Description
Solve a Jeffcott rotor eigenvalues and plot Campbell Diagram.
Solution Number
SOL 107
Features Used
Massless shaft flexible rotor with lumped mass from disk - Jeffcott Rotor
Rigid supported bearing - Pinned at bearing location
Campbell Diagram - Spin speed vs Natural Frequency
Reference
1. Friswell, M. I., Penny, J. E. T., Garvey, S. D. and Lees, A. W. (2010). Dynamics of rotating
machines. Cambridge University Press, P114~116.
Main Index
CHAPTER 3.2 193
Jeffcott Rotor Eigenvalues
The finite element models used in the analysis is shown in the following figure:
Rotor define:
Main Index
194 MSC Nastran Verification Guide
Results
The Campbell Diagram from MSC Nastran's output is shown as:
Main Index
CHAPTER 3.2 195
Jeffcott Rotor Eigenvalues
The comparison of MSC Nastran results and analytical solution from [1]:
Conclusion
The Jeffcott rotor eigenvalues results from Nastran SOL 107 match the analytical results. The natural
frequencies' maximum error is 0.39% at 0 rpm and 0.51% at 900 rpm.
Main Index
196 MSC Nastran Verification Guide
Problem Description
Solve a Jeffcott rotor with multiple disks for complex eigenvalues.
Solution Number
SOL 107
Features Used
Massless shaft flexible rotor with lumped mass from disk - Jeffcott Rotor
Rigid supported bearing - Pinned at bearing location.
Reference
1. Friswell, M. I., Penny, J. E. T., Garvey, S. D. and Lees, A. W. (2010). Dynamics of rotating
machines. Cambridge University Press, P205-206.
Main Index
CHAPTER 3.3 197
Complex Eigenvalue solution for Jeffcott rotor with multiple disks
The finite element model used in the analysis is shown in the following figure.
Main Index
198 MSC Nastran Verification Guide
Rotor define
Campbell define:
Results
The comparison of MSC Nastran results and analytical solution from [1]:
Speed,
Mode Type rpm Freq, Hz Freq, Hz Freq, Hz Freq, Hz Freq, Hz Freq, Hz
SOL107 0 13.65 13.65 43.26 43.26 113.53 113.53
Analytical [1] 0 13.79 13.79 43.66 43.66 114.08 114.08
Error 1.02% 1.02% 0.92% 0.92% 0.48% 0.48%
SOL107 4000 13.45 13.83 39.75 46.44 95.24 130.63
Analytical [1] 4000 13.59 13.97 40.07 46.90 95.52 131.63
Error 1.03% 1.00% 0.80% 0.80% 0.29% 0.76%
Main Index
CHAPTER 3.3 199
Complex Eigenvalue solution for Jeffcott rotor with multiple disks
Conclusion
The Jeffcott rotor eigenvalues results from Nastran SOL 107 match the analytical results. The natural
frequencies' maximum error is 1.02% at 0 rpm and 1.03% at 4000 rpm.
Main Index
200 MSC Nastran Verification Guide
Complex Eigenvalue
3.4 Solution for Jeffcott Rotor
with Overhung Disk and
Anisotropic Bearings
Problem Description
Solve a Jeffcott rotor with overhung disk supported by anisotropic bearings for complex eigenvalues.
Solution Number
SOL 107
Features Used
Massless shaft flexible rotor with lumped mass from disk - Jeffcott Rotor
Rigid anisotropic bearing - Pinned at bearing location.
Reference
1. Friswell, M. I., Penny, J. E. T., Garvey, S. D. and Lees, A. W. (2010). Dynamics of rotating
machines. Cambridge University Press, P278-279.
Main Index
CHAPTER 3.4 201
Complex Eigenvalue Solution for Jeffcott Rotor with Overhung Disk and Anisotropic Bearings
The finite element model used in the analysis is shown in the following figure.
Several different variants of this system are considered in which the difference line in the bearing (and bearing
support) properties as shown in below table.
Main Index
202 MSC Nastran Verification Guide
Main Index
CHAPTER 3.4 203
Complex Eigenvalue Solution for Jeffcott Rotor with Overhung Disk and Anisotropic Bearings
Rotor define
Campbell define:
Results
The comparison of MSC Nastran results and analytical solution from [1]:
Main Index
204 MSC Nastran Verification Guide
1000 SOL107 2 3 -0.020 ± 42.13j -0.028 ± 46.30j -0.625 ± 391.9j -1.352 ± 406.3j
1000 Analytical [1] 2 3 -0.020 ± 42.26j -0.028 ± 46.40j -0.619 ± 392.0j -1.358 ± 406.6j
2000 SOL107 2 3 -0.016 ± 40.10j -0.033 ± 48.33j -0.201 ± 378.0j -1.595 ± 410.1j
2000 Analytical [1] 2 3 -0.017 ± 40.22j -0.033 ± 48.47j -0.190 ± 377.9j -1.605 ± 410.3j
3000 SOL107 2 3 -0.137 ± 38.12j -0.038 ± 50.37j -0.0001 ± 357.8j -1.774 ± 412.8j
3000 Analytical [1] 2 3 -0.139 ± 38.23j -0.039 ± 50.51j -0.0004 ± 357.3j -1.786 ± 413.0j
Conclusion
The Jeffcott rotor eigenvalues results from Nastran SOL 107 match the analytical results for different bearing
stiffness & damping settings.
Main Index
Section 4: SOL 200 Optimization
Problems
Main Index
Main Index
4.1 Optimization of a Three Bar
Truss
Problem Description
Optimization of a three bar truss.
Solution Number
SOL 200
Features Used
Sizing optimization with design constraints on displacements and stresses in the simple structure
under two static load conditions.
Reference
1. DOT-Design Optimization Tools, User's Manual ,Vanderplaats Research and Development,
Colorado Springs, CO, 1999.
Main Index
208 MSC Nastran Verification Guide
10” 10”
1 2 3
X
3
1
A3 10”
A1
2
Results
The table below shows a comparison between the MSC Nastran result and the result in the reference.
Main Index
CHAPTER 4.1 209
Optimization of a Three Bar Truss
Conclusion
Slight differences are observed between the MSC Nastran result and the Referenced Results, which is typical
for the nonlinear programming tasks involved with optimization. MSC Nastran employs sensitivities to
make the optimization search process more efficient so that a converged result is achieved in 4 design cycles
(5 function evaluations) while the reference uses gradients computed using finite differences, and therefore,
many more function evaluations.
Main Index
210 MSC Nastran Verification Guide
Problem Description
Optimization of a cantilevered structure with a constraint on the first natural frequency.
Solution Number
SOL 200
Features Used
Sizing optimization with a design constraint on a natural frequency.
Reference
1. M.J. Turner "Design of Minimum Mass Structures with Specified Natural Frequencies," AIAA
Journal , Vol 5, No. 3, March 1967.
Main Index
CHAPTER 4.2 211
Design with Constraint on Minimum Frequency
Z
A1 A2 A3
6” t1 t2 t3 X
Results
The referenced analysis is from 1967 and used innovative techniques of the time to attack the design task.
The webs of the structure were taken to be piecewise linear with the design variables the thickness at four
spanwise stations. MSC Nastran considered all the elements to be constant thickness so there are three design
variables for the web and another three for the caps (top and bottom caps are taken to be equal). The table
below averages the web thickness from the two ends to provide a direct comparison.
The table below shows a comparison between the MSC Nastran result and the result in the reference.
Main Index
212 MSC Nastran Verification Guide
Conclusion
It is seen that the two methods produce the same basic design but there are differences between them. There
are significant differences in the modeling techniques used and the algorithms employed, so the strong
qualitative agreement is considered acceptable.
Main Index
CHAPTER 4.3 213
Shape Design of a Culvert
Problem Description
Optimization of a Culvert Structure using Shape Variables to provide a least weight design with design
constraints on element stresses.
Solution Number
SOL 200
Features Used
Shape optimization with constraints on the stress in the elements.
Reference
1. Belegundu, A.D. and Rajan, S.P., "Shape Optimal Design Using Isoparametric Elements,”
Proceedings, 29th AIAA/ASME/ASEE/AMS/ASC Structures, Dynamics and Materials Conference,
pp. 696-701, Williamsburg, VA, April 1988.
Main Index
214 MSC Nastran Verification Guide
Results
The final design, including stress contours, is shown in Figure 4.3-2. An exact comparison with the reference
cannot be made because
a. the number of elements used to mesh the surface differs between the referenced paper and the
current analysis and
b. the reference does not contain data (such as grid locations) that can be compared directly.
Main Index
CHAPTER 4.3 215
Shape Design of a Culvert
However there is good qualitative agreement in that the reference indicates the final design has reduced the
weight by 19% and MSC Nastran shows an 18.9% reduction.
Main Index
216 MSC Nastran Verification Guide
Problem Description
A classic problem of topology optimization that seeks the optimal topology design of a simply supported
beam under a point load.
Solution Number
SOL 200
Features Used
Topology optimization with constraints on the weight of the structure and minimal compliance as
the design objective
Reference
1. Sigmund Ole, Morphology-Based Black and White Filters for Topology Optimization. Structural
and Multidisciplinary Optimization, 33(4-5) 401-424.https://doi.org/10.1007/s-158-006-0087-x.
Main Index
CHAPTER 4.4 217
Topology Design of the MBB Beam
Results
The final design is shown in Figure 4.4-2. The minimum thickness, which applies to the members generated
by the design, is very important in constructing this figure as it prevents the possibility of very small
members. An exact comparison with the reference cannot be made because of differences in dimensions in
the two studies and wide number of variations applied to the design task. Nonetheless, many of the designs
of the reference show a similar final design with the result of Figure 5g in the reference very similar to the
one shown here.
Main Index
218 MSC Nastran Verification Guide
Figure 4.4-2 MBB Beam Topology Design with “Minimum Member Size”
Main Index
Section 5: SOL 700 Explicit Dynamics
Problems
Main Index
Main Index
5.1 Impulsively Loaded Strip
Problem Description
An aluminum strip clamped at its left and right edges is subjected to an impulsive initial velocity (vz = -132
m/sec) over its central portion (one-fourth of the length).
The purpose is to investigate the sensitivity of the response to Poisson's ratio ( =0.01, = 0.3) and to check
results against data available from experiments performed at the Air Freight Flight Dynamics Laboratory.
This problem is also presented in the Demonstration Problem Manual- Explicit Nonlinear (SOL 700).
Solution Number
SOL 700
Features Used
Explicit dynamics.
Elastic-plastic material.
Impact (initial velocity)
Reference
1. Balmer, H. A. and Witmer, E. A. “Theoretical-experimental Correlation of Large Dynamic and
Permanent Deformation of Impulsively Loaded Simple Structure,” 1964, Air Force Flight Dynamics
Laboratory report FDRTDR-64-108.
Main Index
222 MSC Nastran Verification Guide
To better represent the actual experimental conditions, the initial velocity distribution is slightly modified so
as to have a smooth transition to the remainder of the strip. The TIC1 entry is used to input initial velocities.
Using the PARAM, INISTEP entry, the initial time step is set to 0.1e-6 sec according to the COURANT
criterion.
length L = 0.254m
width W = 0.0305m
thickness t = 0.00318m
density 2791kg/m3
Young’s modulus E = 7.17 1010N/m2
Poisson’s ratio = 0.3
yield stress y = 2.854 108N/m2
hardening modulus Eh = 0N/m2
The ENDTIME Case Control entry is used to follow the dynamic behavior in the range up to 0.001 seconds.
The Case Control entries TYPE, SAVE, GPOUT, GRIDS, and TIMES are used to build the z displacement
time history of Grid Point 31 (in the plane of symmetry) by saving results every 0.01 msec. The experimental
results are available for every 0.1 msec.
Main Index
CHAPTER 5.1 223
Impulsively Loaded Strip
Results
Figure 5.1-2 Vertical Z-displacement Time History of the (Free) Corner Grid Point (31) of the Half Strip
Conclusion
As viewed in the above figure, the analytical results agree well with the test results.
Main Index
224 MSC Nastran Verification Guide
Problem Description
The effect of a detonation on the environment can be simulated by assuming that the detonated material can
be idealized by a sphere of hot gas with a homogeneous density and specific internal energy. This approach
is suited for problems in which the processes inside of the explosive material are not to be investigated.
In this problem, the propagation of a blast wave will be simulated starting from the initial shock front radius
R0=0.05 m at the time t = 0 sec until it reaches a radius of R = 10 R0. Both, the gas in the sphere and the
surrounding environment behave as an ideal gas (Gas Constant = 1.4).
Solution Number
SOL 700
Features Used
Transient explicit dynamic analysis.
8-node hex elements for Euler
Two regions (high and low pressure) are defined with transient initial values
Reference
1. Baker, W. E., Explosions in Air, University of Texas Press, 1973, Austin and London.
where l denotes the smallest element dimension, c0 the initial speed of sound, and S is a safety factor Dytran’s
default value is 2/3).
Main Index
CHAPTER 5.2 225
Basic Blast Analysis
With l = 10R 0 20 = 0.025 m and c 0 = – 1 e 1 Equation (5.2-1) yields for the time step
t = 1.96E-7 second. Therefore, the value 1.96E-7 sec is used for the INISTEP parameter.
Figure 5.2-1 Radial Pressure Distribution for Several Time Steps and Theoretical Maximum Pressure
Figure 5.2-2 SOL 700 CTETRA Model - Impact of Cylinder - Quarter Model
Main Index
226 MSC Nastran Verification Guide
doc\SOL700_input_files\Fluid_Dynamics\Basic_Blast
Results
Figure 5.2-3 shows pressure profiles in the elements 1 to 20, which are located along the edge y = z = 0 m of
the control volume. The theoretical value of the shock front pressure is also included in this plot. Though
this mesh is very coarse, the analysis results is a fairly good approximation of the theoretical values. Note that
at t = 0, element 2 is not completely inside of the sphere of high energy gas. Therefore, its pressure is below
that of element 1 at t = 0 sec.
Conclusion
The contour plot of the pressure in Figure 5.2-3 shows an almost spherical shape of the shock front. Deviations
are due to the fact that mass transport takes place along element faces only. Thus, the relief of pressure is
hampered for the elements on the diagonal.
Main Index
CHAPTER 5.3 227
JWL Blast Analysis
Problem Description
This example demonstrates detonation of high explosives.
Solution Number
Transient explicit dynamic analysis (SOL 700)
Features Used
Detonation of High Explosives (HE)
Hydrodynamic Theory of Steady-state Plane Detonation
Steady-state Detonation Modeling with SOL 700
At the start time, all the elements are filled with explosive material. Therefore, they all reference a EOSJWL
equation of state defined by a DMAT entry.
The specific chemical energy q0 is assigned as the initial specific internal energy of the explosive by using the
TICEL entry.
Main Index
228 MSC Nastran Verification Guide
The Chapman-Jouguet detonation velocity Dcj, the ignition point, and the ignition time are specified by the
DETSPH entry. The ignition time is taken as the start time of the analysis and the ignition point is the center
of the left face of the mesh. See Figure 5.3-1.
The duration of the analysis is set to 60 sec necessary to burn the slab (50 cm) at a detonation speed of 7980
m/sec.
Edits of pressure profiles are requested every 10 sec.
Boundary conditions: Two regions (high and low pressure) are defined with transient initial values
Element type: 8-node hex elements for Euler
FE results: Pressure distribution plot
doc\SOL700_input_files\Fluid_Dynamics\JWL_Blast\ep3_3.dat
Results
The solution of the detonation analysis is shown in Figure 5.3-2.
The peak pressure of each profile corresponding to the pressure behind the detonation front is depicted in
Figure 5.3-3 in terms of the Chapman-Jouguet pressure fraction.
From the hydrodynamic theory of steady-state plane detonation for a -law gas, the CJ pressure is
Main Index
CHAPTER 5.3 229
JWL Blast Analysis
P cj = 2 c j – 1 q 0 0
The code needs about 60 elements to build up the detonation front at which time the pressure reaches
approximately 0.85 of Pcj. Subsequently, the front propagates with only a small increase of pressure.
Conclusion
These results are acceptable if the fact is taken into account that SOL 700 is a first-order code that smears
the shock front over a number of elements (always conserving momentum and energy). The result of that
will be a reduction in peak pressure.
Main Index
230 MSC Nastran Verification Guide
Problem Description
The Taylor “bar test” is an important laboratory test in the science of ballistics. It enables to determine an
average value of the dynamic yield stress of a material. It consists of accelerating a cylindrical bar (the velocity
being parallel to the axis of symmetry), and then let the bar hit a rigid target. As a result, the bar shortens and
the impact side expands radially acquiring a mushroom like shape.
Much research work has been done on these impact tests. In this example, the experimental work done by
Johnson/Cook is validated with a number of simulations in SOL 700 with different element types. On top
of this work, the results are compared against the theoretical solution developed by Taylor.
Johnson and Cook developed a material model that represents a constitutive model for materials subjected
to large strains, high strain rates and high temperatures. The SOL 700 implementation is validated for a set
of constitutive constants presented in Reference 1.
Solution Number
SOL 700
Features Used
A list of the functionalities used in this problem.
References
1. Johnson, G.R. and Hook, W.H. “A constitutive model and data for metals subjected to large strains,
high strain rates and high temperatures”, April 1983, 7th Ballistic Symposium, The Hague, The
Netherlands
2. Wilkins, L.M and Guinan M.W. “Impact of cylinders on a rigid boundary”, August 1972, Lawrence
Livermore Laboratory
Main Index
CHAPTER 5.4 231
Taylor Test (a rod impacted against a rigid wall)
Theoretical - Taylor
In Figure 5.4-9, Figure 5.4-10, and Figure 5.4-11 the theoretical reduction in length according to Taylor is
plotted:
Lf /L0 = 0.78
This theoretical reduction is based on an assumption of constant yield of 700 MPa.
In this Taylor test, the process stops when all the kinetic energy has converted to plastic work. It is interesting
to note that the governing parameters of the process give a non-dimensional quantity, therefore, the result
depends on this combination and not on the specific value of each parameter.
Main Index
232 MSC Nastran Verification Guide
Figure 5.4-2 SOL 700 CHEXA Model - Impact of Cylinder - Quarter Model
Figure 5.4-3 SOL 700 CTETRA Model - Impact of Cylinder - Quarter Model
Main Index
CHAPTER 5.4 233
Taylor Test (a rod impacted against a rigid wall)
Results
The deformed shape is compared in Figure 5.4-4 through .
Figure 5.4-4 SOL 700 Simulation - Johnson Cook Material Model with CHEXA - One Point Gauss
Main Index
234 MSC Nastran Verification Guide
Figure 5.4-5 SOL 700 Simulation - von Mises Material Model with CHEXA - One Point Gauss
Figure 5.4-6 SOL 700 Simulation - von Mises Material Model with CHEXA - 2x2x2
Figure 5.4-7 SOL 700 Simulation - von Mises Material Model with CTETRA - Old Element
Main Index
CHAPTER 5.4 235
Taylor Test (a rod impacted against a rigid wall)
Figure 5.4-8 SOL 700 Simulation - von Mises Material Model with CTETRA - New Element
Main Index
236 MSC Nastran Verification Guide
Conclusion
In Figure 5.4-11, a detailed section of Figure 5.4-10 is shown. The results all show an oscillating behavior that is
common in all transient dynamic events and is physical. Furthermore, the new 2x2x2 CHEXA element is
deviating from the experimental result. and further developments are targeted to improve the behavior of the
new CHEXA element in impact events involving metal plasticity. The 2x2x2 CHEXA element has
specifically been implemented and tested to reduce the hour-glassing phenomenon in rubber and foam
modeling.
Main Index
Section 6: NAFEMS Problems
Main Index
Main Index
Linear Elastic
Main Index
242 MSC Nastran Verification Guide
Problem Description
Find the tangential edge stress at D (refer below figure) of elliptical membrane which is uniformly loaded with
outward pressure.
Solution Number
SOL 101
Features Used
Linear Static Analysis.
Reference
1. Test LE1 from NAFEMS Publication NNB, Rev. 3, NAFEMS Linear Benchmarks, 5 Oct 1990.
Loading
Uniform outward pressure of 10MPa at outer edge BC, Inner curved edge AD unloaded.
Boundary Condition
Edge AB, symmetry about Y axis, e.g. zero x displacement
Edge CD, symmetry about X axis, e.g. zero y displacement.
Main Index
CHAPTER 6.1 243
LE1: Elliptic Membrane
Engineering Data
E = 210E+03 MPa
= 0.3
P = 10 Mpa
Main Index
244 MSC Nastran Verification Guide
Results
Theoretical Results
Tangential edge stress (yy) at D: 92.7 MPa
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.2 245
LE2: Cylindrical Shell Patch Test
Problem Description
Find the outer surface tangential (- ) stress at E (refer below figure) of cylindrical shell with two load cases
as follow.
1. Uniform normal edge moment, on DC of 1.0 kNm/m
2. Uniform outward normal pressure, at mid-surface ABCD, of 0.6 MPa and Tangential outward
normal pressure, on edge DC, of 60.0 MPa.
Solution Number
SOL 101
Features Used
Linear Static Analysis.
Reference
1. Test LE2 from NAFEMS report TSBM, Publication NNB, Rev. 3, NAFEMS Linear Benchmarks,
5 Oct 1990.
Main Index
246 MSC Nastran Verification Guide
Loading
Boundary Condition
Edge AB, all translations and rotations zero
Edge AD and edge BC, symmetry about r plane e.g. Z translations and normal rotations all zero.
Engineering Data
E = 210E+03 MPa and = 0.3
Main Index
CHAPTER 6.2 247
LE2: Cylindrical Shell Patch Test
Results
Theoretical Results
Outer Convex surface tangential stress at point E: 60 MPa
MSC Nastran results vs theoretical results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
248 MSC Nastran Verification Guide
Problem Description
Find the displacement in x- direction of point A of hemisphere shown in Figure 6.3-1.
Solution Number
SOL 101
Features Used
Linear Static Analysis.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE3. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Uniform normal pressure of 1 MPa on the upper surface of the plate.
Boundary Condition
Edge AB, all translations and rotations zero
Edge AD and edge BC, symmetry about r plane, e.g. Z translations and normal rotations all zero.
Engineering Data
E = 68.25E+03 MPa and = 0.3
Main Index
CHAPTER 6.3 249
LE3: Hemisphere-Point Loads
Results
Theoretical Results
X displacement at point A: 0.185 m
MSC Nastran results vs theoretical results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
250 MSC Nastran Verification Guide
Problem Description
A z-section cantilever beam is subjected to a torsion load. Find the axial stress (X-X) at mid-surface at Point
A shown in Figure 6.4-1.
Solution Number
SOL 101
Features Used
Linear Static Analysis.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE5. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Torque of 1.2 MNm applied at end x = 10 by two uniformly distributed edge shears, S = 0.6 at each flange.
Engineering Data
E = 210E+03 MPa
= 0.3
Boundary Condition
Edge x = 0, all displacements are zero.
Main Index
CHAPTER 6.4 251
LE5: Z-section Cantilever
Results
Theoretical Results
Axial (x-x) stress at mid-surface, point A: -108 MPa (Compression)
MSC Nastran results vs theoretical results are as follows:
Main Index
252 MSC Nastran Verification Guide
Figure 6.4-2 Axial Stress Result in X- direction of Z-section Cantilever Beam by Nastran
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.5 253
LE6: Skew Plate Normal Pressure
Problem Description
A skew plate is subjected uniform normal pressure in the vertical z- direction. Find the maximum principal
at plate center.
Solution Number
SOL 400
Features Used
Linear Static Analysis.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE6. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Normal pressure of -0.7KPa in vertical Z-direction.
Engineering Data
E = 210E+03 MPa and = 0.3
Boundary Condition
Simple supports (no Z-displacement) for all edges AB, BC, CD, DA.
Main Index
254 MSC Nastran Verification Guide
Results
Theoretical Results
Maximum principal stress on the lower surface at the plate center: 0.802 MPa
MSC Nastran results vs theoretical results are as follows:
Conclusion
MSC Nastran result Maximum principal stress at Center of plate for decks with different elements
CQUAD4, CQUAD8, CQUAD8 with advance and CQUADR are comparable with Theoretical results.
Main Index
CHAPTER 6.6 255
LE7: Axisymmetric Cylinder/Sphere-Pressure
Problem Description
Axisymmetric cylinder-sphere shell is subjected uniform internal pressure. Find the axial stress yy on the
outer surface of the upper cylinder at D.
Solution Number
SOL 400
Features Used
Linear Static Axisymmetric.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE7. Glasgow: NAFEMS, Rev. 3, 1990.
Boundary Condition
Point A, zero radial displacement and zero rotation. Point F, zero y-displacement.
Loading
Uniform internal pressure of 1 MPa.
Engineering Data
E = 210E+03 MPa and = 0.3
Main Index
256 MSC Nastran Verification Guide
Results
Theoretical Results
Axial (yy) stress on the outer surface of the upper cylinder, at point D: 25.86 MPa
Nastran result is calculated using interpolation formula.
MSC Nastran results vs theoretical results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.7 257
LE8: Axisymmetric Shell Pressure
Problem Description
Axisymmetric shell is subjected uniform internal pressure. Find the hoop stress on the outer surface of the
shell at D.
Solution Number
SOL 400
Features Used
Linear Static Axisymmetric.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE8. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Uniform internal pressure of 1 MPa.
Boundary Condition
Point A, zero y displacement. Point E, zero radial displacement and zero rotation.
Main Index
258 MSC Nastran Verification Guide
Results
Theoretical Solution
Hoop stress (tt) stress on the outer surface of shell, at point D: 94.55 MPa
Nastran result is calculated using extrapolation formula:
83.97 + 110.14 -
-----------------------------------------
Hoop Stress at node 33: = 97.055 MPa
2
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.8 259
LE9: Axisymmetric Branched Shell-Pressure
Problem Description
Axisymmetric branched shell is subjected uniform internal pressure over edge BCD shown below in figure.
Find the axial stress yy on the outer surface of the upper cylinder at C.
Solution Number
SOL 400
Features Used
Linear Static Axisymmetric.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE9. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Uniform internal pressure of 1 MPa over edge BCD
Boundary Condition
Point A, fully fixed
Main Index
260 MSC Nastran Verification Guide
Results
Theoretical Solution
Axial (yy) stress on the outer surface of the upper cylinder, at point C: -319.9 MPa
Nastran result is calculated using extrapolation formula as follows:
0.0234375
Axial Stress at node 33: 139.515 + ------------------------- 248.565 – 139.515 = 303.09 MPa
0.015625
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.9 261
LE10: Thick Plate Pressure
Problem Description
A elliptical thick plate is subjected uniform normal pressure on the upper surface of the plate. Find the direct
stress yy at Point D shown in Figure 6.9-1.
Solution Number
SOL 101
Features Used
Linear Static Analysis.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE10. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Uniform normal pressure of 1 MPa on the upper surface of the Plate
Boundary Condition
Face DCD'C' zero y-displacement
Face ABA'B' zero x- displacement
Face BCB'C' x and y displacements fixed, z displacements fixed along mid-plane
Main Index
262 MSC Nastran Verification Guide
Results
Theoretical Solution
Normal (yy) stress at point D: 5.38 MPa
Result comparison of theoretical solution and MSC Nastran results are as follows:
Main Index
CHAPTER 6.9 263
LE10: Thick Plate Pressure
Figure 6.9-2 a) Model with CHEXA element with 8 node b) Model with CHEXA element with 20 node
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
264 MSC Nastran Verification Guide
Problem Description
A solid spherical taper cylinder is subjected to temperature loading. Find the direct stress zz at Point A shown
in below figure.
Solution Number
SOL 101
Features Used
Linear Static Analysis.
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
LE11. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Boundary Condition
Symmetry on x-z plane i.e., zero Y-displacement
Symmetry on y-z plane i.e., zero X- displacement
Face on x-y plane zero Z- displacement
Face HIH’I’ zero Z- displacement
Main Index
CHAPTER 6.10 265
LE11: Solid Cylinder/ Taper/ Sphere-Temperature
Results
Theoretical Solution
Direct stress (zz) at point A: -105 MPa
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
266 MSC Nastran Verification Guide
Main Index
CHAPTER 6.11 267
T1: Membrane with Hot Spot
Problem Description
Find the direct stress in y- direction at point D in membrane as shown below figure which is subjected to
thermal load (hot-spot) at center.
Solution Number
SOL 400
Features Used
Heat Transfer
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
T1. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Within hot-spot ( 0 r 1.0 mm )
Boundary Condition
On quarter model,
at y=0, symmetric about x, e.g. y displacement zero
Main Index
268 MSC Nastran Verification Guide
Main Index
CHAPTER 6.11 269
T1: Membrane with Hot Spot
Results
Theoretical Solution
Direct stress at Point D in Y direction: 50 MPa
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
270 MSC Nastran Verification Guide
Heat Transfer
Main Index
CHAPTER 6.12 271
T2: One Dimensional Heat Transfer with Radiation
Problem Description
Find the material temperature at end of the bar when it is subjected to thermal load at other end.
Solution Number
SOL 400
Features Used
Heat Transfer
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
T2. Glasgow: NAFEMS, Rev. 3, 1990.
At boundary B
Emissivity = 0.98
Stefan-Boltzman Constant = 5.67 x 10-8 Wm2 0K-4
Loading
Zero internal heat generation
Boundary Condition
At A, prescribed temperature of 1000 0K
Main Index
272 MSC Nastran Verification Guide
Results
Theoretical Solution
Material temperature at B: 927 oK
Result comparison of theoretical solution and MSC Nastran results are as follows:
Temperature at B in K
Theory MSC Nastran % Error
927 912.7 -1.54%
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.13 273
T3: One Dimensional Transient Heat Transfer
Problem Description
Find the Material Temperature at point C of the bar when it is subjected to transient thermal load at B. (Refer
Figure 6.13-1)
Solution Number
SOL 400
Features Used
Transient Heat Conduction
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
T3. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Zero internal heat generation
Boundary Condition
At time t = 0, All temperature = 0 0C
At time t > 0, All temperature = 0 0C
At B, temperature = nt- 0C
100 sin -----
40
Main Index
274 MSC Nastran Verification Guide
Results
Theoretical Solution
Material temperature at C: 36.60 oC
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.14 275
T4: Two Dimensional Heat Transfer with Convection
Problem Description
Natural Convection heat transfer occurs in one edge insulted plate which is subjected to thermal loading at
one edge. Find the temperature at point E. (Refer Figure 6.14-1)
Solution Number
SOL 400
Features Used
Heat Conduction
Reference
1. NAFEMS Finite Element Methods & Standards, The Standard NAFEMS Benchmarks, Test No.
T4. Glasgow: NAFEMS, Rev. 3, 1990.
Loading
Zero internal heat generation
Boundary Condition
Edge AB, temperature = 100 0C
Edge DA, Zero heat flux
Edges BC, CD, convection to ambient temperature = 0 0C
Main Index
276 MSC Nastran Verification Guide
Results
Theoretical Solution
Material temperature at E: 18.3 oC
Result comparison of theoretical solution and MSC Nastran results are as follows:
Conclusion
MSC Nastran results are comparable with theoretical results.
Main Index
CHAPTER 6.14 277
T4: Two Dimensional Heat Transfer with Convection
Free Vibration
Main Index
278 MSC Nastran Verification Guide
Problem Description
NAFEMS Benchmark FV2 - pin-ended double cross normal modes
Solution Number
SOL 103
Features Used
Normal modes
Beam and bar elements (Euler and considering shear flexibility)
Coupling between flexural and extensional behavior
Repeated eigenvalues
Reference
1. Problem FV2 from The Standard NAFEMS Benchmarks, Oct 5, 1990.
Main Index
CHAPTER 6.15 279
FV2- Pin-Ended Double Cross: In-Plane Vibration
Boundary Conditions:
• Out-of plane motion is constrained at all nodes.
• Nodes A through H are constrained for in-plane translation.
Material Properties
E = 200 109N/m2
8000kg/m3
Main Index
280 MSC Nastran Verification Guide
Results
Natural frequencies (Hz):
Conclusion
Although not specifically mentioned, the NAFEMS values are obviously based on Euler beam theory and
ignore the shear flexibility. If Euler beams are used, the results compare well with the NAFEMS expected
values. If shear flexibility is included, then the natural frequencies are slightly lower.
Main Index
CHAPTER 6.16 281
FV4 - Cantilever with Off-Centre Point Masses
Problem Description
Calculation of modes of a cantilever beam with offset masses at the free end.
Solution Number
SOL 103 - normal modes
Features Used
BAR elements (Euler beam and inclusion of shear flexibility)
3-noded Beam elements
CONM2 with offset
Normal modes
Coupling between bending and torsion
Close eigenvalues
Reference
1. NAFEMS Finite Element Methods & Standards. Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 4. Glasgow: NAFEMS, Nov.,
1987.
Main Index
282 MSC Nastran Verification Guide
Results
Natural frequencies (Hz) of the first 6 modes:
Conclusion
While not an exact match, the calculated natural frequencies agree well with the expected values.
Main Index
CHAPTER 6.17 283
FV5 - Deep Simply-supported Beam
Problem Description
Calculation of modes of a simply supported beam
Solution Number
SOL 103 - normal modes
Features Used
BEAM elements (Timoshenko beam formulation)
2- and 3-noded Beam elements
Normal modes - Lanczos method
Shear deformation and torsion
Reference
1. NAFEMS Finite Element Methods & Standards. Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 5. Glasgow: NAFEMS, Nov.,
1987.
Main Index
284 MSC Nastran Verification Guide
Results
Natural frequencies (Hz) of the first 9 modes
2-Noded Beams
NAFEMS using lumped 2-Noded Beams using
Mode Value mass Coupled Mass 3-Noded Euler Beams
1 and 2 42.649 43.111 43.226 38.635
3 77.542 77.204 77.842 77.524
4 125.0 124.487 125.515 125.002
5 and 6 148.31 149.393 154.519 172.888
7 233.1 224.056 241.242 232.811
8 and 9 284.55 269.578 305.392 253.929
Conclusion
While not an exact match, the calculated natural frequencies agree well with the expected values. While the
bending modes for the 3-noded beam are different than those of the 2-noded beams, they agree well with the
modes one obtains if the models using 2-noded beams are refined to have 10 elements.
Main Index
CHAPTER 6.18 285
FV12- Free Thin Square Plate
Problem Description
Out of plane free vibration of a square plate with in-plane motion constrained.
Solution Number
SOL 103 - normal modes
Features Used
Rigid-body modes (3)
Repeated eigenvalues
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 12. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 8 × 8
Main Index
286 MSC Nastran Verification Guide
Boundary Conditions
x- and y-translation and z-rotation are constrained to 0.0
Material Properties
Young’s modulus, E = 200 109N/m2
Poisson's ratio, 0.3
Results
Natural frequencies (Hz) of modes 4 through 10 (modes 1-3 are rigid-body modes).
Main Index
CHAPTER 6.18 287
FV12- Free Thin Square Plate
QUAD4 results
QUAD8 results
Conclusion
In all but one case (QUAD8, coupled mass, mode 5) the results are closer to the theoretical values than the
expected results. It should be noted that a 45 degree skew angle is well beyond the recommended maximum
skew of 30 degrees, but the elements still performed well.
Main Index
288 MSC Nastran Verification Guide
Problem Description
Solution to find the first 6 modes of a clamped flat plate, which is skewed 45 degrees.
Solution Number
SOL 103 normal modes
Features Used
Lanczos eigenvalue solution
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 15. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 12 × 12
Main Index
CHAPTER 6.19 289
FV15- Clamped Thin Rhombic Plate
Boundary Conditions
Clamped edges
Main Index
290 MSC Nastran Verification Guide
Results
QUAD4 results
QUAD8 results
Conclusion
In all but one case (QUAD8, coupled mass, mode 5) the results are closer to the theoretical values than the
expected results. It should be noted that a 45 degree skew angle is well beyond the recommended maximum
skew of 30 degrees, but the elements still performed well.
Main Index
CHAPTER 6.20 291
FV16- Cantilevered Thin Plate
Problem Description
Normal modes of a cantilevered thin plate modeled using QUAD8 elements.
Solution Number
SOL 103 - normal modes using Lanczos method.
Features Used
Normal modes calculation using QUAD8 elements, distorted in-plane.
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 16. Glasgow: NAFEMS, Nov.,
1987.
Main Index
292 MSC Nastran Verification Guide
Coordinates (m)
Node No. x y
1 4.0 4.0
2 2.25 2.25
3 4.75 2.5
4 7.25 2.75
5 7.5 4.75
6 7.75 7.25
7 5.25 7.25
8 2.25 7.25
9 2.5 4.75
Coordinates (m)
Node No. x y
1 4.0 4.0
Main Index
CHAPTER 6.20 293
FV16- Cantilevered Thin Plate
Main Index
294 MSC Nastran Verification Guide
Results
Test 1 results:
Test 2 results:
Main Index
CHAPTER 6.20 295
FV16- Cantilevered Thin Plate
Test 3 results:
Test 4 results:
Conclusion
Considering the number of elements in the models, the program gave excellent results. With square elements,
the results show excellent correlation, considering the mesh and the results for distorted elements are
acceptable, considering the mesh and element distortions. Using the recommendation of 5 GRID points per
half sine-wave of the deformed shape, Test 1 and 2 might be expected to predict mode 1 and possibly mode
2 well and Test 3 and 4 might be expected to predict the first mode well.
Main Index
296 MSC Nastran Verification Guide
Problem Description
Solution to find the first 6 modes of a clamped flat plate, which is skewed 45 degrees.
Solution Number
SOL 103 - normal modes
Features Used
Lanczos Eigenvalue Solution
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 22. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 10 × 10
Boundary conditions
Clamped edges.
Main Index
CHAPTER 6.21 297
FV22 - Clamped Thick Rhombic Plate
Results
QUAD4
QUAD8
Conclusion
Distortion of elements to 45 degrees skew angle is well beyond the recommended maximum value of 30
degrees. However, the elements still performed well in this example. As the mesh is not that refined, the
results are expected to be less accurate after the first mode or two.
Main Index
298 MSC Nastran Verification Guide
Problem Description
Solution to find the first 6 modes of a tapered membrane plate.
Solution Number
SOL 103-normal modes
Features Used
Lanczos eigenvalue solution
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 32. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 16 × 8
Main Index
CHAPTER 6.22 299
FV32 - Cantilevered Tapered Membrane
Boundary Conditions
Clamped edges
Main Index
300 MSC Nastran Verification Guide
Results
QUAD4 results
QUAD8 results
Conclusion
Distortion of elements is a way to test their accuracy. In this case, the results are acceptably close to the
predicted results.
Main Index
CHAPTER 6.23 301
FV41 - Axisymmetric Vibration
Problem Description
Free modes of a cylinder using axisymmetric elements.
Solution Number
SOL 103 - normal modes
Features Used
Axisymmetric elements
Normal modes - Lanczos method
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 41. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 16 × 3
Main Index
302 MSC Nastran Verification Guide
Results
Main Index
CHAPTER 6.23 303
FV41 - Axisymmetric Vibration
4-noded CQUADX
4-noded QUADX
NAFEMS NAFEMS Expected Using Lumped 4-noded QUADX
Mode Theoretical Value Value Mass Using Coupled Mass
2 243.53 244.01 243.23 244.01
3 377.41 379.41 374.47 379.42
4 394.11 395.41 389.47 395.44
5 397.72 397.94 389.49 401.38
6 405.28 421.87 396.10 421.89
8-noded CQUADX
8-noded QUADX
NAFEMS NAFEMS Expected Using Lumped 8-noded QUADX
Mode Theoretical Value Value Mass Using Coupled Mass
2 243.53 243.50 243.24 243.57
3 377.41 377.46 356.49 377.46
4 394.11 394.28 356.88 394.30
5 397.72 397.94 375.85 397.97
6 405.28 406.41 393.65 406.44
Conclusion
While the results using lumped mass are not bad, the results using coupled mass fall within 1% of the
expected values for all modes.
Main Index
304 MSC Nastran Verification Guide
Problem Description
Radial modes of a hollow sphere using axisymmetric elements.
Solution Number
SOL 103 - normal modes
Features Used
Axisymmetric Elements
Cylindrical Coordinate System
MPC Equations
Normal Modes - Lanczos Method
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C., Selected Benchmarks for Natural Frequency Analysis, Test No. 42. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 10 × 1
= 5°
Constraints:
z' displacements are constrained to 0.0 and GRIDs at the same radius are constrained to have the same radial
displacement.
Main Index
CHAPTER 6.24 305
FV42 - Thick Hollow Sphere
Results
Natural frequencies (hz) of the first 5 modes
4-noded CQUADX
4-noded QUADX
NAFEMS NAFEMS Expected Using Lumped 4-noded QUADX
Mode Theoretical Value Value Mass Using Coupled Mass
1 369.91 370.64 370.47 370.64
2 838.03 841.2 833.49 841.20
3 1451.2 1473.1 1423.77 1473.06
4 2117.0 2192.2 2033.83 2192.20
5 2795.8 2975.7 2050.78 2975.67
8-noded CQUADX
Main Index
306 MSC Nastran Verification Guide
Conclusion
While the results using lumped mass are not bad, the results using coupled mass on the 4-noded elements
matches the baseline values within a fraction of a percent. The 8-noded elements match the first 3 modes
almost exactly, but fall off on the last two modes, indicating a need for a more refined mesh.
Main Index
CHAPTER 6.25 307
FV52 - Simply Supported Solid Square Plate
Problem Description
Calculate the first 10 modes of a plate which is supported in the Z-direction on its edges.
Solution Number
SOL 103 normal modes
Features Used
Solid elements
Normal modes - Lanczos method
Reference
1. NAFEMS Finite Element Methods & Standards. Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 51. Glasgow: NAFEMS, Nov.,
1987.
L.O.E. 8 × 8 × 3
Main Index
308 MSC Nastran Verification Guide
Constraints
z' displacements are constrained to 0.0 along the bottom edges
Results
Natural frequencies (Hz) of modes 4-10. Modes 1-3 are rigid-body modes.
8-noded HEXA
20 noded HEXA
Conclusion
In most cases, the values from MSC Nastran were closer to the theoretical values than the NAFEMS expected
values. As expected with a coarse mesh model, the results using coupled mass were closer to the theoretical
values than those using lumped mass.
Main Index
CHAPTER 6.26 309
FV73 - Cantilevered Thin Square Plate
Problem Description
Calculate the first 6 modes of a plate which is simply supported along the Y-axis.
Solution Number
SOL 103 normal modes
Features Used
8-noded QUAD elements (NAFEMS lists 8-noded semi-loof elements, which do not exist in MSC
Nastran)
Normal modes - Lanczos method
Reference
1. NAFEMS Finite Element Methods & Standards. Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 73. Glasgow: NAFEMS, Nov.,
1987.
Main Index
310 MSC Nastran Verification Guide
Constraints
X, Y, and Z translation and Y rotation are constrained to 0.0 on the GRID points along the Y-axis
Results
Natural frequencies (Hz) of modes 1-6.
Main Index
CHAPTER 6.26 311
FV73 - Cantilevered Thin Square Plate
NAFEMS Values
Lumped Mass
Coupled Mass
Conclusion
While the results using lumped mass are good, the results using coupled mass is closer to the reference
solution than the NAFEMS expected values in all cases.
Main Index
312 MSC Nastran Verification Guide
Forced Vibration
Main Index
CHAPTER 6.27 313
Test 5 - Deep Simply Supported Beam
Problem Description
Normal modes of a simply supported beam model.
Solution Number
SOL 103 - normal modes using Lanczos method
Features Used
Normal modes calculation using 3-noded BEAM3 elements
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 5. Glasgow: NAFEMS, Nov.,
1987.
Main Index
314 MSC Nastran Verification Guide
Results
NAFEMS
Theoretical NAFEMS Iso BEAM3
Mode Value Value Results
1 and 2 42.650 42.657 43.185
3 71.20 71.202 71.209
4 125.00 125.00 125.002
5 and 6 148.15 148.71 154.424
7 213.61 213.72 213.848
8 and 9 283.47 286.91 307.523
Conclusion
The normal recommendation is to have 8 elements per sine wave of the deformed shape. With only 5
elements in the model, the first mode is the only mode which can be expected to be accurate. Modes 1
through 4 are the first bending, axial, and torsional modes and are predicted well using this model. The other
modes are not predicted as well, but mesh refinement would improve the accuracy of these modes.
Main Index
CHAPTER 6.28 315
Test 5H - Deep Simply Supported Beam Harmonic Forced Response
Problem Description
Forced harmonic vibration response of a simply supported beam model.
Solution Number
SOL 108 and 111 - Direct and Modal Frequency response analysis
Features Used
Frequency response calculation using 3-noded BEAM3 elements
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 5H. Glasgow: NAFEMS, Nov.,
1987.
Main Index
316 MSC Nastran Verification Guide
Iso-parametric: 5 elements
Results
Conclusion
The maximum response occurs at the first natural frequency of the system, which for the 3-noded beam
model is 43.185hz. The peak values are within engineering accuracy of the predicted values.
Main Index
CHAPTER 6.29 317
Test 5P - Deep Simply Supported Beam Periodic Forced Vibration Response
Problem Description
Forced periodic vibration response of a simply supported beam model.
Solution Number
SOL 108 and 111 - Direct and Modal Frequency response analysis
Features Used
Frequency response calculation using 3-noded BEAM3 elements
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 5P. Glasgow: NAFEMS, Nov.,
1987.
Main Index
318 MSC Nastran Verification Guide
Although this appears to be a transient solution, it is really a frequency response solution, using Fourier
transforms and inverse Fourier transforms.
Results
Main Index
CHAPTER 6.29 319
Test 5P - Deep Simply Supported Beam Periodic Forced Vibration Response
Conclusion
As residual vectors were used, the modal and direct methods gave the same results (to the number of
significant digits shown), which were within engineering accuracy of the reference solution.
Main Index
320 MSC Nastran Verification Guide
Problem Description
Forced transient vibration response of a simply supported beam model.
Solution Number
SOL 109 and 112 - Direct and Modal transient response analysis
Features Used
Transient response calculation using 3-noded BEAM3 elements
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 5T. Glasgow: NAFEMS, Nov.,
1987.
Main Index
CHAPTER 6.30 321
Test 5T-Deep Simply Supported Beam Transient Forced Vibration Response
F0 = 106 N/m
Results
Main Index
322 MSC Nastran Verification Guide
Conclusion
As residual vectors were used, the modal and direct methods gave the same results (to the number of
significant digits shown), which were within engineering accuracy of the reference solution.
Main Index
CHAPTER 6.31 323
Test 5R-Deep Simply Supported Beam Random Forced Response
Problem Description
Forced random vibration response of a simply supported beam model.
Solution Number
SOL 108 and 111 - Direct and Modal Frequency response analysis
Features Used
Frequency response calculation using BEAM elements
Random response analysis
Reference
1. NAFEMS Finite Element Methods & Standards, Abbassian, F., Dawswell, D. J., and Knowles, N.
C.Selected Benchmarks for Natural Frequency Analysis, Test No. 5R. Glasgow: NAFEMS, Nov.,
1987.
Main Index
324 MSC Nastran Verification Guide
Iso-parametric: 5 elements
Main Index
CHAPTER 6.31 325
Test 5R-Deep Simply Supported Beam Random Forced Response
Results
Main Index
326 MSC Nastran Verification Guide
Conclusion
The maximum response occurs at the first natural frequency of the system, which for this model is 43.166hz.
The peak values are within engineering accuracy of the predicted values.
Main Index
CHAPTER 6.32 327
Test 13: Simply Supported Thin Square Plate: Frequency Extraction
Problem Description
Modal analysis of a thin square plate for extraction of frequencies.
Solution Number
SOL 103
Features Used
Modal analysis using QUAD4 and QUAD8 shell elements
Reference
1. Test 13 from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Steady state harmonic forcing function F = F0 sint
F0 = 100 N/m2 over whole plate with = 2f and f = 4.16 Hz.
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
328 MSC Nastran Verification Guide
Main Index
CHAPTER 6.32 329
Test 13: Simply Supported Thin Square Plate: Frequency Extraction
Results
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
330 MSC Nastran Verification Guide
Problem Description
Forced harmonic vibration response of a thin square plate with modal and direct solutions.
Solution Number
SOL 108 and SOL 111
Features Used
Harmonic response calculation using QUAD4 and QUAD8 shell elements
Reference
1. Test 13H from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Steady state harmonic forcing function F = F0 sint
F0 = 100 N/m2 over whole plate with = 2f and f = 4.16 Hz.
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
CHAPTER 6.33 331
Test 13H: Simply Supported Thin Square Plate: Harmonic Forced Vibration
Main Index
332 MSC Nastran Verification Guide
Results
Frequency f (Hz)
MSC Nastran
NAFEMS Values % Error
QUAD4 (Modal) 2.377 2.3318 -1.90%
Main Index
CHAPTER 6.33 333
Test 13H: Simply Supported Thin Square Plate: Harmonic Forced Vibration
Frequency f (Hz)
MSC Nastran
NAFEMS Values % Error
QUAD4 (Direct) 2.377 2.3318 -1.90%
QUAD8 (Modal) 2.377 2.3757 -0.05%
QUAD8 (Direct) 2.377 2.3757 -0.05%
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
334 MSC Nastran Verification Guide
Problem Description
Forced periodic vibration response of a thin square plate with modal and direct solutions.
Solution Number
SOL 108 and SOL 111
Features Used
Periodic response calculation using QUAD4 and QUAD8 shell elements
Reference
1. Test 13P from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Steady state periodic forcing function F = F0 (sint - sin3t)
F0 = 100 N/m2 over whole plate with = 2f and f = 1.2 Hz.
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
CHAPTER 6.34 335
Test 13P: Simply Supported Thin Square Plate: Periodic Forced Vibration
Main Index
336 MSC Nastran Verification Guide
Results
Main Index
CHAPTER 6.34 337
Test 13P: Simply Supported Thin Square Plate: Periodic Forced Vibration
Conclusion
The displacement and stress results from the direct solution with QUAD4 elements are shown in Figure 6.34-2
and Figure 6.34-3. The results from MSC Nastran are compared with the NAFEMS results in the above table.
Main Index
338 MSC Nastran Verification Guide
Problem Description
Forced transient vibration response of a thin square plate with modal and direct solutions.
Solution Number
SOL 109 and SOL 112
Features Used
Transient response calculation using QUAD4 and QUAD8 shell elements.
Reference
1. Test 13T from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Suddenly applied step load F0 = 100 N/m2 over whole plate.
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
CHAPTER 6.35 339
Test 13T: Simply Supported Thin Square Plate: Transient Forced Vibration
Main Index
340 MSC Nastran Verification Guide
Main Index
CHAPTER 6.35 341
Test 13T: Simply Supported Thin Square Plate: Transient Forced Vibration
Results
Conclusion
The displacement and stress results from the direct solution with QUAD4 elements are shown in Figure 6.35-2
and Figure 6.35-3. The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
342 MSC Nastran Verification Guide
Problem Description
Forced random vibration response of a thin square plate with modal and direct solutions.
Solution Number
SOL 108 and SOL 111
Features Used
Random response calculation using QUAD4 and QUAD8 shell elements
Reference
1. Test 13R from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Random forcing with uniform power spectral density (of force) PSD = (100 N/m2)2 / Hz over whole plate.
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
CHAPTER 6.36 343
Test 13R: Simply Supported Thin Square Plate: Random Forced Vibration
Main Index
344 MSC Nastran Verification Guide
Results
Frequency f (Hz)
MSC Nastran
NAFEMS Values % Error
QUAD4 (Modal) 2.377 2.3318 -1.90%
QUAD4 (Direct) 2.377 2.3318 -1.90%
QUAD8 (Modal) 2.377 2.3757 -0.05%
QUAD8 (Direct) 2.377 2.3757 -0.05%
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
CHAPTER 6.37 345
Test 21: Simply Supported Thick Square Plate: Frequency Extraction
Problem Description
Modal analysis of a thick square plate for extraction of frequencies.
Solution Number
SOL 103
Features Used
Modal analysis using QUAD4 and QUAD8 shell elements.
Reference
1. Test 21 from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
346 MSC Nastran Verification Guide
Main Index
CHAPTER 6.37 347
Test 21: Simply Supported Thick Square Plate: Frequency Extraction
Results
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above table.
Main Index
348 MSC Nastran Verification Guide
Problem Description
Forced harmonic vibration response of a thick square plate with modal and direct solutions.
Solution Number
SOL 108 and SOL 111
Features Used
Harmonic response calculation using QUAD4 and QUAD8 shell elements.
Reference
1. Test 21H from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Steady state harmonic forcing function F = F0 sint
F0 = 106 N/m2 over whole plate with =2f and f = 0 to 78.17 Hz.
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
CHAPTER 6.38 349
Test 21H: Simply Supported Thick Square Plate: Harmonic Forced Vibration
Main Index
350 MSC Nastran Verification Guide
Results
Frequency f (Hz)
MSC Nastran
NAFEMS Values % Error
QUAD4 (Modal) 45.9 45.2339 -1.47%
QUAD4 (Direct) 45.9 45.2339 -1.47%
QUAD8 (Modal) 45.9 46.1462 0.54%
QUAD8 (Direct) 45.9 46.1462 0.54%
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
CHAPTER 6.39 351
Test 21P: Simply Supported Thick Square Plate: Periodic Forced Vibration
Problem Description
Forced periodic vibration response of a thick square plate with modal and direct solutions.
Solution Number
SOL 108 and SOL 111
Features Used
Harmonic response calculation using QUAD4 and QUAD8 shell elements.
Reference
1. Test 21P from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2 = 0.3 = 8000 kg/m3
Loading
Steady state harmonic forcing function F = F0 (sint - sin3t)
Main Index
352 MSC Nastran Verification Guide
Boundary Conditions
Ux = Uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m.
Main Index
CHAPTER 6.39 353
Test 21P: Simply Supported Thick Square Plate: Periodic Forced Vibration
Results
Main Index
354 MSC Nastran Verification Guide
Results
Conclusion
The displacement and stress results from the direct solution with QUAD4 elements are shown in Figure 6.39-2
and Figure 6.39-3. The results from MSC Nastran are compared with the NAFEMS results in the above table.
Main Index
CHAPTER 6.40 355
Test 21T: Simply Supported Thick Square Plate: Transient Forced Vibration
Problem Description
Forced transient vibration response of a thick square plate with modal and direct solutions.
Solution Number
SOL 109 and 112
Features Used
Random transient calculation using QUAD4 and QUAD8 shell elements
Reference
1. Test 21T from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Material Properties
E = 200 x 109 N/m2
= 0.3
= 8000 kg/m3
Loading
Suddenly applied step load F0 = 106 N/m2 over whole plate
Boundary Condition
ux = uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Main Index
356 MSC Nastran Verification Guide
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m
Main Index
CHAPTER 6.40 357
Test 21T: Simply Supported Thick Square Plate: Transient Forced Vibration
Results
Main Index
358 MSC Nastran Verification Guide
Main Index
CHAPTER 6.40 359
Test 21T: Simply Supported Thick Square Plate: Transient Forced Vibration
Conclusion
The displacement and stress results from the direct solution with QUAD4 elements are shown in Figure 6.40-2
and Figure 6.40-3. The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
360 MSC Nastran Verification Guide
Problem Description
Forced random vibration response of a thick square plate with modal and direct solutions.
Solution Number
SOL 108 and 111
Features Used
Random response calculation using QUAD4 and QUAD8 shell elements
Reference
1. Test 21R from NAFEMS Selected Benchmarks for Forced Vibration, R0016, March 1993.
Loading
Random forcing with uniform power spectral density (of force) PSD = (106 N/m2)2 / Hz over the whole
plate.
Boundary Condition
ux = uy = Rz = 0.0 at all nodes and Uz = 0 along all 4 edges
Main Index
CHAPTER 6.41 361
Test 21R: Simply Supported Thick Square Plate: Random Forced Vibration
Rx = 0.0 along edges x=0 and x=10m and Ry = 0.0 along edges y=0 and y=10m
Main Index
362 MSC Nastran Verification Guide
Results
Frequency f (Hz)
MSC Nastran
NAFEMS Values % Error
QUAD4 (Modal) 45.9 45.2239 -1.47%
QUAD4 (Direct) 45.9 45.2239 -1.47%
QUAD8 (Modal) 45.9 46.1462 0.54%
QUAD8 (Direct) 45.9 46.1462 0.54%
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above tables.
Main Index
CHAPTER 6.41 363
Test 21R: Simply Supported Thick Square Plate: Random Forced Vibration
Nonlinear
Main Index
364 MSC Nastran Verification Guide
Problem Description
This is a NAFEMS nonlinear benchmark to demonstrate a geometric nonlinearity solution procedure test of
snap-back under displacement control. It is modeled using the modified Riks-Ramm arc length procedure
available in MSC Nastran.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Modified Riks/Ramm arc length method.
Reference
1. Test NL4 from NAFEMS Publication NNB, Rev. 1, NAFEMS Non-Linear Benchmarks, October
1989 from a previous report Benchmark Tests for Solution Procedures for Geometric Non-Linearity
by Crisfield, Duxbury & Hunt. NAFEMS Report SPGNL, October 1987.
Material Properties
AE = 5.0x107, L = 2500, L = 25, K1 = 1.5, K2 = AE/L(1 + 2)1/2 = 19999.0, K3 = 0.25, K4 = 1.0
Loading
Incremental load P=10000.0 is applied to node A in the X-direction using arc length method.
Boundary Condition
The displacements uy = 0.0 at node A and B, along with the displacement ux = 0.0 at node C.
Main Index
CHAPTER 6.42 365
NL4: Snap-back Under Displacement Control
Main Index
366 MSC Nastran Verification Guide
Results
Main Index
CHAPTER 6.42 367
NL4: Snap-back Under Displacement Control
Main Index
368 MSC Nastran Verification Guide
Conclusion
The results from MSC Nastran are compared with NAFEMS results at loads P= 649.9, 1300.0, 1949.0,
2599.0, 3243.0 and -1099.0 in the above figures. The largest deviation from the NAFEMS solution is less
than 1% at all these reference points.
Main Index
CHAPTER 6.43 369
NL6: Straight Cantilever with Axial End Point Load
Problem Description
A cantilever beam with an axial compressive load combines bending and membrane deformation with
bifurcation of initially straight elements using the modified Riks-Ramm arc length procedure available in
MSC Nastran.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Modified Riks/Ramm arc length method.
Reference
1. Test NL6 from NAFEMS Publication NNB, Rev. 1, NAFEMS Non-Linear Benchmarks, October
1989 as well as a previous report Finite Element Benchmarks for 2D beams and Axisymmetric shells
involving Geometric Non-Linearity by P Lyons and S Holsgrove NAFEMS Report FEBNLGBAS,
March 1989.
Material Properties
The material is elastic with a Young's modulus of E = 210 x 109 N/m2 and a Poisson's ratio of 0.0.
Loading
The loading at end A is applied in increments up to a maximum value of PL2/EI = 22.493.
Boundary Condition
The displacements ux = uy = 0.0 at node 1 and ux = 0.0 at nodes 2 and 3 of the end B.
Main Index
370 MSC Nastran Verification Guide
Figure 6.43-1 Straight Cantilever Beam with Axial End Point Load
Results
Main Index
CHAPTER 6.43 371
NL6: Straight Cantilever with Axial End Point Load
Conclusion
The results from MSC Nastran are matching well with NAFEMS results as shown in the above figures.
Main Index
372 MSC Nastran Verification Guide
Problem Description
This example demonstrates the ability of MSC Nastran to model a beam column structure composed of
initially straight beams that will suddenly snap through and back. Modified Riks/Ramm & Crisfield arc
length methods are used to trace the solution path to this unstable problem.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Modified Riks/Ramm and Crisfield arc length methods.
Reference
1. Finite Element Benchmarks for 2D beams and Axisymmetric shells involving Geometric Non-
Linearity by P Lyons and S Holsgrove NAFEMS Report FEBNLGBAS, March 1989.
Material Properties
The material is elastic with a Young's modulus of E = 71.74 x 109 N/m2 and a Poisson's ratio of 0.0.
Loading
A concentrated load P at point C with an ultimate value of 50000 N is applied incrementally.
Boundary Condition
The beams are pinned at points B and D.
Main Index
CHAPTER 6.44 373
NL7: Lee's Frame Buckling Problem
Main Index
374 MSC Nastran Verification Guide
Results
Main Index
CHAPTER 6.44 375
NL7: Lee's Frame Buckling Problem
Conclusion
The results from MSC Nastran are matching well with NAFEMS results as shown in the above figures for
both modified Riks/Ramm and Crisfield methods.
Main Index
376 MSC Nastran Verification Guide
Fracture Mechanics
Main Index
CHAPTER 6.45 377
Test 1.1: Center Cracked Plate in Tension
Problem Description
The following test case describe finite element analysis of a centre cracked plate tension specimen.
Normalized stress intensity factor for the test case is calculated using virtual crack extension method. The
results are compared to the analytical solution [1].
Solution Number
Sol 400
Features Used
Implicit Nonlinear
Reference
1. Test 1.1 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus, E = 20700 and Poisson's ratio = 0.3
Geometrical properties
a/b = 0.5; h/b = 1.0; b = 20
Loading
Uniform stress of 100 unit is applied on top face as shown in Figure 6.45-1.
Main Index
378 MSC Nastran Verification Guide
Boundary Condition
Symmetry boundary condition as shown in Figure 6.45-1.
uy = 0 on x axis
ux = 0 on y axis
Results
K 0.5
Reference solution [1] of the normalized stress intensity factor is ------I = 1.325 where K o = a
Ko
Main Index
CHAPTER 6.45 379
Test 1.1: Center Cracked Plate in Tension
0.5
K o = a = 560.5
K
------I = 1.3364
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
380 MSC Nastran Verification Guide
Problem Description
The following test case describe finite element analysis of a center cracked plate subjected to quadratic
temperature distribution. Normalized stress intensity factor for the test case is calculated using virtual crack
extension method. The results are compared to the analytical solution [1].
Solution Number
Sol 400
Features Used
Implicit Nonlinear
Reference
1. Test 1.2 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus E = 20700 and a Poisson's ratio = 0.3 and coefficient
of thermal expansion =1.35 × 10-5
Geometrical properties
a/b = 0.1; h/b = 2.5; b = 100
Loading
Quadratic thermal distribution
2
T = 0.01 x
Main Index
CHAPTER 6.46 381
Test 1.2: Center Cracked Plate with Thermal Load
Boundary Condition
Symmetry boundary condition as shown in Figure 6.45-1.
uy = 0 on x axis
ux = 0 on y axis
Main Index
382 MSC Nastran Verification Guide
Results
K 2 2
------I = 1.00 ET b----- a----- 0.5
Reference solution [1] of the normalized stress intensity factor is where K o = -----------
2 3
– a
Ko c 2
K0 = (a)0.5 =514.3
K
------I = 1.008
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
CHAPTER 6.47 383
Test 2.1: Single Edge Cracked Plate in Tension
Problem Description
The following test case describe finite element analysis of a Single Edge Cracked Plate subjected to uniform
stress. Normalized stress intensity factor for the test case is calculated using virtual crack extension method.
The results are compared to the analytical solution [1].
Solution Number
Sol 400
Features Used
Implicit Nonlinear
Reference
1. Test case 2.1 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus of E = 20700 and a Poisson's ratio = 0.3
Geometrical properties
a/b = 0.5; h/b = 0.5; b = 20
Loading
Uniform stress of 100 unit is applied on top face as shown in Figure 6.47-1.
Main Index
384 MSC Nastran Verification Guide
Boundary Condition
Symmetry boundary condition as shown in Figure 6.47-1.
uy = 0 on x axis
ux = 0 at point A
Results
K 0.5
Reference solution [1] of the normalized stress intensity factor is ------I = 3.00 where K o = a
Ko
0.5
K o = a = 560.5
Main Index
CHAPTER 6.47 385
Test 2.1: Single Edge Cracked Plate in Tension
K
------I = 3.00
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
386 MSC Nastran Verification Guide
Problem Description
The following test case describe finite element analysis of an angle crack plate subjected to uniform stress.
Normalized stress intensity factor for the test case is calculated using virtual crack extension method. The
results are compared to the analytical solution [1].
Solution Number
Sol 400
Features Used
Implicit Nonlinear
Reference
1. Test case 3.1 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus of E = 20700 and a Poisson's ratio = 0.3
Geometrical properties
a/b = 0.5; h/b = 1.25; b = 50
Loading
Uniform stress of 100 unit is applied on top face as shown in Figure 6.48-1.
Boundary Condition
Boundary condition is shown in Figure 6.48-1.
Main Index
CHAPTER 6.48 387
Test 3: Angle Crack Embedded in a Plate
uy = 0 on x axis at y = 0
ux = 0 at point A
Results
K K
Reference solution [1] of the normalized stress intensity factor is ------I = 0.190 ; ------II- = 0.405 where
Ko Ko
0.5
K o = a
Main Index
388 MSC Nastran Verification Guide
K II = RK I = 354.21
0.5
K o = a = 886.2
K K II
------I = 0.186 , ------- = 0.400
Ko Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
CHAPTER 6.49 389
Test 4: Cracks at a Hole in a Plate
Problem Description
The following test case describe finite element analysis of cracks at a hole in a plate and subjected to uniform
stress. Normalized stress intensity factor for the test case is calculated using virtual crack extension method.
The results are compared to the analytical solution [1].
Solution Number
Sol 400
Features Used
Implicit Nonlinear
Reference
1. Test case 4.1 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus E = 207000 and a Poisson's ratio = 0.3
Geometrical properties
a/b = 0.3; R/b = 0.25; b = 10
Loading
Uniform stress of 100 unit is applied on top face as shown in Figure 6.45-1.
Main Index
390 MSC Nastran Verification Guide
Boundary Condition
Boundary condition as shown in Figure 6.45-1.
uy = 0 on x axis
ux = 0 on y axis
Main Index
CHAPTER 6.49 391
Test 4: Cracks at a Hole in a Plate
Results
K 0.5
Reference solution [1] of the normalized stress intensity factor is ------I = 1.050 where K o = a
Ko
0.5
K o = a = 307.0
K
------I = 1.068
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
392 MSC Nastran Verification Guide
Problem Description
The following test case describe finite element analysis of axisymmetric crack in a bar and subjected to
uniform stress. Normalized stress intensity factor for the test case is calculated using virtual crack extension
method. The results are compared to the analytical solution [1].
Solution Number
Sol 400
Features Used
Implicit Nonlinear
Reference
1. Test case 5 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus, E = 207000 and Poisson's ratio, = 0.3
Geometrical properties
b/R = 0.5; h/R = 1.5; R = 20
Loading
Uniform stress of 100 unit is applied on top face as shown in Figure 6.45-1.
Main Index
CHAPTER 6.50 393
Test 5: Axisymmetric Crack in a Bar
Boundary Condition
Boundary condition as shown in Figure 6.45-1.
uy = 0 on x axis
ux = 0 on y axis
Main Index
394 MSC Nastran Verification Guide
Results
K P 0.5
Reference solution [1] of the normalized stress intensity factor is ------I = 0.475 where K o = --------2- b
Ko b
P 0.5
K o = --------2- b = 2242.0
b
K
------I = 0.4842
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
CHAPTER 6.51 395
Test 6: Compact Tension Specimen
Problem Description
The following test case describe finite element analysis of compact tension specimen subjected to a point
load. Normalized stress intensity factor for the test case is calculated using virtual crack extension method.
The results are compared to the analytical solution [1].
Solution Number
SOL 400
Features Used
Implicit Nonlinear
Reference
1. Test case 6 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus E = 207000 and a Poisson's ratio = 0.3
Geometrical properties
a/W=0.5; W= 50.0; t = 1.0.
Loading
Point load of 100 unit is applied on top face as shown in Figure 6.45-1.
Main Index
396 MSC Nastran Verification Guide
Boundary Condition
Boundary condition as shown in Figure 6.45-1.
uy = 0 on x axis
ux = 0 at point A.
Results
K P-
Reference solution [1] of the normalized stress intensity factor is ------I = 9.659 where K o = -----------
Ko tw
0.5
P-
K o = -----------
0.5
= 141.42
tw
Main Index
CHAPTER 6.51 397
Test 6: Compact Tension Specimen
K
------I = 9.671
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
398 MSC Nastran Verification Guide
Problem Description
The following test case describe finite element analysis of T joint subjected to uniform stress at the top edge.
Normalized stress intensity factor for the test case is calculated using virtual crack extension method. The
results are compared to the analytical solution [1].
Solution Number
SOL 400
Features Used
Implicit Nonlinear
Reference
1. Test 7.1 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus, E = 207000 and Poisson's ratio, = 0.3
Geometrical properties
a/T=0.1; W/T= 0.5; t/T=1.0; L/T = 12.0; =45°; T=50.
Loading
Point load of 100 unit is applied on top face as shown in Figure 6.45-1.
Boundary Condition
Boundary condition as shown in Figure 6.45-1.
Main Index
CHAPTER 6.52 399
Test 7: T-Joint Weld Attachment
uy = 0 on base
ux = 0 on base and top.
Main Index
400 MSC Nastran Verification Guide
Results
K 0.5
Reference solution [1] of the normalized stress intensity factor is ------I = 1.317 where K o = b
Ko
0.5
K o = b = 396.3
K
------I = 1.335
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
CHAPTER 6.53 401
Test 8: V-Notch Specimen in Tension
Problem Description
The following test case describe finite element analysis of V-Notch cracked plate subjected to uniform stress
at the top edge. Normalized stress intensity factor for the test case is calculated using virtual crack extension
method. The results are compared to the analytical solution [1].
Solution Number
SOL 400
Features Used
Implicit Nonlinear
Reference
1. Test case 8.1 from NAFEMS Publication, Ref: R0020, 2D Test Cases in Linear Elastic Fracture
Mechanics, Jan,1992.
Material Properties
The material is linear elastic with a Young's modulus, E = 207000 and Poisson's ratio, = 0.3
Geometrical properties
a/d=0.2; d/W= 0.1; h/W=1.0; =45°; W=250.
Loading
Uniform stress of 100 unit is applied on top face as shown in Figure 6.45-1.
Main Index
402 MSC Nastran Verification Guide
Boundary Condition
Boundary condition as shown in Figure 6.45-1.
uy = 0 on X axis.
ux = 0 on Y axis.
Results
K 0.5
Reference solution [1] of the normalized stress intensity factor is ------I = 2.740 where K o = a
Ko
Main Index
CHAPTER 6.53 403
Test 8: V-Notch Specimen in Tension
0.5
K o = a = 396.3
K
------I = 2.81
Ko
Conclusion
The results from MSC Nastran are in a good agreement with the reference solution.
Main Index
404 MSC Nastran Verification Guide
Creep
Main Index
CHAPTER 6.54 405
Test 1A: 2-D Plane Stress - Uniaxial Load, Secondary Creep
Problem Description
This NAFEMS creep test involves secondary creep of a two dimensional thin plate subjected to uniaxial load.
Solution Number
SOL 400
Features Used
Secondary Creep Analysis with Plane Stress Elements
Reference
1. Test 1A from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n t A = 3.125 × 10-14 per hour ( in N/mm2) n=5
Loading
Prescribed tensile stress 1 = 200 N/mm2 on line BC.
Boundary Condition
Ux = 0.0 at line AD
Uy = 0.0 on midpoint of line AD.
Main Index
406 MSC Nastran Verification Guide
Figure 6.54-1 Plane Stress Model under Uniaxial Load and Secondary Creep
Main Index
CHAPTER 6.54 407
Test 1A: 2-D Plane Stress - Uniaxial Load, Secondary Creep
Results
Conclusion
C C
The reference solutions from NAFEMS are defined as and YY = ZZ = – 0.005 t .
C C
X X = e ff = 0.01t
Main Index
408 MSC Nastran Verification Guide
The X and Y components of creep strain values from MSC Nastran are compared to these reference results
from NAFEMS in Figure 6.54-2 and also in the following table:
C
XX C
YY
Time NAFEMS MSC Nastran % Error NAFEMS MSC Nastran % Error
0.1 0.001 0.001 0.0% -0.0005 -0.0005 0.0%
1.0 0.01 0.01 0.0% -0.005 -0.005 0.0%
5.0 0.05 0.05 0.0% -0.025 -0.025 0.0%
10.0 0.1 0.1 0.0% -0.05 -0.05 0.0%
50.0 0.5 0.5 0.0% -0.25 -0.25 0.0%
100.0 1.0 1.0 0.0% -0.50 -0.50 0.0%
500.0 5.0 5.0 0.0% -2.5 -2.5 0.0%
1000.0 10.0 10.0 0.0% -5.0 -5.0 0.0%
Main Index
CHAPTER 6.55 409
Test 1B: 2-D Plane Stress - Uniaxial Displacement, Secondary Creep
Problem Description
This NAFEMS creep test involves secondary creep of a two dimensional thin plate subjected to uniaxial
displacement.
Solution Number
SOL 400
Features Used
Secondary Creep Analysis with Plane Stress Elements
Reference
1. Test 1B from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n t A = 3.125 × 10-14 per hour ( in N/mm2) n=5
Loading
None.
Boundary Condition
Ux = 0.0 at line AD
Main Index
410 MSC Nastran Verification Guide
Figure 6.55-1 Plane Stress Model under Uniaxial Displacement and Secondary Creep
Main Index
CHAPTER 6.55 411
Test 1B: 2-D Plane Stress - Uniaxial Displacement, Secondary Creep
Results
Conclusion
The X component of stress values from MSC Nastran are compared to the reference results from NAFEMS
in Figure 6.54-2 and also in the following table.
Main Index
412 MSC Nastran Verification Guide
Main Index
CHAPTER 6.56 413
Test 2A: 2-D Plane Stress - Biaxial Load, Secondary Creep
Problem Description
This NAFEMS creep test involves secondary creep of a two dimensional thin plate subjected to biaxial load.
Solution Number
SOL 400
Features Used
Secondary Creep Analysis with Plane Stress Elements
Reference
1. Test 2A from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n t A = 3.125 × 10-14 per hour ( in N/mm2) n=5
Loading
Prescribed tensile stress 1 = 200 N/mm2 on line BC and 2 = 200 N/mm2 on line CD
Boundary Condition
Ux = 0.0 at line AD
Uy = 0.0 on midpoint of line AB
Main Index
414 MSC Nastran Verification Guide
Figure 6.56-1 Plane Stress Model under Biaxial Load and Secondary Creep
Results
Main Index
CHAPTER 6.56 415
Test 2A: 2-D Plane Stress - Biaxial Load, Secondary Creep
Conclusion
C C C
The reference solutions from NAFEMS are defined as X X = e ff = 0.005t and e ff = 0.01t .
The X and Y components of creep strain values from MSC Nastran are compared to these reference results
from NAFEMS in Figure 6.54-2 and also in the following table:
C
XX C
e ff
Time NAFEMS MSC Nastran % Error NAFEMS MSC Nastran % Error
0.1 0.0005 0.0005 0.0% 0.001 0.001 0.0%
1.0 0.005 0.005 0.0% 0.01 0.01 0.0%
5.0 0.025 0.025 0.0% 0.05 0.05 0.0%
10.0 0.05 0.05 0.0% 0.10 0.10 0.0%
50.0 0.25 0.25 0.0% 0.50 0.50 0.0%
100.0 0.50 0.50 0.0% 1.00 1.00 0.0%
500.0 2.50 2.50 0.0% 5.00 5.00 0.0%
1000.0 5.00 5.00 0.0% 10.00 10.00 0.0%
Main Index
416 MSC Nastran Verification Guide
Problem Description
This NAFEMS creep test involves secondary creep of a two dimensional thin plate subjected to biaxial
displacement.
Solution Number
SOL 400
Features Used
Secondary Creep Analysis with Plane Stress Elements
Reference
1. Test 2B from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n t A = 3.125 × 10-14 per hour ( in N/mm2) n=5
Loading
None.
Boundary Condition
Ux = 0.0 at line AD
Uy = 0.0 on midpoint of line AD
Main Index
CHAPTER 6.57 417
Test 2B: 2-D Plane Stress - Biaxial Displacement, Secondary Creep
Figure 6.57-1 Plane Stress Model under Biaxial Displacement and Secondary Creep
Main Index
418 MSC Nastran Verification Guide
Results
Conclusion
The X component of stress values from MSC Nastran are compared to the reference results from NAFEMS
in Figure 6.54-2 and also in the following table.
Main Index
CHAPTER 6.57 419
Test 2B: 2-D Plane Stress - Biaxial Displacement, Secondary Creep
Main Index
420 MSC Nastran Verification Guide
Problem Description
This NAFEMS creep test involves secondary creep of a two dimensional thin plate subjected to biaxial
displacement.
Solution Number
SOL 400
Features Used
Secondary Creep Analysis with Plane Stress Elements
Reference
1. Test 3B from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n t A = 3.125 × 10-14 per hour ( in N/mm2) n=5
Loading
None.
Boundary Condition
Ux = 0.0 at line AD
Uy = 0.0 on midpoint of line AB
Main Index
CHAPTER 6.58 421
Test 3B: 2-D Plane Stress - Biaxial (-ve) Displacement, Secondary Creep
Figure 6.58-1 Plane Stress Model under Biaxial (Negative) Displacement and Secondary Creep
Main Index
422 MSC Nastran Verification Guide
Results
Conclusion
The X component of stress values from MSC Nastran are compared to the reference results from NAFEMS
in Figure 6.54-2 and also in the following table.
XX
Main Index
CHAPTER 6.58 423
Test 3B: 2-D Plane Stress - Biaxial (-ve) Displacement, Secondary Creep
XX
Main Index
424 MSC Nastran Verification Guide
Problem Description
This NAFEMS creep test involves primary creep of a two dimensional thin plate subjected to uniaxial load.
Solution Number
SOL 400
Features Used
Primary Creep Analysis with Plane Stress Elements
Reference
1. Test 8A from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n tm A = 3.125 × 10-14 per hour ( in N/mm2) n=5 m = 0.5
Loading
Prescribed tensile stress 1 = 200 N/mm2 on line BC.
Boundary Condition
Ux = 0.0 at line AD
Uy = 0.0 on midpoint of line AD.
Main Index
CHAPTER 6.59 425
Test 8A: 2D Plane Stress - Uniaxial Load, Primary Creep
Figure 6.59-1 Plane Stress Model under Uniaxial Load and Primary Creep
Results
Main Index
426 MSC Nastran Verification Guide
Conclusion
C C 0.5
The reference solutions from NAFEMS are defined as X X = e ff = 0.01t and
C C 0.5
YY = ZZ = – 0.005 t .
The X & Y components of creep strain values from MSC Nastran are compared to these reference results
from NAFEMS in Figure 6.54-2 and also in the following table:
C
XX C
YY
Time NAFEMS MSC Nastran % Error NAFEMS MSC Nastran % Error
0.1 0.00316 0.003278 3.74% -0.00158 -0.001639 3.74%
1.0 0.01 0.010141 1.41% -0.005 -0.005071 1.42%
5.0 0.0224 0.022530 0.58% -0.0112 -0.011265 0.58%
10.0 0.0316 0.031806 0.65% -0.0158 -0.015903 0.65%
50.0 0.0707 0.070914 0.30% -0.0354 -0.035458 0.16%
100.0 0.1 0.100224 0.22% -0.05 -0.050112 0.22%
500.0 0.2236 0.223871 0.12% -0.112 -0.111935 -0.06%
1000.0 0.3162 0.316511 0.10% -0.158 -0.158251 0.16%
Main Index
CHAPTER 6.60 427
Test 8B: 2-D Plane Stress - Uniaxial Displacement, Primary Creep
Problem Description
This NAFEMS creep test involve primary creep of a two dimensional thin plate subjected to uniaxial
displacement.
Solution Number
SOL 400
Features Used
Primary Creep Analysis with Plane Stress Elements
Reference
1. Test 8B from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n tm A = 3.125 × 10-14 per hour ( in N/mm2) n=5 m = 0.5
Loading
None.
Boundary Condition
Ux = 0.0 at line AD and Uy = 0.0 on midpoint of line AD
Main Index
428 MSC Nastran Verification Guide
Figure 6.60-1 Plane Stress Model under Uniaxial Displacement and Primary Creep
Results
Main Index
CHAPTER 6.60 429
Test 8B: 2-D Plane Stress - Uniaxial Displacement, Primary Creep
Conclusion
The X component of stress values from MSC Nastran are compared to the reference results from NAFEMS
in Figure 6.54-2 and also in the following table:
Main Index
430 MSC Nastran Verification Guide
Problem Description
This NAFEMS creep test involve primary creep of a two dimensional thin plate subjected to stepped load.
Solution Number
SOL 400
Features Used
Primary Creep Analysis with Plane Stress Elements
Reference
1. Test 8C from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n tm A = 3.125 × 10-14 per hour ( in N/mm2) n=5 m = 0.5
Loading
Prescribed tensile stress 1 on line BC
1 = 200 N/mm2 for t = 0 to 100 hoursand 1 = 250 N/mm2 for t > 100 hours
Boundary Condition
Ux = 0.0 at line AD
Main Index
CHAPTER 6.61 431
Test 8C: 2D Plane Stress - Stepped Load Primary Creep
Figure 6.61-1 Plane Stress Model under Stepped Load and Primary Creep
Main Index
432 MSC Nastran Verification Guide
Results
Conclusion
The X component of creep strain values from MSC Nastran are compared to the reference results from
NAFEMS in Figure 6.54-2 and also in the following table:
Main Index
CHAPTER 6.62 433
Test 9A: 2-D Plane Stress - Biaxial Load, Primary Creep
Problem Description
This NAFEMS creep test involve primary creep of a two dimensional thin plate subjected to biaxial load.
Solution Number
SOL 400
Features Used
Primary Creep Analysis with Plane Stress Elements
Reference
1. Test 9A from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n tm A = 3.125 × 10-14 per hour ( in N/mm2) n=5 m = 0.5
Loading
Prescribed tensile stress 1 = 200 N/mm2 on line BC and 2 = 200 N/mm2 on line CD.
Boundary Condition
Ux = 0.0 at line AD
Uy = 0.0 at line AB
Main Index
434 MSC Nastran Verification Guide
Figure 6.62-1 Plane Stress Model under Stepped Load and Primary Creep
Main Index
CHAPTER 6.62 435
Test 9A: 2-D Plane Stress - Biaxial Load, Primary Creep
Results
Conclusion
C C 0.5 C 0.5
The reference solutions from NAFEMS are defined as X X = YY = 0.005t and e ff = 0.01t .
The X component and effective creep strain values from MSC Nastran are compared to these reference results
from NAFEMS in Figure 6.54-2 and also in the following table:
C
XX C
e ff
Time NAFEMS MSC Nastran % Error NAFEMS MSC Nastran % Error
0.1 0.00158 0.001642 3.93% 0.00316 0.003280 3.80%
1.0 0.005 0.005076 1.53% 0.01 0.010200 2.00%
5.0 0.01112 0.011273 1.38% 0.0224 0.022500 0.45%
10.0 0.0158 0.015912 0.71% 0.0316 0.031800 0.63%
50.0 0.0354 0.035467 0.19% 0.0707 0.071325 0.88%
100.0 0.05 0.050121 0.24% 0.1 0.100273 0.27%
500.0 0.1118 0.111942 0.13% 0.2236 0.223990 0.17%
1000.0 0.1581 0.158268 0.11% 0.3162 0.317000 0.25%
Main Index
436 MSC Nastran Verification Guide
Problem Description
This NAFEMS creep test involve primary creep of a two dimensional thin plate subjected to biaxial
displacement.
Solution Number
SOL 400
Features Used
Primary Creep Analysis with Plane Stress Elements
Reference
1. Test 9B from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n tm A = 3.125 × 10-14 per hour ( in N/mm2) n=5 m = 0.5
Loading
None.
Boundary Condition
Ux = 0.0 at line AD and Uy = 0.0 at line AB
Biaxial displacement: u1 = 0.1 mm on line BC and u2 = 0.1 mm on line CD.
Main Index
CHAPTER 6.63 437
Test 9B: 2-D Plane Stress - Biaxial Displacement, Primary Creep
Figure 6.63-1 Plane Stress Model under Biaxial Displacement and Primary Creep
Main Index
438 MSC Nastran Verification Guide
Results
Conclusion
The X component of stress values from MSC Nastran are compared to the reference results from NAFEMS
in Figure 6.63-2 and also in the following table:
Main Index
CHAPTER 6.64 439
Test 9C: 2-D Plane Stress - Biaxial Stepped Load, Primary Creep
Problem Description
This NAFEMS creep test involve primary creep of a two dimensional thin plate subjected to biaxial stepped
load.
Solution Number
SOL 400
Features Used
Primary Creep Analysis with Plane Stress Elements
Reference
1. Test 9C from NAFEMS Publication Ref: R0027, NAFEMS Fundamental Tests of Creep Behavior,
June 1993.
Material Properties
Young’s Modulus E = 200 × 103 N/mm2 Poisson's ratio = 0.3
Creep Law: = A n tm A = 3.125 × 10-14 per hour ( in N/mm2) n=5 m = 0.5
Loading
Prescribed tensile stress 1 on line BC and 2 on line CD.
Main Index
440 MSC Nastran Verification Guide
Boundary Condition
Ux = 0.0 at line AD and Uy = 0.0 at line AB
Figure 6.64-1 Plane Stress Model under Biaxial Stepped Load and Primary Creep
Main Index
CHAPTER 6.64 441
Test 9C: 2-D Plane Stress - Biaxial Stepped Load, Primary Creep
Results
Conclusion
The X component of creep strain values from MSC Nastran are compared to the reference results from
NAFEMS in Figure 6.63-2 and also in the following table:
Main Index
442 MSC Nastran Verification Guide
Composite
Main Index
CHAPTER 6.65 443
R0031(1): Laminated Strip under Three-point Bending
Problem Description
This NAFEMS composite benchmark is a simply supported 7 layer symmetric strip with a central line load.
A quarter model, with appropriate boundary conditions, is used and predictions for stress and displacement
are obtained using composite shell elements of Nastran.
Solution Number
Sol 101
Features Used
Linear Static Analysis
Layered Composite Elements
Reference
1. Test R0031/1 from NAFEMS publication R0031, Composites Benchmarks, February 1995.
Material Properties
E1 = 100GPa 12 = 0.4 G12 = 3GPa
E2 = 5GPa 13 = 0.3 G13 = 2GPa
E3 = 5GPa 23 = 0.3 G23 = 2GPa
Loading
Line load of 10 N/mm at point C.
Boundary Condition
Uz = 0.0 at point A
Main Index
444 MSC Nastran Verification Guide
Results
MSC Nastran
NAFEMS Values % Error
11 at E (MPa) 683.9 682.7 -0.18%
13 at D (MPa) -4.1 -4.0908 -0.22%
uz at E (mm) -1.06 -1.0625 0.24%
Conclusion
The results from MSC Nastran are compared with the NAFEMS results in the above table. The value of 11
at point E are obtained by extrapolating stress values at element centroid of two elements adjacent to point
E. More variants of this model with other composite element types are demonstrated with SOL400 in
Chapter 7: Laminated Strip under Three-point Bending of MSC Nastran Implicit Nonlinear Demonstration Guide.
Main Index
CHAPTER 6.66 445
R0031(2): Wrapped Thick Cylinder under Pressure and Thermal Loading
Problem Description
This NAFEMS composite benchmark is a long thick cylinder made from isotropic material onto which
external hoop windings of orthotropic material have been added. This problem demonstrates the ability to
model pressure and thermal loading for composite laminated material.
Solution Number
Sol 101
Features Used
Linear Static Analysis
Laminated Composites
Reference
1. Test R0031/2 from NAFEMS publication R0031, Composites Benchmarks, February 1995.
Material Properties
Main Index
446 MSC Nastran Verification Guide
Loading
Case-1: Internal Pressure of 200 MPa
Case-2: Internal Pressure of 200 MPa and temperature rise of 130°C. The temperature loading simulates the
pre-stress due to tension in the hoop windings.
Boundary Condition
Axial displacement Uz = 0 at z = 0.
Figure 6.66-1 Wrapped Thick Cylinder under Pressure and Thermal Loading
Main Index
CHAPTER 6.66 447
R0031(2): Wrapped Thick Cylinder under Pressure and Thermal Loading
Results
MSC Nastran
Hoop Stress
(MPa) NAFEMS Values % Error
Case 1 at radius = 24 1483 1414.2 -4.64%
at radius = 26 822 874.8 6.42%
Case 2 at radius = 24 1309 1235.9 -5.58%
at radius = 26 994 1052.9 5.93%
Conclusion
The hoop stress values in the inner and outer cylinders for the two load cases at r = 24 and r = 26 from MSC
Nastran are compared with the NAFEMS results in the above table. SOL400 variant of this model is
demonstrated in Chapter 8: Wrapped Thick Cylinder under Pressure and Thermal Loading of MSC Nastran Implicit
Nonlinear Demonstration Guide.
Main Index
448 MSC Nastran Verification Guide
Problem Description
This NAFEMS composite benchmark is a simply supported square sandwich plate subjected to a normal
pressure loading. The orthotropic materials of the sandwich plate are defined in such a way that the facing
sheets carry the bending and the shear is carried by the core.
Solution Number
Sol 101
Features Used
Linear Static Analysis
Laminated Composites
Reference
1. Test R0031/3 from NAFEMS publication R0031, Composites Benchmarks, February 1995.
Material Properties
E1 = 10×105 psi E2 = 4×106 psi 12 = 0.3 and
G12 = 1.875×106 psi G13 = 1.875×106 psi G23 = 1.875×106 psi for face sheets
E1 = 10 psi E2 = 10 psi 12 = 0.0 and
G12 = 10 psi G13 = 3 ×104 psi G23 = 1.2×104 psi for core
Loading
Plate is subjected to uniform pressure of 100 psi.
Main Index
CHAPTER 6.67 449
R0031 (3): Three-layer Sandwich Shell under Normal Pressure Loading
Boundary Condition
The plate is simply supported at all four edges and fixed at four corners. Symmetric boundary conditions are
applied at the symmetric planes of one quarter model.
Figure 6.67-1 Wrapped Thick Cylinder under Pressure and Thermal Loading
Results
MSC Nastran
NAFEMS Values % Error
Uz at C (in) -0.123 -0.1226 -0.33%
11 at C (kpsi) 34.45 34.03 -1.22%
22 at C (kpsi) 13.93 13.29 -4.59%
12 at E (kpsi) -5.07 -5.04 -0.59%
Conclusion
The Z-displacement at point C and stress values at points C and E from MSC Nastran are compared with
the NAFEMS results in the above table. The variants of this model with other composite element types are
demonstrated with SOL400 in Chapter 9: Three-layer Sandwich Shell under Normal Pressure Loading of MSC
Nastran Implicit Nonlinear Demonstration Guide.
Main Index
450 MSC Nastran Verification Guide
Contact
Main Index
CHAPTER 6.68 451
R0094 (1): Benchmark 1: 2D Cylinder Roller Contact
Problem Description
In this NAFEMS contact benchmark, a steel cylinder is pressed into an aluminum block. The solution with
a 2-D plane strain approximation of this problem from Nastran is compared to an analytical solution for a
frictionless case.
Solution Number
SOL 400
Features Used
Nonlinear static analysis.
Advancing contact area with curved contact surfaces
Deformable to deformable contact
Reference
1. Benchmark 1 from NAFEMS publication R0094, Advanced Finite Element Contact Benchmarks,
2006.
Material Properties
Ecylinder = 210 kN/mm2
cylinder = 0.3
Eblock = 70 kN/mm2
block = 0.3
Loading
Vertical point load F = 35kN
Main Index
452 MSC Nastran Verification Guide
Boundary Condition
Symmetric displacement constraints along vertical symmetric line
Bottom surface of the block is fixed (Ux = Uy = 0.0)
Touching contact between cylinder and block with segment to segment contact
Main Index
CHAPTER 6.68 453
R0094 (1): Benchmark 1: 2D Cylinder Roller Contact
Results
Conclusion
The contact pressure results along the touching surfaces of steel cylinder from MSC Nastran are compared
with the NAFEMS results in the above plot. The length of contact zone and maximum contact pressure
values from MSC Nastran are compared with NAFEMS results in the above table. More details of this
contact benchmark model with and without friction are demonstrated using node to segment contact in 2-
D Cylindrical Roller Contact of MSC Nastran Implicit Nonlinear Demonstration Problems Guide.
Main Index
454 MSC Nastran Verification Guide
Problem Description
In this NAFEMS contact benchmark, an axisymmetric steel punch is compressed on an aluminum
foundation. The solution with a 2-D axisymmetric model of this problem from Nastran is compared to
NAFEMS solution for a frictionless case.
Solution Number
SOL 400
Features Used
Nonlinear static analysis.
Axisymmetric contact
Deformable to deformable contact
Reference
1. Benchmark 2 from NAFEMS publication R0094, Advanced Finite Element Contact Benchmarks,
2006.
Main Index
CHAPTER 6.69 455
R0094 (2): Benchmark 2: 3D Punch (Rounded Edges)
Material Properties
Epunch = 210 kN/mm2
punch = 0.3
Efoundation = 70 kN/mm2
foundation = 0.3
Loading
Uniform pressure in punch P = 100 N/mm2
Boundary Condition
Bottom surface of the foundation is fixed (Ux = Uy = 0.0)
Touching contact between punch and foundation with node to segment contact
Main Index
456 MSC Nastran Verification Guide
Results
Conclusion
The radial displacement results along the touched surfaces of aluminum foundation from MSC Nastran are
compared with the NAFEMS results in the above plot for the cases with and without friction. More details
of this contact benchmark model with 3D solid elements are demonstrated in Chapter 2: 3-D Punch (Rounded
Edges) Contact of MSC Nastran Implicit Nonlinear Demonstration Guide.
Main Index
CHAPTER 6.70 457
R0094 (3): Benchmark 3: 3D Sheet Metal Forming
Problem Description
This NAFEMS contact benchmark is an approximation of Numisheet 2002 - Benchmark B problem.
Simulations are carried out with MSC Nastran with 2D plane strain and 3D shell models to find the angles
before and after spring back.
Solution Number
SOL 400
Features Used
Nonlinear static analysis.
Large sliding contact around circular surface
Deformable to deformable contact
Reference
1. Benchmark 3 from NAFEMS publication R0094, Advanced Finite Element Contact Benchmarks,
2006.
Main Index
458 MSC Nastran Verification Guide
Material Properties
E = 70.5 kN/mm2
= 0.342
Initial yield stress 0 = 194 N/mm2
n
Hollomon hardening: = k
Loading
Prescribed vertical displacement of punch = -28.5 mm
Boundary Condition
Symmetric displacement restraints (half symmetry)
Rigid surface for Die is fixed
Touching contact of sheet with punch and die using node to segment contact
Main Index
CHAPTER 6.70 459
R0094 (3): Benchmark 3: 3D Sheet Metal Forming
Results
Main Index
460 MSC Nastran Verification Guide
Conclusion
The force-displacement curves of punch from MSC Nastran for 2D plane strain and 3D shell models are
compared with the NAFEMS results in the above plots for the cases with and without friction. More details
of this contact benchmark model are presented in Chapter 3: 3-D Sheet Metal Forming of MSC Nastran Implicit
Nonlinear Demonstration Guide.
Main Index
CHAPTER 6.71 461
R0094 (4): Benchmark 4: 3D Loaded Pin
Problem Description
This NAFEMS contact benchmark evaluates the performance of contact algorithms at the curved
boundaries. A cylindrical pin is located in the cylindrical hole of a strip, with identical diameters of the hole
and the pin. Two equal point loads are applied to the center of the pin, resulting to a loss of contact at one
side of the pin and a localized contact area on the other side. A 2D plane strain representation of the model
is used to evaluate the contact normal stress in the contact surfaces of strip and pin.
Solution Number
SOL 400
Features Used
Nonlinear Static Analysis.
Receding Contact Area with Curved Contact Surfaces
Deformable to Deformable Contact with Friction
Reference
1. Benchmark 4 from NAFEMS publication R0094, Advanced Finite Element Contact Benchmarks,
2006.
Main Index
462 MSC Nastran Verification Guide
Material Properties
Epin = 210 kN/mm2
pin = 0.3
Esheet = 70 kN/mm2
sheet = 0.3
Loading
Point load of 100kN is applied at the center of pin.
Boundary Condition
Symmetric displacement restraints (half symmetry).
Left side of sheet is fixed.
Touching contact of strip with pin (Node to Segment contact).
Main Index
CHAPTER 6.71 463
R0094 (4): Benchmark 4: 3D Loaded Pin
Results
Main Index
464 MSC Nastran Verification Guide
Conclusion
The results from MSC Nastran for the variation of contact normal stress along the contact surfaces of pin and
strip are compared with the results from Marc, which is also used with the node to segment contact in the
NAFEMS publication for this benchmark. Segment to segment contact will give more smooth variation of
the contact normal stress in both the pin and strip surfaces. Further details of this contact benchmark model
are presented in Chapter 4: 3-D Loaded Pin with Friction of MSC Nastran Implicit Nonlinear Demonstration
Problems Guide.
Main Index
CHAPTER 6.72 465
R0094 (5): Benchmark 5: 3D Steel Roller on Rubber
Problem Description
In this NAFEMS contact benchmark, a steel cylindrical roller rotates on a rubber base. It tests the ability of
contact algorithms to model the movement of two surfaces when they are in contact, as the rubber will stretch
compared to the steel cylinder.
Solution Number
SOL 400
Features Used
Nonlinear Static Analysis.
Rolling Contact with Friction
Deformable and Rigid Contact Bodies
Reference
1. Benchmark 5 from NAFEMS publication R0094, Advanced Finite Element Contact Benchmarks,
2006.
Material Properties
Esteel = 210 kN/mm2
steel = 0.3
C10, rubber = 10 kN/mm2 with Neo-Hookean material description
Loading
Move bottom surface of rubber 3 mm up during time period 0-1 second and held at this position from 1-2
second. Prescribed rotation of steel roller by 360 degrees during time period 1-2 second.
Main Index
466 MSC Nastran Verification Guide
Boundary Condition
No rotation of roller and sheet x-displacement fixed during time period 0-1 second.
Sheet is free to move in horizontal direction during time period 1-2 second.
Touching contact of roller and sheet and rotating rigid body is glued to roller (Segment to Segment contact).
Main Index
CHAPTER 6.72 467
R0094 (5): Benchmark 5: 3D Steel Roller on Rubber
Results
Main Index
468 MSC Nastran Verification Guide
Conclusion
The results from MSC Nastran for the variation of vertical force on roller are compared with the results from
Marc, which is also used in the NAFEMS publication for this benchmark. Further details of this contact
benchmark model are presented in Chapter 5: Steel Roller on Rubber of MSC Nastran Implicit Nonlinear
Demonstration Problems Guide.
Main Index
CHAPTER 6.72 469
R0094 (5): Benchmark 5: 3D Steel Roller on Rubber
Geometric Nonlinear
Main Index
470 MSC Nastran Verification Guide
Problem Description
This NAFEMS geometric nonlinear benchmark on Z-shaped cantilever under an end load is used to test
beam element's formulation for large displacement and rotations. The model is subjected to significant stress
stiffening accompanied by stress reversal.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
3D Large rotation - large displacement with membrane and bending actions.
Tension stiffening with change in sign of bending moment at section A.
Reference
1. Test 3DNLG-1 from NAFEMS Publication 'A Review of Benchmark Problems for Geometric Non-
Linear Behaviour of 3D Beams & Shells (Summary)' by Prinja. N & Clegg. R, R0024, August 1993.
Material Properties
E = 2.0x107 = 0.3
Loading
Total load P = 4000.0 is applied at the free end in the Z-direction in 3 different steps.
Main Index
CHAPTER 6.73 471
3DNLG1: Elastic Large Deflection Response of A Z-Shaped Cantilever under End Load
Boundary Condition
All degrees of freedom are constrained at the built-in end.
Results
Tip Z-Displacement
Number of Elements
Target
Load Value 30 60 90 120 150
104.5 80.4 80.690 80.686 80.685 80.685 80.685
1263 133.1 133.146 133.126 133.122 133.121 133.120
4000 143.4 143.534 143.482 143.472 143.469 143.467
Main Index
472 MSC Nastran Verification Guide
Conclusion
The results for tip Z-displacement and moment at section A-A from MSC Nastran are compared with
NAFEMS results at loads P = 104.5, 1263.0 and 4000.0 in the above tables. The largest deviation of tip Z-
displacement from the NAFEMS solution is less than 0.4% and largest deviation of moment at section A-A
from the NAFEMS solution is less than 3.8%.The results of MSC Nastran is also found to be accurate even
for the course mesh of 30 elements for this model.
Main Index
CHAPTER 6.74 473
3DNLG2: Elastic Large Deflection Response of a Pear-Shaped Cylinder under End Shortening
Problem Description
This NAFEMS benchmark involves initial softening in the load defection path due to the collapse of the flat
plate sections. This is accompanied by stress redistribution and subsequent stiffening to a maximum load,
which causes buckling of curved plate sections.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis with elastic large defection and large rotation
Initial softening behavior followed by stiffening to final buckling
Interaction of bending, membrane and shear effects.
Reference
1. Test 3DNLG-2 from NAFEMS Publication 'A Review of Benchmark Problems for Geometric Non-
Linear Behavior of 3D Beams & Shells (Summary)' by Prinja. N & Clegg. R, R0024, August 1993.
Material Properties
E = 1.0×107 and = 0.3
Main Index
474 MSC Nastran Verification Guide
Loading
Uniform negative z-displacement at nodes initially lying on plane z =0.4. Due to limitation of arc length
method's implementation in MSC Nastran for enforced displacement, the end shortening is applied with a
maximum load of 2600 in -z direction.
Boundary Condition
Ux = y = z = 0.0 on plane x = 0.0, Uz = x = y = 0.0 on plane z = 0.0 and Ux = Uy = 0.0 on plane z = 0.4.
Main Index
CHAPTER 6.74 475
3DNLG2: Elastic Large Deflection Response of a Pear-Shaped Cylinder under End Shortening
Results
Conclusion
The results from MSC Nastran are compared with NAFEMS results at loads 249.26, 1285.44, 2255.87,
2437 in the above figure. The largest deviation from the NAFEMS solution around 3.3% at all these
reference points. The collapse load predicted by MSC Nastran is 2511.97 which is 3.1% higher compared
to the NAFEMS estimated collapse load of 2437.
Main Index
476 MSC Nastran Verification Guide
Problem Description
This NAFEMS benchmark on buckling of imperfect flat plate demonstrates shell element's formulation for
modeling shear behavior.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Large rotation - large displacement response with in-plane shear
Initial imperfection.
Reference
1. Test 3DNLG-6 from NAFEMS Publication 'A Review of Benchmark Problems for Geometric Non-
Linear Behavior of 3D Beams & Shells (Summary)' by Prinja. N & Clegg. R, R0024, August 1993.
Material Properties
E = 6.4x106 = 0.3
Loading
Maximum shear load/edge = 2.12x108 is applied in all four edges using modified Riks arc length method.
Boundary Condition
Uz = 0.0 in all four edges. Ux = Uy = z = 0.0 at central node.
Main Index
CHAPTER 6.75 477
3DNLG6: Buckling of a Flat Plate When Subjected To In-Plane Shear
Initial Imperfection
The initial imperfection is defined by the following equations:
z = 0.5 t f(x y) f(L 2 L 2)
where;
mx ny
f(x y) = a mn sin ----------
L
- sin ---------
L
m n
a11=1.0a22 = 0.2897
a13 = -0.0706a31 = -0.0691
a13 = 0.0384a42 = 0.0032
These imperfection amplitudes calculated at all grid points are used to define imperfection shape of the
model. The bulk data entry 'IMPGEOM' is used to define this imperfection geometry in MSC Nastran
input file.
Main Index
478 MSC Nastran Verification Guide
Results
Conclusion
The results from MSC Nastran are compared with NAFEMS results at edge loads 0.5442E8, 0.7585E8,
1.074E8, 1.3521E8, 1.7718E8 and 2.1117E8 in the above figure. The largest deviation from the NAFEMS
solution is less than 1.0% at all these reference points.
Main Index
CHAPTER 6.76 479
3DNLG7: Elastic Large Deflection Response of a Hinged Spherical Shell Under Pressure Loading
Problem Description
This NAFEMS benchmark involves snap through buckling response of hinged spherical shell under pressure
load. It demonstrates the ability of Riks arc length method to trace the positive and negative tangents in the
load-displacement equilibrium path.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Snap through buckling with arc length method.
Large displacement - large rotation response.
Reference
1. Test 3DNLG-7 from NAFEMS Publication 'A Review of Benchmark Problems for Geometric Non-
Linear Behaviour of 3D Beams & Shells (Summary)' by Prinja. N & Clegg. R, R0024, August 1993.
Material Properties
E = 69 = 0.3
Loading
Maximum distributed pressure of 0.1 applied normal to shell surface using Riks arc length method.
Main Index
480 MSC Nastran Verification Guide
Boundary Condition
Ux = Uy = Uz = 0.0 on all four edges.
Main Index
CHAPTER 6.76 481
3DNLG7: Elastic Large Deflection Response of a Hinged Spherical Shell Under Pressure Loading
Results
Conclusion
The results from MSC Nastran are compared with NAFEMS results at loads 0.0624, 0.0468, 0.02885 and
0.1 in the above figure. The largest deviations from the NAFEMS solution are around 1.83% and -0.93%
at limit points 1 and 2 respectively.
Main Index
482 MSC Nastran Verification Guide
Problem Description
This NAFEMS benchmark is a good test case for membrane locking of doubly curved shell and it is an
extension of linear NAFEMS benchmark on pinched hemispherical shell to incorporate the effect of large
displacement and large rotation.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Large rigid body rotation with stress redistribution.
Membrane, bending and twisting structural actions.
Reference
1. Test 3DNLG-9 from NAFEMS Publication 'A Review of Benchmark Problems for Geometric Non-
Linear Behavior of 3D Beams & Shells (Summary)' by Prinja. N & Clegg. R, R0024, August 1993.
Material Properties
E = 6.825x107 = 0.3
Loading
Maximum load of Px =.-100 at point A and Py = 100 at point B
Main Index
CHAPTER 6.77 483
3DNLG9: Large Elastic Deflection of a Pinched Hemispherical Shell
Boundary Condition
Ux = y = z = 0.0 on plane x = 0, Uy = x = z = 0.0 on plane y = 0, Uz = 0 at point A to prevent rigid
body motion.
Main Index
484 MSC Nastran Verification Guide
Results
Conclusion
The results from MSC Nastran are compared with NAFEMS results in the above figure. The maximum
deviation of 2.1% and 1.2% are seen at points A and B respectively.
Main Index
CHAPTER 6.78 485
3DNLG10: Elastic-Plastic Behavior of a Stiffened Cylindrical Panel under Compressive End Load
3DNLG10: Elastic-Plastic
6.78 Behavior of a Stiffened
Cylindrical Panel under
Compressive End Load
Problem Description
This NAFEMS benchmark illustrate the snap through behavior of stiffened shell panel shell with interaction
between geometric and material nonlinearities.
Solution Number
SOL 400
Features Used
Implicit nonlinear analysis.
Snap through with elasto-plastic large deflection behavior.
Initial imperfection.
Reference
1. Test 3DNLG-10 from NAFEMS Publication 'A Review of Benchmark Problems for Geometric
Non-Linear Behaviour of 3D Beams & Shells (Summary)' by Prinja. N & Clegg. R, R0024, August
1993.
Material Properties
E = 2.1x105, = 0.3, y=350
Loading
Compressive load of 12875 distributed evenly at end A.
Main Index
486 MSC Nastran Verification Guide
Boundary Condition
Ucirc = radial = z = 0.0 at edges 1 and 2
Ucirc = Uradial = radial = 0.0 at end A.
Uz = circ = radial = 0.0 for panel at Z = L/2
Uz = X = Y = 0.0 for stiffener at Z = L/2
Initial Imperfection
The initial imperfection for cylindrical panel and stiffener are defined by the following equations:
2
Z
R = 0.569 sin ------------ for cylindrical panel (cylindrical system -RZ )
7200
where;
with -9 ≤ ≤ 9 and 0 ≤ Z ≤ 200
x y
and z = – 0.3 sin --------- sin ------
400 20
Main Index
CHAPTER 6.78 487
3DNLG10: Elastic-Plastic Behavior of a Stiffened Cylindrical Panel under Compressive End Load
Results
Conclusion
The results from MSC Nastran are compared with NAFEMS results at loads 23355, 15785 and 18122 in
the above figure. As shown in the following table, the load at limit points 1 and 2 from MSC Nastran show
the deviation of less than 2.0% compared to the NAFEMS solution.
Main Index
488 MSC Nastran Verification Guide
Main Index