!!!FEA NL3: Large deformation of a built-in plate

! Created by Gergely Gyebroszki, 2018.05.02.
! All units in SI

!!Parameters

E=120E6
NU=0.3
A=40
B=30
T=0.4
F=40E3

!!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,SHELL181
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 = 120E6
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
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.,TT = 0.4,1,0.0,3
SECOFFSETSECOFFSET, Location, OFFSET1, OFFSET2, CG-Y, CG-Z, SH-Y, SH-Z Defines the section offset for cross sections.,MID

!!MODELLING

KK, NPT, X, Y, Z Defines a keypoint.,! Auto-numbered at origin
KK, NPT, X, Y, Z Defines a keypoint.,,40
KK, NPT, X, Y, Z Defines a keypoint.,,40,30
KK, NPT, X, Y, Z Defines a keypoint.,,,30
AA, P1, P2, P3, P4, P5, P6, P7, P8, P9, P10, P11, P12, P13, P14, P15, P16, P17, P18 Defines an area by connecting keypoints.,1,2,3,4
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,ALL,,,20
AMESHAMESH, NA1, NA2, NINC Generates nodes and area elements within areas.,ALL

! Built-in end
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,0,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: Static

/SOL/SOLU Enters the solution processor.
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,0! 0: Static
NLGEOMNLGEOM, Key Includes large-deflection effects in a static or full transient analysis.,1! Large displacement
OUTRESOUTRES, Item, Freq, Cname, -- , NSVAR, DSUBres Controls the solution data written to the database.,ALL,1! Write every substep
NEQITNEQIT, NEQIT, FORCEkey Specifies the maximum number of equilibrium iterations for nonlinear analyses.,150! Number of eq. iterations

! Force load
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,S,LOC,X,40,40! Select nodes with X=40
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,R,LOC,Y,0,0! Narrow selection to the single node with Y=0
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,ALL,FZ,-FF = 40E3/5 ! Apply FZ
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
! Another way to select a node is by using the NODE(x,y,z) function
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,NODE(40,0,0),FZ,-FF = 40E3*2/5
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,2
! Load step 3
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,NODE(40,0,0),FZ,-FF = 40E3*3/5
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,3
! Load step 4
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,NODE(40,0,0),FZ,-FF = 40E3*4/5
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,4
! Load step 5
FF, NODE, Lab, VALUE, VALUE2, NEND, NINC Specifies force loads at nodes.,NODE(40,0,0),FZ,-FF = 40E3
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,5
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,1,5,1
FINISHFINISH Exits normally from a processor.

!!Postprocessing

/POST1/POST1 Enters the database results postprocessor.
NA=NODE(40,0,0)
NB=NODE(40,30,0)
*DO*DO, Par, IVAL, FVAL, INC Defines the beginning of a do-loop.,I,1,5 ! Start a cycle, looping on the five load steps
  SETSET, Lstep, Sbstep, Fact, KIMG, TIME, ANGLE, NSET, ORDER Defines the data set to be read from the results file.,I
  /OUTPUT/OUTPUT, Fname, Ext, --, Loc Redirects text output to a file or to the screen.,ADISP,TXT,,APPEND! Create an output file, ADISP
  *VWRITE*VWRITE, Par1, Par2, Par3, Par4, Par5, Par6, Par7, Par8, Par9, Par10, Par11, Par12, Par13, Par14, Par15, Par16, Par17, Par18, Par19 Writes data to a file in a formatted sequence.,UX(NA),UY(NA),UZ(NA)! Write the displacements at A into ADISP
  (E16.8,5X,E16.8,5X,E16.8,5X)! Specify the number format
  /OUTPUT/OUTPUT, Fname, Ext, --, Loc Redirects text output to a file or to the screen.! Finalize output file
  /OUTPUT/OUTPUT, Fname, Ext, --, Loc Redirects text output to a file or to the screen.,BDISP,TXT,,APPEND! Create an output file, BDISP
  *VWRITE*VWRITE, Par1, Par2, Par3, Par4, Par5, Par6, Par7, Par8, Par9, Par10, Par11, Par12, Par13, Par14, Par15, Par16, Par17, Par18, Par19 Writes data to a file in a formatted sequence.,UX(NB),UY(NB),UZ(NB)! Write the displacements at B into BDISP
  (E16.8,5X,E16.8,5X,E16.8,5X)! Specify the number format
  /OUTPUT/OUTPUT, Fname, Ext, --, Loc Redirects text output to a file or to the screen.! Finalize output file
*ENDDO*ENDDO Ends a do-loop and starts the looping action.! End of cycle

! Click here to view ADISP.TXT
! Click here to view ADISP.TXT
! Click here to view Displacement plot of Node A
! Click here to view Displacement plot of Node B
PLNSOLPLNSOL, Item, Comp, KUND, Fact, FileID Displays results as continuous contours.,U,Z,0,1.0! Plot nodal solution UZ Click to view plot