Girkmann Problem Using Axisymmetric Shell Elements
Girkmann Problem Using Axisymmetric Shell Elements
Summary 1178
Introduction 1179
Problem Statement 1179
Modeling Details 1179
Results and Discussions 1182
Conclusion 1184
Input File(s) 1184
Video 1185
References 1185
1178 MD Demonstration Problems
CHAPTER 62
Summary
Title Chapter 62: Girkmann Problem using Axisymmetric Shell Elements
Contact Features • Use of axisymmetric shell elements
• Use of MPCs
Geometry
Introduction
This example showcases the capability of CAXISYM elements. The problem demonstrated here was first discussed
by PitkÄaranta et al. (References). The FEA solution of the Girkmann problem was obtained and compared with the
solutions obtained by classical methods.
Problem Statement
A spherical shell of thickness h = 0.06m , crown radius Rc = 15.00m is connected to a stiffening ring at the meridional
angle = 2 9 40 . The dimensions of the ring are: a = 0.60m and b = 0.50m . The radius of the midsurface of the
sperical shell is Rm = Rc sin .
The notation is shown in Figure 62-1. The z axis is the axis of rotational symmetry.
Consider gravity loading only. The equivalent (homogenized) unit weight of the material comprised of the shell and
the cladding is 32.69 KN\m3 . Assume that uniform normal pressure PAB is acting at the base AB of the stiffening ring.
The resultant of PAB equals the weight of the structure. Assume that the stiffening ring is weightless. The goals of the
computations are as follows:
• To find the shear force at the interface between spherical shell and stiffening ring.
• To find the maximum bending moment.
• To find the meridional angle at which the maximum bending moment occurs.
Modeling Details
The spherical shell was modeled with axisymmetric shell element and stiffening ring with axisymmetric solid
elements as shown in Figure 62-2.
1180 MD Demonstration Problems
CHAPTER 62
Element Modeling
The spherical shell was modeled with 2-node axisymmetric shell elements CAXISYM. The properties for CAXISYM
were entered as given below:
PAXISYM 2 1 0.06 prop_axi
The CAXISYM elements were defined as given below:
CAXISYM 769 2 830 831
CAXISYM 770 2 831 832
CAXISYM 771 2 832 833
CAXISYM 772 2 833 834
CAXISYM 773 2 834 835
The solid stiffener was modeled with axisymmetric solid element CQUADX. First, the properties PLPLANE for the
CQUADX were defined:
PLPLANE 1 2 prop_axi
PSHLN2 1 2 +
+ C4 AXSOLID L prop_axi
The CQUADX elements were defined as given below,
CQUADX 1 1 1 10 11 2
CQUADX 2 1 2 11 12 3
CQUADX 3 1 3 12 13 4
…
Material Properties
The material properties were modeled by using MAT1 entry as given below:
MAT1 12.059+10 0.0 3335.71 0.0 mat_axi_
MAT1 22.059+10 0.0 mat_axi_
As Nastran applies gravity over the entire model by default, two material cards were defined. The material without
density was used for solid stiffener modeled with CQUADX so that the gravity load does not get applied on it.
CHAPTER 62 1181
Girkmann Problem using Axisymmetric Shell Elements
Analysis Parameters
Subcase corresponding to the loadcase that contains above loads and the boundary conditions were defined and
analysis was carried out by SOL 400 by using following parameters:
$# NASTRAN Control Section
$# File Management Section
NASTRAN SYSTEM (316)=19
ASSIGN OUTPUT2 = 'girk_fo_s400.op2', UNIT = 12, FORM = FORMATTED
$# Executive Control Section
SOL 400
CEND
ECHO = NONE
$# Case Control Section
SUBCASE 1
STEP 1
$ Subcase name: DefaultLoadCase
$LBCSET SUBCASE1 DefaultLbcSet
TITLE=DefaultLoadCase
SET 10 = 769 THRU 1123
SET 20 = 830 THRU 1185
MPC = 33
SPC = 1
LOAD = 6
DISPLACEMENT (SORT1, PLOT, PRINT, REAL) =ALL
NLSTRESS (SORT1, PLOT, PRINT, REAL, VONMISES, CENTRE) =ALL
GPFORCE (PLOT, PRINT) =20
ANALYSIS = LNSTAT
The f06 file generated out of the analysis was studied carefully and the STRESS-XY for the axisymmetric shell
element CAXISYM 1123 were noted as
CHAPTER 62 1183
Girkmann Problem using Axisymmetric Shell Elements
The grid point force balance for the interface element CAXISYM 1123 is as given below,:
The bending stress values for all CAXISYM elements are calculated like this and the node where the max. bending
stress occurs was noted. At that particular node, the bending moment would be maximum. Figure 62-5 shows the
variation of bending moment/thickness (N.m/m) with meridional angle.
The maximum bending moment of 255.126 N.m/m occurs at the meridional angle of 38.15° as shown in Figure 62-5.
Conclusion
The following table shows the comparison of the results obtained from FEA modeling with that of theoretical one,
Input File(s)
Files Description
nug_62.bdf MD Nastran input for Girkmann problem using axisymmetric shell and solid elements.
CHAPTER 62 1185
Girkmann Problem using Axisymmetric Shell Elements
Video
Click on the image or caption below to view a streaming video of this problem; it lasts approximately 20 minutes and
explains how the steps are performed.
References
PitkÄaranta, J., Babuska, I. and Szabo, B. The Problem of Verification with Reference to the Girkmann Problem.
IACM Expressions. 24, January 2009, 14-15.