Transient Dynamic Analysis in Ansys
Transient Dynamic Analysis in Ansys
Page 1 of 55
where:
[M] = mass matrix
[C] = damping matrix
[K] = stiffness matrix
= nodal acceleration vector
= nodal velocity vector
{u} = nodal displacement vector
{F(t)} = load vector
At any given time, t, these equations can be thought of as a set of "static" equilibrium equations that
also take into account inertia forces ([M] ) and damping forces ([C] ). The ANSYS program uses
the Newmark time integration method to solve these equations at discrete timepoints. The time
increment between successive timepoints is called the integration time step.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 2 of 55
Graphical User Interface (GUI) to build and solve models no matter what type of analysis you are
doing.
Section 5.8, "Sample Transient Dynamic Analysis (Command or Batch Method)," shows you the
sequence of commands you would issue (either manually or while running ANSYS as a batch job) to
perform an example transient dynamic analysis. Section 5.7, "Sample Transient Dynamic Analysis
(GUI Method)," shows you how to execute the same sample analysis using menu choices from the
ANSYS GUI. (To learn how to use the commands and GUI selections for building models, read the
ANSYS Modeling and Meshing Guide.)
For detailed, alphabetized descriptions of the ANSYS commands, see the ANSYS Commands
Reference.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 3 of 55
Note-If you do not want to include any nonlinearities, you should consider using one of the other
methods because the full method is also the most expensive method of the three.
The advantages of the full method are:
It is easy to use, because you don't have to worry about choosing master degrees of freedom or
mode shapes.
It allows all types of nonlinearities.
It uses full matrices, so no mass matrix approximation is involved.
All displacements and stresses are calculated in a single pass.
It accepts all types of loads: nodal forces, imposed (non-zero) displacements (although not
recommended), and element loads (pressures and temperatures).
It allows effective use of solid-model loads.
The main disadvantage of the full method is that it is more expensive than either of the other methods.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 4 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 5 of 55
If you are interested in wave propagation effects (for example, a bar dropped exactly on its
end), the mesh should be fine enough to resolve the wave. A rule of thumb is to have at least 20
elements per wavelength along the direction of the wave.
New Analysis
Analysis Type:
Transient Dynamic
ANTYPE
Solution Method
Mass Matrix
Formulation
Large Deformation
Effects
Stress Stiffening
Effects
SSTIF
Newton-Raphson
Option
NROPT
Equation Solver
EQSLV
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 6 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 7 of 55
Program-chosen (default)
Full
Modified
Initial Stiffness
The first load step you apply is usually to establish initial conditions. You then specify the loads
and load step options for the subsequent transient load steps. For each load step, you need to
specify both load values and time values, along with other load step options such as whether to
step or ramp the loads, use automatic time stepping, etc. You then write each load step to a file
and solve all load steps together.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 8 of 55
Nonzero initial displacement and nonzero initial velocity-This is similar to the above case,
except that the imposed displacements are actual values instead of "small" values. For example,
if = 1.0 and = 2.5, you would apply a displacement of 1.0 over a time interval of 0.4:
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 9 of 55
...
TIMINT,OFF
! Time integration effects off
D,ALL,UY,1.0
! Initial displacement = 1.0
TIME,.4 ! Initial velocity = 1.0/0.4 = 2.5
LSWRITE ! Write load data to load step file (Jobname.S01)
DDELE,ALL,UY
! Remove imposed displacements
TIMINT,ON
! Time integration effects on
...
Nonzero initial displacement and zero initial velocity-This requires the use of two substeps
[NSUBST,2] with a step change in imposed displacements [KBC,1]. Without the step change
(or with just one substep), the imposed displacements would vary directly with time, leading to
a nonzero initial velocity. The example below shows how to apply = 1.0 and = 0.0:
...
TIMINT,OFF
! Time integration effects off for static solution
D,ALL,UY,1.0
! Initial displacement = 1.0
TIME,.001
! Small time interval
NSUBST,2
! Two substeps
KBC,1
! Stepped loads
LSWRITE ! Write load data to load step file (Jobname.S01)
! Transient solution
TIMINT,ON
! Time-integration effects on for transient solution
TIME,...
! Realistic time interval
DDELE,ALL,UY
! Remove displacement constraints
KBC,0
! Ramped loads (if appropriate)
! Continue with normal transient solution procedures
...
See the ANSYS Commands Reference for a discussion of the ACEL, TIME, NSUBST, KBC,
LSWRITE, DDELE and KBC commands.
Apply Loads for the Transient Loading Portion
Table 5-2 summarizes the loads applicable to a transient dynamic analysis. Except for inertia loads,
you can define loads either on the solid model (keypoints, lines, and areas) or on the finite element
model (nodes and elements). For a general discussion of solid-model loads versus finite-element
loads, see Chapter 2 of the ANSYS Basic Analysis Procedures Guide.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 10 of 55
You can also apply complex boundary conditions by defining a 1-dimensional table (TABLE type
array parameter). See "Applying Loads Using TABLE Type Array Parameters" in Section 2.5.2.3.
Table 5-2 Loads applicable in a transient dynamic analysis
Cmd
Family
Load Type
Category
Constraints D
GUI Path
Main Menu>Solution>-LoadsApply>-Structural- Displacement
SF
Pressure (PRES)
Surface
Loads
Temperature (TEMP)
Fluence (FLUE)
Moisture Content (MSCN)
Body Loads BF
Inertia
Loads
Load Type
Solid
Model Entity
or FE
Displacement
Solid
Keypoints DK
Model
Solid
Lines
Model
DL
Solid
Areas
Model
DA
DDELE
Finite
Elem
Nodes
Apply
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
Delete
List
DLIST
Operate
Apply
Settings
DSCALE DSYM
DCUM
12/22/2015
Page 11 of 55
Solid
Keypoints FK
Model
Finite
Elem
FDELE
Solid
Lines
Model
SFL
Solid
Areas
Model
SFA
Finite
Elem
SF
SFDELE
Force
Pressure
Finite
Elem
Temperature
Fluence
Inertia
Nodes
Nodes
Elements SFE
FLIST
SFLIST
FSCALE FCUM
SFSCALE
SFGRAD
SFCUM
SFGRAD
SFBEAM
SFEDELE SFELIST SFSCALE
SFFUN
SFCUM
Solid
Keypoints BFK
Model
Solid
Lines
Model
BFL
Solid
Areas
Model
BFA
Solid
Volumes BFV
Model
Finite
Elem
Nodes
BFDELE
Finite
Elem
Elements BFE
BF
BFSCALE BFCUM
ACEL
OMEGA
DOMEGA
CGLOC
CGOMGA
DCGOM IRLF
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
BFLIST
12/22/2015
Page 12 of 55
Command
GUI Path
TIMINT
Dynamics Options
Time Integration
Effects
Transient Integration
TINTP
Parameters
Damping
General Options
Time
TIME
Stepped or Ramped
Loads
KBC
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 13 of 55
Automatic Time
Stepping
NSUBST
DELTIM
AUTOTS
Nonlinear Options
Max. No. of
NEQIT
Equilibrium Iterations
Convergence
Tolerances
CNVTOL
Predictor-Corrector
Option
PRED
LNSRCH
Creep Criteria
CRPLIM
Solution Termination
NCNV
Options
OUTPR
Extrapolation of
Results
ERESX
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 14 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 15 of 55
recommend that you turn on automatic time stepping, with upper and lower limits for the
integration time step. These limits, specified using DELTIM or NSUBST, help to limit the
range of variation of the time step; see Section 5.12.2, "Automatic Time Stepping" for more
information. The default is ON when SOLCONTROL is ON. Otherwise, it is OFF.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
LSWRITE
TIME,...
Loads ...
KBC,...
LSWRITE
TIME,...
Loads ...
KBC,...
LSWRITE
Etc.
!
!
!
!
!
!
!
!
!
Page 16 of 55
See the ANSYS Commands Reference for a discussion of the TIME, KBC and LSWRITE
commands.
4. Save the load configuration for this step to the load step file.
Command(s):
LSWRITE
GUI:
Main Menu>Solution>Write LS File
Repeat steps 3. and 4 for each corner of the load-versus-time curve. You may also want to have
an additional load step that extends past the last time point on the curve to capture the response
of the structure after the transient loading.
5. Save a back-up copy of the database to a named file. You can then retrieve your model by reentering the ANSYS program and issuing RESUME.
Command(s):
SAVE
GUI:
Utility Menu>File>Save as
6. Start the transient solution. For additional ways to create and solve multiple load steps (the array
parameter method and the multiple SOLVE method), see Section 3.10 of the ANSYS Basic Analysis
Procedures Guide.
Command(s):
LSSOLVE
GUI:
Main Menu>Solution>-Solve-From LS Files
7. Leave SOLUTION.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 17 of 55
Command(s):
FINISH
GUI:
Close the Solution menu.
5.5.4.1 Postprocessors
You can review these results using either POST26, the time-history postprocessor, or POST1, the
general postprocessor.
POST26 is used to review results at specific points in the model as functions of time.
POST1 is used to review results over the entire model at specific time-points.
Some typical postprocessing operations for a transient dynamic analysis are explained below. For a
complete description of all postprocessing functions, see Section 4.1 of the ANSYS Basic Analysis
Procedures Guide.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 18 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 19 of 55
GUI:
Main Menu>General Postproc>-Read Results-By Time/Freq
Note-If you specify a time for which no results are available, the results that are stored will be a
linear interpolation between the two nearest time points.
3. Perform the necessary POST1 operations. Typical POST1 operations are explained below.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 20 of 55
FSUM
GUI:
Main Menu>General Postproc>Nodal Calcs>Total Force Sum
You can also check the total force and total moment at each selected node. For a body in
equilibrium, the total load is zero at all nodes except where an applied load or reaction load
exists.
Command(s):
NFORCE
GUI:
Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node
The FORCE command (Main Menu>General Postproc>Options for Outp) dictates which
component of the forces is being reviewed:
Total (default)
Static component
Damping component
Inertia component
For a body in equilibrium, the total load (using all FORCE components) is zero at all nodes
except where an applied load or reaction load exists.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 21 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 22 of 55
different shell thicknesses, or other discontinuities meet. To avoid the smearing effect, use selecting
(described in Chapter 7 of the ANSYS Basic Analysis Procedures Guide) to select elements of the
same material, same shell thickness, etc. before issuing PLNSOL. Alternatively, use PowerGraphics
with the AVRES command (Main Menu>General Postproc>Options for Outp) to not average
results across different materials and/or different shell thicknesses.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
----FINISH
Page 23 of 55
! Generate model
See the ANSYS Commands Reference for a discussion of the ANTYPE, TRNOPT, ALPHAD,
BETAD, KBC, TIME, AUTOTS, DELTIM, OUTRES, LSWRITE, LSSOLVE, SOLU, NSOL,
ESOL, RFORCE, PLVAR, PRVAR, PLDISP, PRRSOL, PLNSOL and PRERR commands.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 24 of 55
should consider using this method if you don't want to include nonlinearities (other than simple nodeto-node contact) in the analysis.
The procedure for a reduced transient dynamic analysis consists of five main steps:
1. Build the model.
2. Obtain the reduced solution.
3. Review the results of the reduced solution.
4. Expand the solution (expansion pass).
5. Review the results of the expanded solution.
Of these, the first step is the same as for the full method, except that no nonlinearities are allowed
(other than simple node-to-node contact, which is specified in the form of a gap condition instead of
an element type). Details of the other steps are explained below.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 25 of 55
M
MGEN
TOTAL
MLIST
MDELE
GUI:
Main Menu>Solution>Master DOFs>Define/Copy/Program Selected
Main Menu>Solution>Master DOF>List All
Main Menu>Solution>Master DOFs>Delete
4. Define gap conditions, if any.
Command(s):
GP
GUI:
Main Menu>Solution>Dynamic Gap Cond>Define
You can also list the defined gaps and delete gaps.
Command(s):
GPLIST
GPDELE
GUI:
Main Menu>Solutions>Dynamic Gap Cond>List All
Main Menu>Solutions>Dynamic Gap Cond>Delete
Gap Conditions
Gap conditions can only be defined between two master nodes or between master nodes and
ground, as shown in the following figure.
Figure 5-2 Examples of gap conditions
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 26 of 55
Gap conditions are similar to gap elements and are specified between surfaces that are expected
to contact (impact) each other during the transient. The ANSYS program accounts for the gap
force, which develops when the gap closes, by using an equivalent nodal load vector.
Some guidelines to define gap conditions are presented below:
Use enough gap conditions to obtain a smooth contact stress distribution between the contacting
surfaces.
Define a reasonable gap stiffness. If the stiffness is too low, the contacting surfaces may
overlap too much. If the stiffness is too high, a very small time step will be required during
impact. A general recommendation is to specify a gap stiffness that is one or two orders of
magnitude higher than the adjacent element stiffness. You can estimate the adjacent element
stiffness using AE/L, where A is the contributing area around the gap condition, E is the elastic
modulus of the softer material at the interface, and L is the depth of the first layer of elements at
the interface.
5. Apply initial conditions to the model. The following loading restrictions apply in a reduced
transient dynamic analysis:
Only displacements, forces, and translational accelerations (such as gravity) are valid.
Acceleration loading is not allowed if the model contains any master DOF at any nodes with
rotated nodal coordinate systems.
Forces and non-zero displacements must be applied only at master DOF.
As mentioned for the full method, multiple load steps are usually required to specify the load
history in a transient analysis. The first load step is used to establish initial conditions, and
second and subsequent load steps are used for the transient loading, as explained next.
Establish initial conditions. The only initial condition that may be explicitly established is the
initial displacement ( ); that is, initial velocity and acceleration must be zero ( = 0, = 0).
Displacements cannot be deleted in subsequent load steps, therefore they cannot be used to
specify an initial velocity. In a reduced transient analysis, a static solution is always performed
as the first solution, using the loads given, to determine .
Specify load step options for the first load step. The following options are available for the first
load step.
Table 5-5 Options for the first load step
Option
Dynamics Options
Transient Integration
Parameters
Damping
TINTP
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 27 of 55
General Options
Integration Time Step
DELTIM
OUTPR
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 28 of 55
Use this option to output the displacement solution at the master DOF.
6. Write the first load step to a load step file (Jobname.S01).
Command(s):
LSWRITE
GUI:
Main Menu>Solution>-Solve-Write LS File
7. Specify loads and load step options for the transient loading portion, writing each load step to a
load step file [LSWRITE].
The following load step options are valid for the transient load steps:
General Options
Time (specifies the time at the end of the load step) [TIME]
Stepped [KBC,1] or ramped loads [KBC]
Output Controls
Printed output [OUTPR]
Reduced displacement file [OUTRES]
The only valid label on these commands is NSOL (nodal solution). The default for OUTRES is
to write the solution for every fourth time-point to the reduced displacement file (unless there
are gap conditions defined, in which case the default is to write every solution).
8. Save a back-up copy of the database to a named file.
Command(s):
SAVE
GUI:
Utility Menu>File>Save as
9. Start the transient solution. For additional ways to create and solve multiple load steps (the array
parameter method and the multiple SOLVE method), see Section 3.10.3 of the ANSYS Basic Analysis
Procedures Guide.
Command(s):
LSSOLVE
GUI:
Main Menu>Solution>-Solve-From LS Files
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 29 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 30 of 55
Expansion Pass
On/Off
No. of Solutions to be
Main Menu>Solution>-Load Step OptsNUMEXP
Expanded
ExpansionPass>Range of Solu's
Single Solution to
Expand
EXPSOL
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 31 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 32 of 55
You can review these results using POST1. (If you expanded solutions at several time-points, you can
also use POST26 to obtain graphs of stress versus time, strain versus time, etc.) The procedure to use
POST1 (or POST26) is the same as described for the full method.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 33 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 34 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 35 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 36 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 37 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 38 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 39 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 40 of 55
NUMVAR,0
FILE,file,rdsp
NSOL,2,2,U,Y,NSOL
! Define the variables
PLVAR,2 ! Graph the variables
PRVAR,2 ! List the variables
FINISH
/POST1
SET,FIRST
PLDISP,1
FINISH
! Read in results
! Display deformed and undeformed shape
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 41 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 42 of 55
Specify displacement constraints, if any. These constraints will be ignored if they are specified
in the mode superposition transient solution instead of in the modal solution.
If you need to apply element loads (pressures, temperatures, accelerations, etc.) in the transient
dynamic analysis, you must specify them in the modal analysis. The loads are ignored for the
modal solution, but a load vector will be calculated and written to the mode shape file
(Jobname.MODE). You can then use this load vector for the transient solution.
The modes need not be expanded for the mode superposition solution. (If you want to review
mode shapes from a reduced modal solution, however, you must expand the mode shapes.)
The model data (e.g., nodal rotations) should not be changed between the modal and transient
analyses.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 43 of 55
Command(s):
GP
GUI:
Main Menu>Solution>Dynamic Gap Cond>Define
4. Apply loads to the model. The following loading restrictions apply in a mode superposition
transient dynamic analysis:
Only forces, translational accelerations, and a load vector created in the modal analysis are
valid. Imposed non-zero displacements are ignored. Use the LVSCALE command (Main
Menu>Solution>-Loads-Apply> Load Vector>For Mode Super) to apply the load vector
from the modal solution.
If mode shapes from a reduced modal solution are being used, forces may be applied only at
master DOF.
Multiple load steps are usually required to specify the load history in a transient analysis. The
first load step is used to establish initial conditions, and second and subsequent load steps are
used for the transient loading, as explained next.
Establish initial conditions. The only initial condition that may be explicitly established is the
initial displacement. A static solution using the mode superposition method is always
performed as the first solution with the given loads. If PowerDynamics was used for the modal
solution, no loads or displacements are allowed (that is, only u=0 is valid as the initial
condition). For pseudo-static analyses, the mode superposition method may yield poor results at
TIME=0.
The following load step options are available for the first load step:
Table 5-7 Options for the first load step
Option
Dynamics Options
Transient Integration
Parameters
TINTP
Load Vector
LVSCALE
Damping
ALPHAD
BETAD
Main Menu>Solution>-Load Step Opts-Time/Frequenc>
DMPRAT
Damping
MP,
DAMP
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 44 of 55
General Options
Integration Time Step
DELTIM
OUTPR
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 45 of 55
assumed to be 1/(20f), where f is the highest frequency chosen for the solution. The DELTIM
command is valid only in the first load step and is ignored in subsequent load steps.
Note-If you do issue the TIME command in the first load step, it will be ignored. The first
solution is always a static solution at TIME = 0.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 46 of 55
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 47 of 55
The mode superposition transient solution is written to the reduced displacement file, Jobname.RDSP,
regardless of whether the subspace or reduced method was used for the modal solution. You will
therefore need to expand the solution if you are interested in stress results.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
LSWRITE
------TIME,...
KBC,...
OUTRES,...
--LSWRITE
SAVE
LSSOLVE
FINISH
Page 48 of 55
See the ANSYS Commands Reference for a discussion of the ANTYPE, MODOPT, M, TOTAL,
ACEL, TRNOPT, LVSCALE, MDAMP, DELTIM, TIME, KBC, OUTRES, LSSOLVE, FILE,
SOLU, NSOL, PLVAR, PRVAR, EXPASS, NUMEXP, OUTRES, PLDISP, PRRSOL, PLNSOL
and PRERR commands.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 49 of 55
2. For all subsequent load steps, turn time integration effects on [TIMINT,ON], and proceed using
the full transient dynamic analysis procedures described previously. Once all load steps are written to
files [LSWRITE], you can initiate the multiple load step solution [LSSOLVE].
Note-The static prestress solution must be done as a separate solution if initial conditions are to
be defined with the IC command (Main Menu>Solution>-Loads-Apply> Initial
Condit'n>Define). The IC command is valid only in the first load step.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 50 of 55
2. Resolve the applied load-versus-time curve(s). The time step should be small enough to "follow"
the loading function. The response tends to lag the applied loads, especially for stepped loads, as
shown in Figure 5-5. Stepped loads require a small ITS at the time of the step change so that the step
change can be closely followed. ITS values as small as 1/180f may be needed to follow stepped loads.
Figure 5-5 Transient input versus transient response
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 51 of 55
3. Resolve the contact frequency. In problems involving contact (impact), the time step should be
small enough to capture the momentum transfer between the two contacting surfaces. Otherwise, an
apparent energy loss will occur and the impact will not be perfectly elastic. The integration time step
can be determined from the contact frequency (fc) as:
where k is the gap stiffness, m is the effective mass acting at the gap, and N is the number of
points per cycle. To minimize the energy loss, at least thirty points per cycle of (N=30) are
needed. Larger values of N may be required if acceleration results are needed. For the reduced
and mode superposition methods, N must be at least 7 to ensure stability.
You can use fewer than thirty points per cycle during impact if the contact period and contact
mass are much less than the overall transient time and system mass, because the effect of any
energy loss on the total response would be small.
4. Resolve the wave propagation. If you are interested in wave propagation effects, the time step
should be small enough to capture the wave as it travels through the elements. See "Build the Model"
for a discussion of element size.
5. Resolve the nonlinearities. For most nonlinear problems, a time step that satisfies the preceding
four guidelines is sufficient to resolve the nonlinearities. There are a few exceptions, however: if the
structure tends to stiffen under the loading (for example, large deflection problems that change from
bending to membrane load-carrying behavior), the higher frequency modes that are excited will have
to be resolved.
After calculating the time step using the appropriate guidelines, use the minimum value for your
analysis. By using automatic time stepping, you can let the ANSYS program decide when to increase
or decrease the time step during the solution. Automatic time stepping is discussed next.
Caution: Avoid using exceedingly small time steps, especially when establishing initial conditions.
Exceedingly small numbers can cause numerical difficulties. Based on a problem time scale of unity,
for example, time steps smaller than 10-10 could cause numerical difficulties.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 52 of 55
resource savings. Also, the number of times that you might have to rerun the analysis (adjusting the
time step size, nonlinearities, etc.) is greatly reduced. If nonlinearities are present, automatic time
stepping gives the added advantage of incrementing the loads appropriately and retreating to the
previous converged solution (bisection) if convergence is not obtained. You can activate automatic
time stepping with the AUTOTS command. (For more information on automatic time stepping in the
context of nonlinearities, see Chapter 8.)
Although it seems like a good idea to activate automatic time stepping for all analyses, there are some
cases where it may not be beneficial (and may even be harmful):
Problems that have only localized dynamic behavior (for example, turbine blade and hub
assemblies), where the low-frequency energy content of part of the system may dominate the
high-frequency areas
Problems that are constantly excited (for example, seismic loading), where the time step tends
to change continually as different frequencies are excited
Kinematics (rigid-body motion) problems, where the rigid-body contribution to the response
frequency term may dominate
5.12.3 Damping
Damping is present in most systems and should be specified in a dynamic analysis. Five forms of
damping are available in the ANSYS program:
Only the constant damping ratio and modal damping are available in the ANSYS/LinearPlus program.
You can specify more than one form of damping in a model. The program will formulate the damping
matrix [C] as the sum of all the specified forms of damping.
Table 5-8 shows the types of damping available for different structural analyses.
Table 5-8 Damping for different analysis types
Element
Modal
Damping Damping3
Analysis
Type
Constant
Damping
Ratio
Static
N.A.
N.A.
N.A.
N.A.
N.A.
Modal
No5
No5
No5
No
No
Undamped Yes
Yes
No
No
Yes
(COMBIN7,
[DMPRAT] [MDAMP] etc.)
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 53 of 55
Damped
Harmonic
Yes
Yes
Yes
No
Yes
Reduced
Yes
Yes
Yes
No
Yes
Mode Sup
Yes
Yes4
Yes
Yes
No
Yes
Yes
No
No
Yes
Reduced
Yes
Yes
No
No
Yes
Mode Sup
Yes
Yes4
Yes
Yes
No
Spectrum
SPRS,
MPRS2
Yes1
Yes
Yes
Yes
No
DDAM2
Yes1
Yes
Yes
Yes
No
PSD
Yes
No
Yes
Yes
No
Buckling
N.A.
N.A.
N.A.
N.A.
N.A.
Yes
No
No
Yes
Full
Transient
Full
Substructure Yes
Alpha damping and Beta damping are used to define Rayleigh damping constants and . The
damping matrix [C] is calculated by using these constants to multiply the mass matrix [M] and
stiffness matrix [K]:
The ALPHAD and BETAD commands are used to specify and , respectively, as decimal numbers.
The values of and are not generally known directly, but are calculated from modal damping ratios,
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 54 of 55
. is the ratio of actual damping to critical damping for a particular mode of vibration, i. If
natural circular frequency of mode i, and satisfy the relation
is the
In many practical structural problems, alpha damping (or mass damping) may be ignored ( =0). In
such cases, you can evaluate from known values of and , as
Only one value of can be input in a load step, so choose the most dominant frequency active in that
load step to calculate .
To specify both and for a given damping ratio , it is commonly assumed that the sum of the and
terms is nearly constant over a range of frequencies (see Figure 5-6). Therefore, given and a
frequency range to , two simultaneous equations can be solved for and .
Figure 5-6 Rayleigh damping
Alpha damping can lead to undesirable results if an artificially large mass has been introduced into the
model. One common example is when an artificially large mass is added to the base of a structure to
facilitate acceleration spectrum input. (You can use the large mass to convert an acceleration
spectrum to a force spectrum.) The alpha damping coefficient, which is multiplied by the mass matrix,
will produce artificially large damping forces in such a system, leading to inaccuracies in the
spectrum input, as well as in the system response.
Beta damping and material damping can lead to undesirable results in a nonlinear analysis. These
damping coefficients are multiplied by the stiffness matrix, which is constantly changing in a
nonlinear analysis. The resulting change in damping can sometimes be opposite to the actual change
in damping that can occur in physical structures. For example, whereas physical systems that
experience softening due to plastic response will usually experience a corresponding increase in
damping, an ANSYS model that has beta damping will experience a decrease in damping as plastic
softening response develops.
Material-dependent damping allows you to specify beta damping ( ) as a material property
[MP,DAMP]. Note, however, that MP,DAMP in a spectrum analysis [ANTYPE,SPECTR] specifies
a material-dependent damping ratio , not . Also note that for multi-material elements such as
SOLID46, SOLID65, SHELL91, and SHELL99, can only be specified for the element as a whole,
not for each material in the element. In these cases, is determined from the material pointer for the
element (set with the MAT command), rather than the material pointed to by any real constant MAT
for the element.
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015
Page 55 of 55
The constant damping ratio is the simplest way of specifying damping in the structure. It represents
the ratio of actual damping to critical damping, and is specified as a decimal number with the
DMPRAT command. DMPRAT is available only for spectrum, harmonic response, and mode
superposition transient dynamic analyses.
Modal damping gives you the ability to specify different damping ratios for different modes of
vibration. It is specified with the MDAMP command and is available only for the spectrum and mode
superposition method of solution (transient dynamic and harmonic response analyses).
Element damping involves using element types having viscous damping characteristics, such as
COMBIN7, COMBIN14, COMBIN37, COMBIN40, etc.
For more details about damping, see the ANSYS Theory Reference.
Go to the beginning of this chapter
http://mostreal.sk/html/guide_55/g-str/GSTR5.htm
12/22/2015