NAFEMS LE11: Thick Cylinder Thermal Stress
Theory and Physics
Overview
Teacher! Today's topic is NAFEMS LE11: Thermal Stress in a Thick-Walled Cylinder, right? What is it about?
The NAFEMS LE11 benchmark. A steady-state temperature distribution is applied to a thick-walled cylinder. It verifies the accuracy of thermal stress. The reference solution is σ_yy = −105 MPa.
Governing Equations
So, if you cut corners on the thermal stress of a thick-walled cylinder, you'll face consequences later. I'll keep that in mind!
Discretization Method
How do you actually solve these equations on a computer?
We use spatial discretization by the Finite Element Method (FEM). We assemble the element stiffness matrices and construct the global stiffness equation.
We transform it into the weak form (variational form) and use formulation by the Galerkin method using test functions and shape functions. The choice of element type (low-order elements vs. higher-order elements, full integration vs. reduced integration) directly affects the trade-off between solution accuracy and computational cost.
Matrix Solution Algorithms
What exactly are matrix solution algorithms?
Solve the system of equations using direct methods (LU decomposition, Cholesky decomposition) or iterative methods (CG method, GMRES method). For large-scale problems, preconditioned iterative methods are effective.
| Solver | Classification | Memory Usage | Applicable Scale |
|---|---|---|---|
| LU Decomposition | Direct Method | O(n²) | Small to Medium Scale |
| Cholesky Decomposition | Direct Method (Symmetric Positive Definite) | O(n²) | Small to Medium Scale |
| PCG Method | Iterative Method | O(n) | Large Scale |
| GMRES Method | Iterative Method | O(n·m) | Large Scale / Non-symmetric |
| AMG Preconditioner | Preprocessing | O(n) | Very Large Scale |
So, if you cut corners on the finite element method part, you'll face consequences later. I'll keep that in mind!
Implementation in Commercial Tools
So, what software can be used to do NAFEMS LE11: Thermal Stress in a Thick-Walled Cylinder?
| Tool Name | Developer/Current Status | Main File Formats |
|---|---|---|
| MSC Nastran / NX Nastran | MSC Nastran (Hexagon), NX Nastran (Siemens Digital Industries Software) | .bdf, .dat, .f06, .op2, .pch |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| Ansys Mechanical (formerly ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| COMSOL Multiphysics | COMSOL AB | .mph |
Vendor Lineage and Product Integration History
Do the origins of each software have dramatic stories?
MSC Nastran / NX Nastran
Next is the story about MSC Nastran, right? What's it about?
Developed in the 1960s as NASA Structural Analysis (NASTRAN). Commercialized by MSC Software, later UGS (now Siemens) derived NX Nastran. MSC was acquired by Hexagon AB in 2017.
Current Affiliation: MSC Nastran (Hexagon), NX Nastran (Siemens Digital Industries Software)
Abaqus FEA (SIMULIA)
What exactly is Abaqus FEA?
Developed in 1978 by HKS (Hibbitt, Karlsson & Sorensen). Acquired by Dassault Systèmes in 2005 and integrated into the SIMULIA brand.
Current Affiliation: Dassault Systèmes SIMULIA
Wait, wait, structural analysis... so, can it also be used for cases like this?
Ansys Mechanical (formerly ANSYS Structural)
Tell me about "Ansys Mechanical"!
Developed in 1970 by Swanson Analysis Systems Inc. (SASI). Based on APDL (Ansys Parametric Design Language).
Current Affiliation: Ansys Inc.
Wow, the story about structural analysis is super interesting! Tell me more.
File Formats and Interoperability
Are there any points to note when transferring data between different software?
| Format | Extension | Type | Overview |
|---|---|---|---|
| STEP | .stp/.step | Neutral CAD | ISO 10303 compliant 3D CAD data exchange format. Supports geometry + PMI. |
| IGES | .igs/.iges | Neutral CAD | Early CAD data exchange standard. Has issues with surface data compatibility. Transition to STEP is progressing. |
| VTK | .vtk/.vtu | Visualization | Visualization Toolkit format. Used by ParaView, etc. |