FEA D7: Impact loading on a beam

! Created by Gergely Gyebroszki, 2018.04.18.

Units in SI: [N],[N/m^2=Pa],[N/m]

H=0.025
E=2E11
NU=0.32
RHO=7850
A=6.25E-4
IZ=32.552E-9
FRCY=450
T1=0.0001
T2=0.0002
T3=0.5
SUBSTEPS1=20
SUBSTEPS2=20
SUBSTEPS3=200
DAMPING=0.0005

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,BEAM3! BEAM3 element type is not available in the GUI
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 = 2E11
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,1,,NUNU = 0.32
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DENS,1,,RHORHO = 7850
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,1,AA = 6.25E-4,IZIZ = 32.552E-9,HH = 0.025

MODELLING: Geometry and meshing

KK, NPT, X, Y, Z Defines a keypoint.,,-0.5
KK, NPT, X, Y, Z Defines a keypoint.,,0.5
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,,,30
LMESHLMESH, NL1, NL2, NINC Generates nodes and line elements along lines.,ALL

! Select nodes at X=-0.5 and also select at X=0.5
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,-0.5
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,AA = 6.25E-4,LOC,X,0.5
! We can use D on ALL nodes when applying constraints, but don't forget to select everything afterwards
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: Time history is given by 3 time points

/SOL/SOLU Enters the solution processor.
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,4! 4: Transient

! Time point at T1
TRNOPTTRNOPT, Method, MAXMODE, --, MINMODE, MCout, TINTOPT, VAout Specifies transient analysis options.,FULL
LUMPMLUMPM, Key Specifies a lumped mass matrix formulation.,0
AUTOTSAUTOTS, Key Specifies whether to use automatic time stepping or load stepping.,0! Turn off auto time stepping
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,SUBSTEPS1SUBSTEPS1 = 20,0,0
OUTRESOUTRES, Item, Freq, Cname, -- , NSVAR, DSUBres Controls the solution data written to the database.,ALL,1
TIMETIME, TIME Sets the time for a load step.,T1T1 = 0.0001 ! Time when Fy reaches its maximum value
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,0
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,ALL,FY,-FRCYFRCY = 450 ! Max value of Fy
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,1
! Also create a version of this load step with damping
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,DAMPINGDAMPING = 0.0005
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,4
! Time point at T2
TIMETIME, TIME Sets the time for a load step.,T2T2 = 0.0002 ! Time when Fy reaches zero
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,SUBSTEPS2SUBSTEPS2 = 20,0,0
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,0
FDELEFDELE, NODE, Lab, NEND, NINC, Lkey Deletes force loads on nodes.,ALL,FY
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,2
! Also create a version of this load step with damping
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,DAMPINGDAMPING = 0.0005
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,5
! Time point at T3
TIMETIME, TIME Sets the time for a load step.,T3T3 = 0.5 ! Ending time of simulation
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,0
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,SUBSTEPS3SUBSTEPS3 = 200,0,0
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,3
! Also create a version of this load step with damping
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,DAMPINGDAMPING = 0.0005
TIMETIME, TIME Sets the time for a load step.,T3T3 = 0.5 ! Ending time of simulation/2! Stop at 0.25 sec instead of 0.5
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,6
! Solve from Load Step files
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,1,3,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(0,0,0),U,Y,UY_CENTER
/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
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,Case_A_No_Damping,bmp ! Click here to view image
FINISHFINISH Exits normally from a processor.

CASE B: With Damping

/SOL/SOLU Enters the solution processor.
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,4,6,1
FINISHFINISH Exits normally from a processor.
/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(0,0,0),U,Y,UY_CENTER
/COLOR/COLOR, Lab, Clab, N1, N2, NINC Specifies the color mapping for various items.,CURVE,BMAG,1! Change the color of Curve 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,Case_B_With_Damping,bmp ! Click here to view image
FINISHFINISH Exits normally from a processor.