FEA D5: Dynamic analysis of a 4-bar linkage

! Created by Gergely Gyebroszki, 2018.04.19.

Units in mm/MPa

Parameters

AREA=0.063
IZZ=0.00525
IYY=2.084e-5
TKZ=0.063
K1=10e6
K2=10e6
K3=10e6
K4=1.5528e-4
E=10e6
NU=0.3
RHO=2.588e-4
DAMPING1=0.000139
DAMPING2=0.000278

SETUP

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,BEAM4
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,2,COMBIN7
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,1,AREAAREA = 0.063,IZZIZZ = 0.00525,IYYIYY = 2.084e-5,TKZTKZ = 0.063
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,2,K1K1 = 10e6,K2K2 = 10e6,K3K3 = 10e6,K4K4 = 1.5528e-4
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 = 10e6
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 = 2.588e-4
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DMPR,1,,DAMPING1DAMPING1 = 0.000139
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EX,2,,EE = 10e6
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,2,,NUNU = 0.3
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DENS,2,,RHORHO = 2.588e-4
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DMPR,2,,DAMPING2DAMPING2 = 0.000278

MODELLING

NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,0
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,2.5
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,2.5

! New local coord system
LOCALLOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2 Defines a local coordinate system by a location and orientation.,11,0,2.5,0,0,66.024,,,1,1
CSYSCSYS, KCN Activates a previously defined coordinate system.,11
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,2.5
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,5.5
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,8.5
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,11.0
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,11.0
CSYSCSYS, KCN Activates a previously defined coordinate system.,0
*AFUN*AFUN, Lab Specifies units for angular functions in parameter expressions.,DEG
! Second local CS
LOCALLOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2 Defines a local coordinate system by a location and orientation.,12,0,2.5+11*cos(66.024),11*sin(66.024),0,-73.21,,,1,1
CSYSCSYS, KCN Activates a previously defined coordinate system.,12
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,2.5
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,5.25
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,8.0
NN, NODE, X, Y, Z, THXY, THYZ, THZX Defines a node.,,10.5
! Create BEAM4 elements 1-5
MATMAT, MAT Sets the element material attribute pointer.,1
TYPETYPE, ITYPE Sets the element type attribute pointer.,1
REALREAL, NSET Sets the element real constant set attribute pointer.,1
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,1,2
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,3,4
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,4,5
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,5,6
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,6,7
MATMAT, MAT Sets the element material attribute pointer.,2
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,8,9
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,9,10
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,10,11
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,11,12
! Create Joint elements
TYPETYPE, ITYPE Sets the element type attribute pointer.,2
REALREAL, NSET Sets the element real constant set attribute pointer.,2
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,2,3
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,7,8
! Constraints
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,1,UX
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,1,UY
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,12,UX
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,12,UY
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,ALL,UZ
! Gravity
ACELACEL, ACEL_X, ACEL_Y, ACEL_Z Specifies the linear acceleration of the global Cartesian reference frame for the analysis.,209.4,324.7,0
FINISHFINISH Exits normally from a processor.

SOLUTION

/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
NLGEOMNLGEOM, Key Includes large-deflection effects in a static or full transient analysis.,1 ! Large displacement static
AUTOTSAUTOTS, Key Specifies whether to use automatic time stepping or load stepping.,0 ! Auto stepping: OFF
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,1 ! Number of substeps
NEQITNEQIT, NEQIT, FORCEkey Specifies the maximum number of equilibrium iterations for nonlinear analyses.,50 ! Number of Equilibrium iterations
TIMETIME, TIME Sets the time for a load step.,0.001! Time at the end of step
TIMINTTIMINT, Key, Lab Turns on transient effects.,0 ! Transient effect: OFF
OUTRESOUTRES, Item, Freq, Cname, -- , NSVAR, DSUBres Controls the solution data written to the database.,ALL,1 ! Write every substep
! Rotz = -0.041888 on left linkage (-1.2°)
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,1,ROTZ,-0.041888
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,1
! Load step 2
TIMETIME, TIME Sets the time for a load step.,0.6! Time at the end of step
NSUBSTNSUBST, NSBSTP, NSBMX, NSBMN, Carry Specifies the number of substeps to be taken this load step.,99 ! Number of substeps
TIMINTTIMINT, Key, Lab Turns on transient effects.,1 ! Transient effect: ON
KBCKBC, KEY Specifies ramped or stepped loading within a load step.,0.0 ! Ramped loading
NEQITNEQIT, NEQIT, FORCEkey Specifies the maximum number of equilibrium iterations for nonlinear analyses.,30 ! Number of Equilibrium iterations
! Rotz = -25.13274 on left linkage (-720°)
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,1,ROTZ,-25.13274
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

TIME HISTORY POSTPROCESSING

/POST26/POST26 Enters the time-history results postprocessor.
/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
/XRANGE/XRANGE, XMIN, XMAX Specifies a linear abscissa (X) scale range.,0,0.6

! Displacement of Node 2
NSOLNSOL, NVAR, NODE, Item, Comp, Name, SECTOR Specifies nodal data to be stored from the results file.,2,2,U,Y,UY_NODE2
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,UY_Node2,bmp ! Click here to view image
! Reaction moment on Node 1
RFORCERFORCE, NVAR, NODE, Item, Comp, Name Specifies the total reaction force data to be stored.,3,1,M,Z,MZ_NODE1
/COLOR/COLOR, Lab, Clab, N1, N2, NINC Specifies the color mapping for various items.,CURVE,RED,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,MZ_Node3,bmp ! Click here to view image
! Animate over time
ANTIMEANTIME, NFRAM, DELAY, NCYCL, AUTOCNTRKY, RSLTDAT, MIN, MAX Produces a sequential contour animation over a range of time.,100,0.2,,1,2,0,0.6 ! Click here to view animation