FEA H2: EXCENTRICALLY LOADED U-SECTION BEAM

Initial setup

! Isometric view
/VIEW/VIEW, WN, XV, YV, ZV Defines the viewing direction for the display.,1,1,1,1
/ANG/ANGLE, WN, THETA, Axis, KINCR Rotates the display about an axis.,1
/REP/REPLOT, Label Automatically reissues the last display command for convenience.,FAST
! 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
! Constant parameters
! Units: [N],[Nmm],[N/mm^2=MPa],[N/mm]
ELAST=2E5
NU=0.3
RHO=7.8e-9
T=10
F=3000

SETUP: element types, materials, layup

KEYWKEYW, Keyword, KEY Sets a keyword used by the GUI for context filtering (GUI).,PR_STRUC,1! Preferences / Structural
/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,SHELL181

!
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,,ELASTELAST = 2E5
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,1,,NUNU = 0.3
!
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,,RHORHO = 7.8e-9
!
! Layup definition
SECTSECTYPE, SECID, Type, Subtype, Name, REFINEKEY Associates section type information with a section ID number.,1,shell,,
SECDATASECDATA, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8, VAL9, VAL10, VAL11, VAL12 Describes the geometry of a section.,t,1,0.0,3
SECOFFSETSECOFFSET, Location, OFFSET1, OFFSET2, CG-Y, CG-Z, SH-Y, SH-Z Defines the section offset for cross sections.,MID
SECCONTROLSECCONTROL, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8, VAL9, VAL10, VAL11, VAL12, VAL13 Supplements or overrides default section properties.
!
! Define keypoints
KK, NPT, X, Y, Z Defines a keypoint.,1,0,0,0,
KK, NPT, X, Y, Z Defines a keypoint.,2,0,0,1000,
KK, NPT, X, Y, Z Defines a keypoint.,3,100,0,0
KK, NPT, X, Y, Z Defines a keypoint.,4,0,200,0
KK, NPT, X, Y, Z Defines a keypoint.,5,100,200,0,
!
! Create lines
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,5,4
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,4,1
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,1,3
! Linene along the beam axis
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,1,2
! Create the surface model by extruding the contour of the U-section
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
! Connect the extruded areas
AGLUEAGLUE, NA1, NA2, NA3, NA4, NA5, NA6, NA7, NA8, NA9 Generates new areas by "gluing" areas.,ALL
!
! Set the number of elements on the lines
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,1,,,10
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,2,,,20
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,3,,,10
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,5,,,10
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,6,,,60
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,7,,,60
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,8,,,20
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,9,,,60
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,10,,,10
!
! Mesh the areas using mapped meshing algorithm
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
! Reset mesh option
MSHKEYMSHKEY, KEY Specifies whether free meshing or mapped meshing should be used to mesh a model.,0
!
! Select nodes at z = 0
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Z,0,0
!
! Constrain all DoFs at this location
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,ALL,0
!
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
! select nodes at x = 0, z = 1000
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,0,0
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,R,LOC,Z,1000,1000
! Apply the calculated force at every node at this location
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,ALL,FY,-FF = 3000/21
! Select everything
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
FINISHFINISH Exits normally from a processor.
!
! Solve the finite element equation K.U = F
/SOL/SOLU Enters the solution processor.
SOLVESOLVE, Action Starts a solution.
! Post processing
/POST1/POST1 Enters the database results postprocessor.
! Save the displaced structure as bmp
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,X,0,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,disp_x,bmp
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,Y,0,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,disp_y,bmp
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,Z,0,1.0
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,disp_z,bmp
! Plott the stress values at the longitudinal edge of the beam
! x = 0, y = 0
PATHPATH, NAME, nPts, nSets, nDiv Defines a path name and establishes parameters for the path.,2,2,30,100
PPATHPPATH, POINT, NODE, X, Y, Z, CS Defines a path by picking or defining nodes, or locations on the currently active working plane, or by entering specific coordinate locations.,1,,0,0,0
PPATHPPATH, POINT, NODE, X, Y, Z, CS Defines a path by picking or defining nodes, or locations on the currently active working plane, or by entering specific coordinate locations.,2,,0,0,1000
PDEFPDEF, Lab, Item, Comp, Avglab Interpolates an item onto a path.,SX,S,X,AVG
PDEFPDEF, Lab, Item, Comp, Avglab Interpolates an item onto a path.,SY,S,Y,AVG
PDEFPDEF, Lab, Item, Comp, Avglab Interpolates an item onto a path.,SZ,S,Z,AVG
! Save the stress distribution along the path as bmp
PLPATHPLPATH, Lab1, Lab2, Lab3, Lab4, Lab5, Lab6 Displays path items on a graph.,SX
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,sx,bmp
!
PLPATHPLPATH, Lab1, Lab2, Lab3, Lab4, Lab5, Lab6 Displays path items on a graph.,SY
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,sy,bmp
!
PLPATHPLPATH, Lab1, Lab2, Lab3, Lab4, Lab5, Lab6 Displays path items on a graph.,SZ
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,sz,bmp
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
!
! Modal analysis, first 10 natural frequencies
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,0,0,,OFF
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.,0
!
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
! Save the eigenshapes as bmp files
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,III,1,10,1
  PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,Y,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.