FEA D4: Transient dynamic analysis of a built-in beam

! Created by Gergely Gyebroszki, 2018.04.19.

Units in mm/MPa

S=10
E=2.1E5
NU=0.3
RHO=7.85E-9
P0=3
T1=0.05
T2=1.0
BETA=0.0002
TSTEP=0.005

SETUP: Element types, Materials, Sections

KEYWKEYW, Keyword, KEY Sets a keyword used by the GUI for context filtering (GUI).,PR_STRUC,1
/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
KEYOPTKEYOPT, ITYPE, KNUM, VALUE Sets element key options.,1,3,3! Plane stress with thickness
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,1,SS = 10! Real constant to define thickness
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,,EE = 2.1E5
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
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DENS,1,,RHORHO = 7.85E-9

MODELLING

BLC4BLC4, XCORNER, YCORNER, WIDTH, HEIGHT, DEPTH Creates a rectangular area or block volume by corner points.,0,0,1000,200
BLC4BLC4, XCORNER, YCORNER, WIDTH, HEIGHT, DEPTH Creates a rectangular area or block volume by corner points.,1000,0,300,200

! Set explicit element size of 50mm to achieve same divisions
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,ALL,50
AMESHAMESH, NA1, NA2, NINC Generates nodes and area elements within areas.,ALL
NUMMRGNUMMRG, Label, TOLER, GTOLER, Action, Switch Merges coincident or equivalently defined items.,NODE,0.01,,,LOW
! Built-in end at X=0
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 10,LOC,X,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,ALL
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
FINISHFINISH Exits normally from a processor.

SOLUTION: Transient

/SOL/SOLU Enters the solution processor.

! Load step 1
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,4 ! Analysis type: Transient
TRNOPTTRNOPT, Method, MAXMODE, --, MINMODE, MCout, TINTOPT, VAout Specifies transient analysis options.,FULL! Solution method: FULL
LUMPMLUMPM, Key Specifies a lumped mass matrix formulation.,0 ! Lumped mass approximation: OFF
TIMETIME, TIME Sets the time for a load step.,T1 ! Time at the end
AUTOTSAUTOTS, Key Specifies whether to use automatic time stepping or load stepping.,1 ! Auto time-stepping: ON
DELTIMDELTIM, DTIME, DTMIN, DTMAX, Carry Specifies the time step sizes to be used for the current load step.,TSTEP ! Time step size
OUTRESOUTRES, Item, Freq, Cname, -- , NSVAR, DSUBres Controls the solution data written to the database.,ALL,1 ! Write every substep
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,BETA ! Damping
KBCKBC, KEY Specifies ramped or stepped loading within a load step.,0 ! Ramped loading (instead of stepped)
! Pressure load on right-top line
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,SS = 10,LOC,Y,200,200! Select lines with Y in (200,200)
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,R,LOC,X,1000,1300 ! Reduce selection to lines with X in (1000,1300)
SFLSFL, Line, Lab, VALI, VALJ, VAL2I, VAL2J Specifies surface loads on lines of an area.,ALL,PRES,P0P0 = 3
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,1
! Load step 2 (only the time changes)
TIMETIME, TIME Sets the time for a load step.,T2T2 = 1.0
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,2
! Solve from LS file 1-2
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,1,2,1
FINISHFINISH Exits normally from a processor.

TIME HISTORY POSTPROCESSING

/POST26/POST26 Enters the time-history results postprocessor.
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,2,NODE(1300,100,0),U,Y,UY_A
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,3,NODE(1300,100,0),V,Y,VY_A
/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
/YRANGE/YRANGE, YMIN, YMAX, NUM Specifies a linear ordinate (Y) scale range.,-4,-3.7,1
PLVARPLVAR, NVAR1, NVAR2, NVAR3, NVAR4, NVAR5, NVAR6, NVAR7, NVAR8, NVAR9, NVAR10 Displays up to ten variables in the form of a graph.,2
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,Disp_vs_time_A,bmp
/YRANGE/YRANGE, YMIN, YMAX, NUM Specifies a linear ordinate (Y) scale range.,-140,60,1
PLVARPLVAR, NVAR1, NVAR2, NVAR3, NVAR4, NVAR5, NVAR6, NVAR7, NVAR8, NVAR9, NVAR10 Displays up to ten variables in the form of a graph.,3
/IMAGE/IMAGE, Label, Fname, Ext, -- Allows graphics data to be captured and saved.,SAVE,Veloc_vs_time_A,bmp
FINISHFINISH Exits normally from a processor.