! FEA T1: ANALYSIS OF AN AXISYMMETRIC DRUM

Initial setup

! White canvas
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,100,100,100,0
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,0,0,0,15
! Define constant parameters
! Units: [N],[Nmm],[N/mm^2=MPa],[N/mm]
D1 = 40
D2 = 50
D3 = 60
D4 = 100
D5 = 160
R1 = 20
R2 = 15
H1 = 40
H2 = 100
TH = 5
ELAST = 2.1E5
NU = 0.33
RHO = 7.86E-9
DELTA = 0.005
OM = 200

SETUP: element types, materials,

/PREP7/PREP7 Enters the model creation preprocessor.
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,1,PLANE182

! set axysimmetric option
KEYOPTKEYOPT, ITYPE, KNUM, VALUE Sets element key options.,1,3,1
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EX,1,,ELAST
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,1,,NU
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DENS,1,,rho
! create keypoints
KK, NPT, X, Y, Z Defines a keypoint.,1,D1/2,H1,0
KK, NPT, X, Y, Z Defines a keypoint.,2,D2/2,H1,0
KK, NPT, X, Y, Z Defines a keypoint.,3,D1/2,0,0
KK, NPT, X, Y, Z Defines a keypoint.,4,D1/2,TH,0
KK, NPT, X, Y, Z Defines a keypoint.,5,D3/2,0,0
KK, NPT, X, Y, Z Defines a keypoint.,6,D5/2-TH-R1,0,0
KK, NPT, X, Y, Z Defines a keypoint.,7,D5/2-TH-R1,R1 R1 = 20+ TH,0
KK, NPT, X, Y, Z Defines a keypoint.,8,D5/2,R1 R1 = 20+ TH,0
KK, NPT, X, Y, Z Defines a keypoint.,9,D5/2,H2-TH-R2,0
KK, NPT, X, Y, Z Defines a keypoint.,10,D5/2-TH-R2,H2,0
KK, NPT, X, Y, Z Defines a keypoint.,11,D4/2,H2,0
KK, NPT, X, Y, Z Defines a keypoint.,12,D5/2-TH-R2,H2-TH-R2,0
! create lines
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,4,3
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,3,5
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,5,6
LARCLARC, P1, P2, PC, RAD Defines a circular arc.,6,8,7,R1+TH
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,8,9
LARCLARC, P1, P2, PC, RAD Defines a circular arc.,9,10,12,R2+TH
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,10,11
ADRAGADRAG, NL1, NL2, NL3, NL4, NL5, NL6, NLP1, NLP2, NLP3, NLP4, NLP5, NLP6 Generates areas by dragging a line pattern along a path.,1,,,,,,2,3,4,5,6,7
AA, P1, P2, P3, P4, P5, P6, P7, P8, P9, P10, P11, P12, P13, P14, P15, P16, P17, P18 Defines an area by connecting keypoints.,1,2,13,4
AGLUEAGLUE, NA1, NA2, NA3, NA4, NA5, NA6, NA7, NA8, NA9 Generates new areas by "gluing" areas.,ALL
! set element size on specific lines for mapped mesh
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,1,,,8
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,III,8,23,3
  LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,III,,,8
*ENDDO*ENDDO Ends a do-loop and starts the looping action.
!*
! proportional element size
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,III,9,24,3
  LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,III,TH/8
*ENDDO*ENDDO Ends a do-loop and starts the looping action.
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,28,TH/8
! create mapped mesh on every area
MSHKEYMSHKEY, KEY Specifies whether free meshing or mapped meshing should be used to mesh a model.,1
AMESHAMESH, NA1, NA2, NINC Generates nodes and area elements within areas.,ALL
DKDK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6 Defines DOF constraints at keypoints.,3,UY,0

load case 1

DLDL, LINE, AREA, Lab, Value1, Value2 Defines DOF constraints on lines.,1,,UX,DELTA
DLDL, LINE, AREA, Lab, Value1, Value2 Defines DOF constraints on lines.,28,,UX,DELTA
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,1

load case 2

! apply global angular velocity about the Y axis
DLDELEDLDELE, LINE, Lab Deletes DOF constraints on a line.,1,UX
DLDELEDLDELE, LINE, Lab Deletes DOF constraints on a line.,28,UX
OMEGAOMEGA, OMEGX, OMEGY, OMEGZ Specifies the rotational velocity of the structure.,0,OM,0
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,2
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,1,2,1
FINISHFINISH Exits normally from a processor.
/POST1/POST1 Enters the database results postprocessor.
LCDEFLCDEF, LCNO, LSTEP, SBSTEP, KIMG Creates a load case from a set of results on a results file.,1,1,Last,
LCDEFLCDEF, LCNO, LSTEP, SBSTEP, KIMG Creates a load case from a set of results on a results file.,2,2,Last
! 1st load case
LCASELCASE, LCNO Reads a load case into the database.,1
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,S,INT,2,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,S_INT_LCASE_1,bmp
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,X,2,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,DISP_X_LCASE_1,bmp
! 2nd load case
LCASELCASE, LCNO Reads a load case into the database.,2
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,S,INT,2,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,S_INT_LCASE_2,bmp
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,X,2,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,DISP_X_LCASE_2,bmp
! 1st and 2nd load case
LCASELCASE, LCNO Reads a load case into the database.,1
LCOPERLCOPER, Oper, LCASE1, Oper2, LCASE2 Performs load case operations.,ADD,2
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,X,2,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,DISP_X_LCASE_12,bmp
FINISHFINISH Exits normally from a processor.
! modal analysis with prestress effect
/PREP7/PREP7 Enters the model creation preprocessor.
DLDELEDLDELE, LINE, Lab Deletes DOF constraints on a line.,1,UX
DLDELEDLDELE, LINE, Lab Deletes DOF constraints on a line.,28,UX
OMEGAOMEGA, OMEGX, OMEGY, OMEGZ Specifies the rotational velocity of the structure.,0,200,0
FINISHFINISH Exits normally from a processor.
! calculate initial stress state
/SOL/SOLU Enters the solution processor.
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,0
PSTRESPSTRES, Key Specifies whether prestress effects are calculated or included.,1
SOLVESOLVE, Action Starts a solution.
FINISHFINISH Exits normally from a processor.
! solve modal analysis using prestress effect
/SOL/SOLU Enters the solution processor.
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,2
MODOPTMODOPT, Method, NMODE, FREQB, FREQE, Cpxmod, Nrmkey, ModType, BlockSize, --, --, Scalekey Specifies modal analysis options.,LANB,10
EQSLVEQSLV, Lab, TOLER, MULT, --, KeepFile Specifies the type of equation solver.,SPAR
MXPANDMXPAND, NMODE, FREQB, FREQE, Elcalc, SIGNIF, MSUPkey, ModeSelMethod Specifies the number of modes to expand and write for a modal or buckling analysis.,10,,,1
LUMPMLUMPM, Key Specifies a lumped mass matrix formulation.,0
PSTRESPSTRES, Key Specifies whether prestress effects are calculated or included.,1
MODOPTMODOPT, Method, NMODE, FREQB, FREQE, Cpxmod, Nrmkey, ModType, BlockSize, --, --, Scalekey Specifies modal analysis options.,LANB,10,0,0,,OFF
SOLVESOLVE, Action Starts a solution.
FINISHFINISH Exits normally from a processor.
/POST1/POST1 Enters the database results postprocessor.
SETSET, Lstep, Sbstep, Fact, KIMG, TIME, ANGLE, NSET, ORDER Defines the data set to be read from the results file.,FIRST
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,III,1,10
  PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,SUM,2,1.0
  /IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,strcat(chrval(III),'_eigenshape'),bmp
  SETSET, Lstep, Sbstep, Fact, KIMG, TIME, ANGLE, NSET, ORDER Defines the data set to be read from the results file.,NEXT
*ENDDO*ENDDO Ends a do-loop and starts the looping action.