Abaqus tutorial
1 like4,326 views
This document provides instructions for creating a finite element model and analysis in ABAQUS. It describes how to generate nodes and elements, define section properties, apply loads and boundary conditions, and output results for a nonlinear beam-column model with spring, damper, and user-defined elements. Key steps include mesh generation, defining nonlinear beam properties, applying gravity loads, performing a displacement-controlled pushover analysis, and conducting a time-history dynamic analysis under base excitation. Variables for the user-defined element subroutine are also outlined.
Abaqus tutorial
- 1. ABAQUS TutorialABAQUS Tutorial by Shiby Shi--Yu XuYu Xu
- 2. Creating the MeshCreating the Mesh **NODENODE, NSET=, NSET=WholeModelWholeModel **10, 0.0, 0.00, 0.0**10, 0.0, 0.00, 0.0 20, 0.0, 0.00, 0.020, 0.0, 0.00, 0.0 30 28 28 30, 0.0,3.048, 0.030, 0.0,3.048, 0.0 ** node#,X,Y,Z** node#,X,Y,Z **NGENNGEN, NSET=, NSET=BENT_COLBENT_COL 28 26 26 24 20, 30, 220, 30, 2 ** 1** 1stst node#, 2node#, 2ndnd node#, increment in the #node#, increment in the # *NSET, NSET=Enforced_node*NSET, NSET=Enforced_node 24 22 24 22 Y,, 3030 **ELEMENTELEMENT,TYPE=,TYPE=B31B31 20, 20, 2220, 20, 22 20 20 X Y Z, ,, , ** ele #, 1** ele #, 1stst node, 2node, 2ndnd nodenode **ELGENELGEN, ELSET=, ELSET=BENTCOL_M_phiBENTCOL_M_phi 20,5, 2, 220,5, 2, 2, , ,, , , ** Master ele #, # of eles to be defined, Increment in node #, Increment in ele #** Master ele #, # of eles to be defined, Increment in node #, Increment in ele #
- 4. Defining the Nonlinear BeamDefining the Nonlinear Beam--Column ElementColumn Element **BEAM GENERAL SECTIONBEAM GENERAL SECTION, DENSITY=2405.0,, DENSITY=2405.0, SECTION=NONLINEARSECTION=NONLINEAR AA TT T hT h Set SECTION=GENERAL (default) for linear beam; Set SECTION=NONLINEAR GENERAL for NL beam. n1GENERALGENERAL, ELSET=, ELSET=BENTCOL_M_phiBENTCOL_M_phi 1.1678, 0.1085, 0.0, 0.1085, 0.21691.1678, 0.1085, 0.0, 0.1085, 0.2169 0.0, 0.0,0.0, 0.0, --1.01.0 **AXIALAXIAL, LINEAR, LINEAR A, I11, I12, I22, J Direction cosines of t1 axis n1 t1,, 2.899E10,2.899E10, **M1M1 0.0000E+00,0.0000E+00, 0.0000E0.0000E--0000 2 7937E+052 7937E+05 8 2008E8 2008E 0505 Axial stiffness t2 Y M 2.7937E+05,2.7937E+05, 8.2008E8.2008E--0505 5.5857E+05,5.5857E+05, 1.6406E1.6406E--0404 ………… **M2M2 X Z φ ………… **TORQUETORQUE, LINEAR, LINEAR 2.244E9,2.244E9, **TRANSVERSE SHEAR STIFFNESSTRANSVERSE SHEAR STIFFNESS Torsion stiffness φ **TRANSVERSE SHEAR STIFFNESSTRANSVERSE SHEAR STIFFNESS 1.075E10, 1.075E101.075E10, 1.075E10 Shear stiffness K23, K13
- 5. Spring and Dashpot ElementsSpring and Dashpot Elements 30 28 **ELEMENTELEMENT, TYPE=, TYPE=SPRING2SPRING2, ELSET=S2, ELSET=S2 201, 10, 20201, 10, 20 **SPRINGSPRING ELSET=S2ELSET=S2 28 26 28 26 130 SPRINGSPRING, ELSET S2, ELSET S2 2, 22, 2 1.E14,1.E14, 24 24 22 128 126 128 126 230 **ELEMENTELEMENT, TYPE=, TYPE=DASHPOT2DASHPOT2, ELSET=C1, ELSET=C1 300, 10, 20300, 10, 20 301 110 120301 110 120 22 20 10 20 126 124 124 122 230 228 228 226301, 110, 120301, 110, 120 302, 210, 220302, 210, 220 303, 310, 320303, 310, 320 **DASHPOTDASHPOT ELSET=C1ELSET=C1 10 122 120 122 120 226 224 224 DASHPOTDASHPOT, ELSET=C1, ELSET=C1 1, 11, 1 1.5589E7,1.5589E7, 110 222 220 222 220 220 210
- 6. UserUser--DefinedDefined Element (UEL)Element (UEL) **USER ELEMENTUSER ELEMENT, NODES=2, TYPE=, NODES=2, TYPE=U1U1,, PROPERTIESPROPERTIES=3,=3, I PROPERTIESI PROPERTIES=1,=1, COORDINATES=3, VARIABLES=14COORDINATES=3, VARIABLES=14 1,2,31,2,3 **ELEMENTELEMENT, TYPE=U1, ELSET=NLSPRING, TYPE=U1, ELSET=NLSPRING 102, 1, 2102, 1, 2 **UEL PROPERTYUEL PROPERTY, ELSET=NLSPRING, ELSET=NLSPRING 1 84E+9 1 84E+8 1 79676E+7 11 84E+9 1 84E+8 1 79676E+7 1 k1 k2QD 1.84E+9, 1.84E+8, 1.79676E+7, 11.84E+9, 1.84E+8, 1.79676E+7, 1 PROPERTIESPROPERTIES I PROPERTIESI PROPERTIES NODES=2: number of nodes connected to the UELNODES=2: number of nodes connected to the UEL TYPE=TYPE=U1U1: label of the UEL: label of the UEL PROPERTIESPROPERTIES=3: number of floating point number parameters required by the UEL=3: number of floating point number parameters required by the UEL I PROPERTIESI PROPERTIES=1: number of integer number parameters required by the UEL=1: number of integer number parameters required by the UEL COORDINATES=3: maximum number of active DOFs at each nodeCOORDINATES=3: maximum number of active DOFs at each node VARIABLES=14: number of solutionVARIABLES=14: number of solution--dependent state variables used by the UELdependent state variables used by the UEL
- 7. MiscellaneousMiscellaneous **DAMPINGDAMPING, ALPHA=0.6641, BETA=0.001884, ALPHA=0.6641, BETA=0.001884 **ELEMENTELEMENT, TYPE=, TYPE=MASSMASS, ELSET=M1, ELSET=M1 801 30801 30801, 30801, 30 **MASSMASS, ELSET=M1, ELSET=M1 389571.389571. **AMPLITUDEAMPLITUDE NAME=linear pushover DEFINITION=TABULARNAME=linear pushover DEFINITION=TABULARAMPLITUDEAMPLITUDE, NAME linear_pushover, DEFINITION TABULAR, NAME linear_pushover, DEFINITION TABULAR 0.00, 0.00, 10.00, 1.00,0.00, 0.00, 10.00, 1.00, Amp. 1 t 10 **AMPLITUDEAMPLITUDE, NAME=accel_ctrl, INPUT=ElCentro033x_ABAQUS.txt, NAME=accel_ctrl, INPUT=ElCentro033x_ABAQUS.txt **BASELINE CORRECTIONBASELINE CORRECTION
- 8. Data OutputData Output S TS T S T C lTS T C lT 30 28 **NSETNSET, NSET=ColTOP, NSET=ColTOP 30,30, **NSETNSET, NSET=ColBOTTOM, NSET=ColBOTTOM 20,20, 28 26 28 26 ,, **NSETNSET, NSET=FNDNbase, NSET=FNDNbase 20,20, **NSETNSET, NSET=Nopt1, NSET=Nopt1 C lTOP C lBOTTOM FNDNbC lTOP C lBOTTOM FNDNb 24 24 22 ColTOP, ColBOTTOM, FNDNbaseColTOP, ColBOTTOM, FNDNbase **ELSETELSET, ELSET=ELopt1, ELSET=ELopt1 20,20, **OUTPUTOUTPUT, FIELD, FREQ=1, FIELD, FREQ=1 22 20 20 **NODE FILENODE FILE, FREQUENCY=1, NSET =Nopt1, FREQUENCY=1, NSET =Nopt1 UU **NODE FILENODE FILE, FREQUENCY=1, NSET =FNDNbase, FREQUENCY=1, NSET =FNDNbase RFRF Output to .fil file Must compose Fortran subroutine to read. **OUTPUTOUTPUT, HISTORY, FREQUENCY=1, HISTORY, FREQUENCY=1 **NODE OUTPUTNODE OUTPUT, NSET=Nopt1, NSET=Nopt1 U1, U2, U3U1, U2, U3 **NODE OUTPUTNODE OUTPUT, NSET=FNDNbase, NSET=FNDNbase Output to .obd file RF1, RF2, RF3, RM1, RM2, RM3RF1, RF2, RF3, RM1, RM2, RM3 **ELEMENT OUTPUTELEMENT OUTPUT, ELSET=ELopt1, ELSET=ELopt1 SF1, SF2, SF3, SM1, SM2, SM3SF1, SF2, SF3, SM1, SM2, SM3 p Can be post-processed in CAE.
- 9. Applying Gravity LoadApplying Gravity Load 30 28 **NSETNSET, NSET=ColTOP, NSET=ColTOP 30,30, **NSETNSET NSET=FNDNbaseNSET=FNDNbase 28 26 28 26 NSETNSET, NSET FNDNbase, NSET FNDNbase 20,20, **STEPSTEP 24 24 22 STEPSTEP STEP 0STEP 0 -- GRAVITY LOADGRAVITY LOAD **STATICSTATIC 1 0E1 0E--10 10 010 10 0 **Applying gravity load in 10 secs. Initial time increment duration of step 22 20 20 1.0E1.0E--10, 10.0,10, 10.0, **BOUNDARYBOUNDARY COLTOP, 3, 5COLTOP, 3, 5 FNDNbase 1 6FNDNbase 1 6 Initial time increment, duration of step Fix DOFs 3 ~ 5 at node set “COLTOP” Fix DOFs 1 ~ 6 at node set “FNDNbase” X Y FNDNbase , 1, 6FNDNbase , 1, 6 **DLOADDLOAD , GRAV, 9.81, 0,, GRAV, 9.81, 0, --1,01,0 **END STEPEND STEP Fix DOFs 1 6 at node set FNDNbase Apply gravity load in -Y direction X Z **END STEPEND STEP
- 10. DisplacementDisplacement--Controlled PushoverControlled Pushover A **NSETNSET, NSET=Enforced_node, NSET=Enforced_node 30 28 28 Amp. 1 3030 **STEPSTEP, INC=8000, INC=8000 28 26 26 24Maximum # of increments to run. t 10 STEP 1STEP 1 -- Static PushoverStatic Pushover **STATICSTATIC, DIRECT, DIRECT 0.05, 10.0,0.05, 10.0, 24 22 22Constant time increment Initial time increment, duration of step ** 0.05, 10.0, 1.E** 0.05, 10.0, 1.E--15, 0.0515, 0.05 **BOUNDARYBOUNDARY,OP=NEW,OP=NEW COLTOP, 3, 5COLTOP, 3, 5 20 20Same, same, min increment, max increment Fix DOFs 3 ~ 5 at node set “COLTOP” FNDNbase , 1, 6FNDNbase , 1, 6 **BOUNDARYBOUNDARY,OP=NEW, TYPE=,OP=NEW, TYPE=DISPLACEMENTDISPLACEMENT, AMPLITUDE=linear_pushover, AMPLITUDE=linear_pushover Enforced_node, 1, 1, 0.4Enforced_node, 1, 1, 0.4 Fix DOFs 1 ~ 6 at node set “FNDNbase” Apply displacement in DOF 1 with SF=0.4 **END STEPEND STEP
- 11. TimeTime--History AnalysisHistory Analysis **NSETNSET NSET=CNSET=CololTOPTOP 30 28**NSETNSET, NSET=C, NSET=CololTOPTOP 30,30, **NSETNSET, NSET=ColBOTTOM, NSET=ColBOTTOM 2020 28 26 28 26 20,20, **NSETNSET, NSET=FNDNbase, NSET=FNDNbase 110,0, 24 24 22 **STEPSTEP, INC=10000, INC=10000 STEP 1STEP 1--11 -- EARTHQUAKE EXCITATION with gravity loadEARTHQUAKE EXCITATION with gravity load **DYNAMICDYNAMIC, DIRECT, DIRECT 22 20 10 20 Maximum # of increments to run. Constant time increment,, 0.004, 25.0,0.004, 25.0, **BOUNDARYBOUNDARY,OP=NEW,OP=NEW COL, 3, 5COL, 3, 5 10 Initial time increment, duration of step Fix DOFs 3 ~ 5 at node set “COL”, ,, , ColTOP, 3, 5ColTOP, 3, 5 FNDNbase, 2, 6FNDNbase, 2, 6 **BOUNDARYBOUNDARY,OP=NEW, TYPE=,OP=NEW, TYPE=ACCELERATIONACCELERATION, AMPLITUDE=accel ctrl, AMPLITUDE=accel ctrl Fix DOFs 3 ~ 5 at node set “COLTOP” Fix DOFs 2 ~ 6 at node set “FNDNbase” , ,, , , _, _ FNDNbase, 1, 1, 9.81FNDNbase, 1, 1, 9.81 **END STEPEND STEP Apply acceleration in DOF 1
- 12. Running the AnalysisRunning the Analysis Change to ABAQUS working directory firstChange to ABAQUS working directory firstChange to ABAQUS working directory first.Change to ABAQUS working directory first. Under the command prompt, type:Under the command prompt, type: The string after theThe string after the ““jobjob”” parameter is theparameter is the script filenamescript filename (w/o extension).(w/o extension). The string after theThe string after the ““useruser”” parameter is the userparameter is the user--defined subroutine (optional).defined subroutine (optional).
- 13. Variables Need to Be DefinedVariables Need to Be Defined RHSRHS ( ) ( ) ( ) extMu x Cu x ku x f RHS+ + − = RHSRHS Contribution of the UEL to the RightContribution of the UEL to the Right--HandHand--Side of the equation of motion.Side of the equation of motion. Size of the RHS array is NRHS.Size of the RHS array is NRHS. AMATRXAMATRX Contribution of the UEL to the Jacobian of the system.Contribution of the UEL to the Jacobian of the system. SVARSSVARS UserUser--defined solutiondefined solution--dependent state variables associated with the UEL.dependent state variables associated with the UEL. Size of the SVARS array is NSVARS.Size of the SVARS array is NSVARS. ENERGYENERGY Energy quantities associated with the UEL.Energy quantities associated with the UEL. ENERGY(1): Kinetic energyENERGY(1): Kinetic energyENERGY(1): Kinetic energy.ENERGY(1): Kinetic energy. ENERGY(2): Elastic strain energy.ENERGY(2): Elastic strain energy. ENERGY(4): Plastic dissipation.ENERGY(4): Plastic dissipation. ENERGY(5): Viscous dissipation.ENERGY(5): Viscous dissipation.( ) p( ) p ENERGY(8): Incremental work done by loads applied within the user element.ENERGY(8): Incremental work done by loads applied within the user element.
- 14. PassedPassed--in Variables (1)in Variables (1) NDOFELNDOFEL Number of DOFs associated with the UEL.Number of DOFs associated with the UEL. NNODENNODE Number of nodes in the UEL.Number of nodes in the UEL. PROPSPROPS UserUser--defined floating point number parameters. Size of the PROPS array is NPROPS.defined floating point number parameters. Size of the PROPS array is NPROPS. JPROPSJPROPS UserUser--defined integer parameters. Size of the JPROPS array is NJPROPS.defined integer parameters. Size of the JPROPS array is NJPROPS. COORDSCOORDS C di t f th d f th UEL ( d f d fi ti )C di t f th d f th UEL ( d f d fi ti )Coordinates of the nodes of the UEL (undeformed configuration).Coordinates of the nodes of the UEL (undeformed configuration). U, DU, V, AU, DU, V, A Current estimates of the solution variables at the end of the increment.Current estimates of the solution variables at the end of the increment. U: Total displacementU: Total displacementU: Total displacement.U: Total displacement. DU: Incremental displacement.DU: Incremental displacement. V: Velocities.V: Velocities. A: Accelerations.A: Accelerations.
- 15. PassedPassed--in Variables (2)in Variables (2) JTYPEJTYPE UserUser--assigned integer number for a specific type of UEL.assigned integer number for a specific type of UEL. (i.e., *ELEMENT, TYPE=U1(i.e., *ELEMENT, TYPE=U1)))) TIMETIME(1)(1) Current step time.Current step time. TIMETIME(2)(2) Current total time.Current total time. DTIMEDTIME Time increment.Time increment. KSTEPKSTEPKSTEPKSTEP Current step number.Current step number. KINCKINC Current increment numberCurrent increment numberCurrent increment number.Current increment number. JELEMJELEM UserUser--assigned element number.assigned element number.
- 16. PassedPassed--in Variables (3)in Variables (3) PARAMSPARAMSPARAMSPARAMS Constants used in the HilberConstants used in the Hilber--HughesHughes--Taylor time integration scheme (i.e., The αTaylor time integration scheme (i.e., The α--Method).Method). PARAMS(1)=α ; PARAMS(2)=β ; PARAMS(3)=γPARAMS(1)=α ; PARAMS(2)=β ; PARAMS(3)=γ LFLAGSLFLAGSLFLAGSLFLAGS An array defining the tasks that the UEL should do when being called.An array defining the tasks that the UEL should do when being called. LFLAGS(1): Defines the procedure type.LFLAGS(1): Defines the procedure type. 1=Static, automatic incrementation1=Static, automatic incrementation,, 2=Static, direct incrementation2=Static, direct incrementation 11=Implicit dynamic, half11=Implicit dynamic, half--step residual tolerance givenstep residual tolerance given 12=Implicit dynamic, fixed time increments12=Implicit dynamic, fixed time increments 41=Eigenvalue frequency extraction41=Eigenvalue frequency extractiong q yg q y LFLAGS(3)=1: Normal implicit time incrementation procedure. Define the RHS and AMATRX arrays.LFLAGS(3)=1: Normal implicit time incrementation procedure. Define the RHS and AMATRX arrays. LFLAGS(3)=2: Define the current stiffness matrix only.LFLAGS(3)=2: Define the current stiffness matrix only. LFLAGS(3)=3: Define the current damping matrix only.LFLAGS(3)=3: Define the current damping matrix only. LFLAGS(3)=4: Define the current mass matrix only.LFLAGS(3)=4: Define the current mass matrix only.( ) y( ) y LFLAGS(3)=5: Define the current residual vector only.LFLAGS(3)=5: Define the current residual vector only. LFLAGS(3)=6: Define the current mass matrix and residual vector.LFLAGS(3)=6: Define the current mass matrix and residual vector. LFLAGS(3)=100: Define perturbation quantities for output.LFLAGS(3)=100: Define perturbation quantities for output. LFLAGS(4)=0: General step.LFLAGS(4)=0: General step.( ) p( ) p LFLAGS(4)=1: Linear perturbation step.LFLAGS(4)=1: Linear perturbation step.
- 17. Compiling the UELCompiling the UEL UseUse ““makemake”” command:command:
- 18. ExampleExample *USER ELEMENT, NODES=2, TYPE=U1, PROPERTIES=4, COORDINATES=3, VARIABLES=12 1, 2, 3 T T 1*ELEMENT, TYPE=U1 101, 101, 102 *ELGEN, ELSET=UTRUSS 101, 5 *UEL PROPERTY, ELSET=UTRUSS 0.002, 2.1E11, 0.3, 7200. k k2QD k1