!!!FEA NL4: IMPACT OF A SQUARE BLOCK WITH AN ELASTIC 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
/ESHAPE/ESHAPE, SCALE, KEY Displays elements with shapes determined from the real constants or section definition.,1.0
! Define constant parameters
! Units: [N],[Nm],[N/m^2=Pa],[N/m],[kg=N/(m/s^2)]
ELAST=70E9
NU=0.3
RHO=2700
H=0.025
V=0.025
L=2.5
BEAMH=28.0624e-3
BEAMW=11.4031e-3
delT=0.005
ELEMSIZE=0.15

!!SETUP: element types, real constants, materials, sections

! plane stress with thickness
KEYOPTKEYOPT, ITYPE, KNUM, VALUE Sets element key options.,2,3,3
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,3,CONTA175
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,4,TARGE169
! Define crossection properties
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.,BEAMWBEAMW = 11.4031e-3,BEAMHBEAMH = 28.0624e-3
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,1,3.2e-4,2.1e-8,0.025
! Thickness of the plane element
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,2,HH = 0.025
! Penalty stiffness factor of the contact element
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,3,2000
! Define 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 = 70E9
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.,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 = 2700

!!Modeling

!!Meshing

! Mesh lines
TYPETYPE, ITYPE Sets the element type attribute pointer.,1
MATMAT, MAT Sets the element material attribute pointer.,1
SECNSECNUM, SECID Sets the element section attribute pointer.,1
LMESHLMESH, NL1, NL2, NINC Generates nodes and line elements along lines.,1
! Mesh area with plane elements using mapped algorithm
TYPETYPE, ITYPE Sets the element type attribute pointer.,2
REALREAL, NSET Sets the element real constant set attribute pointer.,2
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.,1
! Click here to see current state.
FINISHFINISH Exits normally from a processor.

!!Transient analysis definition

! Fix the ends of the beam
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,NODE(0,0,0),ALL,0
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,NODE(L,0,0),ALL,0
! Fix the nodes of rectangle
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Y,2.5,2.65
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
! Define gravitational acceleration
ACELACEL, ACEL_X, ACEL_Y, ACEL_Z Specifies the linear acceleration of the global Cartesian reference frame for the analysis.,0,9.81,0
NLGEOMNLGEOM, Key Includes large-deflection effects in a static or full transient analysis.,1
DELTIMDELTIM, DTIME, DTMIN, DTMAX, Carry Specifies the time step sizes to be used for the current load step.,0.0001
ALPHADALPHAD, VALUE Defines the mass matrix multiplier for damping.,1.5
BETADBETAD, VALUE Defines the stiffness matrix multiplier for damping.,6e-4
KBCKBC, KEY Specifies ramped or stepped loading within a load step.,1.0
TIMETIME, TIME Sets the time for a load step.,0.0002
TIMINTTIMINT, Key, Lab Turns on transient effects.,0
CNVTOLCNVTOL, Lab, VALUE, TOLER, NORM, MINREF Sets convergence values for nonlinear analyses.,F,,0.0001,2
SOLVESOLVE, Action Starts a solution.
! Delete constraints on the nodes of the rectangle
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,Y,2.5,2.65
DDELEDDELE, NODE, Lab, NEND, NINC, Rkey Deletes degree-of-freedom constraints.,ALL,ALL
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
FINISHFINISH Exits normally from a processor.
/PREP7/PREP7 Enters the model creation preprocessor.
TYPETYPE, ITYPE Sets the element type attribute pointer.,3
TSHAPTSHAP, Shape Defines simple 2-D and 3-D geometric surfaces for target segment elements.,LINE
TYPETYPE, ITYPE Sets the element type attribute pointer.,4
TSHAPTSHAP, Shape Defines simple 2-D and 3-D geometric surfaces for target segment elements.,LINE
! Create conta175 elements on the sides of the rectangle
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,LINE,,2,5,1
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
TYPETYPE, ITYPE Sets the element type attribute pointer.,3
REALREAL, NSET Sets the element real constant set attribute pointer.,3
ESURFESURF, XNODE, Tlab, Shape Generates elements overlaid on the free faces of selected nodes.
! Create targe169 elements on the beam
LSELLSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Selects a subset of lines.,S,LINE,,1
NSLLNSLL, Type, NKEY Selects those nodes associated with the selected lines.,S,1
TYPETYPE, ITYPE Sets the element type attribute pointer.,4
REALREAL, NSET Sets the element real constant set attribute pointer.,3
ESURFESURF, XNODE, Tlab, Shape Generates elements overlaid on the free faces of selected nodes.
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.,ALL
GPLOTGPLOT Controls general plotting.
FINISHFINISH Exits normally from a processor.
/SOL/SOLU Enters the solution processor.
DELTIMDELTIM, DTIME, DTMIN, DTMAX, Carry Specifies the time step sizes to be used for the current load step.,0.02,0.0002,0.02
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
KBCKBC, KEY Specifies ramped or stepped loading within a load step.,1.0
PREDPRED, Sskey, --, Lskey Activates a predictor in a nonlinear analysis.,ON,,ON
TIMETIME, TIME Sets the time for a load step.,3.5
TIMINTTIMINT, Key, Lab Turns on transient effects.,1
CNVTOLCNVTOL, Lab, VALUE, TOLER, NORM, MINREF Sets convergence values for nonlinear analyses.,,-1,1
SOLVESOLVE, Action Starts a solution.
FINISHFINISH Exits normally from a processor.

!!Time History Postprocessing

! Displacements
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,2,NODE(L/2,0,0),U,Y,MID
/AXLAB/AXLAB, Axis, Lab Labels the X and Y axes on graph displays.,X,'t [s]'
/AXLAB/AXLAB, Axis, Lab Labels the X and Y axes on graph displays.,Y,'uy [m]'
/XRANGE/XRANGE, XMIN, XMAX Specifies a linear abscissa (X) scale range.,0,3.5
PLVARPLVAR, NVAR1, NVAR2, NVAR3, NVAR4, NVAR5, NVAR6, NVAR7, NVAR8, NVAR9, NVAR10 Displays up to ten variables in the form of a graph.,2! Click here to see the y displacement of the middle of the beam.
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,3,NODE(1.4,2.5,0),U,Y,BL
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,4,NODE(1.55,2.65,0),U,Y,BR
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,5,NODE(1.4,2.65,0),U,Y,TL
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,6,NODE(1.55,2.65,0),U,Y,TR
/YRANGE/YRANGE, YMIN, YMAX, NUM Specifies a linear ordinate (Y) scale range.,-2.8,0
PLVARPLVAR, NVAR1, NVAR2, NVAR3, NVAR4, NVAR5, NVAR6, NVAR7, NVAR8, NVAR9, NVAR10 Displays up to ten variables in the form of a graph.,3,4,5,6! Click here to see the y displacement of the corner of the rectangles.
! Reaction forces and moments
RFORCERFORCE, NVAR, NODE, Item, Comp, Name Specifies the total reaction force data to be stored.,7,NODE(0,0,0),F,Y,FY_0
RFORCERFORCE, NVAR, NODE, Item, Comp, Name Specifies the total reaction force data to be stored.,8,NODE(0,0,0),M,Z,MZ_0
RFORCERFORCE, NVAR, NODE, Item, Comp, Name Specifies the total reaction force data to be stored.,9,NODE(L,0,0),F,Y,FY_L
RFORCERFORCE, NVAR, NODE, Item, Comp, Name Specifies the total reaction force data to be stored.,10,NODE(L,0,0),M,Z,MZ_L
/YRANGE/YRANGE, YMIN, YMAX, NUM Specifies a linear ordinate (Y) scale range.,DEFAULT
/AXLAB/AXLAB, Axis, Lab Labels the X and Y axes on graph displays.,Y,'Fy [N]'
PLVARPLVAR, NVAR1, NVAR2, NVAR3, NVAR4, NVAR5, NVAR6, NVAR7, NVAR8, NVAR9, NVAR10 Displays up to ten variables in the form of a graph.,7,9! Click here to see the y component of reaction forces at the fixed ends.
/AXLAB/AXLAB, Axis, Lab Labels the X and Y axes on graph displays.,Y,'Mz [Nm]'
PLVARPLVAR, NVAR1, NVAR2, NVAR3, NVAR4, NVAR5, NVAR6, NVAR7, NVAR8, NVAR9, NVAR10 Displays up to ten variables in the form of a graph.,8,10! Click here to see the z component of reaction moments at the fixed ends.
FINISHFINISH Exits normally from a processor.

!!Postprocessing

! Uncomment the following do block to save images
!*DO,TIME,0,3.5,delT
!SET,,,1,,TIME
!PLDISP,0
!/image,save,strcat('disp_',chrval(200*TIME)),bmp
!*ENDDO
! Click here to see the animation of the whole motion.

! This webpage was generated with APDL highlighter (version: 0.1-beta) on: 09.05.2018 11:20:14, from APDL script: script.txt