tbpt,,0.006422389,31.73! set the axles labels for the stress-strain curve plot
/axlab,x,log strain (n/m^2)
/axlab,y,true stress (n/m^2)tbpl,kinh,1! plot and verify the material stress-strain curve
! define a rectangle which is the axisymmetric cross section of the plate.
! the rectangle has a length equal to the radius of the plate and a height equal
! to the thickness of the plate
rect,,radius,,thick
! select the left and right bounding lines of the created rectangle and set
! the line division to 8 (8 elements through the thickness of the plate)
flst,5,2,4,orde,2
fitem,5,2
fitem,5,4
cm,_y,line
lsel, , , ,p51x
cm,_y1,line
cmsel,,_y
lesize,_y1, , ,8,1,
cmdel,_y
cmdel,_y1
! select the top and bottom bounding lines of the created rectangle and set
! the line division to 40 (40 elements through the radius of the plate)flst,5,2,4,orde,2
fitem,5,1
fitem,5,3
cm,_y,line
lsel, , , ,p51x
cm,_y1,line
cmsel,,_y
lesize,_y1, , ,40,1,
cmdel,_y
cmdel,_y1
cm,_y,area
asel, , , , 1
cm,_y1,area
chkmsh,'area'
cmsel,s,_y
amesh,all
cmdel,_y
cmdel,_y1
cmdel,_y2
fini
/solve
nlgeom,on! turn on geometric nonlinearity
! get the node numbers for the nodes located at the top
! of the axis of symmetry and at bottom right of the modelntop = node(0,thick,0)
nright = node(radius,0,0)
! activate the monitoring of the displacement and reaction force histories
! during the analysis. the histories are written out to the monitor file ratch.mntr
monitor,1,ntop,uy
monitor,2,nright,fy
outres,all,all! output all the results for all substeps to theRelease 15.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential informationNonlinear Structural Analysis