MSC Nastran 2022.1 SOL 400 Getting Started Guide
MSC Nastran 2022.1 SOL 400 Getting Started 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: [email protected]
Email: [email protected]
Worldwide Web
www.mscsoftware.com, www.hexagon.com
Support
https://simcompanion.hexagon.com
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 2022 Hexagon AB and/or its subsidiaries. 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.
The Hexagon logo, Hexagon, MSC Software logo, MSC, Dytran, Marc, MSC Nastran, Patran, 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
March 11, 2022
NA:V2022.1:Z:Z:Z:DC-SOL400GS-PDF
Documentation Feedback
At Hexagon Manufacturing Intelligence, 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:
[email protected].
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 Hexagon Manufacturing Intelligence 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
SOL 400 Getting Started Guide
Content
Contents
Preface
Using this Manual. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Prerequisites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Goals of this Manual. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Contents of This Guide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
List of MSC Nastran Guides . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Using other Manuals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
Patran Documentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Typographical Conventions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Accessing MSC Nastran Manuals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Downloading the PDF Documentation Files. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Navigating the PDF Files. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Printing the PDF Files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
Training and Internet Resources . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Technical Support . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Visit SimCompanion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Help Us Help You . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1 Nonlinear Analysis
Linear vs. Non-linear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
Nonlinear Analysis: Basics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
When to Consider Nonlinear Analysis? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Causes of Nonlinearity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
Nonlinear Characteristics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
General Classification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
General Recommendations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
Applications for Nonlinear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Industry Uses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Limitations of Nonlinear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
Performing Nonlinear Static Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
Main Index
2 MSC Nastran SOL 400 Getting Started Guide
3 Capabilities
Nonlinear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Applications for Nonlinear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Static Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Linear Analysis. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Linear Perturbation Analysis. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Linear Static. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
Normal Modes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
Direct and Modal Linear Transient Dynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
Direct and Modal Frequency Response. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Direct and Modal Complex Eigenvalue . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
General Nonlinear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Nonlinear Static . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Nonlinear Transient Dynamic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Creep . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Viscoelastic . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Heat Transfer Procedures. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Nonlinear Transient Response Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Nonlinear Transient Response Analysis Interface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Time Step Definition. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
Coupled Thermal-Mechanical . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Thermal Contact. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Benefits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Main Index
Contents 3
5 Elements
Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
Element Classes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
0-D. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
1-D Elements – Not Numerically Integrated . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
1-D Elements that are Numerically Integrated. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
Beam Element Considerations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
Large Displacement/Large Strain . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Planar Continuum Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Axisymmetric Shell Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3-D Membrane, Plate, and Shell Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
Shells . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
Shear Panel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
3-D Solid Shell Element . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
3-D Volumetric Solid Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
Composite Solid Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
Interface Elements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
Automatic Property Mapping . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
6 Materials
Material Model Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
Material Property Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
SOL 400 Material Entries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
Linear Elastic Behavior . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
Element Selection for Incompressible Materials . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
Main Index
4 MSC Nastran SOL 400 Getting Started Guide
7 Contact
Introduction to Contact . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
Contact Types . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
Touching Contact . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
Glued Contact . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
Cohesive Contact . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
Contact Definition Method. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 123
Contact Bodies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
Deformable Contact Bodies (3D/2D/1D) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125
Creating a Deformable Body. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126
Rigid Contact Bodies (3D/2D) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128
Control of Rigid Body . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
Contact in MSC Nastran. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
Defining Contact Interactions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131
Contact Settings for a SOL400 Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
Node-to-segment (NTS) Contact. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
Segment-to-Segment (STS) contact . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
Contact Considerations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
Contact Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
Main Index
Contents 5
8 Constraints
Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Constraints . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Single-Point Constraints (SPC and SPC1) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Enforced Motion Constraints (SPCD and SPCR) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
Applying Constraints . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
Single Point Constraint . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
Multipoint Constraints (MPC). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
Static Loads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
Using Patran to Apply Loads and Boundary Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
CWELD/CFAST/CSEAM Element Enhancements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
Benefits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143
Description of Features. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143
Enhanced Search Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
Additional Information. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
Connector Stiffness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
Detailed Projection Algorithm for Best Possible Projection. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 150
9 Boundary Conditions
Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 153
Zero and Enforced Displacements. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 153
Enforced Motion Constraints (SPCD and SPCR) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 153
Fixed Direction Grid Point Forces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 153
p-Element Loads and Constraints . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
Thermal Loads (TEMP and TEMPD) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
Inertial and Dynamic Loads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156
Gravity and Centrifugal Force . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156
Initial Stress and Initial Plastic Strain Mapping from Previous Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
Main Index
6 MSC Nastran SOL 400 Getting Started Guide
11 Trouble Shooting
Overview. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
Review Fatal Error Message . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182
Review the .sts File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 183
Review the .f06 File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 184
Review NLSTEP parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
Review Nonlinear Iteration Diagnostics. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
Request for More Diagnostics. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 187
Review the Intermediate Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 188
SOL 400 Analysis Messages . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 188
Reviewing Convergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 190
Determining if Your Simulation has failed to Converge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 191
Review Plot Results from Converged Increments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 192
Reviewing the MSC Analysis Manager . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 193
Tips for Starting with Nonlinear Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 194
Output Messages . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 194
Main Index
Preface
Preface
Using this Manual 2
List of MSC Nastran Guides 3
Using other Manuals 4
Typographical Conventions 5
Accessing MSC Nastran Manuals 6
Training and Internet Resources 7
Technical Support 8
Main Index
2 MSC Nastran SOL 400 Getting Started Guide
Using this Manual
Prerequisites
No prior experience with commercial finite element software is assumed and no finite element-specific
university coursework is required. It is assumed that you have a bachelor’s degree in any of the fields relevant
to structural analysis: mechanical engineering, civil engineering, engineering mechanics, or the equivalent.
Main Index
CHAPTER 3
Preface
8. Chapter 8: Constraints:This chapter discusses the procedures for applying loads and constraints which
can be used in SOL 400 models.
9. Chapter 9: Boundary Conditions: This chapter describes the types of loads and boundary conditions which
can be applied using SOL 400
10. Chapter 10: Iteration Control in Nonlinear Analysis: This chapter provides a comprehensive FE solution for
multi-physics problems such as structure analysis, thermal analysis, as well coupled analysis.
11. Chapter 11: Trouble Shooting: This chapter provides some trouble shooting tips that you can review,
when you receive non-convergence error messages (when a job does not converge).
12. Appendix A: Frequently Asked Questions: This Appendix contains a list of frequently asked question
regarding SOL 400 and the answers to those questions.
Main Index
4 MSC Nastran SOL 400 Getting Started Guide
Using other Manuals
Rotordynamics
Implicit Nonlinear (SOL 400)
Explicit Nonlinear (SOL 700)
Aeroelastic Analysis
User Defined Services
Non Linear (SOL 600)
High Performance Computing
DEMATD
You may find any of these documents from Hexagon at:
https://simcompanion.hexagon.com/customers/s/article/MSC-Nastran-Support-Home-Page
Main Index
CHAPTER 5
Preface
This guide contains many highlighted links (in blue) to other MSC Nastran documents and all the
documents were delivered together as a collection. If you keep the collection together the links between
documents will work.
Two ways of working with links are as follows:
Use alt to return back to the window your curser is in.
Open the other linked to documents in a new window from an Adobe Reader.
Edit Preferences Documents Open.
a. Choose
b. Cross-document links in the same window.
c. Uncheck the checkbox.
d. Select OK.
Patran Documentation
Three key books from the Patran library may be of assistance in running MSC Nastran Nonlinear:
Patran User’s Guide: This introductory guide gives you the essential information you need to
immediately begin using Patran for MSC Nastran Nonlinear projects. Understanding and using the
information in this guide requires no prior experience with CAE or finite element analysis. •
Patran Reference Manual: This a counterpart to the MSC Nastran Reference Guide, this manual
provides complete descriptions of basic functions in Patran, geometry modeling, finite element
modeling, material models, element properties, loads and boundary conditions, analysis, and results.
MSC Nastran Preference Guide: This gives specific information that relates to using Patran with
MSC Nastran as the intended analysis code. All application forms and required input are tailored to
MSC Nastran.
Typographical Conventions
The section provides a brief overview of the typographical conventions used in the document to help the user
better follow the MSC Nastran documentation.
This section describes some syntax that will help you in understanding text in the various chapters and thus
in facilitating your learning process. It contains stylistic conventions to denote user action, to emphasize
particular aspects of a MSC Nastran run or to signal other differences within the text.
Courier New Represents command-line options of MSC Nastran and results from f04/f06 files.
Example: memorymax=16gb
Arial font To represent elements.
Example: RBE3 and RSPLINE are interpolation elements and are not rigid.
Main Index
6 MSC Nastran SOL 400 Getting Started Guide
Accessing MSC Nastran Manuals
Example: smp=16
Bold Text Represents items in the text that we want to emphasize.
Example: dmp=4
Italic Text Represents references to manuals/documents.
Note: Since there is no user interface in MSC Nastran, we shall use bold font to emphasize.
Main Index
CHAPTER 7
Preface
If the page is too large to fit on your paper size, you can reduce it by doing the following:
1. Select the File -> Print.
2. Under Page Scaling, choose the Shrink to Printable Area option.
The SimCompanion link above gives you access to the wealth of resources for Hexagon 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 Hexagon, allows you to manage your support cases, and participate in our
discussion forums.
http://www.mscsoftware.com/msc-training
Main Index
8 MSC Nastran SOL 400 Getting Started Guide
Technical Support
The MSC-Training link above will point you to schedule and description of MSC Seminars. Following
courses are recommended for beginning 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
results. This seminar provides the foundation required for intermediate and advanced MSC Nastran
applications.
Technical Support
If you encounter difficulties while using MSC Nastran, first refer to the section(s) of the manual containing
information you are trying to use or the type of problem you are trying to solve.
Visit SimCompanion
The product documentation is available in SimCompanion (http://simcompanion.hexagon.com). The
SimCompanion gives you access to a wealth of resources for Hexagon products. You will find various
information such as:
Product documentations
Knowledge base articles
Product error lists (fixed and known issues for each release)
SimAcademy webinars
Product and support contact information
SimCompanion 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 Hexagon, allows you to manage
your support cases, and participate in our discussion forums.
Main Index
Chapter 1: Nonlinear Analysis
1 Nonlinear Analysis
Linear vs. Non-linear 10
Nonlinear Analysis: Basics 10
When to Consider Nonlinear Analysis? 11
Causes of Nonlinearity 11
General Classification 12
General Recommendations 13
Applications for Nonlinear Analysis 14
Limitations of Nonlinear Analysis 14
Performing Nonlinear Static Analysis 15
Main Index
10 MSC Nastran: SOL 400 Getting Started Manual
Linear vs. Non-linear
P = K P u u
As the structure displaces due to loading, the stiffness changes, and as the stiffness changes the structure’s
response changes. As a result, nonlinear problems require incremental solution schemes that divide the
problem up into steps calculating the displacement, then updating the stiffness. Each step uses the results
from the previous step as a starting point. As a result the stiffness matrix must be generated and inverted many
times during the analysis adding time and costs to the analysis.In addition, because the response is not
proportional to the loads, each load case must be solved separately and the principle of superposition is not
applicable.
Main Index
MSC Nastran: SOL 400 Getting Started Manual 11
Nonlinear Analysis
Causes of Nonlinearity
Nonlinear response can be caused by any of several characteristics of a structure:
Large deformations and strains. (Geometric Nonlinearity)
Material behavior (Material Nonlinearity)
The effect of contact (Contact Nonlinearity)
Boundary condition nonlinearities
In reality many structures exhibit combinations of these various nonlinearities. Most of the time either
material behavior is not linear in the operating conditions, or geometry of the structure itself keeps it from
responding linearly.
Due to cost or weight advantage of nonmetals (polymers, woods, composites etc.) over metals, nonmetals are
replacing metals for variety of applications, which have nonlinear load to response characteristics, even under
mild loading conditions. The structures are also optimized to make most of its strength, pushing the load
level so close to the strength of the material, that it starts behaving nonlinearly. In order to accurately predict
the strength of the structures in these circumstances, it is necessary to perform nonlinear analysis.
Main Index
12 MSC Nastran: SOL 400 Getting Started Manual
Nonlinear Characteristics
Nonlinear Characteristics
Modeling for nonlinear analysis is not exempted from the guidelines for good modeling practice pertaining
to linear analysis, which are summarized as follows:
The analyst should have some insight into the behavior of the structure to be modeled; otherwise, a
simple model should be the starting point.
Substructuring should be considered for the modularity of the model and/or synergism between
projects and agencies involved. The structure represented by a substructure is always linear.
The size of the model should be determined based on the purpose of the analysis, the trade-off
between accuracy and efficiency, and the scheduled deadline.
Prior contemplation of the geometric modeling will increase efficiency in the long run. Factors to be
considered include selection of coordinate systems, symmetric considerations for simplification, and
systematic numbering of nodal points and elements for easy classification of locality.
Discretization should be based on the anticipated stress gradient, i.e., a finer mesh in the area of
stress concentrations.
Element types and the mesh size should be judiciously chosen. For example, avoid highly distorted
and/or stretched elements (with high aspect ratio).
The model should be verified prior to the analysis by some visual means, such as plots and graphic
displays.
General Classification
Nonlinear problems of any type require iterative solution methods and incremental loading to obtain
(converge to) a solution, and are generally far more computationally difficult than linear problems. Nonlinear
problems are classified into the following three broad categories:
Geometric nonlinearity: In the structures where stiffness is dependent on the displacement which
they may undergo are termed geometrically nonlinear. If a continuum body undergoes large finite
deformations, the strain-displacement relations become nonlinear. For structural mechanics
problems, under large deformations, the stiffness changes with deformation thus making the
problem nonlinear. (Buckling problems are also nonlinear).
Geometric nonlinearity accounts for phenomena such as:
• The stiffening of a loaded clamped plate,
• Buckling or snap-through behavior in slender structures or components.
Without taking these geometric effects into account, a computer simulation can fail to predict the
real structural behavior.
Material nonlinearity: It refers to the ability for a material to exhibit a nonlinear stress-strain
(constitutive) response. This is one of the most common forms of nonlinearity. Material
nonlinearity is often, but not always, characterized by a gradual weakening of the structural response
as an increasing force is applied, due to some form of internal decomposition.
Main Index
MSC Nastran: SOL 400 Getting Started Manual 13
Nonlinear Analysis
Elasto-plastic, hyperelastic, crushing, and cracking are good examples, but this can also include
temperature and time-dependent effects such as visco-elasticity or visco-plasticity (creep). For
thermal problems, a temperature dependent thermal conductivity will produce nonlinear equations.
This is required to predict plastic strains in metallic parts, cracking or crushing of concrete, or
extreme deformation of plastic or rubber materials. It is nonlinear constitutive relation.
Boundary/Contacts nonlinearity: When considering either highly flexible components, or
structural assemblies comprising multiple components, progressive displacement gives rise to the
possibility of either self or component-to-component contact. This characterizes to a specific class of
geometrically nonlinear effects known collectively as boundary condition or contact nonlinearity.
In boundary condition nonlinearity the stiffness of the structure or assembly may change
considerably when two or more parts either contact or separate from initial contact. Examples
include bolted connections, toothed gears, and different forms of sealing or closing mechanisms.
Problems involving contact mechanics normally include a boundary condition that depends on the
deformation thereby producing a nonlinear formulation. Thermal problems involving melting or
freezing (phase change) also include such nonlinear boundary conditions.
They are required to predict change in status and/or sliding friction between assembly parts. It is on-
constant displacement BCs.
General Recommendations
With these points in mind, additional recommendations are imperative for nonlinear analysis:
Identify the type of nonlinearity and localize the nonlinear region for computational efficiency. If
unsure, perform a linear analysis to help understand the problem.
Segregate the linear region by using superelements and/or linear elements if possible. Notice that the
potentially nonlinear elements can be used as linear elements.
The nonlinear region usually requires a finer mesh. Use a finer mesh if severe element distortions or
stress concentrations are anticipated.
The step/subcase structure should be utilized properly to divide the load or time history for
conveniences in data recovery, and database storage control, not to mention changing constraints
and loading paths.
The load or time for the subcase with NLSTEP should then be further divided into increments, for
the purpose of convergence control. Automatic adaptive load incrementation, such as is available via
NLSTEP, is the recommended method.
Many options are available in solution methods to be specified on the NLSTEP data entries. The
defaults should typically be used on all options to gain experience before experimenting with other
options.
Normal rotation for the drilling degree of freedom of shell elements restrained by the default value of
100 on the K6ROT parameter when the geometric nonlinearity is involved. In rare cases it may be
necessary to adjust his value. This can help with convergence, but may also affect the results.
Understand the basic theory of plasticity, creep, or rubber elasticity before using these capabilities.
Main Index
14 MSC Nastran: SOL 400 Getting Started Manual
Applications for Nonlinear Analysis
The time step size for a transient response analysis should be carefully considered based on the
highest natural frequency of interest because it has significant effects on the efficiency as well as the
accuracy of the solution. The automated procedures used by NLSTEP is adequate for this purpose.
Industry Uses
Aerospace and Defense: Landing gear, Wing structures, Fuselage, Seals and hoses, Sheet metal
forming
Automotive: Power train, Tire, Seals and gaskets, Exhaust systems, Brakes, Suspension, Gear
contact, welding, joints and connectors
Electronics: Soldering, Welding, Drop tests, Sealing, Switches and connectors
Energy: Wind turbine blades, Composite blade failure, Gear trains, Packers, blow out preventers,
Seals and gaskets, Pipes and casings, Weldments, Drill bits
Heavy Equipment and Machinery: Gears, Steering yokes, Belts, Hoses, Metal forming, Hose
crimping, Wire crimping, Curing, Welding, Extrusion
Medical: Stents, Catheters, Pacemakers, Dental and knee implants, Prosthetics, Muscle and tissue,
Hospital equipment like beds, wheel chairs
Rail: Tip-over stability study, Structural components, Welding, Joints and connectors
Shipbuilding: Structural analysis, Riveting, Bolts, Welding, Sealing
Main Index
MSC Nastran: SOL 400 Getting Started Manual 15
Nonlinear Analysis
The initial state of stress (e.g. residual stresses from welding, temperature, or prestressing of
reinforcement and cables) may be extremely important for the overall response
Main Index
16 MSC Nastran: SOL 400 Getting Started Manual
Performing Nonlinear Static Analysis
If the FE analysis is being done to virtually validate any test done at lab, then it is a good practice to
visit the testing facility and take important measurements on test fixtures and loading devices. These
measurements will help apply loads and boundary conditions in FE model to in the same way as the
part or test specimen is subjected to at the time of testing. For example, measurements of loading
device would help placing the loads on certain nodes or elements. Similarly the fixture dimensions
would dictate constraints locations (nodes) and its degrees of freedom.
4. Specify nonlinear analysis control parameters: The basic controlling parameters for nonlinear static
analysis are initial increment, minimum and maximum increment, maximum numbers of iterations,
the interval at which results file are to be output and convergence criteria for iterations (acceptable
residual load).
5. Run the analysis: The FE model is now ready to be run.The analysis run command may have options
to specify solver version, memory size, and number of CPUs to better control execution.
6. Review and interpret results: It is highly recommended that the analysis results should be carefully
reviewed and checked for accuracy before making any conclusions based on simulations. There are
many ways the FEA results can be checked, some of them are:
• Observe for unexpected movements in the animation.
• Compare the reaction forces against applied forces.
• Check if stresses and strains are as per material properties supplied to the FE model
• Check interacting surfaces in the contact set for any malfunction
• Make quick hand calculations by simplifying the problem and compare it with the FE results.
While reporting the FE results, you should always share all the assumptions made while building the FE
model. The FEA results should always be verified using engineering judgment and past results with similar
FE model.
If there are unrealistically high or low results, you should review the model. It is recommended to maintain
consistency in the model to increase accuracy of the results, especially when analysis is being performed to
evaluate effect of changing certain parameter in the model.
Main Index
MSC Nastran: SOL 400 Getting Started Manual 17
Nonlinear Analysis
Main Index
Chapter 2: Overview of SOL 400
What’s New in SOL 400 since 2016 Release 19
Introduction to SOL 400 21
History of SOL 400 22
Advantages of SOL400 22
Capabilities of SOL400 23
Converting Nastran Linear to Nastran Nonlinear 24
Analysis Types 25
Analysis Procedures 27
Main Index
CHAPTER 2 19
Overview of SOL 400
2021.1
Accelerated separation check for node-to-segment contact
2021
Brake squeal enhancement
Support specification of contact status in linear perturbation analysis with touching contact
Segment-to-segment contact default settings
SOL 400 algorithm improvements
Element strain energy (ESE) and kinetic energy (EKE) with linear perturbation analysis
External Superelement MONPNT1 and MONPNT3 Data Recovery
Addition of Lossy Compression for Nastran HDF5 Matrix and NLOUT Outputs (NH5RDB)
2020
Geometric perturbation
Nonlinear buckling analysis
New MATVE format
Check Input errors when using Advanced Nonlinear Elements
Support for NLLOAD output in HDF5
2019FP1
Linear perturbation buckling analysis
Surface contact
Linear contact in segment-to-segment contact
2019
Model contact check
Support for SEGANGL on BCPARA (and BCAUTOP)
Support for N3DSUM on NLOPRM, NLDBG to provide the simplified debug output
Output of advanced composites
Monitor points with NLSTAT and NLTRAN
2018
Nonlinear CFAST/PFAST
Main Index
20 MSC Nastran: SOL 400 Getting Started Manual
What’s New in SOL 400 since 2016 Release
2017
Solder creep material model
Cohesive contact
Output the resultant contact force and moment for each contact pair
THRU capability for BCTABL1 entry
2016
Support advance elements in Superelements and NLRESTART
Contact in small deformation simulations
Segment-to-segment contact enhancements
Beam contact
Interference fit
Clearance (initial gap)
Contact separation improvement
Support connector elements (CWELD/CFAST/CSEAM) in SOL 400
RC Network heat transfer analysis
Main Index
CHAPTER 2 21
Overview of SOL 400
In MSC Nastrans advanced integrated nonlinear solution (SOL 400), all the MSC Nastran infrastructure is
available. This is the recommended default solution for solving nonlinear problems. SOL 400 is a powerful,
easy to use tool for simulating manufacturing processes and component designs. An extensive finite element
library for building your simulation model, and a set of solution procedures for the nonlinear analysis, which
can handle very large matrix equations, are available in both solution sequences of MSC Nastran implicit
nonlinear.
There are two forms of this combination.
The first form is used in SOL 400 where the algorithms of Marc are embedded completely in MSC
Nastran to form a completely integrated MSC Nastran solution. We recommend this as the default
solution for solving MSC Nastran nonlinear problems.
The second method is SOL 600 where MSC Nastran preprocesses the data and calls the Marc
solver. In SOL 400, all MSC Nastran infrastructure is available, while in SOL 600 only selected
MSC Nastran infrastructure capabilities have been integrated.
This is why MSC recommends that SOL 400 be the default solution method for solving nonlinear
problems. An extensive finite element library for building your simulation model, and a set of
solution procedures for the nonlinear analysis, which can handle very large matrix equations, are
available in both solution sequences of MSC Nastran implicit nonlinear.
Main Index
22 MSC Nastran: SOL 400 Getting Started Manual
History of SOL 400
Advantages of SOL400
The advantages of SOL 400 are as follows:
Allows for analysis chaining which are:
• Automatically chains together sequences of analyses with output state of one used as input state
for another
• Model complete processes in a single simulation through analysis chaining.
Utilizes native Nastran elements, where no translation is required.
Follower force distributed loads.
Temperature applied to nodes. Temperature can be applied as a load in a structural analysis. The
reference temperature is user definable.
Inertial body forces, acceleration and velocity can be applied in the global coordinate system.
Combines static and transient into one analysis
• Pre-stress, transient, steady-state analysis chaining
• Thermal-structural analysis chaining
• Multiple, independent loadcases in one run
• Linear perturbation
Use general contact capability: Solid-to-solid, surface-to-surface, edge-to-edge and beam-to-beam,
etc.
Not restricted to small-strain element limitations
Main Index
CHAPTER 2 23
Overview of SOL 400
Capabilities of SOL400
It is an advanced nonlinear solution process that combines capabilities of multiple solution sequences and
software components into a common solution.
SOL 400 solves linear and nonlinear (material, contact and/or geometric) static, heat transfer,
modal (vibration), buckling, and transient dynamic structural finite element problems.
Eigenvalue solutions are available in SOL 400 for solving linear or nonlinear modal analyses and
linear buckling analysis using either Lanczos or inverse power sweep methods of iteration. Through
the use of parameters you can control the convergence of the eigenvalues, and the modes to retain.
Forces can be applied to nodes in any coordinate system.
Constrained nodal displacements (zero displacements at specified DOF). Enforced nodal
displacements (nonzero displacements at specified DOF in the nodal coordinate system). Both
constrained and enforced displacements can be specified as relative or absolute.
Defines contact between two bodies by selecting the contacting bodies and defining the contact
interaction properties, where gluing and un-gluing properties are provided. Enforced motion or
velocity of rigid contacts surfaces is available.
Supports the following elements types:
• Scalar elements
• Beams
• Shells
• 2D plane strain
• 2D plane stress
• Axisymmetric
• 3D solids
• Lower-order elements
• Higher-order elements
RBE elements and multi-point constraint equations are supported to tie specific nodes or DOFs
with each other. Special MPC entities are supported, (e.g. rigid links) which can be used to tie two
nodes together or equate the motion of two degrees of freedom. Both small and large rotations are
supported.
SOL 400 supports both temperature independent and dependent, isotropic, orthotropic, and
anisotropic material properties. They can be defined for elastic, elastic-plastic, hyper-elastic, hypo-
elastic, visco-elastic, and creep constitutive models.
Main Index
24 MSC Nastran: SOL 400 Getting Started Manual
Converting Nastran Linear to Nastran Nonlinear
SOL 400 also supports cohesive material, gasket material, and thermo-mechanical shape memory
material. Nonlinear elastic-plastic materials can be defined by specifying piecewise linear stress-strain
curves, which may be temperature and / or rate dependent.
Physical properties can be associated with SOL 400 elements such as the cross-sectional properties of
the beam element, the area of the beam and rod elements, the thickness of shell, plane stress, plane
strain, and membrane elements, spring parameters, and point masses among others.
Fracture mechanics capabilities include VCCT crack propagation and cohesive zone interface and
closure analysis, and a large number of failure index criterion for analyzing delamination of
composite elements.
Laminated composite solid and shell elements are supported in SOL 400 through the PCOMP,
PCOMPG, and PCOMPLS entries of the materials capability. Each layer has its own material, thickness,
and orientation and represents linear or nonlinear material behavior, failure index calculations are
also supported. Fast integration techniques are available with the PCOMPF entry. Equivalent material
models may be incorporated using PSHELL.
Analysis jobs can consist of (possibly) complicated loading histories (such as would occur in a multi-
step manufacturing process). A single SOL 400 analysis (sub-case) may consist of multiple steps that
specify the loading sequence.
Geometric Nonlinearity
• Large displacement RBEs, Bush, and Connectors
• Large displacement beam and shell offsets
Material Nonlinearity
• Extensive nonlinear material library
• Large strain elements
• User subroutines and digimat support
Boundary Nonlinearity
• 3D contact: touching and glued
Adaptive load incrementing (NLSTEP)
Multidiscipline Analysis: chained, linear perturbation, coupled, multiple independent sub-cases
Main Index
CHAPTER 2 25
Overview of SOL 400
Analysis Types
A large class of stress analysis problems can be solved with SOL 400. A fundamental division of stress
problems is whether a static, transient dynamics, or perturbation analysis is to be performed. In a dynamic
response, the inertia effects are important. SOL 400 allows complete flexibility in making this distinction, so
that the same analysis may contain several static, dynamic, and pertubation phases.
An important aspect of the flexibility is the manner in which MSC Nastran SOL 400 allows you to step
through the loading history to be analyzed. This is accomplished by defining the job steps for the analysis.
A basic concept in MSC Nastran SOL 400 is the division of the problem history into job steps. A clear
distinction is made in MSC Nastran SOL 400 between linear analysis and nonlinear analysis procedures.
Loading conditions are defined differently for the two cases, time measures are different, and the results
should be interpreted differently. A step is any convenient phase of the history – a thermal transient, a creep
period, a dynamic transient, etc. In its simplest form, a step is just a static analysis of a load change from one
magnitude to another.
An analysis step during which the response may be nonlinear is called a general analysis step. An analysis step
during which the response may only be linear is called a linear perturbation analysis step. Since MSC Nastran
SOL 400 treats such linear analysis as a linear perturbation about a preloaded, predeformed state, its
capability for doing linear analysis is rather more general than that of a purely linear analysis program.
Analysis types for steps/subcases in SOL 400 include the following:
Type Description
Linear Static Static stress analysis is used when inertia effects can be neglected. During a linear
static step, the model’s response is defined by the linear elastic stiffness at the base
state, the state of deformation and stress at the beginning of the step.
Nonlinear Static Nonlinear static analysis requires the solution of nonlinear equilibrium equations,
for which SOL 400 uses Newton’s method. Many problems involve history
dependent response, so that the solution is usually obtained as a series of
increments, with iteration within each increment to obtain equilibrium. For most
cases, the automatic incrementation provided by SOL 400 is preferred, although
direct user control is also provided for those cases where you have experience with
a particular problem. This includes creep, viscoelastic and body approach
Normal Modes This solution type uses eigenvalue techniques to extract the frequencies of the
current system,
Modal Transient The Duhamel Integral method integrates all of the equations of motion through
Dynamic time. The accuracy is based upon the number of modes extracted.
Nonlinear This solution type is used when nonlinear dynamic response is being studied. For
Transient Dynamic most cases, the automatic incrementation provided is preferred, although direct
user control is also provided for those cases where you have experience with a
particular problem. The Generalized-alpha method has been presented as an
unconditionally stable (for linear systems), second-order algorithm that allows
user-controllable numerical dissipation.
Main Index
26 MSC Nastran: SOL 400 Getting Started Manual
Analysis Types
Direct Frequency Frequency response analysis is a method used to compute structural response to
Response steady-state oscillatory (such as rotating machinery) excitation. In frequency
response analysis the excitation is explicitly defined in the frequency domain. The
direct method solves the coupled equations of motion in terms of forcing
frequency
Modal Frequency The modal method utilizes the mode shapes of the structure to reduce and
Response uncouple the equations of motion. The solution for a particular forcing frequency
is obtained through the summation of the individual modal responses.
Direct Complex Complex eigenvalue analysis is used to compute the damped modes of structures
Eigenvalue and assess the stability of systems modeled with transfer functions. Direct
Complex eigenvalue analysis solves the coupled equations of motion.
Modal Complex The modal method utilizes the mode shapes of the structure to reduce and
Eigenvalue uncouple the equations of motion.
Steady State Heat This solution type is for heat transfer problems where the temperature field can be
Transfer found for the current contact and deformation states of the bodies being studied.
(e.g., it is assumed that there are no changes in the structure during the heat
transfer analysis). For cases where there are mechanical changes in the structure a
coupled thermal-structural solution must be performed.
Transient Heat This solution type is for transient heat transfer problems where the temperature
Transfer field can be found based on the current state of stress and deformation of the
bodies being studied (e.g., it is assumed that there are no changes in the structure
during the heat transfer analysis). For cases where there are mechanical changes in
the structure a coupled thermal-structural solution must be performed. For all
transient heat transfer cases, the time increments may be specified directly, or will
be selected automatically based on a user prescribed maximum nodal temperature
change in a step. Automatic time incrementation is generally recommended.
A complete description of the available analysis types is provided in Chapter 3: Nonlinearity and Analysis
Types.
In each step, choose the solution type. This defines the type of analysis to be performed during the step:
dynamic stress analysis, eigenvalue buckling, transient heat transfer analysis, etc. The procedure choice may
be changed from step to step in any meaningful way, so that you have great flexibility in performing analyses.
Since the state of the model (stresses, strains, temperatures, etc.) is updated throughout all nonlinear analysis
steps, the effects of previous history are always included in the response in each new step. Thus, for example,
if natural frequency extraction is performed after a geometrically nonlinear static analysis step, the preload
stiffness will be included. Superposition cannot be applied in nonlinear problems. In general, a different
loading sequence (reordering of the steps) requires a complete new analysis.
In a nonlinear static analysis, you first determine the total value of loading to be applied at a particular stage
of the analysis. This loading value is selected with the LOAD case control command specifying a load set ID
that exists in the bulk data. In this case, the step functions as a type of landmark in the loading history. It may
be an expected point or a point at which the nature of the loading changes (for example, first applying an
internal pressure loading and then an axial loading on a cylinder). The steps is a major partition of the loading
Main Index
CHAPTER 2 27
Overview of SOL 400
history. The loading history should be divided into subcases since this provides you with more control over
the solution and restart strategy.
Analysis Procedures
This section describes the practical steps involved in setting up and running MSC Nastran SOL 400 jobs,
including the use of Patran. The Patran UI is set up to guide you through the process of setting up the SOL
400 analysis, including
The job information for the executive section (select the solution type and solution parameters),
The subcases and steps (create subcases and steps and use subcase parameters to specify the required
step/subcase control information)
The job submission (use Analyse --> Entire Model -> Full Run)
Monitoring the job while it is running and once it is done (user Monitor -> Job -> View sts to
monitor and Monitor -> Job -> View f06 to debug).
It covers the Patran user interface (UI) and its capabilities for setting up, submitting, and monitoring the job.
It shows the user how to specify the analysis type, but does not go into any detail about the analysis types of
the individual steps being set up. See Nonlinearity and Analysis Types for the details of analysis types.
The STEP is the SOL 400 mechanism for associating loads and boundary conditions, output requests, and
various other parameters to be used during part of a complete run. Each step can be designated with one of
the analysis types listed below. For each Analysis Type, define the Solution Parameters and Output
Requests; these collectively constitute the Analysis Procedures.
In MSC Nastran, case control commands provide the loads and constraints, and load incrementation
method, and controls the program after the initial elastic analysis. Case control commands also include
blocks which allow changes in the initial model specifications. Case control commands can also specify print-
out and postprocessing options.
Each set of load sets must begin with a SUBCASE/STEP command and be terminated by another
SUBCASE/STEP or a BEGIN BULK command. If there is only one load case, the SUBCASE/STEP entry is
not required. The SUBCASE option requests that the program perform another increment or series of
increments. The input format for these options is described in MSC Nastran Quick Reference Guide.
Main Index
28 MSC Nastran: SOL 400 Getting Started Manual
Analysis Procedures
not required.
The SUBCASE option requests that the program perform another increment or series of increments. The
input format for these options is described in MSC Nastran Quick Reference Guide.
In SOL 400, all SUBCASEs are independent from each other. Sub-cases are independent, multidiscipline
load cases. Like SOL 101, many independent subcases can be performed in one single run. Each STEP can
be performed with another type of analysis, i.e. statics, transient, normal modes, and frequency response, etc.
Main Index
Chapter 3: Capabilities
3 Capabilities
Applications for Nonlinear Analysis 30
Nonlinear Analysis 30
Static Analysis 30
Linear Analysis 31
General Nonlinear Analysis 34
Nonlinear Transient Response Analysis 36
Coupled Thermal-Mechanical 38
Thermal Contact 38
Main Index
30 MSC Nastran: SOL 400 Getting Started Manual
Nonlinear Analysis
Nonlinear Analysis
However, we know that in many structures the deflections and the stresses do not change proportionately
with the loads. In these problems the structure’s response depends upon its current state and the equilibrium
equations reflect the fact that the stiffness of the structure is dependent on both u and P .
P = K P u u
As the structure displaces due to loading, the stiffness changes, and as the stiffness changes the structure’s
response changes. As a result, nonlinear problems require incremental solution schemes that divide the
problem up into steps calculating the displacement, then updating the stiffness. Each step uses the results
from the previous step as a starting point. As a result the stiffness matrix must be generated and inverted many
times during the analysis adding time and costs to the analysis.In addition, because the response is not
proportional to the loads, each load case must be solved separately and the principle of superposition is not
applicable.
Static Analysis
Static stress analysis is used when inertia effects can be neglected. The problem may still have a real time scale,
for example when the material has a viscoplastic response, such as rate dependent yield. The analysis may be
linear or nonlinear. Nonlinearity may arise from large displacement effects, material nonlinearity and
boundary nonlinearity (such as contact and friction).
Linear static analysis involves the specification of load cases and appropriate boundary conditions. Solutions
may be combined in a postprocessing mode.
Nonlinear static analysis requires the solution of nonlinear equilibrium equations, for which the program uses
Full Newton-Raphson or Modified Newton-Raphson. Many problems involve history dependent response,
so that the solution is usually obtained as a series of increments, with iteration within each increment to
obtain equilibrium. Increments must sometimes be kept small (in the sense that rotation and strain
increments must be small) to assure correct modeling of history dependent effects, but most commonly the
choice of increment size is a matter of computational efficiency – if the increments are too large, more
iteration will be required. Each solution method has a finite radius of convergence, which means that too
large an increment can prevent any solution from being obtained because the initial state is too far away from
the equilibrium state that is being sought – it is outside the radius of convergence. Thus, there is an
Main Index
CHAPTER 3 31
Capabilities
algorithmic restriction on the increment size. For most cases, the automatic incrementation scheme is
preferred, because it will select increment sizes based on these considerations. Direct user control of
increment size is also provided because there are cases when you have considerable experience with his
particular problem and can therefore select a more economic approach. A complete discussion of the
numerical methods used to solver nonlinear static problems is included in this chapter; to get a converged
solution, see Chapter 6: Setting Up, Monitoring, and Debugging the Analysis.
Geometrically nonlinear static problems frequently involve buckling or collapse behavior, where the load-
displacement response shows a negative stiffness, and the structure must release strain energy to remain in
equilibrium. Several approaches are possible in such cases. One is to treat the buckling response dynamically,
thus actually modeling the kinetic response with inertia effects included as the structure snaps. This is easily
accomplished by using the restart option to terminate the static procedure and switch to a dynamic procedure
when the static solution goes unstable. In some simple cases, displacement control can provide a solution,
even when the conjugate load (the reaction force) is decreasing as the displacement increases. More generally,
static equilibrium states during the unstable phase of the response can be found by using the modified Riks
method. This method is for cases where the loading is proportional – that is, where the load magnitudes are
governed by a single scalar parameter. The method obtains equilibrium solutions by controlling the path
length along the load-displacement curve within each increment (rather than controlling the load or
displacement increment), so that the load magnitude becomes an unknown of the system. The method can
provide solutions even in cases of complex, unstable response. The Riks method cannot be used in contact,
heat transfer, coupled, or enforced motion.
Linear Analysis
In a linear analysis, we implicitly assume that the deflections and strains are very small and the stresses are
smaller than the material yield stresses. Consequently, there is assumed to be a linear relationship between
the applied loads and the response of the structure. The stiffness can be considered to remain constant (i.e.,
independent of the displacements and forces) and the finite element equilibrium equations
P = Ku
are linear where the stiffness matrix K is independent of both u , the generalized displacement vector, and
P , the generalized force vector. This linearity implies that any increase or decrease in the load will produce
proportional increase or decrease in displacements, strains, and stresses. Because of the linear relationship,
you only need to calculate the stiffness of the structure once. From this stiffness representation, you can find
the structure’s response to other applied loads by multiplying the load vectors by the decomposed stiffness
matrix. Linear static problems are solved in a single step. In addition, solutions can be combined using the
principle of superposition.
A linear analysis is the simplest and most cost effective type of analysis to perform. Because linear analysis is
simple and inexpensive to perform and often gives satisfactory results, it is the most commonly used
structural analysis. Nonlinearities due to material, geometry, or boundary conditions are not included in this
type of analysis. The behavior of an isotropic, linear, elastic material can be defined by two material constants:
Young’s modulus, and Poisson’s ratio.
Main Index
32 MSC Nastran: SOL 400 Getting Started Manual
Linear Analysis
In actuality, linear analysis is merely an approximation to the true behavior of a structure. In some cases the
approximation is very close to the true behavior, in other cases linear analysis may provide highly inaccurate
results.
Linear analysis is obtained by considering the response in the step as the linear perturbation response about
the base state. The base state is the current state of the model at the end of the last general, nonlinear analysis
step prior to the linear perturbation step. Thus, the concept of linear analysis in MSC Nastran SOL 400 is
rather general. A simple example of the value of this generalization is the natural frequencies of a violin string
under increasing tension. In this case, geometrically nonlinear analysis of the string can be done in several
steps, in each of which the tension is increased. At the end of each such step, the frequencies can be extracted
in a linear perturbation analysis step.
Load magnitudes (including the magnitudes of prescribed boundary conditions), during a linear perturbation
analysis step, are defined as the magnitudes of the load perturbations only.
Likewise, the value of any solution variable is output as the perturbation value only – the value of the variable
in the base state is not included.
During a linear perturbation analysis step, the model’s response is defined by its linear elastic stiffness at the
base state. Plasticity and other inelastic effects are ignored. For hyperelastic materials, the tangent elastic
moduli in the base state are used. Contact conditions cannot change during a perturbation analysis – they
remain as they are defined in the base state. Frictional slipping is not allowed during perturbation analyses –
all points in contact are assumed to be sticking if friction is present. If geometric nonlinearity is included in
the general, nonlinear analysis upon which the linear perturbation study is being based, stress stiffening or
softening effects and (pressure and other follower force) load stiffness effects are included in the linear
perturbation analysis. In this case, perturbation stresses and strains are defined relative to the base state
configuration. The effects of temperature and field variable perturbations are ignored for materials that are
dependent on temperature and field variables. However, temperature perturbations will produce
perturbations of thermal strain.
Some procedures are purely linear perturbation procedures. These are:
Linear Statics: ANALYSIS = STATICS
Bifurcation Buckling: ANALYSIS = BUCK
Natural Frequency: ANALYSIS = MODES
Modal Linear Transient: ANALYSIS = MTRAN
Modal Complex Eigenvalue: ANALYSIS = MCEIG
Modal linear transient analysis and linear static analysis are done in the time domain. The step time of linear
perturbations is never accumulated into the total time. For linear static perturbations the step time always
begins at zero for each new step.
Main Index
CHAPTER 3 33
Capabilities
not be interrupted to perform perturbation analyses: before performing the perturbation analysis, MSC
Nastran SOL 400 requires that the structure be brought into static equilibrium.
In SOL 400, the ANALYSIS case control command may be used to define a linear perturbation analysis STEP.
Pertubation analysis implies a linearized solution about a nonlinear, preloaded state. The preloaded state may
be either a nonlinear static, nonlinear transient analysis, or a thermo-mechanically coupled analysis. When
performing a perturbation analysis in SOL 400 with large displacement activated:
The preload will be performed including large displacements and other material nonlinearities.
The stress obtained will be included in the initial stress stiffness of the perturbation step.
Linear Static
Static stress analysis is used when inertia effects can be neglected. The problem may still have a real time scale,
for example when the material has a viscoplastic response, such as rate dependent yield. The analysis may be
linear or nonlinear.
In linear static analysis, ANALYSIS = STATICS
Linear static analysis involves the specification of load cases and appropriate boundary conditions.
Traditionally, linear static analysis is performed using SOL 101.
Normal Modes
This solution type uses eigenvalue techniques to extract the frequencies of the current system. The stiffness
determined at the end of the previous step is used as the basis for the extraction, so that small vibrations of a
preloaded structure or nonlinearly deformed structure can be modeled.
Main Index
34 MSC Nastran: SOL 400 Getting Started Manual
General Nonlinear Analysis
Nonlinear Static
Nonlinearity may arise from large displacement effects, material nonlinearity and boundary nonlinearity
(such as contact and friction).
ANALYSIS = NLSTATICS
Main Index
CHAPTER 3 35
Capabilities
This behavior requires the solution by a series of increments, with iteration within each increment to obtain
equilibrium. For most cases, the automatic incrementation provided by SOL 400 is preferred, although
direct user control is also provided for those cases where you have experience with a particular problem.
For static analysis, which involves post buckling behavior, where the load-displacement response shows a
negative stiffness, and the structure must release strain energy to remain in equilibrium, an automatic load
incrementation procedure must be used.
For local buckling, a quasi-static damping procedure via the ADAPT option on the NLSTEP bulk data entry
may be used. However, this option only works in conjunction with advanced nonlinear elements. For global
buckling, an arc length based procedure via the ARCLN option on the NLSTEP bulk date entry may be used.
This option only works for non-contact scenarios.
Creep
This analysis procedure performs a transient, static, stress/displacement analysis.
ANALYSIS = NLSTATICS
It is especially provided for the analysis of materials which are described by the MATVP material form. The
time integration method is controlled by (and described under) the NLSTEP bulk data entry.
Viscoelastic
This is especially provided for the time domain analysis of materials which are described by the MATVE
material options.
ANALYSIS = NLSTATICS
The dissipative part of the material behavior is defined through a Prony series representation of the
normalized shear and bulk relaxation moduli. The time integration method is controlled the same as a creep
analysis using (and described under) the NLSTEP bulk data entry.
Main Index
36 MSC Nastran: SOL 400 Getting Started Manual
Nonlinear Transient Response Analysis
The problems can be transient or steady-state, linear or nonlinear. The heat transfer elements allow for heat
storage (specific heat) and heat conduction, and also allow for forced convection caused by fluid flowing
through the mesh. Heat interface elements are also provided, to model the heat transfer across the boundary
layer between a solid and a fluid, or between two closely adjacent solids. Shell-type heat transfer elements are
included, since so many structures are of this type. The second order elements usually give more accurate
results for the same number of nodes in the mesh. Analyses that involve both thermal and mechanical
solutions are referred to in SOL 400 as multi-physics solutions.
Case Control
Each subcase and step defines a time interval starting from the last time step of the previous subcase or
subcase, subdivided into smaller time steps using the NLSTEP entry. The output time is labeled by the
cumulative time, including all previous subcases. The data blocks containing solutions are generated at the
Main Index
CHAPTER 3 37
Capabilities
end of each subcase for storage in the database for output processing and restarts. As such, converged
solutions are apt to be saved at many intermediate steps in case of divergence and more flexible control
becomes possible with multiple subcases. Results from converged increments can be output to *.op2 files
using the NLOPRM case control entry.
The input loading functions may be changed for each subcase or continued by repeating the same DLOAD
request. However, it is recommended that one use the same TLOAD bulk data for all the subcases in order to
maintain the continuity between subcases, because TLOADi data defines the loading history as a function of
cumulative time. Static loads (PLOADi, FORCEi, MOMENTi) may be associated with time-dependent functions
by matching the EXCITEID on the TLOADi entries. Nonlinear forces as functions of displacements or
velocities (NOLINi) may be selected and printed by the case control commands NONLINEAR and NLLOAD,
respectively. Each subcase may have a different time step size, time interval, and iteration control selected by
the NLSTEP request. Case control requests that may not be changed after the first subcase are SPC, MPC,
DMIG, and TF.
Output requests for each subcase are processed independently. Requested output quantities for all the
subcases are appended after the computational process for actual output operation. See Chapter 8 for a
discussion on output requests and see the MSC Nastran Quick Reference Guide for a complete list of output
requests.
Initial conditions (displacement or velocity) can be specified by the bulk data input, TIC, selectable by the
case control command IC. If initial conditions are given, all of the nonlinear element forces and stresses must
be computed to satisfy equilibrium with the prescribed initial displacements. On the other hand, initial
conditions can be generated by applying static analysis for the preload using PARAM,TSTATIC in the first
subcase. Then the transient analysis can be performed in the ensuing subcases. Associated with the adaptive
time stepping method, the PARAM,NDAMP is used to control the stability in the ADAPT method. The NDAMP
parameter represents the numerical damping (a recommended value for usual cases is 0.01), which is often
required to improve the stability and convergence in contact problems.
Main Index
38 MSC Nastran: SOL 400 Getting Started Manual
Coupled Thermal-Mechanical
In nonlinear problems, the mode shapes and frequencies are strong functions of time because of large
displacement effects, so that the above guidelines can be only a coarse approximation. To obtain a more
accurate estimate, repeat the analysis with a significantly different time step (1/5 to 1/10 of the original) and
compare responses.
Coupled Thermal-Mechanical
Traditionally, structural simulations and thermal simulations have been performed independently of one
another in separate analysis codes specialized in solving each physical discipline. In large companies there are
different structural and thermal specialists and departments. These groups interact with each other as the
thermal analysts provide the structural analysts with temperature data for thermal distortion and thermal
stress analysis. For sophisticated systems, the thermal analysis is dependent on the structural deflection, so an
iterative loop is set up to capture all the effects properly. MSC Nastran has had nonlinear structural and heat
transfer capabilities since its inception.
One of the keys to multi-physics simulations is to capture the structural and thermal load path changes caused
by contact. MD Nastran has had structural contact capabilities since its inception. General structural contact
allows mechanical meshes to come in to contact and change the load path. Glued structural contact will weld
mechanical meshes together. Thermal contact is a similar concept but for thermal analysis rather than
structural analysis. In addition, the multi-physics aspect of MSC Nastran allows for coupling thermal and
mechanical contact in the same run. In thermal analysis, thermal contact is available for both steady state and
transient.
For coupled thermal-mechanical analysis the mechanical analysis can be static or dynamic and the thermal
analysis can be steady state or transient. This allows four combinations of coupled thermo mechanical
analysis:
HSTAT-NLSTAT
HSTAT-NLTRAN
HTRAN-NLSTAT
HTRAN-NLTRAN,
where HSTAT, HTRAN, NLSTAT, and NLTRAN stand for steady state heat transfer, transient heat transfer,
structural nonlinear statics, and structural nonlinear transient respectively.
The coupled analysis can account for plasticity and frictional heat coupling and proper updating of interface
conditions when there is relative motion between bodies. The bi-directional coupling is a weakly coupled
approach between thermal and mechanical passes. The available bi-directional coupling schemes are
identified in the figure below and they allow simulation of thermo-mechanical effects associated with large
deformation problems and frictional contact.
Thermal Contact
As described in the multi-physics introduction, one of the keys to multi-physics simulations is to capture the
structural and thermal load path changes caused by contact. MSC Nastran has had structural contact
capabilities since its inception. General structural contact allows mechanical meshes to come in to contact
and change the load path. Glued structural contact will weld mechanical meshes together.
Main Index
CHAPTER 3 39
Capabilities
Thermal contact is a similar concept but for thermal analysis rather than structural analysis. You will be able
to analyze thermal interactions between different contact bodies for the body areas that are in contact and
thermal interactions between contact bodies and the environment for the body areas that are not in contact.
In addition to contact and no contact from the pure mechanical case, there is near contact which allows thermal
interactions between bodies that are getting near to each other, but not yet in real (mechanical) contact.
Thermal contact is available for both steady state and transient analysis.
Benefits
This functionality provides a user-friendly interface to defining thermal contact conditions. The user defines
contact bodies with their thermal properties and contact tables defining the possible contact pairings with
the thermal properties for each pairing. The program automatically identifies the body areas involved in
contact and the body areas exposed to the environment and for each situation, it sets up the appropriate
interface conditions for the thermal interactions.
In a coupled thermal-mechanical analysis, bodies may have relative motions, causing the contact conditions
to change over time. Due to friction forces in the contact interface, heat may be generated adding an extra
heat flux load in the thermal phase of the analysis. All changes in the interface conditions and resulting loads
are updated automatically.
Main Index
Chapter 4: SOL 400 Files
SOL 400 Input File 41
SOL 400 Output File 44
Setting Up a SOL 400 Job 47
Main Index
CHAPTER 4 41
SOL 400 Files
After the generation of the input file is complete, it is submitted for execution as a batch process. Once the
input file has been submitted, you have no additional interaction with MSC Nastran until the job is
complete.
The MSC Nastran Input File, often referred to as the Bulk Data File (.bdf), (or .dat in the MSC Nastran
manuals), is made up of three distinct sections:
Executive Control: It describes the problem or solution type and optional file management.
Case Control: It defines the load history and output requests.
Bulk Data: It defines a detailed model, load and constraint description.
For details, refer to the MSC Nastran Getting Started Guide or the MSC Nastran Quick Reference Guide.
Main Index
42 MSC Nastran: SOL 400 Getting Started Manual
SOL 400 Input File
BEGIN BULK
$ Direct Text Input for Bulk Data
PARAM POST 1
PARAM PRTMAXIM YES
PARAM MRNOECHO 123
BCPARA 0 NLGLUE 1
PARAM LGDISP 1
NLSTEP 1 1.
GENERAL 10 1 10
ADAPT .01 1.-5 .5 4 1.2 0
6 2.-4 Bulk Data Section
MECH PV PFNT
.2
$ Elements and Element Properties for region : solids
PSOLID 1 1 0
PSLDN1 1 1
C8 SOLID L
$ Pset: "solids" will be imported as: "psolid.1"
CHEXA 1 1 1 2 8 7 37
38
44 43
CHEXA 2 1 2 3 9 8 38
39
45 44
Main Index
CHAPTER 4 43
SOL 400 Files
Analysis Forms
Using Patran, you can run a MSC Nastran analysis or you may generate the MSC Nastran Input File to run
externally. For information on generating the MSC Nastran Input file from within Patran, see Chapter 3:
Running an Analysis, Analysis Form in the Patran Interface to MSC Nastran Preference Guide.
SOL 400 uses the Nastran input format which are same as elements, materials, properties, loadings and
formats.
Main Index
44 MSC Nastran: SOL 400 Getting Started Manual
SOL 400 Output File
Main Index
CHAPTER 4 45
SOL 400 Files
While the DISPLACEMENT, SPCFORCES, and STRESS requests are common with most MSC Nastran
solution sequences, the NLSTRESS and BOUTPUT requests are unique to SOL 400 and are required to get
the output unique to nonlinear analyses such as failure indices and contact interaction status, forces and
stresses. Control of these options is available through the Analysis menu job setup interface in Patran (see
Chapter 8 or the Patran MSC Nastran Preference Guide for more information on this topic).
SUBCASE 1
STEP 1
TITLE=This is a default subcase.
ANALYSIS = NLSTATIC
NLSTEP = 1
BCONTACT = ALLBODY
SPC = 2
LOAD = 1
DISPLACEMENT(SORT1,REAL)=ALL
SPCFORCES(SORT1,REAL)=ALL
STRESS(SORT1,REAL,VONMISES,BILIN)=ALL
NLSTRESS(SORT1)=ALL
BOUTPUT(SORT1,REAL)=ALL
MSC recommends to us either the MASTER/DBALL or new OP2 (PARAM,POST,1) output file formats for
SOL 400 analysis. The advantage of using the MASTER/DBALL format is that it contains all of the database
information from the MSC Nastran run and so can be used to retrieve the attributes of the model at a later
time. The disadvantages to this format is that it is platform specific (e.g. not portable across platforms) and
takes up more disk space. The advantages of the new OP2 format is that it takes less disk space than the
MASTER/DBALL and has limited cross-platform portability (see Chapter 2: File Management Statements,
The File Management Section (FMS) of the QRG for more details on portability).
Main Index
46 MSC Nastran: SOL 400 Getting Started Manual
SOL 400 Output File
Control is provided for manipulating the color/range assignment and other attributes for plot tools, and for
controlling, and creating animations of static, and transient results. Results are selected from the database,
and assigned to plot tools using simple forms. Results transformations are provided to derive scalars from
vectors, and tensors as well as to derive vectors from tensors. This allows for a wide variety of visualization
tools to be used with all of the available results.
If the job was created within Patran where a Patran jobname has the same name as the MSC Nastran
jobname, you are only required to use the Results tools and Patran will import or attach the
jobname.xxx file without you having to select it. If you did not create the job in Patran, you can still import
the model and results and postprocess.
Results
Like the enormous amount of data needed to define the simulation model to an analysis code, there is a large
volume of data returned from the simulation analysis. And just as it is virtually impossible to construct a
model with a text editor alone, it is equally difficult to read and interpret the results by hand. Hence, the use
of a postprocessor with a graphical user interface such as Patran, or SimXpert is highly recommended.
The Patran Results application gives you control of powerful graphical capabilities to display results
quantities in a variety of ways:
Deformed structural plots
Color banded fringe plots
Marker plots (vectors, tensors)
Freebody diagrams
Graph (XY) plots
Animations of most of these plot types.
The Results application treats all results quantities in a very flexible and general manner. In addition, for
maximum flexibility results can be:
Sorted
Reported
Filtered
Derived
Deleted
All of the above features help give meaningful insight into results interpretation of engineering problems that
would otherwise be difficult.
The Results application is object oriented, providing postprocessing plots which are created, displayed, and
manipulated to obtain rapid insight into the nature of results data. The imaging is intended to provide
graphics performance sufficient for real time manipulation. Performance will vary depending on hardware,
but consistency of functionality is maintained as much as possible across all supported display devices.
Capabilities for interactive results postprocessing also exist. Advanced visualization capabilities allow creation
of many plot types which can be saved, simultaneously plotted, and interactively manipulated with results
quantities reported at the click of the mouse button to better understand mechanical behavior. Once defined,
Main Index
CHAPTER 4 47
SOL 400 Files
the visualization plots remain in the database for immediate access and provide the means for results
manipulation and review in a consistent and easy to use manner.
For more information, see Chapter 1: Introduction to Results Postprocessing in the Patran Reference Manual.
Solution Type
MSC Nastran SOL 400 can simulate different types of structural and thermal responses. In general, an
structural analysis can be either static or dynamic. Both static, and dynamic analysis may simulate linear
response, or nonlinear response. SOL 400 incorporates the formulations and functionality to simulate
nonlinear static and dynamic structural responses.
Entry Description
ANALYSIS Specifies the analysis procedure to be used for the step or subcase being set up.
Main Index
48 MSC Nastran: SOL 400 Getting Started Manual
Setting Up a SOL 400 Job
Specifying Subcases
Each subcase is designated with the following case control commands:
Entry Description
SUBCASE Delimits and identifies a subcase.
STEP Delimits and identifies a section of a subcase, typically delineating one load step, or
perturbation step in the analysis (see Chapter 4: Solution Strategies for Nonlinear
Analysis, on multi-stepping for a more detailed explanation of how STEPs and
SUBCASEs work together).
Main Index
CHAPTER 4 49
SOL 400 Files
Main Index
Chapter 5: Elements
5 Elements
Introduction 51
Element Classes 51
Automatic Property Mapping 77
Main Index
CHAPTER 5 51
Elements
Introduction
An element is the basic building block of finite element analysis (FEA). There are several basic types of
elements. Which type of element is used for the analysis depends on two factors:
The type of analysis that is going to be performed
The type of object that is to be modeled
An element is a mathematical relation that defines how the degrees of freedom of a node relate to the next.
It also relates how the deflections create stresses. Elements can be:
Lines (trusses or beams)
Areas (2-D or 3-D plates and membranes)
Solids (bricks or tetrahedrals)
The element is used as a mechanism to integrate a physical quantity over a volume, a surface, a curve, or a
point. For the simplest simulation, it provides a transfer function (or impedance) between the degrees of
freedom of one grid and the degrees of freedom of another grid. This is a very general definition which may
be applied to any type of physics including structural analysis, thermal analysis, acoustic analysis, and fluid
dynamics among others.
The selection of the element type and the design of the finite element mesh are important to obtain an
accurate solution. The design of a finite element mesh is done either fully automatically by a mesh generator
available from MSC (Patran, SimXpert) or any other mesh generator available in the market.
In MSC Nastran, there are two aspects of element definition:
1. The definition of the location of the element, by identifying grid points comprising the element. This
is often called the element topology.
2. The definition of the characteristics of the element.
For nonlinear analysis, the topology of the elements, is generally the same for linear analysis (SOL 101) and
for nonlinear analysis (SOL 400), though there are some restrictions. Additionally, there are elements
available in the implicit nonlinear procedure that are not available in the linear solution sequences. For a
detailed description of MSC Nastran element technology, see the MSC Nastran Linear Static Analysis User’s
Guide, MSC Nastran Elements (Ch. 4), MSC Nastran Reference Guide, Structural Elements (Ch. 3) and QRG, Bulk
Data Entries.
Element Classes
In MSC Nastran, the elements are divided into classes based upon their dimensionality and function. In
structural applications such as automotive body, aerospace, and civil engineering, the use of shells, beams,
and rods is prevalent. In applications such as generators, rockets, and pressure vessels, axisymmetric elements
may be used advantageously; while in automotive engines, housing, etc., 3-D solid elements dominate.
MSC Nastran uses Standard and Advanced elements. Standard elements are the native Nastran elements
which are popular for linear static and dynamic analyses. The elements CQUAD, CTRIA, CHEXA,
CPENTA, CTETRA, CBEAM and CBAR are available in both standard and advanced element versions and
Main Index
52 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
these elements provide identical solutions for all linear simulations. The standard elements are defined using
conventional property cards PSHELL, PLPLANE, PSOLID, PBEAM, PBAR, etc.
Advanced elements are activated explicitly by using additional property cards (in addition to conventional
property card) PSHLN1, PSHLN2, PSLDN1, PBEMN1 and PBARN1. While Nastran uses its standard
elements mostly for the linear simulations and also for the limited nonlinear simulations in SOL 106 and
SOL 129, it uses its advanced elements for more advanced nonlinear simulations (large strain material and
contact nonlinear solutions) in SOL 400. Few advanced elements (new type which have no equivalent in
standard element type) are also directly defined using only the conventional property cards (PLCOMP,
PCOMPLS, PCOHE, PAXISYM, etc.).
Nastran also allows the capability of automatic activation of advanced elements when the model uses certain
nonlinear materials or analysis procedures. For more information, refer to See “Automatic Property Mapping”
on page 77.
0-D
0-D elements are a single grid; hence, they do not really have a geometry associated with them. Because of
this, no numerical integration is required.
These elements include:
CELAS1, CELAS2, CELAS3, and CELAS4 (if only one grid identified) – provides a stiffness matrix
CBUSH (if only one grid identified) – provides a stiffness, damping and mass matrix
CDAMP1, CDAMP2, CDAMP3, and CDAMP4 (if only one grid identified) – provides a damping
matrix
CMASS1, CMASS2, CMASS3, and CMASS4 (if only one grid identified) – provides a mass matrix
CONM1 and CONM2 – provides a mass/inertia matrix
Note that only CBUSH allows the ability to change the stiffness due to the deformation, and hence, is more
powerful for nonlinear analysis.
The property options used with these elements are
Main Index
CHAPTER 5 53
Elements
Main Index
54 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
CROD provides a membrane behavior plus twist, but no bending stiffness. The CROD element
is a 2-node element with linear interpolation. One needs to only define the rod cross-
section area.
CBAR provides membrane, bending and torsion behavior. The CBAR is a 2-grid element that
supports linear or nonlinear material behavior. The type of numerical integration along
the length and the numerical integration across the cross section is determined on the
PBARN1.
CBEAM provides general beam behavior. The CEAM is a 2-grid element that will support linear
or nonlinear material behavior. The type of numerical integration along the length and
the numerical integration across the cross section is determined on the PBEMN1.
CBEAM3 provides general beam behavior. The CEAM3 is a 3-grid element with quadratic
interpolation along the length, but only supports linear elastic material.
CBEND is a tube element that may also be used as a curved element. The CBEND is a 2-grid
element that only supports linear elastic material.
Output of Strains
Uniaxial in the truss member.
Output of Stresses
Uniaxial in the truss member.
Main Index
CHAPTER 5 55
Elements
A detailed description of beam cross section and beam orientations is given in the QRG.
Beam cross-section integration – this is substantially different for nonlinear analysis when the
additional property option. In the case of linear elastic material behavior, the cross section of the
beam is integrated to obtain the area, moments of inertia, and the torsional moment. When
nonlinear material models are present, the behavior has to be integrated through the cross section.
Beam offsets – MSC Nastran has two different methods to apply Beam and Shell offsets:
a. Using rigid elements which is the default.
a. Large rotation method activated by the MDLPRM,OFFDEF,LROFF – this is the recommended
approach for a nonlinear analysis.
Beam pin codes – this is the same for linear and nonlinear solutions.
Main Index
56 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
1-D
Auxiliary Large
Number of Property Rotation/ Permits
Number Integration Interpolation Conventional for Large Nonlinear
Element of Grids Points Scheme BEH INT Property Nonlinear Strain Sect Material
CROD 2 1 L ROD L PROD PRODN1 Yes/Yes S Yes
*
CBAR 2 3 LC PROD LC PBAR(L) PBARN1 Yes/No S No
CBAR 2 3 LC PROD LC PBAR(L) PBARN1 Yes/No N Yes
CBAR 2 1 LC PROD LS PBAR(L) PBARN1 Yes/No S No
CBAR 2 1 LC PROD LS PBAR(L) PBARN1 Yes/No N Yes
CBEAM 2 3 LC BEAM LC PBEAM(L) PBEMN1 Yes/No S No
CBEAM 2 3 LC BEAM LC PBEAM(L) PBEMN1 Yes/No N Yes
CBEAM 2 1 LC BEAM LS PBEAM(L) PBEMN1 Yes/No S No
CBEAM 2 1 LC BEAM LS PBEAM(L) PBEMN1 Yes/No N Yes
CBEAM 2 2 LC BEAM LCC PBEAML PBEMN1 Yes/No N Yes
CBEAM 2 2 LC BEAM LCO PBEAML PBEMN1 Yes/No N Yes
* When a PRODN1 is used with a CROD to permit nonlinear material behavior, the element behavior changes such that it no longer supports torsion.
These elements cannot be used with the Hill, Barlat, Linear Mohr-Coulomb, Parabolic Mohr-Coulomb, or
the IMPLICIT CREEP model specified on the MATEP option.
The user can change the Integration along the length when using PBARN1 or PBEMN1. The choices are:
Main Index
CHAPTER 5 57
Elements
Linear/cubic means that linear interpolation of the displacement is used along the axis and cubic
displacement variation normal to the beam axis. This results in linear variation of curvature.
The CBAR, PBAR, and PBARN1 does not support LCO which requires element warping.
The cross-section behavior can either be specified by entering an S or N on the SECT option, where N means
numerically integrated and S means smeared. When a nonlinear material behavior needs to be captured, you
should use the N option.
The resultant quantities are always given with respect to an element axis attached to the 1-D element. When
large displacement is included, this is co-rotated with the element.
3 3
3
6 5
1 2
1 2 1 4 2
3-node triangular element, 6-node element,
1-point integration point 3-point integration point
Main Index
58 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
4 3
4 3
3 4
1 2
1 2 1 2
4-node quadrilateral, 4-node,
4-point integration point 1-point integration point
4 7 3 4 7 3
7 8 9
3 4
8 4 5 6 6 8 6
1 2
1 2 3
1 5 2 1 5 2
8-node quadrilateral, 8-node,
9-point integration point 4-point integration point
Main Index
CHAPTER 5 59
Elements
Plane Stress
Number Auxiliary Large Permits
of Property Rotation/ Nonlinea
Number Integratio Interpolatio Conventiona for Large r
Element of Grids n Points n Scheme BEH INT l Property Nonlinear Strain Material*
CTRIA3 3 1 L PSTRS L PLPLANE PSHLN2 Yes/Yes Yes
CQUAD4 4 4 L PSTRS L PLPLANE PSHLN2 Yes/Yes Yes
CQUAD4 4 1 L PSTRS LRIH PLPLANE PSHLN2 Yes/Yes Yes
CTRIA6 6 7 Q PSTRS Q PLPLANE PSHLN2 Yes/Yes Yes
CQUAD8 8 9 Q PSTRS Q PLPLANE PSHLN2 Yes/Yes Yes
CQUAD8 8 4 Q PSTRS QRI PLPLANE PSHLN2 Yes/Yes Yes
* With exceptions
Output of Strains
1 = xx
2 = yy
3= zz
4 = xy
Output of Stresses
1 = xx
2 = yy
3= zz = 0
4 = xy
These elements cannot be used with the IMPLICIT CREEP model specified on the MATEP option.
Main Index
60 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
Plane Strain
Number Auxiliary Large
of Property Rotation/ Permits
Number Integratio Interpolatio Convention for Large Nonlinear
Element of Grids n Points n Scheme BEH INT al Property Nonlinear Strain Material
CTRIA3 3 1 L PLSTRN L PLPLANE PSHLN2 Yes/Yes Yes
CTRIA3 4 3 L & Cubic IPS L PLPLANE PSHLN2 Yes/Yes Yes
Bubble
CQUAD4 4 4 L PLSTRN L PLPLANE PSHLN2 Yes/Yes Yes
CQUAD4 4 1 L PLSTRN LRIH PLPLANE PSHLN2 Yes/Yes Yes
CTRIA6 6 7 Q PLSTRN Q PLPLANE PSHLN2 Yes/Yes Yes
CQUAD8 8 9 Q PLSTRN Q PLPLANE PSHLN2 Yes/Yes Yes
CQUAD8 8 4 Q PLSTRN QRI PLPLANE PSHLN2 Yes/Yes Yes
Output of Strains
1 = xx
2 = yy
3= zz = 0
4 = xy
Output of Stresses
1 = xx
2 = yy
3= zz
4 = xy
Caution: The conventional CTRIA element with BEH=PLSTRN is known to give very poor results when
used with incompressible or nearly incompressible behavior (including rubber materials, elastic-
plastic or creep), the BEH=IPS should be used.
Main Index
CHAPTER 5 61
Elements
Axisymmetric
Auxiliary Large
Number of Property Rotation/ Permits
Number Integration Interpolation Convention for Large Nonlinear
Element of Grids Points Scheme BEH INT al Property Nonlinear Strain Material
CTRIAX 3 1 L AXISOLID L PLPLANE PSHLN2 Yes/Yes Yes
CTRIAX 4 4 L & Cubic IAX L PLPLANE PSHLN2 Yes/Yes Yes
Bubble
CQUADX 4 4 L AXISOLID L PLPLANE PSHLN2 Yes/Yes Yes
CQUADX 4 1 L AXISOLID LRIH PLPLANE PSHLN2 Yes/Yes Yes
CTRIAX 6 7 Q AXISOLID Q PLPLANE PSHLN2 Yes/Yes Yes
CQUADX 8 9 Q AXISOLID Q PLPLANE PSHLN2 Yes/Yes Yes
CQUADX 8 4 Q AXISOLID QRI PLPLANE PSHLN2 Yes/Yes Yes
Output of Strains
1 = rr
2 = zz
3=
4 = rz
Output of Stresses
1 = rr
2 = zz
3 =
4 = rz
To apply a distributed load on this element one needs to apply the pressure using PLOADX1. Note that the
distributed loads are integrated over one radian. Hence, any point force applied should also be prescribed
over one radian.
Main Index
62 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
Caution: The conventional CTRIA element with BEH=AXISOLID is known to give very poor results when
used with incompressible or nearly incompressible behavior (including rubber materials, elastic-
plastic or creep), the BEH=IAX should be used.
Output of Strains
1 = rr
2 = zz
3=
4 = rz
5 = z
6 = z
Output of Stresses
1 = rr
2 = zz
3 =
4 = rz
5 = z
6 = z
Main Index
CHAPTER 5 63
Elements
Note the displacements for this element are U r , U z , and U which represents the angular displacement
about the symmetry axis measured in radians. This is given as the third degree of freedom, and all torques or
single points constraints on this quantity should be treated as the third degree of freedom.
Axisymmetric Shell
Number of Large Permits
Number of Integration Property for Rotation/ Large Nonlinear
Element Grids Points Nonlinear Strain Material
CAXISYM 2 2 PAXISYM Yes/Yes Yes
CAXISYM 3 3 PAXISYM Yes/Yes Yes
Output of Strains
1 = s = meridional membrane
2 = = circumferential membrane
3 = t = transverse shear strain
Output Of Stresses
1 = s = meridional stress
2 = = circumferential stress
3 = t = transverse shear stress
The degrees of freedom for this element are:
Ur , Uz , Ut
Main Index
64 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
Main Index
CHAPTER 5 65
Elements
Output of Strains
1 = xx
2 = yy
3= xy
Output of Stresses
1 = xx
2 = yy
3 = xy
Shells
There are two classes of advanced shell elements in MSC Nastran SOL 400.
The first group (marked LDK) is thin shell elements based upon Kirchhoff theory. The second group
(marked L, LRIH, and QRI) is thick shell elements that support transverse shears based upon Mindlin
theory. There are the preferred elements for composite simulation.
3-D Shell
Auxiliary Large
Number of Property Rotation/ Permits
Number Integration Interpolation Convention for Large Nonlinear
Element of Grids Points Scheme BEH INT al Property Nonlinear Strain Material*
CTRIA3 3 1 L DCTN LDK PSHELL or PSHLN1 Yes/Yes Yes
PCOMP or
PCOMPG
CQUAD4 4 4 L DCT L PSHELL or PSHLN1 Yes/Yes Yes
PCOMP or
PCOMPG
CQUAD4 4 1 L DCT LRI PSHELL or PSHLN1 Yes/Yes Yes
H PCOMP or
PCOMPG
CQUAD4 4 4 L DCTN LDK PSHELL or PSHLN1 Yes/Yes Yes
PCOMP or
PCOMPG
CQUAD8 8 4 Q DCT QRI PSHELL or PSHLN1 Yes/Yes Yes
PCOMP or
PCOMPG
* With exceptions
Main Index
66 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
The output of strains and stresses are given with respect to a coordinate system attached to the element. This
coordinate system is updated with the deformation if the LGDISP parameter is activated.
Output of Strains
1 = xx
2 = yy
3 = zz
4 = xy
5 = yz for thick shell only
6 = zx for thick shell only
Output of Stresses
1 = xx
2 = yy
3 = zz = 0
4 = xy
5 = yz for thick shell only
6 = zx for thick shell only
These elements cannot be used with the IMPLICIT CREEP model specified on the MATEP option.
Shear Panel
MSC Nastran supports a 4-node shear panel element for linear analysis. Here, linear analysis means small
deformation and linear isotropic elastic material.
Main Index
CHAPTER 5 67
Elements
While SOL 400 has an alternative formulation when PSHEARN is entered, it provides a membrane
formulation to the element. It is not recommended that one uses the PSHEARN option in SOL 400.
Figure 5-1 Connectivity and Integration Points for Composite and Non Composite Materials
Main Index
68 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
Output of Strains
1 = xx
2 = yy
3 = zz
4 = xy
5 = yz
6 = zx
Output of Stresses
1 = xx
2 = yy
3 = zz
4 = xy
5 = yz
6 = zx
For a homogeneous material, use PSLDN1.
Set BEH=SOLID on PSLDN1
For a layered composite material, use PCOMPLS. Note in the Figure 5-1 that for each layer, there are three
layer points, which allows an accurate calculation of the inter laminar shear.
Large Permits
Number of Property for Rotation/ Nonlinear
Element Grids Integration – INT Nonlinear Large Strain Material
CHEXA 8 ASTN PCOMPLS Yes /Yes Yes
1 per layer
Assumed Strain
Main Index
CHAPTER 5 69
Elements
advanced elements include lower- and higher-order tetrahedral, pentahedral, and hexahedral elements with
a variety of integration schemes that have the following geometry.
Output of Strains
1 = xx
2 = yy
3 = zz
4 = xy
5 = yz
6 = zx
Output of Stresses
1 = xx
2 = yy
3 = zz
4 = xy
5 = yz
6 = zx
Main Index
70 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
4 4
10
3
8
7 9
6
1 1
3 2 2
X
4
3
X
3
X
2
X
1
2
4-point integration
6 G6
+6
+4 G15 G14
5 G12
+5
G13
G4 G5
+3
+1 3
G10 G3
2 G11
1 + G9 G8
G7
G1
2 G2
6-node pentahedral, 15-node pentahedral,
6-point integration 21-point integration
Main Index
CHAPTER 5 71
Elements
8 6
5
5
7
5 7
8 2
6 6 7
1 1
3 4
1 8
1 4
2
3 3
2 4
G18
6 G7
13 G6
5
14
16 G19
15 7
G15 G17
8 18 G14
17 G20
G8
G2
19 G10 G5
20 2 G16 G3
9
1 G11 G13
10 G9
12
3 G12
11 G4
4
G1
20-node hexahedral, 20-node hexahedral,
27-point integration 8-point integration
Main Index
72 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
3-D Solid
Auxiliary Large
Number of Property Rotation/ Permits
Number Integration Interpolation Convention for Large Nonlinear
Element of Grids Points Scheme BEH INT al Property Nonlinear Strain Material
CTETRA 4 1 L SOLID L PSOLID PSLDN1 Yes/Yes Yes
CTETRA 4 4 L & CUBIC ISOL L PSOLID PSLDN1 Yes/Yes MATEP
CPENTA 6 6 L SOLID L PSOLID PSLDN1 Yes/Yes Yes
CHEXA 8 8 L SOLID L PSOLID PSLDN1 Yes/Yes Yes
CHEXA 8 1 L SOLID LRIH PSOLID PSLDN1 Yes/Yes Yes
CTETRA 10 9 Q SOLID Q PSOLID PSLDN1 Yes/Yes Yes
CTETRA 10 4 Q SOLID LRIH PSOLID PSLDN1 Yes/Yes Yes
CPENTA 15 21 Q SOLID Q PSOLID PSLDN1 Yes/Yes Yes
CHEXA 20 27 Q SOLID Q PSOLID PSLDN1 Yes/Yes Yes
CHEXA 20 8 Q SOLID QRI PSOLID PSLDN1 Yes/Yes Yes
Caution: The conventional CTETRA element with BEH=SOLID is known to give very poor results when
used with incompressible or nearly incompressible behavior (including rubber materials, elastic-
plastic or creep), the BEH=ISO should be used.
z
ply 2
ply 3
ply 1
ply 1
ply 3
ply 2
ply 2
ply 1 3 3 3
1 2 1 x 2 1 2
Main Index
CHAPTER 5 73
Elements
Large Permits
Number of Integration – Property for Rotation/ Nonlinear
Element Grids INT Nonlinear BEH Large Strain Material
Plane Strain
CQUAD4 4 L 2 per layer PLCOMP COMPS Yes/Yes Yes
CQUAD8 8 Q 2 per layer PLCOMP COMPS Yes/Yes Yes
Axisymmetric
CQUADX 4 L 2 per layer PLCOMP AXCOMP Yes/Yes Yes
CQUADX 8 Q 2 per layer PLCOMP AXCOMP Yes/Yes Yes
Hexahedral
CHEXA 8 L 4 per layer PCOMPLS SLCOMP Yes/Yes Yes
CHEXA 20 Q 4 per layer PCOMPLS SLCOMP Yes/Yes Yes
Gasket Elements
The lower-order continuum composite elements are also used to model gasket materials. When used in this
manner, the number of layers is one and the material is defined through the MATG option. Note that these
elements may be collapsed to a pentahedral so they can model the gasket between two regions modeled with
tetrahedral elements. Care should be exercised to make sure the collapsing occurs in the plane of the gasket
material. It is not necessary that the elements match the mesh in the surrounding material. One may utilize
the contact capability to overcome the mesh incompatibility.
Large
Auxiliary Rotation/ Permits
Number Conventiona Property for Large Nonlinear
Element of Grids Integration – INT l Property Nonlinear BEH Strain Material
Plane Strain
CQUAD4 4 L 2 per layer PLPLANE PSHNL2 COMPS Yes/No MATG
Axisymmetric
CQUADX 4 L 2 per layer PLPLANE PSHNL2 AXCOMP Yes/No MATG
Hexahedral
CHEXA 8 L 4 per layer PSOLID PSLDN1 SLCOMP Yes/No MATG
Interface Elements
There are a series of elements that are used to model the onset and progression of delamination of the
bonding materials using the Cohesive Zone Method. These elements are available for plane strain,
axisymmetric, and three-dimensional behavior. From a meshing perspective, these elements are unique
Main Index
74 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
because one can enter a zero thickness. The interface elements provide two integration schemes. The first uses
the conventional Gaussian integration scheme while the other uses a nodal lumping scheme (Lobatto-Cotes).
The latter scheme may be advantageous when the interface material is relatively stiff compared to the
surrounding material.
Note that the orientation of the element dictates the direction of the interface/delamination. It is not
necessary that the elements match the mesh in the surrounding material. One may utilize the contact
capability to overcome the mesh incompatibility. The material properties are defined using the MCOHE
material model Cohesive Zone Modeling (MCOHE). Note that delamination simulations are highly
nonlinear and one must exercise caution in applying the boundary conditions.
The element is written with respect to a local coordinate system, relating the relative displacement and the
normal and shear traction. The output is given with respect to the local coordinate system. This system is
updated (rotated) when large displacement is used. The element does not include a mass matrixor a geometric
or initial stress stiffness matrix. It also does not support application of distribute loads.
The higher-order elements are not fully quadratic; they are quadratic in the plane of the interface but linear
through the thickness.
The elements are shown below.
v 3
˜1
2 3
˜1 6
2
7
v 5
˜2
˜2
4 4
8
1 1
4-node linear planar interface element 8-node quadratice planar interface element
4 3
1 1 2 2
4 7 3
4 3 8 6
1 2 1 5 2
1 2 Element mid-line
Gaussian integration scheme
Nodal lumping scheme
1 2 3 1 2 3
Gauss Integration Scheme Newton-Cotes/Lobatto Integration Scheme
Main Index
CHAPTER 5 75
Elements
8 8
16
4 4 19
20
v v1 11
˜1 v 7 12 7
˜3 ˜ v3 15
3 ˜ 3
5 v 18
˜2 13
5
v2
1 1
17 ˜ 10
6 9 6
14
2 2
8
4
7
5 3
1
6
2
Element Mid-plane
3
3
4 4
1
2
1
Gauss Integration Scheme
2
Newton-Cotes/Lobatto Integration Scheme
Main Index
76 MSC Nastran: SOL 400 Getting Started Manual
Element Classes
8
4
7
5 3
1
6
2
Element Mid-plane
7
7 8
8
4 9
5 4 9
1 6
2 5
3
1 6
Gauss Integration Scheme 2
3
Newton-Cotes/Lobatto Integration Scheme
Large Permits
Number of Property for Rotation/ Nonlinear
Element Grids Integration – INT Nonlinear Large Strain Material
Plane Strain
CIFQUAD 4 2 PCOHE Yes/No MCHOE
L
CIFQUAD 8 2 PCOHE Yes/No MCHOE
Q/L
Axisymmetric
CIFQDX 4 2 PCOHE Yes/No MCHOE
L
CIFQDX 8 2 PCOHE Yes/No MCHOE
Q/L
Solid
CIFPENT 8 3 PCOHE Yes/No MCHOE
L
Main Index
CHAPTER 5 77
Elements
Large Permits
Number of Property for Rotation/ Nonlinear
Element Grids Integration – INT Nonlinear Large Strain Material
CIFPENT 15 6 PCOHE Yes/No MCHOE
Q/L
CIFHEX 8 4 PCOHE No MCHOE
L
CIFHEX 20 8 PCOHE No MCHOE
Q/L
Output of Strain
The three strain components are given at the element integration points. They are determined by the relative
displacements between the top and bottom face and are given in the local element system:
1 = u top – u bottom
˜ ˜
2 = v top – v bottom
˜ ˜
3 = w top – w bottom for 3-D element only
˜ ˜
Output of Stress
1 = n
2 = s1
3 = s2 for 3-D element only
Main Index
78 MSC Nastran: SOL 400 Getting Started Manual
Automatic Property Mapping
Unsupported
Auxiliary Conventional Features of
Property Property Dimension Material Notes Primary Entry
PBARN1 PBARL 1-D MAT4 Note 4
MATS1 Note 1
MATEP
MATF
MATSMA Note 2
MATVE
MATVP Note 3
PBEMN1 PBEAML 1-D MAT4 Note 4 Tapered Sections
MATS1 Note 1
MATEP
MATF
MATSMA Note 2
MATVE
MATVP Note 3
PRODN1 PROD 1-D MAT4 Note 4 J,C
MATS1 Note 1
MATEP
MATF
MATSMA Note 2
MATVE
MATVP Note 3
PSHEARN PSHEAR 3-D MAT4 Note 4 F1,F2
MAT8
MATS1 Note 1
MATS8 Note 1
MATEP
MATF
MATORT
MATSMA Note 2
MATVE
MATVP Note 3
Main Index
CHAPTER 5 79
Elements
Unsupported
Auxiliary Conventional Features of
Property Property Dimension Material Notes Primary Entry
PSHLN1 PSHELL 3-D MAT4 Note 4 TS/T, nondefault Z1
and Z2, 121/T3
MAT5 Note 4
MATS1 Note 1
MATS8 Note 1
MATEP
MATF
MATORT
MATSMA Note 2
MATVE
MATVP Note 3
PSHLN1 PCOMP/ 3-D MAT4 Note 4 FT, GE, LAM options
PCOMPG other than BLANK
MAT5 Note 4 and SYM, SOUTi
MAT8
MATS1 Note 1
MATS2 Note 1
MATS8 Note 1
MATEP
MATF
MATORT
MATSMA Note 2
MATVE
MATVP Note 3
Main Index
80 MSC Nastran: SOL 400 Getting Started Manual
Automatic Property Mapping
Unsupported
Auxiliary Conventional Features of
Property Property Dimension Material Notes Primary Entry
PSHLN2 PLPLANE 2-D MAT4 Note 4
MAT5 Note 4
MATG Note 5
MATS1 Note 1
MATS3 Note 1
MATS8 Note 1
MATEP
MATF
MATORT
MATSMA Note 6
MATVE
MATVP Note 7
PSLDN1 PSOLID 3-D MAT4 Note 4 IN, ISOP, FCTN
MAT5 Note 4
MATS1 Note 1
MATEP
MATF
MATORT
MATSMA
MATVE
MATVP
Main Index
CHAPTER 5 81
Elements
For shape memory materials, define through MATSMA and BEH=PLSTRS only the
thermo-mechanical model is available.
For creep, define through MATVP and BEH=PLSTRS VALC=0 must be set on
NLMOPTS for explicit formulation.
Main Index
82 MSC Nastran: SOL 400 Getting Started Manual
Automatic Property Mapping
Main Index
Chapter 6: Materials
6 Materials
Material Model Overview 84
Material Property Definitions 84
Linear Elastic Behavior 88
Linear Elastic Materials 89
Viscoelastic 93
Elasto-plastic Behavior 97
Elastoplastic Material Entries 99
Strain Rate Dependent Yield 101
Creep (MATVP, CREEP) 103
Composite (PCOMP or PCOMPG) 107
Cohesive Zone Modeling (MCOHE) 113
Progressive Composite Failure 113
Main Index
84 MSC Nastran Getting Started Guide
Material Model Overview
Linear Deformations are directly proportional to the applied load (i.e., strain is
directly proportional to stress).
Elastic An elastic structure returns to its original, undeformed shape when the load is
removed.
Homogeneous The material is the same throughout — material properties are independent
of location within the material.
Isotropic Material properties do not change with the direction of the material.
Modulus of elasticity E is the constant of proportionality relating stress-to-strain for uniaxial
(Young’s modulus) E behavior in the linear region. The greater the value of E, the stiffer the
material.
Shear modulus (Modulus G is the constant of proportionality relating shear stress to shear strain in the
of rigidity) G linear region.
Poisson’s ratio It is the absolute value of the ratio of lateral linear strain to axial linear strain.
Main Index
Chapter 6: Materials 85
Material Property Definitions
Modeling of nonlinear material behavior is a critical component to obtain the structural response of
structures; especially, when either the loads are large, the temperatures are high, and/or the non-metallic
material is used. MSC Nastran SOL 400 provides a variety of models and, in some cases, multiple numerical
implementation of these models. This latter is due to the historical incorporation of solution sequences 106
and 129. The newer implementations of material models is based upon incorporating technology from Marc
and is the preferred approach. The newer material formulation also requires the use of the newer element
methodology which is often labeled the “advanced” elements. Aspects of these elements are described in
Chapter 11: Element Library. Hence, to use the advanced material models, it is necessary to supplement the
traditional property entries with new property entries (PBARN1, PBEMN1, PRODN1, PSHLN1,
PSHLN2, and PSLDN1). An alternate approach is to use the NLMOPTS bulk data option with the
PROPMAP keyword.
Main Index
86 MSC Nastran Getting Started Guide
Material Property Definitions
Main Index
Chapter 6: Materials 87
Material Property Definitions
Main Index
88 MSC Nastran Getting Started Guide
Linear Elastic Behavior
The following sections describe how to model material behavior in SOL 400. Modeling material behavior
consists of both specifying the constitutive models used to describe the material behavior and defining the
actual material data necessary to represent the material. Directional dependency can be included for materials
other than isotropic materials. Data for the materials can be entered either directly through the input file or
by user subroutines, or material models may be defined in the Patran Materials Application. Each section of
this chapter discusses various options for organizing material data for input. Each section also discusses the
constitutive (stress-strain) relation and graphic representation of the models and includes recommendations
and cautions concerning the use of the models.
Strain
Experiments show that axial elongation is always accompanied by lateral contraction of the bar. The ratio for
a linear elastic material is:
v = (lateral contraction)/(axial elongation)
This is known as Poisson’s ratio. Similarly, the shear modulus (modulus of rigidity) is defined as:
G (shear modulus) = (shear stress)/(shear strain) (6-2)
Main Index
Chapter 6: Materials 89
Linear Elastic Materials
Material Characteristics
Isotropic (MAT1) Isotropic material – No Directional Dependency – most commonly
used material property
Two Dimensional Anisotropic General anisotropic two-dimensional material used with plate and
(MAT2) shell elements
Axi-symmetric / Solid Orthotropic Orthotropic three-dimensional material for use with CTRIAX6,
(MAT3) PSHLN2, or PLCOMP
Two Dimensional Orthotropic Two-dimensional orthotropic stress-strain used with the plate and
(MAT8) shell elements
Three-Dimensional Anisotropic Orthrotropic material property for the CHEXA, CPENTA, and
Material (MATORT) CTETRA solid elements and solid shell and plane strain element
Isotropic Materials
Most linear elastic materials are assumed to be isotropic (their elastic properties are the same in all directions).
For an isotropic material, every plane is a plane of symmetry and every direction is an axis of symmetry. It
can be shown that for an isotropic material:
G = E 21 + v (6-3)
The shear modulus G can be easily calculated if the modulus of elasticity E and Poisson’s ratio v
are known.
Main Index
90 MSC Nastran Getting Started Guide
Linear Elastic Materials
Entry Description
MATBV Defines the material properties for linear isotropic materials.
The material density, used to define the mass of the structure, and the damping value are used in dynamic
loadings, while the expansion coefficient is used to identify the thermal strains.
Orthotropic Materials
An orthotropic material has three mutually orthogonal planes of symmetry. With respect to a coordinate
system parallel to these planes, the constitutive law for this material is given by the following more general
form of Hooke’s Law:
Main Index
Chapter 6: Materials 91
Linear Elastic Materials
11 1 E1 – 12 E 1 – 13 E 1 0 0 0 11
22 – 12 E 1 1 E2 – 23 E 2 0 0 0 22
33 – 13 E 1 – 23 E 2 1 E3 0 0 0 33
=
12 0 0 0 1 G 12 0 0 12
23 0 0 0 0 1 G 23 0 23
13 0 0 0 0 0 1 G 13 13
3-D Orthotropic
Due to symmetry of the compliance matrix, 11 21 = 22 12 , 22 32 = 33 23 , and 33 13 =
11 31 . Using these relations, a general orthotropic material has nine independent constants:
11 , 22 , 33 , 12 , 23 , 31 , G 12 , G 23 , G 31
Note: The inequalities E22 > 23 E33, E11 > 12 E22, and E33 > 31 E11 must be satisfied in
order for the orthotropic material to be stable. This is checked by MSC Nastran
Implicit Nonlinear.
2-D Orthotropic
Orthotropic material models can be used with 2-D elements, such as plane stress, plane strain, and
axisymmetric elements. For example, the orthotropic stress-strain relationship for a plane stress element is:
E1 21 E 1 0
1
C = ------------------------------ 12 E 2 E2 0 (6-4)
1 – 12 21
0 0 1 – 12 21 G
Main Index
92 MSC Nastran Getting Started Guide
Linear Elastic Materials
Entry Description
MAT3 Defines the material properties for linear orthotropic materials used by the
CTRIAX6 element entry.
MAT2 Defines the material property for an orthotropic material for solids and
isoparametric shell elements.
MAT8
MATORT Specifies elastic orthotropic material properties for three-dimensional and plane
strain behavior for linear and nonlinear analyses in MSC Nastran Implicit
Nonlinear in a more general way than MAT2 or MAT8.
Anisotropic Materials
Anisotropic material exhibits different elastic properties in different directions. The significant directions of
the material are labeled as preferred directions, and it is easiest to express the material behavior with respect
to these directions.
Main Index
Chapter 6: Materials 93
Viscoelastic
The stress-strain relationship for an anisotropic linear elastic material can be expressed as
ij = C ijkl kl (6-5)
The values of C ijkl (the stress-strain relation) and the preferred directions (if necessary) must be defined for
an anisotropic material.
Entry Description
MAT2 Defines the material properties for linear anisotropic materials for two-
dimensional elements.
Viscoelastic
Viscoelasticity models rate effects in the deformation of materials. Their behavior becomes time dependent
and typical phenomena associated with this behavior are relaxation and creep. Relaxation is the diminishing
of stress at constant deformation level (figure 1). Creep is the progression of deformation at constant load
level (figure 2)
Main Index
94 MSC Nastran Getting Started Guide
Viscoelastic
The data required to perform this analysis are the static material constants like Young's moduls and Poisson's
ratio or the material constants used in the strain energy function of a hyperelastic material. In a viscoelastic
analysis these material constants correspond to the instantaneous or short term behavior of the material. In
addition viscoelastic material data are needed to describe the rate effects in the material.
Viscoelastic analysis can be performed in the time domain or in the frequency domain.
In time domain viscoelasticity we study the transient response due to arbitrary time varying loads to capture
the rate effects in the deformation process.
In the time domain two models are available:
The hereditary integral model for linear viscoelasticity and the Simo model for nonlinear
viscoelasticity. The Simo model is an extension of the linear model to finite strain viscoelasticity and
leads to similar hereditary integrals as the linear model.
The Bergstrom-Boyce (B-B) model, which is a phenomenological highly nonlinear model.
In frequency domain viscoelastity we study the stationary dynamic response due to harmonic loads to capture
the rate effects in the deformation process. The stationary response due to harmonic loads that are applied
with a certain excitation frequency is a response with the same frequency. Frequency domain viscoelasticity
is always based on a linear perturbation around a static equilibrium state, which can be the undeformed stress-
free state or some other deformed equilibrium state and is therefore limited to small amplitude vibrations
around the static equilibrium state. Frequency domain viscoelasticity leads to the concepts of storage modulus
and loss modulus, which characterize the frequency dependent stiffness and damping properties of a material.
Within the considered amplitude range, the size of the amplitude may still have a substantial effect on the
harmonic response, which is known as the Payne effect or Fletcher-Gent effect. With the Payne effect, storage
modulus and loss modulus of a material also become vibration amplitude dependent, turning the frequency
Main Index
Chapter 6: Materials 95
Viscoelastic
response problem into a nonlinear problem. Without the Payne effect the frequency response problem is
always linear.
In the frequency domain the hereditary integral models for linear viscoelasticity and nonlinear viscoelasticity
(Simo model) are available only. The Bergstrom-Boyce model is not available in the frequency domain.
Because of the presence of damping there is energy dissipation, which in general will be frequency dependent
Viscoelasticity can be applied:
To determine the current state of deformation based on the entire time history of loading.
To include temperature dependencies.
To use in conjunction of isotropic, orthotropic, and anisotropic material for small strain problems.
To use in conjunction with rubber or foam material for large strain problems.
Main Index
96 MSC Nastran Getting Started Guide
Viscoelastic
one at constant temperature in real time. If the shifted time is used, however, the transformed viscoelastic
equilibrium and compatibility equations are not equivalent to the corresponding elastic equations.
In the case where the temperature varies with time, the extended constitutive law implies a nonlinear
dependence of the instantaneous stress state at each material point of the body upon the entire local
temperature history. In other words, the functionals are linear in the strains but nonlinear in the temperature.
The time scale of experimental data is extended for Thermo-Rheologically Simple materials. All characteristic
functions of the material must obey the same property. The shift function is a basic property of the material
and must be determined experimentally. As a consequence of the shifting of the mechanical properties data
parallel to the time axis (see Figure 6-3), the values of the zero and infinite frequency complex moduli do not
change due to shifting. Hence, elastic materials with temperature-dependent characteristics neither belong to
nor are consistent with the above hypothesis for the class of Thermo-Rheologically Simple viscoelastic solids.
T0
f(T1)
T1
T2 f(T2)
GT
ln t
Figure 6-3 Relaxation Modulus vs. Time at Different Temperatures
Narayanaswamy Model
The annealing of flat glass requires that the residual stresses be of an acceptable magnitude, while the
specification for optical glass components usually includes a homogenous refractive index. The design of heat
treated processes (for example, annealing) can be accomplished using the Narayanaswamy model. This allows
you to study the time dependence of physical properties (for example, volumes) of glass subjected to a change
in temperature.
The glass transition is a region of temperature in which molecular rearrangements occur on a scale of minutes
or hours, so that the properties of a liquid change at a rate that is easily observed. Below the glass transition
temperature T g , the material is extremely viscous and a solidus state exists. Above T g , the equilibrium
structure is arrived at easily and the material is in liquidus state. Hence, the glass transition is revealed by a
change in the temperature dependence of some property of a liquid during cooling. If a mechanical stress is
applied to a liquid in the transition region, a time-dependent change in dimensions results due to the
phenomenon of visco-elasticity.
If a liquid in the transition region is subjected to a sudden change in temperature, a time-dependent change
in volume occurs as shown in Figure 6-4. The latter process is called structural relaxation. Hence, structural
relaxation governs the time-dependent response of a liquid to a change of temperature.
Main Index
Chapter 6: Materials 97
Elasto-plastic Behavior
T1
T(t)
T2
t0 t
(a) Step Input for Temperature
V(0,T1)
g(T2-T1)
V(0,T2) l(T2-T1)
V(,T2)
t0 t
(b) Volume Change as Function of Temperature
Entry Description
MATVE Specifies isotropic viscoelastic material properties to be used for quasi-static or
dynamic analysis in SOL 400.
MATTVE Specifies temperature-dependent visco-elastic material properties in terms of
Thermo-Rheologically Simple behavior to be used for quasi-static or transient
dynamic analysis in SOL 400.
Elasto-plastic Behavior
Most materials of engineering interest initially respond elastically. Elastic behavior means that the
deformation is fully recoverable, so that, when the load is removed, the specimen returns to its original shape.
If the load exceeds some limit (the “yield load”), the deformation is no longer fully recoverable. Some parts
of the deformation will remain when the load is removed as, for example, when a paper clip is bent too much,
Main Index
98 MSC Nastran Getting Started Guide
Elasto-plastic Behavior
or when a billet of metal is rolled or forged in a manufacturing process. Plasticity theories model the material’s
mechanical response as it undergoes such nonrecoverable deformation in a ductile fashion. The theories have
been developed most intensively for metals, but they are applied to soils, concrete, rock, and ice. These
materials behave in very different ways (for example, even large values of pure hydrostatic pressure cause very
little inelastic deformation in metals, but quite small hydrostatic pressure may cause a significant, non-
recoverable volume change in a soil sample), but the fundamental concepts of plasticity theories are
sufficiently general that models based on these concepts have been successfully developed for a wide range of
materials. A number of these plasticity modes are available in the SOL 400 material library.
In nonlinear material behavior, the material parameters depend on the state of stress. Up to the proportional
limit, i.e., the point at which linearity in material behavior ceases, the linear elastic formulation for the
behavior can be used. Beyond that point, and especially after the onset of yield, nonlinear formulations are
required. In general, two ingredients are required to ascertain material behavior:
1. an initial yield criterion to determine the state of stress at which yielding is considered to begin
2. mathematical rules to explain the post-yielding behavior
There are two major theories of plastic behavior that address these criterion differently. In the first, called
deformation theory, the plastic strains are uniquely defined by the state of stress. The second theory, called flow
or incremental theory, expresses the increments of plastic strain (irrecoverable strains) as functions of the
current stress, the strain increments, and the stress increments. Incremental theory is more general and can
be adapted in its particulars to fit a variety of material behaviors. The plasticity models in SOL 400 are
“incremental” theories, in which the mechanical strain rate is decomposed into an elastic part and a plastic
(inelastic) part through various assumed flow rules.
The incremental plasticity models are formulated in terms of:
A yield surface, which generalizes the concept of “yield load” into a test function which can be used
to determine if the material will respond purely elastically at a particular state of stress, temperature,
etc.;
A flow rule that defines the inelastic deformation that must occur if the material point is no longer
responding purely elastically;
and some evolution laws that define the hardening – the way in which the yield and/or flow
definitions change as inelastic deformation occurs.
The models also need an elasticity definition, to deal with the recoverable part of the strain models divide
into those that are rate-dependent and those that are rate-independent.
SOL 400 includes the following models of inelastic behavior.
Metal Plasticity (von Mises, Hill, or Barlat)
Pressure-Dependent models – models the behavior of granular (soil and rock) materials or polymers,
in which the yield behavior depends on the equivalent pressure stress.
• Linear Mohr-Coulomb
• Parabolic Mohr-Coulomb
Main Index
Chapter 6: Materials 99
Elastoplastic Material Entries
Entry Description
MATEP Specifies elasto-plastic material properties to be used for large
deformation analysis.
MATTEP Specifies temperature-dependent elasto-plastic material properties to be used
for static, quasi-static, or transient dynamic analysis.
Table 6-4 shows the various input options and criteria available to you for defining elastoplastic behavior.
Drucker-Prager
Parabolic Mohr-
Coulomb
Hardening Slope Isotropic von Mises None
Kinematic Piecewise Linear
Tresca
Combined Cowper-
Mohr-Coulomb
Piecewise Symonds
Drucker-Prager
Linear
None
Main Index
100 MSC Nastran Getting Started Guide
Elastoplastic Material Entries
Stress/Strain Curve – All stress-strain curves are input as piecewise linear. Patran transfers the stress-
strain curve input on the material property field directly to the TABLES1 entry.
The number of linear segments used to define the stress-strain curve may be different from one
material to another. The same strain breakpoints need not be used for all of the different material’s
stress-strain curves. It is recommended to define the stress-strain curves throughout the range of
strains which the analysis is likely to predict. If the analysis predicts a plastic strain greater than the
last point defined by you, SOL 400 continues the analysis after shifting the last strain breakpoint on
that curve to match the predicted value, thereby changing (reducing) the work hardening slope for
the last segment of the curve.
Hardening Slope – The hardening slope and the yield point are required with this Nonlinear Data
Input option.
Perfectly Plastic – Perfect plasticity is described by simply specifying the yield point.
Tables 6-5 through 6-7 provide descriptions for the input data for each of the four types of nonlinear input.
Table 6-5 Isotropic - Stress/Strain Curve or Perfectly Plastic: All Yield Functions
Property Name Description
Stress /Strain Curve or Defines the Cauchy stress vs. logarithmic strain (also called equivalent
Yield Stress tensile stress versus total equivalent strain) by reference to a tabular field.
The field is selected from the Field Definition list. The field is created
using the Fields application. For Perfectly Plastic models, only a Yield
Stress needs to be entered.
Can also be strain rate dependent if Strain Rate Method is Piecewise
Linear. Accepts field of yield stress vs. strain rate.
Coefficient C Visible if Strain Rate Method is Cowper-Symonds.
Inverse Exponent P Visible if Strain Rate Method is Cowper-Symonds.
Alpha When set to Linear Mohr-Coulomb, defines the slope of the yield surface
in square root J2 versus J1 space. This property is required.
Beta When set to Parabolic Mohr-Coulomb, defines the beta parameter in the
equation that defines the parabolic yield surface in square root J2 versus
J1 space. This property is required.
Note: Perfectly Plastic is identical to Stress/Strain except that no hardening rules apply.
Main Index
Chapter 6: Materials 101
Strain Rate Dependent Yield
Table 6-6 Anisotropic/Orthotropic - Stress/Strain Curve or Perfectly Plastic: All Yield Functions
Property Name Description
Stress vs. Strain or Tensile Same as description for Isotropic Elastic-Plastic. If Strain Rate Method is
Yield Stress Piecewise Linear, accepts field of yield stress vs. strain rate.
Or defines an isotropic yield stress. It is a required property when the
Plasticity Type is Perfectly Plastic.
Stress 11/22/33 Yield Defines the ratios of direct yield stresses to the isotropic yield stress in the
Ratios element’s coordinate system.
Stress 12/23/31 Yield Defines the ratios of shear yield stresses to the isotropic shear yield stress (yield
Ratios divided by square root three) in the element’s coordinate system.
Note: Perfectly Plastic is identical to Elastic-Plastic except that no hardening rules apply.
Yield stress variation with strain rate is given using one of four options:
1. The breakpoints and slopes for a piecewise linear approximation to the yield stress strain rate curve
are given. The strain rate breakpoints should be in ascending order, or
2. The Cowper and Symonds model is used. The yield behavior is assumed to be completely determined
by one stress-strain curve and a scale factor depending on the strain rate.
3. The yield stress may be given as a function of the plastic strain, strain, and/or the temperature using
the TABLD3 bulk data entry.
4. The Power Law, Rate Power Law, Johnson Cook model, or Kumar model.
Main Index
102 MSC Nastran Getting Started Guide
Strain Rate Dependent Yield
Note: If multiple material models are used, they must all be expressed as piecewise linear or as Cowper
and Symonds model.
Perfectly Plastic
A material is said to be “perfectly plastic” if, upon the stress state touching the yield surface, an infinitesimal
increase in stress causes an arbitrarily large plastic strain. The uniaxial stress-strain diagram for an elastic-
perfectly plastic material is shown in Figure 6-5. Some materials, such as mild steel, behave in a manner which
is close to perfectly plastic.
xx
.
YS
xx
where y is the yield stress, 0 is the strain corresponding to initial yield stress, p is equivalent plastic
·
strain, p is rate of equivalent plastic strain rate and A , B , m , and n are material parameters.
Main Index
Chapter 6: Materials 103
Creep (MATVP, CREEP)
Johnson-Cook Model
· m
pn
p T – T room
----
-
y = A + B 1 + C ln · 1 – --------------------------------
-
0 T melt – T room
·
where y is the yield stress, p is the equivalent plastic strain, p is the current equivalent plastic strain
·
rate, 0 is strain rate of material characterization and A , B , C , m and n are material parameters. T ,
T room , T melt are, respectively, the current, room, and melting temperatures of the material in absolute
scale.
The following conditions should be noted for the Johnson-Cook model.
T should be between T room and T melt . If T T room , T is set to T room . If T T melt ,
T = T melt – 0.01 .
· · ·
0 0 and 0 . If either condition is violated, the middle term in the above expression is set to
1.0.
Kumar Model
· 1n
p
y = B0 * sin h – 1 ----- e Q nRT
A
·
where y is the yield stress, p is equivalent plastic strain rate, Q is the activation energy, T is temperature,
and A , B 0 , and n are material parameters.
Main Index
104 MSC Nastran Getting Started Guide
Creep (MATVP, CREEP)
Creep Strain
C
Tertiary
Creep
Secondary
Creep
Primary
Creep
Time (t)
Note: Primary Creep:Fast decrease in creep strain rate
Secondary Creep:Slow decrease in creep strain rate
Tertiary Creep:Fast increase in creep strain rate
Figure 6-6 Creep Strain Versus Time (Uniaxial Test at Constant Stress and Temperature)
Main Index
Chapter 6: Materials 105
Creep (MATVP, CREEP)
ee
evp
p vp
= Plastic Element
Inactive if < y
3 i ij
iji = --- --------
2
where ij is the deviatoric stress at the end of the increment and y is the back stress. A is a function of
temperature, time, etc. Creep only occurs if sigma is greater y .
One of three tangent matrices may be formed. The first uses an elastic tangent, which requires more
iterations, but can be computationally efficient because re-assembly might not be required. The second uses
an algorithmic tangent that provides the best behavior for small strain power law creep. The third uses a
secant (approximate) tangent that gives the best behavior for general viscoplastic models.
As an example of the usage of MATVP for defining the creep behavior, see Creep of a Tube (Ch. 34) in the
Demonstration Problems Manual - Implicit Nonlinear.
Main Index
106 MSC Nastran Getting Started Guide
Creep (MATVP, CREEP)
Entry Description
MATVP Specifies viscoplastic or creep material properties to be used for quasi-static
analysis in SOL 400.
Figure 6-8
Main Index
Chapter 6: Materials 107
Composite (PCOMP or PCOMPG)
Main Index
108 MSC Nastran Getting Started Guide
Composite (PCOMP or PCOMPG)
Figure 6-9 identifies the locations of integration points through the thickness of beam and shell elements with
and without a composite formulation.
Note that when the COMPOSITE option is used, as shown in Figure 6-9, the layer points are positioned
midway through each layer. When the COMPOSITE option is not used, the layer points are equidistantly
spaced between the top and bottom surfaces. MSC Nastran Implicit Nonlinear performs a numerical
integration through the thickness. If the COMPOSITE option is used, the trapezoidal method is employed;
otherwise, Simpson’s rule is used.
*
* *
* *
* *
*
Beams or Shells with *
Beams or Shells without
Composite Option Composite Option
Figure 6-9 Integration Points through the Thickness of Beam and Shell Elements
Figure 6-10 shows the location of integration points through the thickness of continuum elements. MSC
Nastran Implicit Nonlinear forms the element stiffness matrix by performing numerical integration based on
the standard isoparametric concept.
Main Index
Chapter 6: Materials 109
Composite (PCOMP or PCOMPG)
* *
* *
* *
* *
Figure 6-10 Integration Points through the Thickness of Continuum Elements
PCOMP
1 2 3 4 5 6 7 8 9 10
PCOMP PID Z0 NSM SB FT TREF GE LAM
MID1 T1 THETA1 SOUT1 MID2 T2 THETA2 SOUT2
MID3 T3 THETA3 SOUT3 etc.
Field Contents
PID Property identification number.
Z0 Distance from the reference plane to the bottom surface.
NSM Nonstructural mass per unit area.
SB Allowable shear stress of the bonding material.
FT Failure theory.
TREF Reference temperature.
LAM “Blank”, “SYM”, “MEM”, “BEND” option.
MIDi Material ID of the various plies. The plies are identified by serially numbering them from
1 at the bottom layer.
Ti Thicknesses of the various plies.
THETAi Orientation angle of the longitudinal direction of each ply with the material axis of the
element.
SOUTi Stress or strain output request.
Main Index
110 MSC Nastran Getting Started Guide
Composite (PCOMP or PCOMPG)
An alternative to the PCOMP entry is the PCOMPG entry. The PCOMPG entry includes a global ply ID,
so it is easier to track the output for the same ply across the model. The format for the PCOMPG bulk data
entry in the QRG is as follows:
PCOMPG
1 2 3 4 5 6 7 8 9 10
PCOMPG PID Z0 NSM SB FT TREF GE LAM
GPLYID1 MID1 T1 THETA1 SOUT1
GPLYID2 MID2 T2 THETA2 SOUT2
Field Contents
PID Property identification number. (0 < Integer < 10000000)
Z0 Distance from the reference plane to the bottom surface. (Real; Default = -0.5 times the
element thickness.)
NSM Nonstructural mass per unit area. (Real)
SB Allowable shear stress of the bonding material (allowable interlaminar shear stress).
Required if FT is also specified. (Real > 0.0)
FT Failure theory. The following theories are allowed (Character or blank. If blank, then no
failure calculation will be performed)
“HILL” for the Hill theory.
“HOFF” for the Hoffman theory.
“TSAI” for the Tsai-Wu theory.
“STRN” for the Maximum Strain theory.
TREF Reference temperature. (Real; Default = 0.0)
GE Damping coefficient. (Real; Default = 0.0)
LAM Laminate Options. (Character or blank, Default = blank).
“Blank” All plies must be specified and all stiffness terms are developed.
“MEM” All plies must be specified, but only membrane terms (MID1 on the
derived PSHELL entry) are computed.
“BEND” All plies must be specified, but only bending terms (MID2 on the derived
PSHELL entry) are computed.
“SMEAR” All plies must be specified, stacking sequence is ignored MID1=MID2 on
the derived PSHELL entry and MID3, MID4, TS/T, and 12I/T**3 terms
are set to zero).
Main Index
Chapter 6: Materials 111
Composite (PCOMP or PCOMPG)
Field Contents
“SMCORE” All plies must be specified, with the last ply specifying core properties and
the previous plies specifying face sheet properties. The stiffness matrix is
computed by placing half the face sheet thicknesses above the core and the
other half below with the result that the laminate is symmetric about the
mid-plane of the core. Stacking sequence is ignored in calculating the face
sheet stiffness.
GPLYIDi User-defined Global (External) Ply ID. The global ply identification number should be
unique with respect to other plies in the entry. The plies are defined in stacking sequence
starting with the bottom layer. (Integer > 0)
MIDi Material ID of the various plies. The plies are identified by serially numbering them from
1 at the bottom layer. The MIDs must refer to MAT1, MAT2, or MAT8 bulk data entries.
(Integer > 0 or blank, except MID1 must be specified.)
Ti Thicknesses of the various plies. (Real or blank, except T1 must be specified.)
THETAi Orientation angle of the longitudinal direction of each ply with the material axis of the
element. (If the material angle on the element connection entry is 0.0, the material axis and
side 1-2 of the element coincide.) The plies are to be numbered serially starting with 1 at
the bottom layer. The bottom layer is defined as the surface with the largest -Z value in the
element coordinate system. (Real; Default = 0.0)
SOUTi Stress or strain output request. (Character: “YES” or “NO”; Default = “NO”)
Using the material properties for each of the lamina with conventional shell elements, MSC Nastran
calculates the equivalent PSHELL and MAT2 entries generated as shown in Figure 6-11. If nonlinear material
behavior including progressive failure occurs, then the advanced elements should be activated using the
PSHLN1 bulk data entry.
Main Index
112 MSC Nastran Getting Started Guide
Composite (PCOMP or PCOMPG)
CQUAD4
PCOMP
or
EQUIV PSHELL
Output if
ECHO = SORT
MID1 MID2 MID3 MID4 (the default)
Main Index
Chapter 6: Materials 113
Cohesive Zone Modeling (MCOHE)
change the ID numbers. Also, if you use the equivalent PSHELL and MAT2s, you will not be able to obtain
the laminae stress or the failure index table. For the failure index, you have a choice of four failure theories:
Hill’s theory, Hoffman’s theory, Tsai-Wu’s theory, and the maximum strain theory.
As an example of the usage of PCOMP, see Laminated Strip under Three-point Bending (Ch. 7) in the
Demonstration Problems Manual - Implicit Nonlinear.
Considering a 3-D interface element, the relative displacement components are given by one normal and two
shear components, expressed with respect to the local element system:
top bottom
vn = u1 – u1
top bottom (6-8)
vs = u2 – u2
top bottom
vt = u3 – u3
Based on the relative displacement components, the effective opening displacement is defined as:
2 2 2
v = vn + vs + vt (6-9)
Main Index
114 MSC Nastran Getting Started Guide
Progressive Composite Failure
different methods for the material degradation as described below. While MSC Nastran allows up to three
failure criteria to be used for failure index calculation, only the first criterion is used for progressive failure.
The material will not heal; the damaged elements keep the degraded properties after unloading.
r i = – 1 – e 1 – F (6-10)
This is done differently for different failure criteria as described below. Six such reduction factors are stored
and updated. They are then used for scaling the respective material modulus according to
new orig
E 11 = r 1 E 11 (6-11)
new orig
E 22 = r 2 E 22 (6-12)
new orig
E 33 = r 3 E 33 (6-13)
new orig
G 12 = r 4 G 12 (6-14)
new orig
G 23 = r 5 G 23 (6-15)
new orig
G 31 = r 6 G 31 (6-16)
The Poisson’s ratios are scaled in the same way as the corresponding shear modulus.
For the maximum stress and maximum strain criteria the reduction factors are calculated separately from each
separate failure index: r 1 is calculated from the first failure index as given by equation 6-11 above, r 2 is
calculated from the second failure index from equation 6-12 etc. Thus, there is no coupling of the different
failure modes for these criteria.
Main Index
Chapter 6: Materials 115
Progressive Composite Failure
Interface
Digimat-CAE/MSC Nastran SOL 400 contains the material library containing the Digimat capabilities and
the required interfaces in order to be linked with the MSC Nastran SOL 400 Finite Element (FE) solver.
Linking MSC Nastran SOL 400 libraries with Digimat-CAE/MSC Nastran SOL 400 gives the user access
to all linear and nonlinear small-strain material models available in Digimat for FE small-strain analyses, just
like any other MSC Nastran SOL 400 material models. It also enables to take into account fiber orientation
computed by an injection molding code. The interface thus allows to model the impact of the injection
process on the structural behavior of composite parts. Note that not all the available Digimat material models
can be used with the Digimat-CAE/MSC Nastran SOL 400 interface.
Main Index
116 MSC Nastran Getting Started Guide
Progressive Composite Failure
Main Index
Chapter 7: Contact
7 Contact
Introduction to Contact 118
Contact Types 118
Contact Definition Method 123
Contact Bodies 124
Defining Contact Interactions 131
Contact in MSC Nastran 131
Contact Settings for a SOL400 Analysis 132
Contact Considerations 134
Contact Parameters 135
Main Index
118 MSC Nastran: SOL 400 Getting Started Manual
Introduction to Contact
Introduction to Contact
The analysis of a finite element model (FEM) with contact bodies, interacting with each other is called
contact analysis. It is the analysis of contact bodies (deformable or rigid) interacting with each other. Many
engineering problems involve contact between two or more components. In these problems a force normal
to the contacting surfaces acts on the two bodies when they touch each other. If there is friction between the
surfaces, shear forces may be created that resist the tangential motion (sliding) of the bodies. The aim of
contact simulations is to identify the areas on the surfaces that are in contact and to calculate the contact
pressures generated.
Contact can be a useful tool when performing FEA simulations. It allows the interaction between multiple
bodies without adding additional elements to the model. By adding the contact condition to a model
correctly it is possible to create a model that more realistically represents reality. However, by adding contact
to a model, it might cause the analysis not to complete because of the added computational complexity and
the number of things the analyst has to take into account.
In a finite element analysis (FEA), contact conditions are a special class of discontinuous constraint, allowing
forces to be transmitted from one part of the model to another. The constraint is discontinuous because it is
applied only when the two surfaces are in contact. When the two surfaces separate, no constraint is applied.
The analysis has to be able to detect when two surfaces are in contact and apply the contact constraints
accordingly. Similarly, the analysis must be able to detect when two surfaces separate and remove the contact
constraints.
By default, MSC Nastran does not assume contact exist between all bodies in a model, simply because of the
computation cost involved in checking for contact between every element and every other element in a
model. This can however be specified, at the expense of the run times, but another reason why this is not
specified by default is because different contact conditions such as varying friction values might exist
throughout the model and therefore the analyst needs a means to specify and control it.
Contact in a model is applied mathematically and unless otherwise specified, MSC Nastran defaults will be
used that are suitable for general contact cases. Default parameters will in most cases be sufficient, however,
it might cause problems under more complex or less common scenarios, or induce unnecessarily longer
solving times if the analyst doesn’t understand these parameters and set them incorrectly.
Contact can be deformable-deformable or rigid-deformable.
Contact Types
Two types of contact interaction exist between two bodies namely glue and touch contact. Although there are
variations in how these contact interactions are applied depending on the software used, these remain the two
main contact interactions.
Touching Contact: In this case, both glue and touching contact can be separate
Main Index
CHAPTER 7 119
Contact
The difference between glue and touching is glue contact does not allow sliding, but touching does.
Touching Contact
In touching contact bodies can separate
Touch contact is a valuable tool to use in FEA models although it comes with the cost of a much
longer analysis time and it does require additional considerations to ensure that the model solves
successfully.
With this method of contact, the parts in the model can move as it is allowed to by the constraints
on the body, however with touch contact active between two bodies the additional condition is that
the one body is not allowed to penetrate the other. With the movement of the body allowed it is
therefore important to consider the position of one body relative to another.
Glued Contact
Glue contact is an ideal option to use in a model with multiple parts in a structure. It allows two bodies to
be fixed at the contact surface without having the elements of both bodies be the same size and attempting
to have the nodes on the two surfaces match.
It can be applied as a linear contact which means that two element faces remain in contact no matter what
happens with the structure around it. As long as the contact condition is detected between two bodies, this
will not cause for a longer analysis time and it should not cause further struggle with the model.
Glue contact can also be applied with certain conditions. For instance, it can be specified that if the stress on
the contact surface exceeds a specified stress, the contact should be released and for the remainder of the
analysis be treated as a broken contact.
MSC Nastran internally generates MPC equations to represent the glued contact. MPCs are generated by
MSC Nastran internally to represent the glued contact. The MPCs generated to model glued contact can be
written during the analysis in a punch (*.pch) file. You can also request for the MPC equations be output
to a punch file so that the you can visualize what is glued and what is not by plotting MPCs.
These MPCs can then be directly imported into the Patran database to provide a visual reference to the glued
contact. The MPC punch file is created with the following Case Control Command:
Main Index
120 MSC Nastran: SOL 400 Getting Started Manual
Contact Types
Main Index
CHAPTER 7 121
Contact
The contact status is checked at the beginning of the step, and those nodes or segments that are in
contact will remain in glued contact for the entire step. The constraints will change due to large
rotations.
Furthermore, if a large tensile force or stress developed over the interface in the current loadcase, no
separation would occur for these regions which are initially in contact. Performing an unglue of
breaking glue would also not be enforced during the step for these regions. This may be successfully
used to model the union of dissimilar meshes, where at a later time one wanted to separate the
bodies (e.g., opening of a door).
When using step glue conventional contact occurs for the nodes/segments of the body which are not
in contact at the beginning of the step. That means when they come into contact, they will glue, but
they may separate within the same step
Main Index
122 MSC Nastran: SOL 400 Getting Started Manual
Contact Types
In SOLs 101 and 400, if contact is initially not true set NLGLUE on BCPARA to 1. For SOL 400 with a
mixture of glued and non-glued bodies, BCPARA,0,NLGLUE,1 must also be used
Cohesive Contact
This is a special case of glued contact. In the modern industries, the product structure like airplane,
automobile, and so on, is becoming more and more complicated with many parts. The assembly process may
be done through a variety of processes such as rivets, bolts, spot weld, seam welds, or adhesives. In the
numerical simulations of the assemblies of many parts into a single structure, it is often too costly to model
each one of these discrete connectors, and the glued contact capability provides an effective and efficient way
to simplify and reduce the computation costs.
While this method is easy to use, however, it often results in too stiff of a structure because effectively the
connection is rigid. To alleviate this, a flexible glued contact capability is available. With cohesive glued
Main Index
CHAPTER 7 123
Contact
contact, it is necessary to model the detailed connectors but provides the stiffness of the connectors in the
glued contact.
There are many cases in which parts are connected using either discrete entities, such as rivets or bolts, or a
spot or seam weld. It is often too expensive to model these discrete entities so the glued approach is often
used to model these connections. While easy to use, this method is effectively a rigid connection and results
in behavior that is too stiff.
To compensate for this undesirable behavior, instead of entering a non-default constant penalty factor for the
segment-to-segment contact algorithm, you can use cohesive contact. It allows you to apply a separate finite
stiffness in the normal and tangential directions.
Main Index
124 MSC Nastran: SOL 400 Getting Started Manual
Contact Bodies
Both methods of applying contact will have the same outcome but with the contact table it is easier to define
contact for a small number of contact bodies, whereas the contact pairs method will be easier to define contact
for a large number of contact bodies.
Both methods have their own merits. Any one of the two methods can be used for any model although it is
recommended to use the contact table for a model with a small number of contact bodies and contact pairs
for a model with a large number of contact bodies. Contact can be deformable-deformable or rigid-
deformable.
Contact Bodies
These types of contact can be applied between various types of bodies in various combinations. Contact
bodies can be created from:
Main Index
CHAPTER 7 125
Contact
Main Index
126 MSC Nastran: SOL 400 Getting Started Manual
Contact Bodies
Main Index
CHAPTER 7 127
Contact
Main Index
128 MSC Nastran: SOL 400 Getting Started Manual
Contact Bodies
Main Index
CHAPTER 7 129
Contact
Main Index
130 MSC Nastran: SOL 400 Getting Started Manual
Contact Bodies
Main Index
CHAPTER 7 131
Contact
Main Index
132 MSC Nastran: SOL 400 Getting Started Manual
Contact Settings for a SOL400 Analysis
Primary and secondary option flag: It is a tool that specifies where the contact is applied to. In this
option select the A=1 for all 3D contact bodies, which will specify that contact should occur on the
outside surface of the contact body, whereas for 1D and 2D contact bodies some other options are
available. These include options such as which side of a shell element should be in contact or
whether the edges of a 1D element should be included.
Coordinate modification: It is a useful tool that will allow MSC Nastran to move nodes (within a
tolerance distance) to avoid some contact problems that might occur. Two options are available here
with the combination of both options available as the third.
In the case of stress free initial contact, nodes are moved in the beginning of the analysis to ensure that the
analysis will start with the nodes in contact. For this option, the nodes will be moved if either the nodes are
already penetrating the contact bodies or a small gap is exist between the two bodies. The nodes will then be
moved (within the tolerance distance) until the nodes are aligned.
Main Index
CHAPTER 7 133
Contact
Main Index
134 MSC Nastran: SOL 400 Getting Started Manual
Contact Considerations
Contact Considerations
When running an analysis that contains contact, whether it is glue or touch contact, there are certain
considerations to be made for the model itself in order to ensure that the contact condition is applied
correctly.
Always start by using a small version of the model in order to test whether the model is set up
correctly. This applies to models with contact especially since touch contact will significantly
increase the solving time, which will result in an unnecessary loss of time if an analysis fails or
produces incorrect results due to an error with the contact condition.
Although it is theoretically possible to have elements on the contact surface that vary greatly in size,
this is not always the case. Some errors have been found in detecting the contact correctly when the
size of the element faces on the contact surface varied too greatly.
This is because you are trying to have contact between one surface with element on the other surface.
For example: fifty contact points with another surface with only five. Although you can find the
contact, it is possible to miss a few of elements when checking the contact. The different sizes of
elements can be further evaluated to find the required element size for the contact surface. However
a good starting point is to have a maximum of five elements on one surface in contact with one
element on the other surface.
A possible way to decrease the size of the elements on the contact surface without increasing the total
amount of elements significantly is to separate the section in contact from the rest of the body and
only increase the amount of elements in that section and apply self-contact to that body. This will
allow the element size to gradually reduce from the one body to the next across the contact surface.
When working with load controller rigid contact bodies, especially with geometries such as a wheel,
it is important to remember that the body will be treated as rigid and therefore the contact surface
will not increase due to deformation. In such cases the deformable contact body might require even
smaller element faces on the contact surface to ensure that the contact condition is detected correctly
and a penetration condition for the rigid body into the deformable body is avoided.
Main Index
CHAPTER 7 135
Contact
Another case is when applying a linear motion to a rigid body for a model where the rigid body is
moved over a deformable body and a load being applied at each step after the body is moved. In
order to select an appropriate distance to move in each step (by choosing the step size) it is
important to select a distance small enough that the rigid body does not move into the deformable
body due to the deformation of the deformable body.
Although having smaller time steps will increase the number of steps to solve for and therefore
increase the total analysis time, having too large steps will either cause the analysis to fail or have an
even longer analysis time in order to solve the penetration condition. It is therefore important to
choose a number of steps for the analysis by considering the distance that the rigid body will move
in each step as well as the size of the elements on the contact face.
Applying contact to a model correctly can therefore allow you to decrease the analysis time by reducing the
total amount of elements. Alternatively it is possible to create a more realistic model by including a realistic
interaction between bodies in your model.
Contact Parameters
Contact parameters defines options for detecting and handling contact
Main Index
136 MSC Nastran: SOL 400 Getting Started Manual
Contact Parameters
Main Index
Chapter 8: Constraints
MSC Nastran Implicit Nonlinear (SOL 600) User’s Guide
8 Constraints
Introduction 138
Constraints 138
Enforced Motion Constraints (SPCD and SPCR) 139
Applying Constraints 139
Multipoint Constraints (MPC) 140
Static Loads 141
CWELD/CFAST/CSEAM Element Enhancements 142
Main Index
138 MSC Nastran: SOL 400 Getting Started Manual
Introduction
Introduction
This chapter describes the procedures for applying loads and constraints which can be used in SOL 400
models. Each type of load or constraint is described on how it is used in finite element modeling and in the
MSC Nastran input. Further information on how the various types of loads and constraints are used can be
found in standard finite element text books, and further information on the MSC Nastran input format can
be found in the MSC Nastran Quick Reference Guide.
Constraints
In static analysis, the rigid body modes must be restrained in order to eliminate singularity of the stiffness
matrix. The required constraints may be supplied with single point constraints, multipoint constraints, or
free body supports. If free body supports are used, the rigid body characteristics will be calculated and checked
for the sufficiency of the supports.
Boundary conditions are imposed in the form of constraints on selected degrees of freedom on the model.
Several degrees of freedom (at least six) are constrained to ground using either SPC bulk data entries or the
PS field of the GRID entry.
Other than single-point constraints, MSC Nastran provides a method for creating linear constraint
relationships between several degrees of freedom known as multiple-point constraints or MPCs.
Main Index
CHAPTER 8 139
Constraints
y
n
CROD 1
x
G
CROD 2
Applying Constraints
Once you have designed a model, constraints are added that forces selected portions of your model to remain
fixed or to move by a specified value. These constraints can be either:
Single Point Constraints
Multipoint Constraints
Main Index
140 MSC Nastran: SOL 400 Getting Started Manual
Multipoint Constraints (MPC)
where u is the value of the prescribed displacement on the degree of freedom i . The case of u = 0 is the
most common case, and is often used as a boundary condition, to “fix” or “ground” the movement of a point
in a certain direction.
Since the value of i is known, one could, in principle, eliminate the specified degree of freedom from the
other degrees of freedom to be solved for as unknowns. This would reduce the size of the system of equations
to be solved, but on the other hand it would take time to perform the elimination, and this approach adds
complexity to the code.
SOL 400 uses a different technique. A number which is large compared to the stiffness coefficients (say, for
discussion, 10 20 ) is added to the diagonal term K 11 of the equation for the degree of freedom to be
20
constrained. Also, if the degree of freedom is to be constrained to a nonzero value u, then u 10 is added
to the right hand side of the modified equation. This modified equation is now:
K i1 1 + + K ii 10 20 i + + K in n = F i + u 10 20 (8-2)
Assuming all K ij to be small with respect to 10 20 , the solution of the system of equations is obtained with
negligible error.
The modified system of equations remain well conditioned. The value used by SOL 400 for the large number
is 10 10 times the largest stiffness coefficient found on the diagonal of the stiffness matrix.
Main Index
CHAPTER 8 141
Constraints
Static Loads
A load is applied to the model only if it is specifically called out in the case control section. If you forget to
request any load in the case control section, the problem will be solved with zero loads applied. There are no
error or warning messages indicating that there is no load being applied. Forgetting to specify a load request
in the case control section can be a common mistake many new users make. An indication of this problem
is when all of the displacements and stresses come out to be zero.
For more information related on different types of loads refer to MSC Nastran Nonlinear (SOL400)
User’s guide.
Main Index
142 MSC Nastran: SOL 400 Getting Started Manual
CWELD/CFAST/CSEAM Element Enhancements
Object Type
Displacement/Velocty/Acceleration
Nodal
Element Uniform
Element Variable
Force • Nodal
Pressure • Element Uniform
Element Variable
Temperature • Nodal
Element Uniform
Element Variable
Inertial Load • Element Uniform
Initial Displacement • Nodal
Initial Velocity • Nodal
Distributed Load • Element Uniform
Element Variable
CID Distributed Load • Element Uniform
Element Variable
Total Load • Element Uniform
Contact • Element Uniform
Crack (VCCT) • Nodal
Initial Plastic Strain • Element Uniform
Initial Stress • Element Uniform
Initial Temperature • Nodal
The Loads and Boundary Conditions application controls which loads and boundaries and contact
information will be created in the MSC Nastran input file.
For more information related on using Patran to apply loads and boundary conditions of loads refer to MSC
Nastran Nonlinear (SOL400) User’s guide.
Main Index
CHAPTER 8 143
Constraints
Benefits
The CWELD/CFAST elements provide the same element output formats in both linear and SOL 400
nonlinear solutions.
In the CFAST/CSEAM/CWELD analysis, the auxiliary points generated are in the solution set.
The CWELD/CFAST algorithm has been improved to find the Best Possible Projection with zero
projection tolerance improvements.
The improved CWELD with options PARTPAT and ELPAT, and CFAST elements do not move GA and
GB if both are supplied by the user, thus maintaining user mesh integrity.
The 3x3 mesh limitation has been removed for the CWELD with options PARTPAT and ELPAT and the
CFAST elements.
There are no required changes in the user element input description of the CFAST/CSEAM/CWELD
elements.
The CFAST/CWELD/CBUSH provides nonlinear force output for SOL 400 ANALYSIS=NLTRAN.
MPC force output is available for the connector element constraints.
Besides global search algorithm control there is now local connector element connectivity control
via the new CONCTL bulk data entry.
A brief summary of connector projection results is output in the F06 file for each connector type.
A new SWLDPRM, CSVOUT, UNITNUM entry will produce a comma separated file useful for reports.
The CSEAM and CWELD (not by default) can now contribute mass to the structure.
Description of Features
Details of the improved CWELD/CFAST algorithm are described below.
Formulation Changes
In the new consistent formulation for the CWELD/CFAST elements for linear analysis, RBE3 elements are
written internally, and the auxiliary points are in the solution set and both are identifiable by the
SWLDPRM, PRTSW entry.
The auxiliary grids generated start with GRID ID 101000001. There are always four auxiliary grids for
patch A and four auxiliary grids for patch B.
The RBE3 elements generated start with 100001002. An RBE3 is generated for each auxiliary point
for each patch A and B tying each patch grid to that auxiliary point. There is a RBE3 generated for
GA tying GA to its patch auxiliary points and a RBE3 generated for GB tying GB to its patch auxiliary
points.
Both linear and nonlinear output is consistent.
The new consistent CWELD/CFAST is selected by default. The old CWELD/CFAST can be actuated by using the
PARAM, OLDWELD, YES entry.
It has been determined in testing that the above formulation changes produce little or no change in solution
results when comparing the old CWELD/CFAST against the new CWELD/CFAST results.
Main Index
144 MSC Nastran: SOL 400 Getting Started Manual
CWELD/CFAST/CSEAM Element Enhancements
d. Use the new CONCTL Bulk Data entry with SWLDPRM, MOVGAB, 1" to allow local control of
specific welds to correct the locations of grids GA and GB.
CONCTL, 83, , MOVGAB, 1
Where: SET3, 83, ELEM, 1345, 2678
Main Index
CHAPTER 8 145
Constraints
b. The tolerance is increased by 0.02 until a projection is found or the PROJTOL value is reached.
c. This can be turned off while computing the auxiliary grid projection onto EIDA/EIDB by setting
PROJTOL= - value where 0.0 . value . 1.0. In this case, the projection calculation starts at tolerance
= PROJTOL|. For the rest of the projection search, the algorithm reverts back to (a) and (b) above.
6. A brief Connector Summary of projection results is always output in the f06 file for each connector
type: FST-ELEM, WLD-PARTPAT, WLD-ELEMID, etc.
7. For linear connectors, MPCFORCE output is available. In nonlinear SOL 400, the RBE3 elements
generated become Lagrange elements if the default RIGID=LAGRANGE is used and are no longer in
the MPC set; hence, there will be NO MPCFORCE output for RIGID=LAGRANGE.
CWELD will not contribute to MASS by default even if its associated MATi entry has a nonzero
density. To react to a nonzero density SWLDPARM,WMASS,1 is required. If mass is computed, the
PARAM,COUPMASS effects the mass calculation.
8. In the improved CFAST and CWELD, GA and GB are not moved and internal coincident grids are
not generated at a new location; thus, two additional restrictions are required.
a. There can be no user-supplied constraints on GA and GB. A fatal message will be issued if there
are any.
b. The CWELD length must be > 10-6. The point to patch option defined by ELEMID or GRIDID will,
however, still create a new GS internally to obtain a minimum required length; i.e., LDMIN .
length/D . LDMAX. For the point-to-patch connection, GS is used as GB. The algorithm will
use the new GS as GB but keep the user-supplied GS unchanged. Since the point-to-patch is
often used to tack two shell corners, the default LDMIN may cause the connector to be unstable.
To avoid this, it is recommended that the user set LDMIN=1.E-6 on the appropriate PWELD entry.
9. CFAST and the CWELD with options PARTPAT and ELPAT with the improved formulation has
removed the restriction that a connector patch cannot span more than three elements. It will now
span over a patch of as many elements that the value of diameter D of the patch encloses and for
which projections can be found.
The following Figure, all element grids contained in the green circle region of say patch A will be
used in a RBE3 connection in addition to the RBE3 connections generated for the four auxiliary
points. The example for this figure can be found at /tpl/connectr14/cei_103.dat.
The CWELD ELEMID option still only connects by design two elements. The diameter is only used to
compute the beam stiffness.
The green circle passes through the four auxiliary points of the patch (the nine digit grid IDs.). The
user-specified diameter D on the PWELD and PFAST entries determine the locations of the four
auxiliary points. The green circle diameter is approximately 1.253D. The element grids shown
outside the green circle belong only to the respective auxiliary points.
Main Index
146 MSC Nastran: SOL 400 Getting Started Manual
CWELD/CFAST/CSEAM Element Enhancements
For higher-order shell elements CQUAD8 or CTRIA6 with no missing midside nodes, the RBE3
relationships use only the midside nodes. If one or more midside node is missing, then the corner
nodes are used.
The diameter D on the PFAST/PWELD entry is used to determine the projection location of these
auxiliary points as well as the stiffness properties of the patch to patch connection.
A single RBE3 then connects the four auxiliary points and the shell grids within the green circle to
the connector grid GA=4000065.
The SWLDPRM, PRTSW entry will list the additional grids located within the green circle under
FMESH SHELL A or B GRIDS where FMESH is the f06 file listing title of the additional shell grids
connected in the RBE3 relationships for Finer MESH.
Table 2-1 shows the grids associated with auxiliary grid 101000023 of Figure 2-1 for its RBE3
generation. The WTi's are weight factors based on patch shape functions. Grids G1, G2, G3 are
selected for RBE3 EID 100001026 because they are the shell grids of the triangular element that
contains the projected auxiliary point.
Table 2-2 shows the grids associated with grid GA=4000065 of Figure 2-1 through
RBE3=100001022. G5 through Gn are the grids contained within the green circle. Grids
G2=4007884, 4007869, 4007815, and 4007830 are NOT included in any of the G5-Gn entries
because they are included in their associated auxiliary point RBE3 elements 100001026 because they
are the shell grids of the triangular element that contains the projected auxiliary point.
Main Index
CHAPTER 8 147
Constraints
10. The CWELD/CFAST/CBUSH has the additional enhancement that for nonlinear transient SOL 400
with ANALYSIS=NLTRAN they will request Element FORCE output for the CWELD/CFAST/CBUSH
elements. FORCE=ALL will request force output for all CWELD/CFAST/CBUSH elements. Element
FORCE output in SOL 400 nonlinear transient analysis is unique to CWELD/CFAST/CBUSH elements.
Other elements will not list force output in SOL 400 nonlinear transient.
If in a SOL 400 nonlinear transient, you have FORCE=ALL and carry this over to say an
ANALYSIS=NLSTAT, you will get Element FORCE of ALL elements capable of force output, such
as CWELD, CFAST, CBEAM, CQUAD4, etc.
11. For user convenience, an additional SWLDPRM command useful for reports creates a comma separated
file of the SWLDPRM, PRTSW, using the command SWLDPRM, CSVOUT, UNITNUM where UNITNUM is
assigned via the:
ASSIGN USERFILE=myfile.csv, UNIT= UNITNUM , FORM= FORMATTED, DELETE,
STATUS=NEW.
Additional Information
1. In SOL 400, for ANALYSIS=NLSTAT or ANALYSIS=NLTRAN, the generated RBE3 constraints
become Lagrange elements and may undergo large rotation. For ANALYSIS=NLTRAN with initial
conditions (IC=n) in case control that cause large initial stresses in the structure at time t=0, the case
control entry RIGID needs to have the value RIGID=LINEAR to insure convergence.
2. For user desiring to postprocess the CFAST/CSEAM/CWELD connectors with their own methods, the
following is useful:
a. The GEOM2 table contains, after module MODGM2, a record ELCORR that correlates the
CFAST/CSEAM/CWELD and its associated RBE3 elements. Also, this module will, for linear
analysis and for nonlinear SOL 400 analysis run with RIGID=LINEAR, place the internal
generated RBE3 into the GEOM4 table.
b. In SOL 400 with RIGID=LAGRANGE (Default), internally generated RBE3 elements go into the
GEOM2 (as do all other user specified rigid elements) not the GEOM4 table.
c. The CWELD/CFAST/CBUSH force output for ANALYSIS=NLTRAN in SOL 400 is OP2 file
output on OEFNLXX data block. If SCR=POST is run, then this force data is also written to the
data base file OEFNL3 and op2 file OEFNL.op2 is also written.
3. The DISPLACEMENT (CONNECTOR=) Case Control Command works in the same fashion for both
the old connector formulation and the new connector formulation.
For the CFAST with option ELEM and the CWELD with option ELPAT, the shell elements connected
on each patch must have same property identification number of PSHELL entry.
If parameter OSWPPT is used to specify the offset for internally generated grid IDs, its value should
be greater than the maximum identification number of GRID entries to avoid conflict IDs.
Main Index
148 MSC Nastran: SOL 400 Getting Started Manual
CWELD/CFAST/CSEAM Element Enhancements
Connector Stiffness
Connector contribution to a structural model's overall stiffness is sensitive to the models mesh size and the
orientation of the connector relative to the mesh. Thus, the discretization process itself may cause, for
example, a model using a fine mesh to be stiffer in torsion than a corresponding model using a coarse mesh.
For production models that correlate well with test, refining the mesh may cause an inherent overall loss of
stiffness due to mesh refinement and hence loss of correlation.
To allow you some control over stiffness, improved connectors (CWELD with ELPAT or PARTPAT or CFAST)
are provided with two options to provide additional connector stiffness. The two options may be used
individually or in combination.
Main Index
CHAPTER 8 149
Constraints
A disadvantage of this method is that as DRATIO is increased using the global command SWLDPRM, DRATIO,
value. So, some connector elements may begin to fail because they may no longer be able to find a patch
projection. To overcome this, the SWLDPRM, NREDIA, can be increased to a value as high as 8 to allow failing
welds to halve their patch diameters up to eight times.
If the SWLDPRM, NREDIA is not an approach the user wishes to pursue, then for these failing elements, the
bulk data entry CONCTL, SETID, ,DRATIO, value can be used to define a set for failing connectors and
set a value of DRATIO for these connectors that allows them to find a projection.
The second stiffening algorithm attempts, based on the diameter of the connector, to determine a measure
of the mesh discretization. This feature is activated by SWLDPRM, SKIN, 1 or CONCTL, SETID, ,SKIN,
1. The default is a 0 which implies no stiffening. There is an associated stiffening factor SWLDPRM,
SCLSKIN with value = 0.10 as default.
Depending on the complexity of the model and the overall mesh size and the number of connectors within
the model and the diameter of the connectors relative to the mesh, the default value tends to stiffen a
structural model from about 0.4% to about 4%. A value of SCLSKIN=10.0 stiffens coarser mesh models by
about 10% to 11% and finer mesh models by about 2% to 6%.
The contribution of the stiffening algorithm to the overall stiffness of the FEM model eventually reaches a
limit. For example, a very large value SCLSKIN=100 increases the stiffness of the models overall by only
about 0.1% to 2% over the stiffness obtained for SCLSKIN=10.
For a correlated structural model evaluated at a specific mesh size, with an aim to refine the mesh for some
portion of this model containing connectors, while leaving other portions containing connectors with an
unmodified mesh, it is recommended that you enter the SKIN, 1 and SCLSKIN, real value on the CONCTL
bulk data entry referring to the connectors within the area of the refined mesh. Different refined mesh areas
within the structural model can have different values of SCLSKIN associated to the specific connectors in
each refined region.
Main Index
150 MSC Nastran: SOL 400 Getting Started Manual
CWELD/CFAST/CSEAM Element Enhancements
For postprocessing of the SKIN option, for the affected shell elements, an updated EPT table is available after
module MODGM2. It contains the PSKNSHL record that correlates the property data of the shells involved and
a list of shell elements for each patch modified.
Main Index
CHAPTER 8 151
Constraints
Main Index
Chapter 9: Boundary Conditions
MSC Nastran Implicit Nonlinear (SOL 600) User’s GuideMS
9 Boundary Conditions
Introduction 153
Zero and Enforced Displacements 153
Fixed Direction Grid Point Forces 153
p-Element Loads and Constraints 154
Thermal Loads (TEMP and TEMPD) 155
Inertial and Dynamic Loads 156
Main Index
CHAPTER 9 153
Boundary Conditions
Introduction
SOL 400 supports the following loads and boundary conditions:
• Constrained nodal displacements (zero displacements at specified degrees of freedom). Enforced
nodal displacements (nonzero displacements at specified degrees of freedom in the nodal
coordinate system). Both constrained and enforced displacements can be specified as relative or
absolute.
• Forces applied to nodes in any coordinate system.
• Follower force distributed loads.
• Temperature applied to nodes. Temperature can be applied as a load in a structural analysis. The
reference temperature is user definable.
• Inertial body forces, acceleration and velocity can be applied in the global coordinate system.
Main Index
154 MSC Nastran: SOL 400 Getting Started Manual
p-Element Loads and Constraints
or nonuniform. Special loading types also exist in various analyses. For example, centrifugal loading exists in
stress analysis, and convection and radiation exist in heat transfer analysis. You can add point loads directly
to the nodal force vector, but equivalent nodal forces first must be calculated by MSC Nastran from
distributed loads and then added to the nodal force vector.
For detailed information on Grid Point Forces and Grid Point Follower Forces refer to MSC Nastran
Nonlinear (SOL400) User’s guide.
Body Loads
Body loads consists of thermal loads and gravity loads.
For thermal loading, the TEMP Bulk Data entry is used for a temperature distribution that is trilinear over
the element, whereas the TEMPF entry is used over a set of elements. The initial or reference temperature of
a body must be supplied using these entries. Gravity loads are defined by the GRAV entry.
Boundary Conditions
The SPC and SPCD entries are used for the point constraints that are allowed only on corner GRID points.
The GMSPC entry is used to define zero constraints for FEEDGEs, FEFACEs, GMCURVs, and GMSURFs.
For nonzero constraints, the GMBC entry is used. In general, for multiple input data for FEEDGE,
FEFACE, GMCURV, and GMSURF entities, the hierarchy set to resolve these conflicts is as follows:
GRIDs
FEEDGEs
GMCURVs
FEFACEs
Main Index
CHAPTER 9 155
Boundary Conditions
GMSURFs
Table 9-1 Bulk Data Entries Used for Temperature Definition on Elements
Elements Temperature Data Bulk Data Entry
CROD, CONROD, Average temperature at ends A and B. TEMPRB
CTUBE
CBAR, CBEAM, CBEND Average temperature and cross-sectional temperature TEMPRB
gradients at ends A and B.
CBEAM3 Temperature field and gradients along the beam TEMPB3
CQUAD4, CTRIA3, Average temperature and gradient in the thickness TEMPP1
CQUAD8, CTRIA6 direction.
Average temperatures specified directly for an element take precedence over the temperatures interpolated
from the element’s connected grid points. Solid elements obtain their temperatures only by interpolation
from connected grid points. Note that interpolated grid point temperatures provide temperature gradients
over the neutral surface of shell elements, whereas the TEMPPi entries do not.
The temperature data and the thermal expansion coefficients are used internally to calculate equivalent forces
and moments acting at the grid points.
The TEMPERATURE (Case) (INIT) and TEMPERATURE (Case) (LOAD) case control commands specify
the initial temperature and applied temperature, respectively. The TEMP(INIT) command must appear
either above the first subcase or inside the first subcase.
Main Index
156 MSC Nastran: SOL 400 Getting Started Manual
Inertial and Dynamic Loads
Centrifugal (RFORCE)
The RFORCE entry is used when you need to apply a force to your structure due to rotational velocity and/or
acceleration. On the RFORCE entry, you input the components of a spin vector that are used internally to
compute centrifugal forces. Each component of the spin vector is multiplied by the same scale factor.
You must select one of two methods for the internal calculation of the loading vector.
Method=1 yields correct results only when there is no coupling in the mass matrix. This occurs when the
lumped mass option is used with or without the ZOFFS option (see the CQUAD4 entry for a description of
ZOFFS). Method=2 yields correct results for lumped or consistent mass matrix only if the ZOFFS option is
not used. The acceleration terms due to the mass offset (X1, X2, X3) on the CONM2 entry are not computed
with method=2. All the possible combinations of mass matrices and offset and the correct method to be used
are shown in Table 9-2.
In addition, for problems with elements that have edge grid points (CQUAD8, CTRIA6, CTRlAX6, CHEXA,
CPENTA, and CTETRA), correct centrifugal loads are produced only if the parameter
PARAM,COUPMASS,x (where x is greater than 1), is included in the input file and Method 2 is used.
Note for PARAM,COUPMASS=-1 (the default) the generation of lumped mass matrices that contains only
translational components for the elements listed above. Notable exceptions to this are the CBAR and CBEAM
elements, both of which will yield rotational and coupling terms in order to preserve the mass center when
element offsets are defined. This offset mass is ‘lumped’ in the sense that it has low matrix rank, and is
‘coupled’ in the sense that there are nonzero off diagonal terms in the mass matrix. The CBEAM element will
Main Index
CHAPTER 9 157
Boundary Conditions
also yield a mass moment of inertia about the local X axis of the element, and if NASTRAN BARMASS > 0,
then this is also true of the CBAR element.
In order to yield a lumped mass matrix containing translational components only for the CBAR and CBEAM
elements, set SYSTEM(414) = 1, along with the default value for PARAM,COUPMASS (-1). The default
value (0) for SYSTEM(414) produces the coupled mass matrices for CBAR and CBEAM.
Initial Stress and Initial Plastic Strain Mapping from Previous Results
MSC Nastran allows you to enter a set of initial stresses through the ISTRESS bulk data entry that simulates
the stress state in the structure at the beginning of an analysis. A typical example is prestress in a tensioned
fabric roof. The set of initial stresses must be self-equilibrating and should not exceed the yield stress of the
material.
MSC Nastran allows you to define the equivalent plastic strain using the IPSTRAIN bulk data entry
throughout the model. This is useful in metal forming analysis in which the previous amount of equivalent
plastic strain is often required. This history dependent variable represents the amount of plastic deformation
that the model was subjected to, and is used in the work (strain) hardening model. This is only used to
determine the value of the strain hardening once plasticity occurs.
Preprocessors such as Patran and SimXpert are very useful in mapping stress states between analyses.
Main Index
Chapter 10: Iteration Control in Nonlinear Analysis
Introduction 159
Nonlinear Characteristics and General Recommendations 160
Starting the Analysis 161
Load Incrementation and Iteration 166
Load Increment Size 167
Convergence Controls 169
Solution Parameters 170
Defining Subcases in Patran 171
Main Index
CHAPTER 10 159
Iteration Control in Nonlinear Analysis
Introduction
MSC Nastran SOL 400 provides a comprehensive FE solution for multi-physics problems such as structure
analysis, thermal analysis, as well coupled analysis. Even though it is originally targeted to deal with various
nonlinearities as geometry nonlinearity, material nonlinearity, as well boundary nonlinearity (contact), it is
also embedded with the powerful capabilities of MSC Nastran on linear analyses.
The FEM is a powerful tool for analyzing complex problems in structural and continuum mechanics. The
analysis of a structure using the FEM has four basic steps:
1. Modeling, in which the structure is subdivided into an assemblage of discrete volumes called finite
elements, and properties are assigned to each element.
2. Evaluation of element characteristics, such as stiffness and mass matrices, followed by assembling the
element characteristic matrices to obtain the assembled or “global” matrices characteristic of the
entire structure. A similar process is followed to obtain the total loads, in vector form, applied to the
structure.
3. Solution of the system equations for displacements, natural frequencies and mode shapes, or buckling
load factors.
4. Calculating other quantities of interest, such as strains, stresses and strain energy. MSC Nastran SOL
400 provides a comprehensive FE solution for multi-physics problems such as structure analysis,
thermal analysis, as well coupled analysis. Even though it is originally targeted to deal with various
nonlinearities as geometry nonlinearity, material nonlinearity, as well boundary nonlinearity
(contact), it is also embedded with the powerful capabilities of MSC Nastran on linear analyses. The
FEM is a powerful tool for analyzing complex problems in structural and continuum mechanics. The
analysis of a structure using the FEM has four basic steps:
SOL 400 uses the finite element displacement method, in which a system of equations is solved to obtain
the displacements at all node points of the structure. Comprehensive presentations of the FEM together with
numerous applications are available in textbooks and the research literature. Solving the equations of finite
element analysis involves the manipulation of large matrices of numbers, which is best done using computers.
While modern computers are extremely fast and have vast amounts of memory, the “bottom line” in finite
element analysis is that it is very easy to discretize a structure to the point that it is considered to be a “large”
problem.
Typically, the geometry of the problem will dictate how fine of a mesh is required to get an acceptable
solution. While modern graphical user interfaces (GUI's) such as Patran can provide some guidance and tools
for evaluating the mesh density, in the end it is up to the user to ensure that the problem is adequately
discretized. The size of FEA problems is typically measured by the number of degrees of freedom (or degrees
of freedom) in the finite element mesh.
For nonlinear problems, the size is even more important because the methods used to solve these nonlinear
equations are usually iterative in nature, meaning that the system of equations must be solved many, many
times to follow the behavior of the structure as it changes. The changes characterized may be the shape (called
large deformation), or the status of the material (metal materials typically yield when a certain stress level is
exceeded), or loading changes (the load orientation may follow the deformation of the structure, or different
parts of the structure may come into contact changing the load path). The purpose of this chapter is to
Main Index
160 MSC Nastran: SOL 400 Getting Started Manual
Nonlinear Characteristics and General Recommendations
describe the numerical methods and procedures required to solve the linear and nonlinear equations used to
perform finite element simulations.
We begin this chapter with some general observations and recommendations on how to approach the solution
to nonlinear FEA problems, and then move through the details of how to get an accurate, efficient solution
using the algorithms and methods available in MSC Nastran SOL 400. Along with the efficiency, we discuss
the concept of “robustness”, which means choosing alogorithms that will tolerate changes in the problem
(sometimes abrupt, such as contact characterizing impact) with out causing the solution algorithms to fail. A
robust numerical method automates changes in the solution parameters, such as load increment size, as
required by the problem yet still provides a reasonably accurate and efficient solution.
The characteristics of the solution methods and strategies in this chapter include: solution iteration methods,
transient analysis effects such as time step selection, and time integration method, arc-length methods for
post-buckling, and convergence measurements and controls. This chapter may be read as a tutorial, but is
probably best used as a reference when considering specific problems encountered while trying to solve a
nonlinear finite element problem.
Main Index
CHAPTER 10 161
Iteration Control in Nonlinear Analysis
geometric nonlinear effects should be significant if the deformed shape of the structure appears distinctive
from the original geometry without amplifying the displacements. There is no distinct limit for large
displacements because geometric nonlinear effects are related to the dimensions of the structure and the
boundary conditions. Strain is nondimensional, and a strain greater than 2% indicates a geometrically
nonlinear problem.
Main Index
162 MSC Nastran: SOL 400 Getting Started Manual
Starting the Analysis
You are guaranteed a solution if the boundary conditions and material properties are set up correctly. The
stiffness matrix is assembled and solved only once in the entire analysis.
Main Index
CHAPTER 10 163
Iteration Control in Nonlinear Analysis
Newton-Raphson Iteration scheme begins: Apply a portion of the total load to start: (1% in this case):
• Assembly of stiffness Matrix.
• Solution of the Stiffness Matrix.
• Check for convergence (this is an important step which is only seen in nonlinear analysis).
If converged, the solution/structure is in equilibrium. Go to step 3.
If not converged, update information and reassemble, resolve stiffness matrix.
Keep iterating till convergence is achieved.
3. (After convergence) Get displacements, strains, stresses.
4. Apply the next increment of load and go to Step 2. Keep doing this until all the load is applied.
5. View Results.
Iterations
In the incremental solution process, the unbalanced forces that occur during a load increment are
reintroduced internally into the solution until the solution has converged. The process of redistributing the
unbalanced force within a load increment is known as an iteration. The iteration is the lowest level of the
solution process. Iterations continue within a load increment until the solution converges or any of the
specified convergence parameters are exceeded. A complete description of the numerical procedure used to
solve the nonlinear problem is given in this chapter.
Main Index
164 MSC Nastran: SOL 400 Getting Started Manual
Starting the Analysis
Input Phase:
Read Input Data
Space Allocation
Data Check
Incremental
Loads
Equivalent Nodal
Load Vector
Matrix Assembly
Iteration Loop
Matrix Solution
Time Step Loop
Stress Recovery
No
Convergence
Yes
Output Phase
Yes Next
Increment
No
Stop
The important point to note is that the total load is applied gradually in steps (or increments) and for each
load step, the solution is arrived at after one or more iterations. If the behavior of the model is generally linear,
few iterations are required to solve that load step. If the model behavior is complex/nonlinear, many iterations
might be required. Each iteration involves an assembly and solution of the stiffness matrix. Hence, nonlinear
problems inherently take longer than linear models (of the same size) to solve. At the end of each iteration, a
check is made to see if the solution has converged. If the convergence check fails, the iteration is repeated with
the new information. This process repeats until convergence is achieved. Following that, the next increment
of load is applied. The load increments are applied until the full load of the model is solved.
Main Index
CHAPTER 10 165
Iteration Control in Nonlinear Analysis
To achieve accurate results, three key points have to be paid attention to:
Iteration method;
Load increment control; and
Convergence criterion.
The three points will be discussed in later sections.
Main Index
166 MSC Nastran: SOL 400 Getting Started Manual
Load Incrementation and Iteration
Normal rotation for the "drilling degree of freedom" of shell elements restrained by the default value
of 100 on the K6ROT parameter when the geometric nonlinearity is involved. In rare cases it may
be necessary to adjust this value. This can help with convergence, but may also affect the results.
Understand the basic theory of plasticity, creep, or rubber elasticity before using these capabilities.
The time step size for a transient response analysis should be carefully considered based on the
highest natural frequency of interest because it has significant effects on the efficiency as well as the
accuracy of the solution. The automated procedures used by NLSTEP is adequate for this purpose.
Main Index
CHAPTER 10 167
Iteration Control in Nonlinear Analysis
Main Index
168 MSC Nastran: SOL 400 Getting Started Manual
Load Increment Size
make other options for mechanical, thermal and coupled analysis. MSC has made every attempt to make
MSC Nastran SOL 400 as robust, efficient, and user friendly as possible; especially when used with either of
the MSC graphical user interfaces (GUI’s) Patran or SimXpert. For many problems the defaults entered are
appropriate to minimize the job setup and obtain accurate results.
The NLSTEP option has a keyword CTRLDEF that automatically sets up the entries for the time stepping
adjustment and convergence tolerance based upon how nonlinear you believe the problem is. This makes it
possible to use smart default based on users’ judgment of the nonlinearity of the model to be analyzed.
Under the keyword CTRLDEF, these three options (QLINEAR, MILDLY, and SEVERELY) adjust the parameter
to provide you the desired results. As the names imply, CTRLDEF should be set to QLINEAR for linear
solutions, MILDLY for mildly nonlinear, and SEVERELY for severely nonlinear behavior.
Main Index
CHAPTER 10 169
Iteration Control in Nonlinear Analysis
For details, refer to chapter 4 of Sol 400 UG (refer to the exact section)
Convergence Controls
Three methods are available for determining if convergence is obtained on any given iteration: residual force,
displacement, and strain energy. You can select one of these three criteria for convergence or you may specify
a combination of residual and displacement. The AND combination signals that both residual and
displacement must be met, while the OR combination specifies that either one can satisfy convergence
criteria. If you are using residual there may be cases in which the force residuals are null in which case is it
necessary to switch over to displacement. An Autoswitching option (on by default) allows for this switching.
The default measure for convergence in SOL 400 is residual which is based on the magnitude of the
maximum residual load compared to the maximum reaction force. This method is appropriate since the
residuals measure the out-of-equilibrium force, which should be minimized. This technique is also
appropriate for Newton methods, where zero-load iterations reduce the residual load. The method has the
additional benefit that convergence can be satisfied without iteration. You have complete control over how
convergence is defined through the Iterations Parameters form in Patran or through the options on the
NLSTEP entry.
The basic procedures are outlined below.
1. Residual checking: Residual checking has one drawback. In some special problems, such as free
thermal expansion, there are no reaction forces. If the value of the residuals and reactions is less than
1.e-6, this test is ignored. If the AUTOSW flag on the NLSTRAT entry is ON, the program automatically
uses displacement checking in this cases.
2. Displacement checking: With this method, convergence is satisfied if the maximum displacement of
the last iteration is small compared to the actual displacement change of the increment. If the value
of the incremental and iterative displacement is less than 1.e-8, this test is ignored. A disadvantage of
this approach is that it results in at least one iteration, regardless of the accuracy of the solution.
3. Strain energy checking: This is similar to displacement testing where a comparison is made between
the strain energy of the latest iteration and the strain energy of the increment. With this method, the
entire model is checked. A disadvantage of this approach is that it results in at least one iteration,
regardless of the accuracy of the solution. The advantage of this method is that it evaluates the global
accuracy as opposed to the local accuracy associated with a single node.
4. Residual or displacement checking: This procedure does convergence checking on both residuals
(Residual 1) and displacements (Procedure 2). Convergence is obtained if one converges.
5. Residual and displacement checking: This procedure does a convergence check on both residuals and
displacements (Procedure 4). Convergence is achieved if both criteria converge simultaneously.
Different problems require different schemes to detect the convergence efficiently and accurately. To do this,
the combinations of residual checking and displacement checking are also available (as mentioned in the last
two steps.
Main Index
170 MSC Nastran: SOL 400 Getting Started Manual
Solution Parameters
Solution Parameters
Solution parameters control a range of functions in the SOL 400 analysis. Functions such as selecting the
solver type, establishing a restart, specifying domain decomposition are all part of the solution parameters.
Solver Options Specifies the solver to be used in numerically inverting the system of
linear equilibrium equations.
Contact Parameters Defines options for detecting and handling contact.
Direct Text Input This subform is used to directly enter entries in the File Management,
Executive Control, Case Control, and bulk data sections of the MSC
Nastran input file.
Restart Parameters Includes a Restart option in the MSC Nastran input file.
Advanced Job Sets alternate versions of the solver and alternate formats for the results
Control file.
Domain Designates that domain decomposition be done manually,
Decomposition semi-automatically, or automatically.
Main Index
CHAPTER 10 171
Iteration Control in Nonlinear Analysis
Main Index
172 MSC Nastran: SOL 400 Getting Started Manual
Defining Subcases in Patran
Subcase Parameters
The subcase parameters represent the settings in MSC Nastran Case Control and Bulk Data Section that take
effect within a subcase and do not affect the analysis in other subcases. Subcase parameters are dependent on
the type of analysis being performed. The set of subcase parameters applicable for each analysis type are
described in the following sections. For more information, see Chapter 3: Solution Methods and Strategies
in Nonlinear Analysis in the MSC Nastran Implicit Nonlinear User’s Guide.
Main Index
CHAPTER 10 173
Iteration Control in Nonlinear Analysis
Entry Description
NLSTEP Defines parameters for automatic load/time stepping used in SOL 400.
NLPARM Nonlinear Static Analysis Parameter Selection.
NLPCI Defines a set of parameters for the arc-length incremental solution strategies in
nonlinear static analysis.
Note: The arc length method cannot be used with contact.
Main Index
174 MSC Nastran: SOL 400 Getting Started Manual
Defining Subcases in Patran
Load Increment Params Defines whether the load increments will be fixed or adapted in each
iteration, and the method by which adaptive load increments will be
determined.
Iteration Parameters Sets forth the iterative procedures that are employed to solve the
equilibrium problem at each load increment.
Contact Table Activates, deactivates, and controls the behavior of contact bodies in
the analysis.
Entry Description
NLSTEP Nonlinear static analysis parameter selection for doing creep analysis.
MATVPMA Defines creep characteristics based on experimental data or known empirical creep law.
This material definition should be used with advanced elements.
CREEP Defines creep characteristics based on experimental data or known empirical creep law.
This material definition should be used with conventional elements.
Main Index
CHAPTER 10 175
Iteration Control in Nonlinear Analysis
Iteration Parameters Sets forth the iterative procedures that are employed to solve the
equilibrium problem at each load increment.
Contact Table Activates, deactivates, and controls the behavior of contact bodies
in the analysis.
Entry Description
EIGRL Defines data needed to perform real eigenvalue (vibration or buckling) analysis with
the Lanczos method (recommended).
EIGR Defines data needed to perform real eigenvalue analysis.
Extraction Method Defines the method to use to extract the real eigenvalues.
Lancozs Parameters
Indicates an estimate of the number of eigenvalues to be
Number of Modes
located.
Main Index
176 MSC Nastran: SOL 400 Getting Started Manual
Defining Subcases in Patran
Entry Description
NLSTEP Parameters for automatic load/time stepping for both static and transient nonlinear analysis
(recommended).
TSTEPNL Traditional way to specify nonlinear dynamic analysis parameter selection, has been replaced
by NLSTEP
Main Index
CHAPTER 10 177
Iteration Control in Nonlinear Analysis
Entry Description
BCMOVE Specifies movement of rigid surfaces.
Main Index
178 MSC Nastran: SOL 400 Getting Started Manual
Defining Subcases in Patran
If ON, specifies that when the first rigid body comes into contact, the rest
Synchronized
stop moving.
Contact Table... Activates, deactivates, and controls the behavior of contact bodies in the
analysis.
Main Index
CHAPTER 10 179
Iteration Control in Nonlinear Analysis
Initial Condition Load Factor... Multiplier applied to loads and constraints used as initial
conditions of complex eigenvalue analysis.
Enable Rotor Dynamics... Toggle this ON to enable Specify Spinning Properties
form which allows user to input rotordynamics
properties.
Contact Table... Toggle Use Contact Table ON to use of contact table
rather than contact pairs in the analysis. Once toggled
ON, the Contact Table button brings up contact table
input form.
Break Squeal... Allows user to input parameters used for break squeal
analysis.
Main Index
180 MSC Nastran: SOL 400 Getting Started Manual
Defining Subcases in Patran
Structural-thermal analysis.
These analysis types employ nonlinear solution algorithms so that nonlinear material properties or boundary
conditions can be included in the model.
MSC Nastran Implicit Nonlinear (SOL 600) User’s Guide
Main Index
Chapter 11: Trouble Shooting
11 Trouble Shooting
Overview 182
Review Fatal Error Message 182
Review the .sts File 183
Review the .f06 File 184
SOL 400 Analysis Messages 188
Reviewing Convergence 190
Review Plot Results from Converged Increments 192
Reviewing the MSC Analysis Manager 193
Tips for Starting with Nonlinear Analysis 194
Review Fatal Error Message 182
Main Index
182 MSC Nastran: SOL 400 Getting Started Manual
Overview
Overview
Sometimes you may receive non-convergence error messages when a job does not converge. In such a
scenario, review the following diagnostics tools to help you:
The fatal error message
Plot results from converged increments
The .sts file
The .f06 file
The .log file
For details, go the end of the .f06 file and search up for %. The nonlinear iteration printout will provide
helpful diagnostics on your job
Main Index
CHAPTER 11 183
Trouble Shooting
The % symbol does not mean percent. It is a convenient way for you to search for the iteration printout. A
few key things to look for are highlighted in the figure.
The first column shows the procedural step, while the second column shows the increment within the step.
Not every increment size is equal, as can be seen in the TIME STEP OF INC column. The third column
(cycl #) indicates the number of attempts made during this increment.
The status file is updated in real time so a running job can be monitored. Review the .sts file. Check for
the following:
Did the job finish normally?
How far did the job get to?
What is the maximum displacement?
• Is the job running away?
• Does the maximum displacement make sense?I
Main Index
184 MSC Nastran: SOL 400 Getting Started Manual
Review the .f06 File
Also review the exit definition. It will inform you if the job terminated normally or abnormally. You can then
check the .f06 file for more information.
Patran provides a convenient way to access the status file using the Monitor option of the Analysis form.
Main Index
CHAPTER 11 185
Trouble Shooting
In the .f06 file, there is another shorter summary of NLSTEP parameters right above this printout. That
one is a subset of this one and is not as useful.
Main Index
186 MSC Nastran: SOL 400 Getting Started Manual
Review the .f06 File
LOAD STEP: Step number minus 1 plus fraction of step, i.e. 0.08 = 8% of first step
NO. INC: Increment number
ITR: Iteration number within the load increment
Error Factors: This must be smaller than the tolerances EPSU, EPSP and EPSW before the
solution will move ahead.
• DISP: Displacement errors
• LOAD: load errors
• WORK: work erroes
CONV RATE: Should be between 0 and 1, bigger than 1 means, the solution will never converge
ITR DIV: Divergence counter, MAXDIV triggers the divergence process
MAT DIV: Divergence counter for element and material routines
NO. BIS: Number of bisections
AVG R_FORCE: Average residual force (forces and moments). Should be small
TOTL WORK: Approximate total work
DISP - AVG MAX AT GRID C: Average and maximum displacements at grid in direction c
LINE_S:FACT NO: Line search factor a and number of line searches (not used for PFNT)
NO. QNV: Number of Quasi Newton Vectors (not used for PFNT)
TOT KUD/ITR: Total number of stiffness updates / iterations
Main Index
CHAPTER 11 187
Trouble Shooting
Main Index
188 MSC Nastran: SOL 400 Getting Started Manual
SOL 400 Analysis Messages
This is useful for tracking down where the worst residuals are located.
For each converged increment an OP2 results file is written at the end of the increment. Allow to check results
before the end of the job. If final OP2 is not written, intermediate results can be checked to help debuging.
Main Index
CHAPTER 11 189
Trouble Shooting
*** JOB DOES NOT CONVERGE AT THE CURRENT TIME STEP OR INCREMENT.
*** SOLUTION DIVERGES FOR SUBCASE m STEP n.
MSC Nastran has only two exit number, i.e.,
• 0: normal exit
• -1: fatal error
MSC Nastran provides many user fatal error messages (UFM). A normal/successful run exit message is as
follows:
For a job that terminates abnormally you are referred to the FATAL ERROR MESSAGE in the jobname.f06
file.
One way to debug convergence issues is to use the NLOPRM bulk data entry, such as NLDBG = ADVDBG,
DBGPOST options and so on, to load increment diagnostics in the.f06 file, which may help debug the
model. Each iteration of each load increment generates a report in the.f06 file giving the convergence and
stiffness update information.
If the job terminates before 100% of the loading is applied this convergence information can be used to
identify possible reasons. The default numerical method used by SOL 400 is the full Newton method with
the load increment size adjusted according to how many iterations are required to achieve convergence. The
incremental load size will be adjusted up if convergence is achieved easily, and will be reduced if more than
the target number of iterations is required to achieve convergence.
Convergence is achieved when the error factors on the selected criteria are below the required values (called
the convergence criteria). In this convergence and stiffness update information, the following values are
given:
LOAD STEP: percentage of total load reached at this incremental load value
NO. INC: increment number of this increment relative to all increments in this step
ITR: number of iterations required to get convergence in this load increment
Under - - - - - - ERROR FACTORS - - - - - Convergence / Error Factors that measure convergence
DISP: displacement vector residual
LOAD: load vector residual
Main Index
190 MSC Nastran: SOL 400 Getting Started Manual
Reviewing Convergence
Reviewing Convergence
From a user’s perspective, when performing a nonlinear simulation, the most difficult thing is to resolve
convergence problems encountered in the analysis.
In MSC Nastran, convergence may refer to:
Convergence in the iterative solver
Convergence in eigenvalue extraction
Convergence in aero-elasticity flutter calculation
Convergence in optimization
Convergence in equilibrium.
In this section, we will only be focusing on convergence equilibrium of the structural system. Figure 11-2 gives
a quick summpar of steps you should undertake.
For structural analysis, a lack of convergence implies that the numerical solution has not reached equilibrium
to the desired level of accuracy. Even if convergence is achieved, it is strongly recommended to determine how
many iterations were required and to consider if this was a reasonable number or an excessive number.
Recall that the number of iterations has a significant influence on the computational costs, so not only is there
the requirement that the solution converges, but that it also converges efficiently.
Main Index
CHAPTER 11 191
Trouble Shooting
Main Index
192 MSC Nastran: SOL 400 Getting Started Manual
Review Plot Results from Converged Increments
Main Index
CHAPTER 11 193
Trouble Shooting
Main Index
194 MSC Nastran: SOL 400 Getting Started Manual
Tips for Starting with Nonlinear Analysis
Files are automatically copied to where they are needed, the analysis is performed and pertinent information
is relayed back to you. When the analysis is complete files are returned/deleted. Time consuming system tasks
are reduced so that more time is available for productive engineering.
The Analysis Manager can also be used to monitor job progress (even non-jobs)
Shows all job and host information
Abort function can be used to stop running jobs
Output Messages
Additional Output
When convergence is a problem it is necessary to obtain as much information as possible on the numerical
process. The NLOPRM case control command is used to obtain additional result information.
In particular:
NLOPRM OUTCTRL=STD,INTERIM, NLDBG=NLBASIC,NRDBG, DBGPOST=LTIME
will provide the additional information to evaluate the performance. If contact bodies are present, it may be
preferable to use:
NLOPRM OUTCTRL=STD,INTERIM, NLDBG=NLBASIC,NRDBG, DVDBG,N3DMED
DBGPOST=LTIME
Main Index
CHAPTER 11 195
Trouble Shooting
Exit 0 – success. The job ran to completion and did everything you asked.
Exit 1 – syntax error in the input file. You should check the input syntax of
the line the error message points to, but it is likely that the actual
error was in the input block prior to where the message points.
Exit 1 with user fatal – typically means no convergence due to rigid body motions or a
EXCESSIVE PIVOT numerically ill-conditioned system. See recommendations for
RATIO in .f06 equilibrium below.
Exit 1 with other user – means the analysis ran into convergence problems part way through
fatal and did not complete. Interpret the error message, look for model
problems
For SOL 400 User http://simcompanion.mscsoftware.com/infocenter/
Fatal Messages (UFMs)
and their
Interpretation see:
If you receive a Nastran input from another source, it is strongly recommended that you read it into either
Patran or SimXpert and display the model to check for completeness. Any.sts file with nonzero total time
lines means there are converged increments. Display the converged increments to see what is going on.
Main Index
196 MSC Nastran: SOL 400 Getting Started Manual
Tips for Starting with Nonlinear Analysis
This exit message or a similar one may indicate excessive element deformation during a particular
load increment in a particular iteration. The way to get around this error is to reduce step size.
However, SOL 400 does that automatically and if the problem still persists, it gives this message. If
the program is unable to recover after several cutbacks, it implies that the last converged solution is
not really a good solution or the deformation is excessive for this mesh. The solution may be to
either use Restart or rerun-the complete analysis with a tighter tolerance.
It may be necessary to change the original mesh to anticipate where the deformation will occur. Also,
it should be recognized that lower-order triangular or tetrahedral elements (with PSLDN1) is the
preferred solution for large strain/distortion models.
User-specified Upper Limit on Number of Iterations: The solution may fail and give this message if
SOL 400 reaches the user-specified upper limit on the number of (Newton-Raphson) iterations
within a load increment. The default is set to 10. One of the first things to try, if you think the
problem is otherwise well posed (e.g., no other problems are known) would be to increase the value
(specified on NLSTEP) to a value of 20.
What is happening is that SOL400 keeps iterating and tries to converge to a solution for that
increment. If that does not happen, SOL 400 will cut back the load (by half ) and re-solve that
increment. Sometimes this is not enough to get convergence, and it will exit with a UFM.
Check the residual values in the ADVDBG section of the .f06 file to confirm this.
Minimum Time Step Size Exceeded/Time-Step Size Too Small: This exit message indicates that
SOL 400 cuts back to a time-step size too small for the analysis to continue. The load stepping
algorithm has a cut-back feature where the load step is automatically reduced (when the time step is
decreased, the factor is calculated internally based upon the minimum time step).
When an increment runs into these problems, it automatically cuts the load-step size and resolve
that increment. If the problem persists, it cut-back the load-step again. This happens until the limit
of the number of cut-backs is reached. This can result in a very small time step. In such a case, SOL
400 stops the analysis with a UFM.
Allowable Number of Bisections Exceeded: SOL 400's automatic load stepping scheme is set up
such that the applied load in an increment scales up (or down) depending on how easy (or difficult)
the solution was in the previous increment.
The degree of difficulty is determined based on the NDESIR parameter of NLSTEP: desired
number of iterations per increment (default = 4). SOL 400 will scale down the step size
until it reaches a lower limit on the step size (default = 0.001% of total time step) and then exit with
a UFM worded similar to the above. This is an indication to you that the analysis encountered some
difficulty at that stage.
Main Index
CHAPTER 11 197
Trouble Shooting
Main Index
Appendix A: Frequently Asked Questions
Main Index
CHAPTER A 199
Frequently Asked Questions
Q1: Is SOL 400 only for nonlinear analysis, or can I also do linear analysis?
SOL 400 can do linear, nonlinear, modal, buckling, and transient structural analysis along with a
host of other analysis types. SOL 400 is designed to be the only solution sequence you need. For a
complete list, see ANALYSIS case control entry documentation in the Quick Reference Guide.
Q2: Will I get the same results from a SOL 400 linear analysis as I do from a SOL 101 analysis?
SOL 400 has additional capabilities that require controls not available in SOL 101. This means the
additional controls for a SOL 400 analysis must be set in the same way as they are for a SOL 101
analysis.
Having said that, SOL 400 control defaults are set in such a way that you should get the same results
from a SOL 400 analysis that you do from a SOL 101 analysis.
Q3: How different is SOL 400 from SOL 600?
In SOL 600, MSC Nastran preprocesses the data and calls the Marc solver. In SOL 400, all MSC
Nastran infrastructure is available, while in SOL 600 only selected MSC Nastran infrastructure
capabilities have been integrated. This is why MSC recommends that SOL 400 be the default
solution method for solving nonlinear problems.
An extensive finite element library for building your simulation model, and a set of solution
procedures for the nonlinear analysis, which can handle very large matrix equations, are available in
both solution sequences of MSC Nastran Implicit Nonlinear.
Q4: Will I get the same results from a SOL 400 nonlinear analysis as I do from a SOL 106 static or SOL 129
transient nonlinear analysis?
There are additional capabilities and additional controls in SOL 400 that must be set in the same
way as they are for SOL 106 or SOL 129 in order to get the same results. Also, SOL 400 has a
complete set of large strain elements that make it unlikely that you will get the same solution unless
you are performing a small deformation analysis.
However, if you set your analysis up to use the same elements as the SOL 106 or 129 solution, and
you limit the conditions applied to the model you should get the same results as a SOL 106 or 129
analysis.
Q5: Will I get the same results from a SOL 400 linear, modal, buckling, or nonlinear analysis as I will from
a SOL 101, 103, 105, and 106 analysis?
SOL 400 encapsulates those solution sequences and should be able to reproduce their results.
Getting the same answers requires equivalence in element types, analysis control settings, and
procedures.
Q6: Is it possible to use the same model for linear and nonlinear analysis. If so what is the best way to add
the additional input required for the nonlinear solution to an existing linear model?
Yes. SOL 400 has been designed to allow you use the same model for linear and nonlinear analysis.
To accomplish this, do the following:
a. Take the linear model.
b. Read it into a graphical pre-and post processor like Patran.
Main Index
200 MSC Nastran: SOL 400 Getting Started Manual
Q10:If I am new to SOL 400 what's the best way to learn it?
There are a few ways of learning:
Main Index
CHAPTER A 201
Frequently Asked Questions
• The MSC documentation system has a complete set of solved SOL 400 problems in a document
called the MSC Nastran Demonstration Problems. The input files from these solved problems are
included in the documentation system.
• Patran has a set of demo problems that can be found in the Analysis menu, under the Run a
Demo tab. After running these demos, Patran can be used to interrogate the model and
investigate how the models are set up.
• The SimCompanion website has a knowledge base with answers to commonly asked questions,
as well as links to the MSC training courses on SOL 400.
Q11:Where can I find example input files for SOL 400 demo problems?
Example problems can be found in the MSC Nastran Demonstration Problems Manual. These
example problems also include input files. There is a test suite of example input files in the MSC
Nastran test problem library. Also, the Run a Demo problems in Patran will leave the input files
that were used to run the example.
Q12:What kind of results are supported in the op2, xdb, HDF5, and/or dball/master?
While you can request SOL 400 output in any valid MSC Nastran output format, only the
MASTER/DBALL and new OP2 (PARAM,POST,1) and HDF5 formats will contain all of the results
data blocks of SOL 400.
Specifically, the nonlinear stress/strain and contact results data blocks will be missing if any other
format is requested. Thus, MSC strongly recommends that you use MASTER/DBALL, new OP2,or
HDF5 output formats for SOL 400 solutions.
Q13:What do I do if I run the analysis and do not get a solution? What if it does not converge?
Check the solution files (*.f06, *.log, *.f04) to look for error messages. Typically there will be
a message telling why the solution did not converge.
There are a couple of answers to this question depending on why the analysis did not get a solution.
Reasons why an analysis does not result in a solution can be grouped into two categories.
• The first group relates to having a correctly formatted input file. These tend to be the most
common, and also the easiest to fix.
• The second group relates to having what might be called a well-posed problem.
This means that our problem does not violate any of the laws of physics, and that there actually
is a valid numerical solution to the problem we are trying to solve. A typical example of this type
of error would be trying to come up with a static solution to a problem that is not statically
determinate. There are methods used to determine if a model is properly constrained, such as
running a modal analysis to look for unconstrained rigid body modes.
Q14:How do I know if I need to do a nonlinear analysis?
While this is not an exact answer as loading conditions will affect this, you should do a nonlinear
analysis in one of the following conditions:
• If your strains are approaching 5%
• The deflection of any node in your model approaches 5% of the smallest dimension.
Main Index
202 MSC Nastran: SOL 400 Getting Started Manual
Q15:How do I know if linear contact will give me the right answer, or if I need to do a nonlinear analysis?
The term linear contact is an oxymoron. The very nature of contact problems are that finite
deformation of the nodes must be tracked to determine if the notes are in contact or not in contact.
This deformation is almost always large enough to invalidate linear contact solutions. The best use
for linear contact is to use it to connect dissimilar meshes or full assembly modeling where no
separation occurs.
Q16:What is the difference between linear buckling analysis and nonlinear buckling analysis?
A linear buckling analysis is based on the un-deformed configuration of the structure.
A nonlinear buckling analysis is based on a deformed shape of the structure. The deformation and
stresses make the structure stiffer (if in tension) and cause the natural frequencies to increase. Hence,
it is often necessary to do a nonlinear buckling analysis to get accurate mode shapes and frequencies.
Q17: What is a stress stiffened model analysis?
Before extracting the eigenvalues from the structural system, the model may be pre-loaded prior to
the modal extraction. Initial stress effects are then included in the stiffness which tend to raise the
eigenvalues. This effect tends to affect the eigenvalues more than the eigenvectors (mode shapes).
Q18:Will RBE's rotate with the rest of the model?
If a large displacement, nonlinear analysis is performed, then yes, the RBEs rotates with the model.
In a linear analysis, the displacement and rotation of the nodes is assumed to be infinitesimally small.
Q19:What are Marc elements? When would I use them? How do I activate them?
Most MSC Nastran elements, also called Advanced nonlinear elements or large strain elements were
originally formulated for linear analysis. At some point, some of those original elements were
modified for use in nonlinear analysis (SOLs 106 and 129), but certainly not all of them.
• Instead of developing new elements and capabilities from scratch for SOL 400, MSC Nastran
decided to copy the well-proven, robust, mature large deformation/large strain analysis
procedures and element formulations from Marc.
• Instead of developing new nonlinear element designations, which would have increased work to
convert a linear SOL 101 model to a SOL 400 nonlinear model, MSC Nastran decided to use the
same element designations and allow you to control the formulation used through secondary
property entries such as PBEMN1, PSHLN1, PSHNL2, and PSLDN1.
While this method gives you direct control over the element formulation, SOL 400 automatically
selects the proper formulation for the problem. This means that if the problem is a large
deformation, large strain problem, SOL 400 automatically uses the large strain element
formulation.
Unless you want to directly control the formulation (such as use reduced integration or some
other special formulation), it is best to allow MSC Nastran to select the formulation used. This
also means that no additional input is required.
Q20:Do shell and beam offsets rotate with the model?
Yes, if MDLPRM,OFFDEF,LROFF is set in the deck when parameter LGDISP>=0.
Q21:Will CBAR elements rotate with the model?
Main Index
CHAPTER A 203
Frequently Asked Questions
Main Index
204 MSC Nastran: SOL 400 Getting Started Manual
It is recommended that you run the example problems and make variations of these problems. When
learning a new MSC Nastran feature, engineers all too often generate a large problem using several
hundred or several thousand degrees of freedom as a test case. This practice has become the norm in
recent years with the advent of graphics preprocessors and automatic meshing.
Rarely is such a large model necessary to learn a new feature; in most cases, it just adds unnecessary
complexity. For this reason most of the examples in this book are small—generally less than 100
degrees of freedom.
To facilitate the use of these example problems, example problems referenced in this guide are
delivered with MSC Nastran system under the MSC_DOC_DIR/doc/linstat with the extensions
of .dat. Copy an example problem to your local directory, so you can see the files created, and you
don't inadvertently create files in the delivery directory.
Q31:How can I see a plot of load vs deflection for the rigid body that loaded my model?
Only available for load controlled rigid bodies, you can plot the results of the control grid of the
rigid body.
Q32:How can I monitor the nonlinear solution to see if it is progressing toward convergence?
Monitor the jobname.sts file with a text editor that automatically updates. As each load
increment converges a new line with relevant information will appear in the .sts file.
For details, see Chapter 6: Setting Up, Monitoring, and Debugging the Analysis.
Q33:Is it possible to see the results of any of the intermediate loading steps before the analysis is complete?
Yes. It is possible to have MSC Nastran write out the results of converged increments to intermediate
OP2 files. See the INTERM option of the OUTCTRL parameter under the NLOPRM entry in the MSC
Nastran Quick Reference Guide.
Q34:How can I find out what the normal and frictional forces or stresses are between the contact bodies?
Use the BOUTPUT case control output request and these values will be in the MASTER/DBALL or new
OP2 files along with the other output requested.
Q35:Is an RBE spider the best way to fix a surface of my structure?
It depends on the constraints on the surface. If all the grids on the surface have the same behavior in
some direction, an RBE is easy to use. It is also easy to glue the area to a rigid surface and use the
rigid surface controls. You may also use SPC/SPC1.
Q36:What are alternatives to an RBE spider for displacement control of part of my structure?
SPC/SPC1/SPCD/SPCR may be combined together for any complicated displacement control.
Glueing the area to a rigid surface and using the rigid surface controls is also easy.
Q37:Will I get the same answer if I apply my loads and constraints using a glued-on rigid surface as I would
using an RBE spider or displacement constraint?
Yes, if both loading methods are done correctly you will get the same results.
Q38:How do a simulate a complex sequence of loading, such as the process of an engine head going through
the sequence of bolt tightening, then heating up?
Main Index
CHAPTER A 205
Frequently Asked Questions
Break the loading up into discrete loadcases. For more information on this process, see Chapter 3:
Nonlinearity and Analysis Types in this manual for more information on this topic.
Q40:What is the simplest way to include contact in my model?
If you create contact bodies and run a SOL 400 nonlinear analysis, the contact interaction is
included by default, but does not include friction between the components. Patran has a tool under
Tools - Modeling that automatically creates contact bodies based on groups, materials, properties,
etc.)
To include friction, specify the friction model and friction coefficients. You can specify a global
model value, a value for each contact body, or for each contact pair—the most localized value, in the
most general to most specific order just provided, will be used.
Q41:What Thermal Solutions are supported in SOL 400?
SOL 400 supports steady state and transient thermal analysis procedures (ANALYSIS = HSTAT and
HTRAN), updated to include thermal contact and coupled thermal-structural interaction, plus the
Sinda RC network approach. Chapter 3: Nonlinearity and Analysis Types has a description of the SOL 400
thermal analysis capabilities, plus a description of the coupled multi-physics capabilities that are
supported.
Q42:How do the RCNS and RCNT RC Network analysis options in SOL 400 work?
These RC analysis options run a subset of the Sinda Solver. Standalone Sinda is a modern equivalent
of traditional Resistor-Capacitor type solver. It is also called Finite Difference Lumped Parameter
Network Solver. Simulation can include various aspects of conduction, convection, radiation and
with optional Fortran access predefined or custom loading function can be specified.
Using Patran with Sinda preference, or in this case from MSC Nastran, an RC network is generated
from traditional shaped elements. MSC Nastran implementation does not include the Fortran
access nor standard/custom loading function. Loading must be manually defined but the same effect
can be achieved. Additionally, MSC Nastran implementation does include automated radiation
refinement techniques available in Sinda.
MSC Nastran Implicit Nonlinear (SOL 600) User’s Guide
Main Index