FEA NL4: DYNAMIC ANALYSIS OF A SIMPLY SUPPORTED BEAM

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
! 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
! Define constant parameters
! Units: [N],[Nmm],[N/mm^2=MPa],[N/mm]
L=1000
A=20
AREA=400
Iy=1.333E4
Iz=1.333E4
ELAST=2.0E5
NU=0.3
RHO=8.0E-9
P0=1
ELEMNUM=20

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,BEAM188
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 = 2.0E5
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 = 8.0E-9
SECTYPESECTYPE, SECID, Type, Subtype, Name, REFINEKEY Associates section type information with a section ID number.,1,BEAM,RECT
SECOFFSETSECOFFSET, Location, OFFSET1, OFFSET2, CG-Y, CG-Z, SH-Y, SH-Z Defines the section offset for cross sections.,CENT
SECDATASECDATA, VAL1, VAL2, VAL3, VAL4, VAL5, VAL6, VAL7, VAL8, VAL9, VAL10, VAL11, VAL12 Describes the geometry of a section.,AA = 20,AA = 20
KK, NPT, X, Y, Z Defines a keypoint.,1,0,0,0
KK, NPT, X, Y, Z Defines a keypoint.,2,1000,0,0
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,1,2
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,1,,,10
LMESHLMESH, NL1, NL2, NINC Generates nodes and line elements along lines.,1
DKDK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6 Defines DOF constraints at keypoints.,1,UX,0
DKDK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6 Defines DOF constraints at keypoints.,1,UY,0
DKDK, KPOI, Lab, VALUE, VALUE2, KEXPND, Lab2, Lab3, Lab4, Lab5, Lab6 Defines DOF constraints at keypoints.,2,UY,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UZ,0
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,4
TRNOPTTRNOPT, Method, MAXMODE, --, MINMODE, MCout, TINTOPT, VAout Specifies transient analysis options.,FULL
LUMPMLUMPM, Key Specifies a lumped mass matrix formulation.,0
DELTIMDELTIM, DTIME, DTMIN, DTMAX, Carry Specifies the time step sizes to be used for the current load step.,0.005
OUTRESOUTRES, Item, Freq, Cname, -- , NSVAR, DSUBres Controls the solution data written to the database.,ALL,ALL
AUTOTSAUTOTS, Key Specifies whether to use automatic time stepping or load stepping.,1
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,0.0002
TIMETIME, TIME Sets the time for a load step.,0.01
SFBEAMSFBEAM, Elem, LKEY, Lab, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST, LENRAT Specifies surface loads on beam and pipe elements.,ALL,2,PRES,p0
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,1
KBCKBC, KEY Specifies ramped or stepped loading within a load step.,0
TIMETIME, TIME Sets the time for a load step.,0.1
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,2
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,1,2,1,
GPLOTGPLOT Controls general plotting.
FINISHFINISH Exits normally from a processor.
/POST26/POST26 Enters the time-history results postprocessor.

! time history post processing comes here
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
SFEDELESFEDELE, ELEM, LKEY, Lab Deletes surface loads from elements.,ALL,2,PRES
! Modal analysis
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.,0
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.,LIST,2
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
! harmonic analysis
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,3
*AFUN*AFUN, Lab Specifies units for angular functions in parameter expressions.,DEG
AMPL=1
PHASE=120
FR=AMPL*COS(PHASE)
FI=AMPL*SIN(PHASE)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,,,3,ELEMNUMELEMNUM = 20 ! has to be an even number and larger than 2 to work+1,1
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,ALL,FY,FRFR = AMPL*COS(PHASE)/(ELEMNUM-1),FIFI = AMPL*SIN(PHASE)/(ELEMNUM-1)
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
HARFRQHARFRQ, FREQB, FREQE, --, LogOpt, FREQARR, Toler Defines the frequency range in a harmonic analysis.,0,500,
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,501
KBCKBC, KEY Specifies ramped or stepped loading within a load step.,1
SOLVESOLVE, Action Starts a solution.
FINISHFINISH Exits normally from a processor.
/POST26/POST26 Enters the time-history results postprocessor.
! time history post processing goes here