Apostila Ansys Training
Apostila Ansys Training
Apostila Ansys Training
3rd Edition
ANSYS Workbench -
Simulation Introduction
Training Manual
ANSYS Workbench Simulation Introduction
Registered Trademarks:
ANSYS is a registered trademark of SAS IP Inc.
All other product names mentioned in this manual are trademarks or registered trademarks of their respective
manufacturers.
Disclaimer Notice:
This document has been reviewed and approved in accordance with the ANSYS, Inc. Documentation Review
and Approval Procedures. This ANSYS Inc. software product (the Program) and program documentation
(Documentation) are furnished by ANSYS, Inc. under an ANSYS Software License Agreement that contains
provisions concerning non-disclosure, copying, length and nature of use, warranties, disclaimers and
remedies, and other provisions. The Program and Documentation may be used or copied only in accordance
with the terms of that License Agreement.
Introduction
Introduction
Welcome! Training Manual
ANSYS, Inc.
On-line documentation
Training manuals from other ANSYS courses
March 29, 2005
Inventory #002215
1-8
Introduction
Technical Support Training Manual
DesignModeler FE Modeler
Thelist here reflects analysis capabilities possible within the Simulation GUI. March 29, 2005
Inventory #002215
Note, however, that the ANSYS license used dictates what functionality is available to the user. 1-12
Introduction
Simulation Overview Training Manual
Thelist here reflects analysis capabilities possible within the Simulation GUI. March 29, 2005
Inventory #002215
Note, however, that the ANSYS license used dictates what functionality is available to the user. 1-13
Introduction
Simulation Overview Training Manual
Depending on what is
highlighted on the Outline
Tree, the available options
in the Context Toolbar and
Details View will change.
The Graphics Window will
also update to reflect the
part of the branch selected.
The Geometry Interface license can be run in reader mode for all
licenses. The Geometry Interface can be run in plug-in mode for the
CAD software listed under Associative
Detailed CAD support listing in Simulation online help
IGES Geometry Interface reader is free
On UNIX (Solaris, HP-UX), UG NX 1.0 is supported
For SolidWorks 2004, Service Pack 4 is required
DesignModeler
Display
Triad
Simulation Basics
Simulation Basics
Chapter Overview Training Manual
Tree Outline
Simulation
Details View Graphics Window Wizard
Short descriptions of each component are covered next March 29, 2005
Inventory #002215
2-4
Simulation Basics
GUI Menus Training Manual
Standard Toolbar
Graphics Toolbar
Named Selection Toolbar
Context Toolbar
Select mode Select entities Select Adjacent Graphics Manipulation Fit All Wireframe Viewports
The Outline Tree shows icons for each branch, along with a
This startup panel can be turned off in Tools > Options > Simulation: Startup March 29, 2005
Inventory #002215
2-13
Simulation Basics
Startup Panel Training Manual
General Preprocessing
General Preprocessing Procedure
Chapter Overview Training Manual
Part Relevance=+100
Part Relevance=-100
Cut Boundaries
Matched Faces
Procedure:
Rotation CS
March 29, 2005
Inventory #002215
3-49
General Preprocessing Procedure
Hex-Dominant Meshing Training Manual
Virtual cells provide the user with another set of tools to aid
K x F
This results in certain assumptions related to the analysis:
[K] is essentially constant
Linear elastic material behavior is assumed
Small deflection theory is used
Some nonlinear boundary conditions may be included
{F} is statically applied
No time-varying forces are considered
No inertial effects (mass, damping) are included
Pressure loading:
Force loading:
Forces can be applied on vertices, edges, or surfaces.
The force will be distributed on all entities. This
means that if a force is applied to two identical
surfaces, each surface will have half of the force
applied. Units are mass*length/time2
A force is defined via vector and magnitude or by
components (in user-defined Coordinate System)
ANSYS License Availability
DesignSpace Entra x
DesignSpace x March 29, 2005
Professional x Inventory #002215
Structural x
Mechanical/Multiphysics x 4-25
Linear Static Structural Analysis
Bearing Load Training Manual
Moment Load:
Remote Load:
Fixed Support:
Given Displacement:
Applies known displacement on vertex, edge, or surface
Allows for imposed translational displacement in x, y, and z (in
user-defined Coordinate System)
Entering 0 means that the direction is constrained.
Leaving the direction blank means that the entity is free to
move in that direction
Frictionless Support:
Cylindrical Constraint:
Applied on cylindrical surfaces
User can specify whether axial, radial, or tangential
components are constrained
Suitable for small-deflection (linear) analysis only
Simply Supported:
Fixed Rotation:
Can be applied on surface, edge, or vertex of surface or line
bodies
Constrains rotations but translations are free
All of the contour and vector plots are usually shown on the
Contact Results:
Retaining Ring
March 29, 2005
Inventory #002215
4-54
Chapter Five
Vibration Analysis
Vibration Analysis
Chapter Overview Training Manual
K M 0
i
2
i
This results in certain assumptions related to the analysis:
[K] and [M] are constant:
Linear elastic material behavior is assumed
Small deflection theory is used, and no nonlinearities included
[C] is not present, so damping is not included
{F} is not present, so no excitation of the structure is assumed
The structure can be unconstrained (rigid-body modes present) or
partially/fully constrained, depending on the physical structure
Mode shapes {} are relative values, not absolute
After setting up the model, one can solve the free vibration
K xo F
Based on the stress state from the static analysis, a stress
stiffness matrix [S] is calculated:
o S
The free vibration with pre-stress analysis is then solved, including
the [S] term
K S i M i 0
March 29, 2005
Inventory #002215
5-18
Vibration Analysis
Procedure w/ Pre-Stress Effects Training Manual
Thermal Analysis
Steady-State Thermal Analysis
Chapter Overview Training Manual
K T T QT
This results in the following assumptions:
No transient effects are considered in a steady-state analysis
[K] can be constant or a function of temperature
Temperature-dependent thermal conductivity can be input for each
material property
{Q} can be constant or a function of temperature
Temperature-dependent film coefficients can be input for
convective boundary conditions
Pinball Radius
Heat Flow:
Heat Flux:
A heat flux can be applied to surfaces only.
Heat flux has units of energy/time/area (i.e., power/area)
A positive value for heat load will add energy to the system.
Also, if multiple loads are present, the effect is cumulative.
ANSYS License Availability
DesignSpace Entra x
DesignSpace x March 29, 2005
Professional x Inventory #002215
Structural
Mechanical/Multiphysics x 6-22
Steady-State Thermal Analysis
Adiabatic Conditions Training Manual
Perfectly Insulated:
Given Temperature:
This imposes a temperature on vertices, edges, or surfaces.
Temperature is the degree of freedom solved for, but this fixes
the temperature on selected entities to a given value.
ANSYS License Availability
DesignSpace Entra x
DesignSpace x March 29, 2005
Professional x Inventory #002215
Structural
Mechanical/Multiphysics x 6-24
Steady-State Thermal Analysis
Thermal Boundary Conditions Training Manual
Convection:
q hATsurface Tambient
temperature:
If any temperature-dependent
convection load is applied, this
will result in a nonlinear solution
since the surface temperature is
solved for, but the film coefficient
h is based on a function of the
surface temperature.
The only exception is if the film
coefficient h is based on a
function of the bulk temperature
only. In Simulation, the bulk
Right mouse click on the table temperature is constant and input
to add or delete values. by the user, so this load will not
be nonlinear.
There are some types of loads that do not have any analogy
There is no thermal equivalent for inertial loads such as rotational
velocity or acceleration
The analogy of convective boundary condition is a foundation
stiffness support in structural terms, similar to a grounded spring
F F
Stable Unstable
K i S i 0
This results in certain assumptions related to the analysis:
[K] and [S] are constant:
Linear elastic material behavior is assumed
Small deflection theory is used, and no nonlinearities included
The response based on loading {F} is a linear function of i
Some additional restrictions:
Nonzero displacement supports or thermal loads are not allowed
Actual _ Load
BucklingLo ad
In the first case the load Multiplier is simply interpreted as the
buckling load.
In the second case the Load Multiplier is interpreted as a safety
factor.
Results Postprocessing
Results Postprocessing
Chapter Overview Training Manual
Displacement Scaling Display Method Contour Settings Outline Display Slice Planes
Proportional Vectors Equal Length Vectors Grid Aligned Element Aligned Line Form Solid Form
Examples
After the solution is complete, one can view the results and
Region of interest
March 29, 2005
Inventory #002215
Model shown is from a sample Mechanical Desktop assembly. 8-44
Results Postprocessing
Convergence & Scoping Example Training Manual
The Geometry Interface license can be run in reader mode for all
licenses. The Geometry Interface can be run in plug-in mode for the
CAD software listed under Associative
Detailed CAD support listing in Simulation online help
IGES Geometry Interface reader is free
One may wish to compare the results of two different load cases, as shown here March 29, 2005
Inventory #002215
9-31
CAD & Parameters
Multiple Design Studies Training Manual
DesignXplorer and
DesignXplorer VT can be used
to examine the effect input
parameters have on output
parameters.
After running the analysis,
through XY graphs, response
surface plots, sensitivity pie
charts, and spider plots, the
user can get a better
understanding of which
parameters affect certain
output and by what degree.
This extends the use of
parameters far beyond what is
possible with the Parameter
Manager.
Harmonic Analysis
Harmonic Analysis
Chapter Overview Training Manual
F Fmax e ej jt
x xmax e j e jt
The use of complex notation is an efficient representation of
the response. Since ejA is simply (cos(A)+jsin(A)), this
represents sinusoidal motion with a phase shift, which is
present because of the imaginary (j=-1) term.
The excitation frequency is the frequency at which the
loading occurs. A force phase shift may be present if
different loads are excited at different phases, and a
displacement phase shift may exist if damping or a force
phase shift is present.
March 29, 2005
Inventory #002215
10-4
Harmonic Analysis
Background on Harmonic Analysis Training Manual
Force Value
Force 1
0
Force 2
F Fmax e j e jt
Fmax cos jFmax sin e jt
F1 jF2 e jt
In this way, a real component F1 and
an imaginary component F2 are used.
The response {x} is analogous to {F}
March 29, 2005
Inventory #002215
Model shown is from a sample SolidWorks assembly. 10-5
Harmonic Analysis
Basics of Harmonic Analysis Training Manual
M jC K x jx F jF
2
1 2 1 2
This results in the following assumptions:
[M], [C], and [K] are constant:
Linear elastic material behavior is assumed
Small deflection theory is used, and no nonlinearities included
Damping [C] should be included. Otherwise, if the excitation
frequency is the same as the natural frequency of the
structure, the response is infinite at resonance.
The loading {F} (and response {x}) is sinusoidal at a given
frequency , although a phase shift may be present
f max f min
2
n
Simulation will solve n frequencies,
starting from .
In the example above, with a
frequency range of 0 10,000 Hz
at 10 intervals, this means that
Simulation will solve for 10
excitation frequencies of 1000,
ANSYS License Availability
DesignSpace Entra
2000, 3000, 4000, 5000, 6000, 7000,
DesignSpace
Professional x
8000, 9000, and 10000 Hz. March 29, 2005
Inventory #002215
Structural x
Mechanical/Multiphysics x 10-15
Harmonic Analysis
Solution Methods Training Manual
M jC K x jx F jF
2
1 2 1 2
For linear systems, one can express the displacements x as a
linear combination of mode shapes i :
n
x yi i
i 1
where yi are modal coordinates (coefficient) for this relation.
For example, one can perform a modal analysis to determine the
natural frequencies i and corresponding mode shapes i.
One can see that as more modes n are included, the approximation
for {x} becomes more accurate.
ANSYS License Availability
DesignSpace Entra
DesignSpace March 29, 2005
Professional x Inventory #002215
Structural x
Mechanical/Multiphysics x 10-17
Harmonic Analysis
Mode Superposition Method Training Manual
M jC K x jx F jF
2
1 2 1 2
In the Full method, this matrix equation is solved for directly in
nodal coordinates, analogous to a linear static analysis except
that complex numbers are used:
K C M jC K
2
xC x1 jx2
FC F1 jF2
ANSYS License
DesignSpace Entra
Availability
K C xC FC
DesignSpace March 29, 2005
Professional Inventory #002215
Structural x
Mechanical/Multiphysics x 10-21
Harmonic Analysis
Full Method Training Manual
Asynchronous Solution
Asynchronous Solution
Chapter Overview Training Manual
General Preprocessing
General Preprocessing Procedure
Global Mesh Controls (ANSYS) Training Manual
M x C x K x F t
where [M] is the mass matrix, [C] is the damping matrix, [K]
is the stiffness matrix, {x} is displacements, and {F} is force
Because this is a static analysis, all time-dependent terms
are removed, leaving the following subset:
K x F
March 29, 2005
Inventory #002215
A4-2
General Preprocessing Procedure
Elements Used Training Manual
Mohr-Coulomb Theory:
c
1
t
Equivalent Stress:
e
1
2
1 2 2 2 3 2 3 1 2
This criterion is commonly used for ductile metals.
When uniaxial tensile tests of specimens are performed to
determine the yield strength and stress-strain relationships,
the engineer needs a way to relate the uniaxial data to the
stress state (tensor). Hence, the equivalent stress is a
commonly used scalar invariant for this purpose.
Stress Intensity:
The stress intensity is twice the value of the maximum shear
stress.
The stress intensity provides the value of the largest
difference between principal stresses
Modal Analysis
General Preprocessing Procedure
Basics of Free Vibration Analysis Training Manual
M x C x K x F t
M x K x 0 March 29, 2005
Inventory #002215
A5-2
General Preprocessing Procedure
Basics of Free Vibration Analysis Training Manual
x i cos it
where i is the mode shape (eigenvector) and i is the
natural circular frequency for mode i.
By substituting this value in the earlier equation, the
following is obtained:
M i cos i t K i cos i t 0
i
2
M K cos t 0
i
2
i i
Noting that the solution i =0 is trivial, i is solved for:
K M 0
i
2
i March 29, 2005
Inventory #002215
A5-3
General Preprocessing Procedure
Requesting Results Training Manual
Thermal Analysis
General Preprocessing Procedure
Basics of Steady-State Heat Transfer Training Manual
C T T K T T Qt , T
where t is time and {T} is temperature, [C] is the specific
heat (thermal capacitance) matrix, [K] is the conductivity
matrix, and {Q} is the heat flow rate load vector.
In a steady-state analysis, all time-dependent terms are
removed. However, nonlinearities can be present:
K S x F
If we consider the analysis to be linear, we can multiply the
load and the stress state by a constant :
K S x F
In a buckling mode, displacements can be large (x+) without
an increase in load, so the following is also true:
If the last two equations are subtracted from each other, the
K S 0
The above equation is what is solved for during a linear
buckling analysis.
The buckling load multiplier is multiplied to the applied loads to
get the critical load for buckling
The buckling mode shape expresses the shape of buckling.
However, the magnitude is not known since is indeterminate.
There are actually many buckling load multipliers and modes,
although the user is usually interested in the first few modes since
these would occur before any higher buckling modes.
Note the similarity of linear buckling equation with the free
vibration equation (Chapter 5). Both are known as eigenvalue
problems which are solved for with similar matrix methods.
K xo F
Based on the stress state from the static analysis, a stress
stiffness matrix [S] is calculated:
o S
The aforementioned eigenvalue problem is then solved to get
the buckling load multiplier i and buckling modes i:
K i S i 0
March 29, 2005
Inventory #002215
A7-4
General Preprocessing Procedure
Requesting Results Training Manual
Harmonic Analysis
General Preprocessing Procedure
Background on Harmonic Analysis Training Manual
2
2 i i j i2 y1i jy2i i F1 jF2
T
Fatigue Module
Fatigue Module
Chapter Overview Training Manual
max
min
Biaxiality Indication B
i
Stress biaxiality contour plot helps to j
determine the state of stress at a location
Biaxiality indication is the ratio of the
smaller to larger principal stress (with
principal stress nearest to 0 ignored).
Hence, locations of uniaxial stress report 0,
pure shear report -1, and biaxial reports 1.
Recall that usually fatigue test data is reflective of a test specimen under
uniaxial stress (although torsional tests would be in pure shear).
The biaxiality indication helps to determine if a location of interest is in a
stress state similar to testing conditions. In this example, the location of March 29, 2005
interest (center) has a value of -1, so it is predominantly in shear. Inventory #002215
A12-31
Fatigue Module
Reviewing Fatigue Results Training Manual
Fatigue Sensitivity:
A fatigue sensitivity chart displays how life,
damage, or safety factor at the critical
location varies with respect to load
Load variation limits can be input (including
negative percentages)
Defaults for chart options available under
Tools menu > Options Simulation:
Fatigue > Sensitivity
Set up a stress analysis (linear, Define fatigue material properties, Specify loading type and
proportional loading) including S-N curve(s) treatment of mean stress effects
Set up a stress analysis (linear, Define fatigue material properties, Specify loading history data and
proportional loading) including S-N curve(s) treatment of mean stress effects
Specify number of bins for Solve and review fatigue results, (e.g., damage matrix, damage contour, March 29, 2005
rainflow cycle counting life contour, etc.) Inventory #002215
A12-39
Fatigue Module
Variable Amplitude Procedure Training Manual
The bin size can range from 10 to 200. The default value is
32, and it can be changed in the Control Panel. March 29, 2005
Inventory #002215
A12-43
Fatigue Module
Specifying Bin Size Training Manual
On the other hand, the bolt load only distributes load on the
compressive side. Hence, to mimic the loading in reverse, the
bolt load needs to be applied in a separate Environment in the
opposite direction.
Shape Finder
Shape Finder
Chapter Overview Training Manual
Any loads and supports may be used with the Shape Finder
For bodies that results are scoped to (see next Chapter), these