ANSYS Mechanical APDL Advanced Tutorials
ANSYS Mechanical APDL Advanced Tutorials
Guide
ANSYS, Inc.
Southpointe
275 Technology Drive
Canonsburg, PA 15317
ansysinfo@ansys.com
http://www.ansys.com
(T) 724-746-3304
(F) 724-514-9494
Release 14.5
October 2012
ANSYS, Inc. is
certified to ISO
9001:2008.
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFIDENTIAL AND PROPRIETARY PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. The software products
and documentation are furnished by ANSYS, Inc., its subsidiaries, or affiliates under a software license agreement
that contains provisions concerning non-disclosure, copying, length and nature of use, compliance with exporting
laws, warranties, disclaimers, limitations of liability, and remedies, and other provisions. The software products
and documentation may be used, disclosed, transferred, or copied only in accordance with the terms and conditions
of that software license agreement.
ANSYS, Inc. is certified to ISO 9001:2008.
Third-Party Software
See the legal information in the product help files for the complete Legal Notice for ANSYS proprietary software
and third-party software. If you are unable to access the Legal Notice, please contact ANSYS, Inc.
Published in the U.S.A.
Table of Contents
1. Probabilistic Design ................................................................................................................................ 1
1.1. Understanding Probabilistic Design .................................................................................................. 1
1.1.1. Traditional (Deterministic) vs. Probabilistic Design Analysis Methods ......................................... 2
1.1.2. Reliability and Quality Issues .................................................................................................... 3
1.2. Probabilistic Design Terminology ...................................................................................................... 3
1.3. Using Probabilistic Design ................................................................................................................. 8
1.3.1. Create the Analysis File ............................................................................................................. 9
1.3.1.1. Example Problem Description ........................................................................................ 10
1.3.1.2. Build the Model Parametrically ....................................................................................... 10
1.3.1.3. Obtain the Solution ....................................................................................................... 11
1.3.1.4. Retrieve Results and Assign as Output Parameters .......................................................... 11
1.3.1.5. Prepare the Analysis File ................................................................................................ 12
1.3.2. Establish Parameters for Probabilistic Design Analysis ............................................................. 12
1.3.3. Enter the PDS and Specify the Analysis File ............................................................................. 13
1.3.4. Declare Random Input Variables ............................................................................................. 14
1.3.5. Visualize Random Input Variables ........................................................................................... 20
1.3.6. Specify Correlations Between Random Variables ..................................................................... 20
1.3.7. Specify Random Output Parameters ....................................................................................... 23
1.3.8. Select a Probabilistic Design Method ...................................................................................... 24
1.3.8.1. Probabilistic Method Determination Wizard ................................................................... 24
1.3.9. Execute Probabilistic Analysis Simulation Loops ...................................................................... 25
1.3.9.1. Probabilistic Design Looping ......................................................................................... 26
1.3.9.2. Serial Analysis Runs ....................................................................................................... 27
1.3.9.3. PDS Parallel Analysis Runs .............................................................................................. 27
1.3.9.3.1. Machine Configurations ........................................................................................ 29
1.3.9.3.1.1. Choosing Slave Machines ............................................................................. 29
1.3.9.3.1.2. Using the Remote Shell Option ..................................................................... 29
1.3.9.3.1.3. Using the Connection Port Option ................................................................ 31
1.3.9.3.1.4. Configuring the Master Machine ................................................................... 33
1.3.9.3.1.5. Host setup using port option ........................................................................ 35
1.3.9.3.1.6. Host and Product selection for a particular analysis ....................................... 35
1.3.9.3.2. Files Needed for Parallel Run ................................................................................. 36
1.3.9.3.3. Controlling Server Processes ................................................................................. 37
1.3.9.3.4. Initiate Parallel Run ............................................................................................... 38
1.3.10. Fit and Use Response Surfaces .............................................................................................. 38
1.3.10.1. About Response Surface Sets ....................................................................................... 39
1.3.10.2. Fitting a Response Surface ........................................................................................... 39
1.3.10.3. Plotting a Response Surface ......................................................................................... 40
1.3.10.4. Printing a Response Surface ......................................................................................... 40
1.3.10.5. Generating Monte Carlo Simulation Samples on the Response Surfaces ........................ 41
1.3.11. Review Results Data ............................................................................................................. 41
1.3.11.1. Viewing Statistics ......................................................................................................... 42
1.3.11.2. Viewing Trends ............................................................................................................ 43
1.3.11.3. Creating Reports .......................................................................................................... 44
1.4. Guidelines for Selecting Probabilistic Design Variables ..................................................................... 44
1.4.1. Choosing and Defining Random Input Variables ..................................................................... 44
1.4.1.1. Random Input Variables for Monte Carlo Simulations ..................................................... 45
1.4.1.2. Random Input Variables for Response Surface Analyses .................................................. 45
1.4.1.3. Choosing a Distribution for a Random Variable ............................................................... 45
1.4.1.3.1. Measured Data ..................................................................................................... 45
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
iii
iv
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
vi
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
vii
viii
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
ix
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
xi
xii
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
List of Figures
1.1. A Beam Under a Snow Load .................................................................................................................... 4
1.2. Probabilistic Design Data Flow ................................................................................................................ 8
1.3. A Beam Under a Snow Load .................................................................................................................. 10
1.4. Histograms for the Snow Height H1 and H2 ........................................................................................... 19
1.5. A Scatter Plot of Snow Height H1 vs. H2 ................................................................................................. 21
1.6. The PDS Method Determination Wizard ................................................................................................. 25
1.7. Graph of X1 and X2 Showing Two Samples with Close Values ................................................................. 52
1.8. Graph of X1 and X2 Showing Good Sample Distribution ........................................................................ 53
1.9. Locations of Sampling Points for Problem with Three Input Variables for CCD ......................................... 57
1.10. Location of Sampling Points for Problem with Three Input Variables for BBM ........................................ 58
1.11. Cumulative Distribution Function of X ................................................................................................. 61
1.12. Sensitivities ........................................................................................................................................ 63
1.13. Range of Scatter .................................................................................................................................. 63
1.14. Effects of Reducing and Shifting Range of Scatter ................................................................................ 66
1.15. The Simple Indeterminate Three-Bar Truss for the Example Problem ..................................................... 70
2.1. Element Plot for Waveguide Example .................................................................................................... 79
2.2. Graph of Phase Angle ............................................................................................................................ 80
2.3. Graph of Magnitude ............................................................................................................................. 81
3.1. Selective Adaptivity .............................................................................................................................. 87
4.1. Rezoning Using a Program-Generated New Mesh .................................................................................. 99
4.2. Rezoning Using a Generic New Mesh Generated by Another Application ............................................... 99
4.3. Rezoning Using Manual Splitting of an Existing Mesh .......................................................................... 100
4.4. Boundary Geometry of a Generic (CDB) New Mesh .............................................................................. 108
4.5. Remeshing Options when Using a Generic (CDB) New Mesh ................................................................ 108
4.6. Splitting of Quadrilateral and Degenerate Linear Elements (PLANE182) ................................................ 111
4.7. Splitting of Quadrilateral, Degenerate and Triangular Quadratic Elements (PLANE183) .......................... 112
4.8. Transition Element Generation Methods ............................................................................................. 113
4.9. /PREP7 Mesh-Control Commands Available in Rezoning ...................................................................... 115
5.1. 2-D to 3-D Process Flow ....................................................................................................................... 138
5.2. 2-D Plane Strain to 3-D Solid Extrusion ................................................................................................. 140
5.3. Axisymmetric to 3-D Solid Extrusion .................................................................................................... 140
6.1. Hydro Rotor -- Model of a Cyclically Symmetric Structure ..................................................................... 143
6.2. A Basic Sector in a Cyclically Symmetric Structure ................................................................................ 145
6.3. Basic Sector Definition ........................................................................................................................ 145
6.4. Full Cyclic Model ................................................................................................................................. 147
6.5. Cyclic Sector ....................................................................................................................................... 147
6.6. Successful Auto Detection with Default FACETOL = 15 Deg ................................................................ 148
6.7. Auto Detection Failure Due to Large Face Tolerance ............................................................................. 149
6.8. Connecting Low and High Edges of Basic and Duplicate Sectors .......................................................... 152
6.9. Process Flow for a Static Cyclic Symmetry Analysis (Cyclic Loading) ...................................................... 156
6.10. Process Flow for a Static Cyclic Symmetry Analysis (Non-Cyclic Loading) ............................................. 156
6.11. Examples of Nodal Diameters (i) ........................................................................................................ 157
6.12. Process Flow for a Stress-Free Modal Cyclic Symmetry Analysis ........................................................... 159
6.13. Process Flow for a Prestressed Modal Cyclic Symmetry Analysis .......................................................... 160
6.14. Process Flow for a Large-Deflection Prestressed Modal Cyclic Symmetry Analysis ............................... 161
6.15. Process Flow for a Linear Buckling Cyclic Symmetry Analysis .............................................................. 162
6.16. Process Flow for a Harmonic Cyclic Symmetry Analysis (Non-Cyclic Loading) ...................................... 163
6.17. Process Flow for a Prestressed Harmonic Cyclic Symmetry Analysis .................................................... 166
6.18. Traveling Wave Animation Example ................................................................................................... 169
6.19. Example Modal Cyclic Symmetry Analysis Results .............................................................................. 177
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
xiii
xiv
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
List of Tables
6.1. Valid Non-Cyclically Symmetric Loads .................................................................................................. 153
6.2. Buckling Cyclic Symmetry: Load Factor Iteration Results ....................................................................... 184
9.1. Loads Applicable in a Substructure Analysis ......................................................................................... 243
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
xv
xvi
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Probabilistic Design
It is neither physically possible nor financially feasible to eliminate the scatter of input parameters
completely. The reduction of scatter is typically associated with higher costs either through better and
more precise manufacturing methods and processes or increased efforts in quality control; hence, accepting the existence of scatter and dealing with it rather than trying to eliminate it makes products
more affordable and production of those products more cost-effective.
To deal with uncertainties and scatter, you can use the Probabilistic Design System (PDS) to answer the
following questions:
If the input variables of a finite element model are subjected to scatter, how large is the scatter of the
output parameters? How robust are the output parameters? Here, output parameters can be any parameter
that the program can calculate. Examples are the temperature, stress, strain, or deflection at a node, the
maximum temperature, stress, strain, or deflection of the model, etc.
If the output is subjected to scatter due to the variation of the input variables, then what is the probability
that a design criterion given for the output parameters is no longer met? How large is the probability that
an unexpected and unwanted event takes place (what is the failure probability)?
Which input variables contribute the most to the scatter of an output parameter and to the failure probability? What are the sensitivities of the output parameter with respect to the input variables?
Probabilistic design can be used to determine the effect of one or more variables on the outcome of
the analysis. In addition to the probabilistic design techniques available, the program offers a set of
strategic tools that can be used to enhance the efficiency of the probabilistic design process. For example,
you can graph the effects of one input variable versus an output parameter, and you can easily add
more samples and additional analysis loops to refine your analysis.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
~16%
~2.3%
~22%
~8%
When a conservative assumption is used, this actually tells you that uncertainty or randomness is involved.
Conservative assumptions are usually expressed in terms of safety factors. Sometimes regulatory bodies
demand safety factors in certain procedural codes. If you are not faced with such restrictions or demands,
then using conservative assumptions and safety factors can lead to inefficient and costly over-design.
You can avoid over-design by using probabilistic methods while still ensuring the safety of the component.
Probabilistic methods even enable you to quantify the safety of the component by providing a probability that the component will survive operating conditions. Quantifying a goal is the necessary first step
toward achieving it. Probabilistic methods can tell you how to achieve your goal.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Probabilistic Design
Figure 1.1: A Beam Under a Snow Load
H2
H1
E
PDS Term
random input variables (RVs)
Description
Quantities that influence the result of an analysis.
In probabilistic design, RVs are often called "drivers" because they
drive the result of an analysis. You must specify the type of statistical
distribution the RVs follow and the parameter values of their distribution functions.
For the beam example, the heights H1 and H2 and the Young's modulus
E are clearly the random input variables. Naturally, the heights H1 and
H2 cannot be negative and more often there will be little snow and
only a few times there will be a lot of snow. Therefore, it might be
appropriate to model the height of the snow as an exponential or a
lognormal distribution, both of which have the bulk of the data at
lower values.
correlation
Two (or more) RVs which are statistically dependent upon each other.
In the beam example, it is unlikely that one side of the roof (beam)
supports a great deal of snow while almost no snow exists on the
other. It is not necessary that H1 and H2 are exactly identical, but with
a lot of snow then H1 and H2 both likely have larger values and with
little snowfall then both would be lower. Therefore, H1 and H2 are
correlated, although the correlation must not be mistaken for a direct
mathematical dependency. In the beam example, no numerical dependency exists but rather a certain trend between the two heights;
with this particular H1 and H2 it is unlikely that their values will be
drastically different at any given point in time.
Note
Mathematical dependencies have some numerical dependence on each other. For example, a true correlation exists
if, when one parameter value doubles, another parameter
value also doubles.
random
output
parameters
(RPs)
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
output parameters. In our beam example, the maximum beam deflection is a random output parameter.
probabilistic
design variables
The random input variables (RVs) and random output parameters (RPs)
are collectively known as probabilistic design variables.
In the Probabilistic Design System (PDS), you must identify which
parameters in the model are RVs and which are RPs.
sample
simulation
The collection of all samples that are required or that you request for
a given probabilistic analysis.
The simulation contains the information used to determine how the
component would behave under real-life conditions (with all the existing uncertainties and scatter); therefore, all samples represent the
simulation of the behavior.
analysis file
loop
or
In each loop, the PDS uses the values of the RVs from one sample and
executes the user-specified analysis. The PDS collects the values for
the RPs following each loop.
simulation
loop
loop file
The probabilistic design loop file (Jobname.LOOP), created automatically by via the analysis file.
The PDS uses the loop file to perform analysis loops.
probabilistic
model
Probabilistic Design
PDS Term
Description
Random input variables (RVs)
Correlations
Random output parameters (RPs)
The selected settings for probabilistic method and its parameters.
If you change any part of the probabilistic model, then you will generate different results for the probabilistic analysis (that is, different
results values and/or a different number of results). For example,
modifying the analysis file may affect the results file. If you add or
take away an RV or change its distribution function, you solve a different probabilistic problem (which again leads to different results). If
you add an RP, you will still solve the same probabilistic problem, but
more results are generated.
probabilistic
design database (PDS
database)
mean value
median
value
The statistical point where 50% of the data is below the median value
and the 50% is above.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
For symmetrical distribution functions (Gaussian, uniform, etc.) the
median value and the mean value are identical, while for nonsymmetrical distributions they are different.
standard
deviation
solution set
response
surface set
remote host
The following figure shows the flow of information during a probabilistic design analysis. The analysis
file must exist as a separate entity, and that the probabilistic design database is not part of the model
database.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Probabilistic Design
Figure 1.2: Probabilistic Design Data Flow
File.DB
ANSYS
Database
File
SAVE
RESUME
ANSYS
/EXIT
(Save
Everything)
Analysis File
(parametrically
defined model)
PDEXE
Model
Database
/CLEAR
Probabilsitic
Design
Database
PDSAVE
PDEXE
File.LOOP
Loop File
PDRESU
File.PDS
Prob. Design
Database
File
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Probabilistic Design
You prefer to perform the initial analysis interactively, and then use the resulting command log as the
basis for the analysis file. In this case, you must edit the log file to make it suitable for probabilistic
design looping. For more information about using the log files, see Using the ANSYS Session and
Command Logs in the Operations Guide.
H2
H1
E
Young's modulus is 20E4.
As mentioned earlier, you can vary virtually any aspect of the design: dimensions, shape, material
property, support placement, applied loads, etc. The only requirement is that the design be defined in
10
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Caution
If you build your model interactively (through the GUI), you will encounter many situations
where data can be input through graphical picking (such as when defining geometric entities).
Because some picking operations do not allow parametric input (and PDS requires parametric
input), you should avoid picking operations. Instead, use menu options that allow direct input
of parameters.
!
!
!
!
!
!
!
!
!
get
get
get
get
get
!
!
!
!
! Leave SOLUTION
This step is not limited to just one analysis. You can, for instance, obtain a thermal solution and then
obtain a stress solution (for thermal stress calculations).
If your solution uses the multiframe restart feature, all changes to the parameter set that are made after
the first load step will be lost in a multiframe restart. To ensure that the correct parameters are used
in a multiframe restart, you must explicitly save (PARSAV) and resume (PARESU) the parameters for
use in the restart. See the Basic Analysis Guide for more information on multiframe restarts.
11
Probabilistic Design
In our beam example, the maximum deflection and the maximum stress of the beam are random output
parameters (RPs). The parameters for these data may be defined as follows:
...
/POST1
SET,FIRST
NSORT,U,Y
*GET,DMAX,SORT,,MIN
!
! Derived data for line
ETABLE,VOLU,VOLU
ETABLE,SMAX_I,NMISC,1
ETABLE,SMAX_J,NMISC,3
!
ESORT,ETAB,SMAX_I,,1
*GET,SMAXI,SORT,,MAX
ESORT,ETAB,SMAX_J,,1
*GET,SMAXJ,SORT,,MAX
SMAX=SMAXI>SMAXJ
FINISH
...
Note
Do not use the /CLEAR command in your analysis file as this will delete the probabilistic
design database during looping. If this happens, the random input variables are no longer
recognized during looping and you will get the same (deterministic) results for all simulation
loops; however, resume the database using the RESUME command as part of your analysis
file. For example, this is helpful if the variations of the random input variables do not require
that meshing is done in every loop (because the mesh is not effected). In this case you can
mesh your model, save the database, and resume the database at the beginning of the
analysis file.
12
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The database does not need to contain model information corresponding to the analysis file
to perform probabilistic design. The model input is automatically read from the analysis file
during probabilistic design looping.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
13
Probabilistic Design
the analysis and the other containing all probabilistic design operations), then you will need to specify
the analysis file using PDANL after entering the PDS (/PDS).
Note
In the analysis file, the /PREP7 and /PDS commands must occur as the first nonblank characters
on a line. (Do not use the $ delimiter on a line containing either of these commands.) This
requirement is necessary for proper loop file construction.
You cannot assign a different analysis file using the PDANL command after a probabilistic
analysis has been performed. This ensures the integrity of the previously generated results
with the specified probabilistic model.
The program cannot restrain you from editing the analysis file or exchanging it with system level commands. If you do so, then it is your responsibility to ensure the integrity of the generated results with
the definitions in the analysis file. If you are not sure that this integrity is maintained or if you know
that it is not, then we recommend that you save the current PDS database via the PDSAVE command
and then clear the probabilistic analysis results from the probabilistic design database using the PDCLR,
POST command. The PDCLR command does not delete the result files that have been generated; it
erases the link to the result files from the database.
In the example of a beam supporting a roof with a snow load you could store the analysis file in a
macro named beam.mac. Here, the analysis is specified with the commands:
...
/PDS
PDANL,beam,mac
...
14
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
You provide values for the mean value and the standard deviation of the random variable x.
Truncated Gaussian (TGAU):
nim
You provide the mean value and the standard deviation of the non-truncated Gaussian distribution
and the truncation limits xmin and xmax.
Lognormal option 1 (LOG1):
You provide values for the mean value and the standard deviation of the random variable x. The
PDS calculates the logarithmic mean and the logarithmic deviation :
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
15
Probabilistic Design
Ax@f
x
You provide values for the logarithmic mean value and the logarithmic deviation . The parameters
and are the mean value and standard deviation of ln(x):
Triangular (TRIA):
vlm
nim
You provide the minimum value xmin, the most likely value limit xmlv and the maximum value xmax.
Uniform (UNIF):
16
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
xmx
x
You provide the lower and the upper limit xmin and xmax of the random variable x.
Exponential (EXPO):
You provide the decay parameter and the shift (or lower limit) xmin of the random variable x.
Beta (BETA):
tDr
You provide the shape parameters r and t and the lower and the upper limit xmin and xmax of the
random variable x.
Gamma (GAMM):
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
17
Probabilistic Design
Ax@f
k
x
rhDm
nim
You provide the Weibull characteristic value xchr , the Weibull exponent m and the minimum value
xmin. Special cases: For xmin = 0 the distribution coincides with a two-parameter Weibull distribution.
The Rayleigh distribution is a special case of the Weibull distribution with = xchr - xmin and m = 2.
You may change the specification of a previously declared random input variable by redefining it. You
may also delete a probabilistic design variable (PDVAR,Name,DEL). The delete option does not delete
the parameter from the database; it simply deactivates the parameter as a probabilistic design variable.
Note
Changing the probabilistic model by changing a random input variable is not allowed after
a probabilistic analysis has been performed. This ensures the integrity of the previously
generated results with the specified probabilistic model. If you need to change one or more
random input variables (for example, because you learned that some specifications were incorrect after running an analysis), then we recommend that you save the current PDS database
(using the PDSAVE command) and then clear the probabilistic analysis results from the
probabilistic design database (using the PDCLR,POST command). The PDCLR command does
not delete the result files that have been generated, it simply removes the link to the results
file from the database.
18
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
6.0E-01
Relative Frequency
5.0E-01
4.0E-01
3.0E-01
2.0E-01
1.0E-01
02
E+
66
02
3.
E+
85
02
2.
E+
2.
03
02
E+
22
1.
4.
07
E+
01
0.0E+00
Snow Height H1
8.0E-01
Relative Frequency
7.0E-01
6.0E-01
5.0E-01
4.0E-01
3.0E-01
2.0E-01
1.0E-01
9.
85
E+
01
2.
95
E+
02
4.
92
E+
02
6.
89
E+
02
8.
86
E+
02
1.
08
E+
03
0.0E+00
Snow Height H2
From these histograms you can conclude that an exponential distribution is suitable to describe the
scatter of the snow height data for H1 and H2. Suppose from the measured data we can evaluate that
the average snow height of H1 is 100 mm and the average snow height of H2 is 200 mm. The parameter
can be directly derived by 1.0 divided by the mean value which leads to 1 = 1/100 = 0.01 for H1,
and 1 = 1/200 = 0.005 for H2. From measurements of the Young's modulus you see that the Young's
modulus follows a Gaussian distribution with a standard deviation of 5%. Given a mean value of 200,000
N/mm2 for the Young's modulus this gives a standard deviation of 10,000 N/mm2. These definitions can
be specified using the following commands:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
19
Probabilistic Design
...
PDVAR,H1,EXPO,0.01
PDVAR,H2,EXPO,0.005
PDVAR,YOUNG,GAUS,200000,10000
...
20
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Snow height H2
500
400
300
200
100
0
100
200
300
400
Snow height H1
Performing a statistical evaluation of the data, we can conclude that the linear correlation coefficient
between the values of H1 and H2 is about 0.8. You can define this correlation using the commands:
...
PDVAR,H1,EXPO,0.01
PDVAR,H2,EXPO,0.005
PDCORR,H1,H2,0.8
...
You may have a more complex correlation where you have a spatial dependency. If so, you can use the
PDCFLD command to calculate a correlation field and store it into an array.
Random fields are random effects with a spatial distribution; the value of a random field not only varies
from simulation to simulation at any given location, but also from location to location. The correlation
field describes the correlation coefficient between two different spatial locations. Random fields can be
either based on element properties (typically material) or nodal properties (typically surface shape
defined by nodal coordinates). Hence, random fields are either associated with the selected nodes or
the selected elements. If a random field is associated with elements, then the correlation coefficients
of the random field are calculated based on the distance of the element centroids.
Note that for correlation fields, the domain distance D({xi} , {xj}) is not the spatial distance |{xi} - {xj}|,
but the length of a path between {xi} and {xj} that always remains inside the finite element domain.
However, exceptions are possible in extreme meshing cases. For elements that share at least one node,
the PDCFLD evaluates the distance by directly connecting the element centroids with a straight line.
If these neighboring elements form a sharp inward corner then it is possible that the domain distance
path lies partly outside the finite element domain, as illustrated below.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
21
Probabilistic Design
After the correlation coefficients have been calculated and stored in a parameter (PDCFLD,ParR), then
you can use the PDCORR command to define the correlations between the elements of the random
field.
Note
When specifying one variable (A) with correlations to two or more other variables (B, C,
etc.), be certain that you consider the relationship implied between the other variables
B and C, etc. If you specify high correlations between A and B and A and C, without
specifying the relationship between B and C, you might receive an error. Specifying a
relatively high correlation between A and B, with only a moderate correlation between
A and C might work because the logical correlation between B and C could still be low
or nonexistent.
Example
The structure illustrated below is modeled with 12 elements. We will evaluate the domain distances of
the element centroids.
Y
4.
1
3.
5
2.
1.
9
10
11
12
0.
0.
1.
2.
3.
4. X
22
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Next, calculate the domain distances and store the results in the array elemdist:
...
/PDS
PDCFLD,elemdist,ELEM,DIST
...
Finally, get all the element domain distances and print them:
...
*GET,numsel,ELEM,0,COUNT
! Get the number of selected elements
!
! Outer loop through all selected elements from first to last
index=0
elem1=0
! Pipe output to file
/OUT,elements,dat
*DO,i,1,numsel
elem1=ELNEXT(elem1)
! get number of next selected element
*IF,elem1,EQ,0,CYCLE
! Leave do loop if no more elements
!
! Inner loop through selected elements from "elem1+1" to last
elem2=elem1
*DO,j,i+1,numsel
elem2=ELNEXT(elem2)
! get number of next selected element
*IF,elem2,EQ,0,CYCLE
! Leave do loop if no more elements
index=index+1
!
! Print out the element distance
*MSG,INFO,elem1,elem2,elemdist(index)
Distance between element %i and %i is %g
*ENDDO
! go to next element for inner loop
*ENDDO
! go to next element for outer loop
...
The print out will show that for the structure illustrated above the "domain distance" between the element
centroids of elements 1 and 9 is 3.8284 and between the element centroids of elements 1 and 12 it is
4.8284. The paths related to these distances are sketched in the illustration with a solid line and a dashed
line respectively. In this example there are 12 elements, thus the array "elemdist" has a length of 12*(121)/2 = 66.
23
Probabilistic Design
GUI: Main Menu> Prob Design> Prob Definitns> Random Output
Note
To use Response Surface Methods, the random output parameters must be smooth and
continuous functions of the involved random input variables. Do not use Response
Surface Methods if this condition is not satisfied.
24
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Based on the information you provide, the wizard will tell you the fastest method for solving your
probabilistic design problem. The wizard will issue the PDMETH command and either the PDLHS or
the PDDOEL command. You will still need to run the analysis, then fit the results to a response surface,
etc. to evaluate your results.
Caution
For security reasons, ANSYS, Inc. recommends that you use parallel processing only within
the firewall of your local area network.
If you use serial processing to perform a probabilistic analysis, you will use only the CPU of the computer
you are working on. If you have access to only one product license, or if you have access to only one
computer, then this is the only way in which you can run a probabilistic analysis. While the simulation
loops are running in serial mode, your session is locked (you cannot perform other tasks in the same
session). If you are running your session in interactive mode, the simulation loops are also performed
in interactive mode. If you are running your session in batch mode, the simulation loops are performed
in batch mode.
If you select the PDS parallel-processing option, you can use other CPUs that you have access to for
running the probabilistic analysis. PDS parallel processing can distribute the necessary jobs in a local
area network. With this option, the simulation loops are sent to CPUs that you can specify, where they
are executed in "server mode." This looks the same as a batch run (in other words, there is no interactive
visualization during the execution of a simulation loop). While the simulation loops are running in parallel mode, your session is locked; however, you can instruct the session to start postprocessing the
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
25
Probabilistic Design
probabilistic results as they are calculated so you can review and visualize the results before all simulation
loops are finished.
In parallel processing, you can monitor the running jobs and the completed results.
When using parallel-processing, for n available licenses, n - 1 will be available for PDS solutions, as one
license is used as a PDS administrator. For example, if you have 8 licenses available, you will be able to
run 7 PDS solutions concurrently.
26
LOOP
X1
X2
X3
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
1
1
1
..
..
1
2
3
..
..
1.619379209e+000
2.237676559e-001
7.931615474e+000
...
...
2.364528435e-001
5.788049712e-001
8.278689033e-001
...
...
1.470789050e+000
1.821263115e+000
2.170793522e+000
...
...
The first line contains the solution label (the parameter Slab is set via the PDEXE command); the
second line contains the headers of the data columns - the iteration number, cycle number, loop
number, and the random variable names. The iteration number and cycle number tell the PDS to which
group (with specific PDS method and settings) the loops belong. Subsequent lines provide specific iteration, cycle, loop, and input sample values for the defined input variables.
The PDS also creates a file where the results are stored. The name of the results file is jobname_Slab.pdrs. Before the job is executed, the file looks like this:
TEST1
ITER CYCL
LOOP ERR
X1
X2
X3
RESULT
In the first line, the PDS enters the solution label. In the second line are the headers for the data columns:
the first four columns are the iteration, cycle number, loop number, and an error flag. The fifth and
subsequent columns are for the random input variable and random output parameter values. If you
run a subsequent analysis (same type of analysis with the same solution label), the cycle is incremented,
the loop count is reset to 1, and the result file is appended.
For example, the content of the result file could look like this:
TEST1
ITER CYCL
1
1
1
1
1
1
..
..
..
..
LOOP ERR
1
0
2
0
3
0
.. ..
.. ..
X1
1.619379209e+000
2.237676559e-001
7.931615474e+000
...
...
X2
2.364528435e-001
5.788049712e-001
8.278689033e-001
...
...
X3
1.470789050e+000
1.821263115e+000
2.170793522e+000
...
...
RESULT
4.162928057e+000
4.744249212e+000
1.149997825e+001
...
...
Note
Loops ending with an error are deemed "not trustworthy", i.e. if the loop lead to an error,
then the calculated results are probably wrong. Those loops will have the error flag in the
fourth column set to "1" instead of "0". Those loops will be excluded from the probabilistic
postprocessing altogether, i.e. the loops will not be used for the response surface fitting and
also the statistical analysis in connection with a Monte Carlo simulation will skip those loops.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
27
Probabilistic Design
1. Configure the remote machines that you want to use for parallel processing. See Machine Configurations (p. 29) for details. You must have account privileges on all remote hosts.
2. Configure the local machine that you want to use for managing the parallel processing. See Configuring the Master Machine (p. 33) for details.
3. Start the parallel execution.
Each of these steps is explained in detail in the following sections.
An understanding of the following terms will be useful in the discussion of the parallel analysis.
Simulation: A simulation is a set of input variables used with the analysis file to produce the output
variables
Parent Process: The parent process is the executable which manages the creation of input parameters,
communication with child processes, and postprocessing of output parameters (i.e. this is where the
session runs, which you started to perform a probabilistic design analysis). There is only one parent
process for a parallel analysis.
Child Process: A child process is a server running on a particular slave machine which processes a
simulation. There can be any number of children for a given parent process.
Master Machine: This is the machine on which you are running the parent process. If this machine
has more than one CPU you may want to also use it as a slave machine.
Slave Machine: This is the name of the machine on which you are running the child process. If the
slave machine has more than one CPUyou may want to run more than one child process.
ANSYS Nanny: This is a program in the parent process that starts a child process on a slave machine;
manages the transfer of files and data between the parent process and child processes and eventually
terminates the child processes on the various slave machines and removes temporary files.
ANSYS Thin Server: This is an application, which must be running on the slave machine, which is
responsible for tranferring files between the master machine and slave machine as well as starting a
session for each child process. This application may be started automatically (remote shell option) or
you can start it manually (connection port option).
ANSYS Server: This is the application using TCP/IP which receives commands and return information
to the parent process (client-server).
Connection Port: This is the port on which the ANSYS Thin Server accepts connections from a Master
Machine.
Communication Port: This is a port on which the ANSYS Thin Server will allow communications from
a Master Machine.
Caution
For security reasons, ANSYS, Inc. recommends that you use parallel processing within your
local network protected by a firewall. There is minimal security when using the ANSYS Thin
Server. The client that connects to the server has all of the permissions allowed by the person
or account starting the ANSYS Thin Server.
28
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
29
Probabilistic Design
On all slave machines, edit/create the .rhosts file in your home directory to include:
MasterMachine1 UserId
MasterMachine2 UserId
MasterMachine3 UserId
MasterMachine4 UserId
each on a separate line. This will allow the user with the account user identification UserId to access
the slave machine from the master machines listed in the file using the remote shell service.
ANSYS in PATH
To be able to run the parallel process the executable and ANSYS Thin Server script must be in your
path. The most effective way to determine this is to do the following.
Determining if the executable is in your path:
On the master machine, use the command rsh SlaveMachine which ansys145 (replace rsh by
remsh if master machine is an HP machine). This should return the string
/ansys145/bin
Here "" is your installation directory If you do not get this then you need to modify your PATH variable
to contain the installation path. The way this is done depends on the operating system being used,
refer to the documentation for your particular O/S environment. Here are some examples for typical
uses:
UNIX
To change the PATH requires you to modify a certain system file depending which shell you are running
under. To find out which shell you are using, issue the UNIX command
echo $shell or echo $SHELL
If the prompt is /bin/csh" or similar, then you are running under c-shell. If the prompt is/bin/ksh
or similar, then you are running under k-shell. If the prompt is /bin/tcsh or similar, then you are
running under tc-shell, which can be treated the same as c-shell.
If you are running under c-shell (or tc-shell) you need include the following line at the end of your
.cshrc file in your home directory:
set path=( .../ansys145/bin $path)
If you are running under k-shell you need include the following line at the end of your .kshrc file in
your home directory:
export PATH={ .../ansys145/bin $PATH }
If you don't have a .cshrc file or .kshrc file in your home directory, then you need to create one
and include the respective commands mentioned above.
PC
To change the PATH variable, go to Control Panel and select System on the Advanced tab. Select
Environment Variables and add the executable to the PATH under System variables.
30
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
31
Probabilistic Design
ANSYS in PATH
See Using the Remote Shell Option (p. 29) for information on how to place the executable in your path.
AnsysClients
The AnsysClients file is read by the ANSYS Thin Server on startup and must be in the directory in
which you are going to run the ANSYS Thin Server. This file must contain a list of each master machine
IP address and the communication port that the ANSYS Thin Server will accept communications from.
The master machine's IP address must be used, but may be specified with a wildcard to allow connection
from a network domain. The communication port is to be a value between 49512 and 65535. The
communication port can also be specified as a range, which allows for that many connections. It is best
to always specify a range for the communication port numbers. For example:
10.3.* 59100-59130
10.3.20.1 59200-59230
10.2.5.55 59300-59330
10.1.1.104 59400-59430
192.1.10.34 59500-59530
Make sure that the port number ranges are not overlapping and are unique for each master machine
or network domain. You should make the range at least as wide as the number of possible users connecting from a master machine.
Start the ANSYS Thin Server
Starting on a PC
On a PC open a command prompt window and go to the directory containing the AnsysClients
file.
Issue the command: ansysts145 connection port
This will start the ANSYS Thin Server using the specific connection port. The value of connection port
is the port on which the master machines will connect to the ANSYS Thin Server. The value of "connection_port" should be between 49512 and 65535. For example:
ansysts145 62000
The ANSYS Thin Server starts without any message and continues running until the command prompt
window is closed or using Ctrl-C to stop the process. The command prompt window may be minimized.
Starting on a UNIX Machine
On a UNIX machine go to the directory containing the AnsysClients file.
Issue the command: ansysts145 connection_port &
The command starts the ANSYS Thin Server using the specific connection port. The value of connection
port is the port on which the master machines will connect to the ANSYS Thin Server. The value of
"connection_port" should be between 49512 and 65535. For example:
ansysts145 62000&
32
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
33
Probabilistic Design
#
# HOST
OS
PORT CPU
MySlaveMachine SGI64 0
CPU TIME
62000
4
LocPORT
I/O Directory
15
59400-59430
/tmp/sdr/pds_runs
This will ensure that the ANSYS Thin Server on the slave machine will be contacted using the same
connection port it has been started with, i.e. 62000 in this case. Also the communication will use the
same communication port numbers the ANSYS Thin Server accepts from the machine "MyMasterMachine"
where you try to connect from.
An example hosts145.ans file looks like this:
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
This file is used to specify those hosts that the ANSYS Nanny may
run children on.
Each host entry is to be on its own line.
several fields which are space delimited.
34
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
ynn gis
oisses xe
Iretsw
PPPPT yn
WWPPSEHHPPS Pretsw
WWIISEHHIIS Iretsw
steilgsse
sFstsoh
QQQQT yn
WWPPSEHHPPS Pretsw
WWIISEHHIIS Iretsw
steilgsse
ynn gis
oisses xe
Pretsw
RRRRT yn
WWPPSEHHPPS Pretsw
WWIISEHHIIS Iretsw
steilgsse
sFstsoh
35
Probabilistic Design
If you are working in batch mode, you must create this file using a text editor. However, for sake of
simplicity it is recommended to let ANSYS create this file using interactive mode and then to proceed
with batch mode operation. This file must reside in the directory where you are running ANSYS PDS.
This file must include the following information:
Remote hosts to be used for this particular parallel run.
Number of processes that can run in parallel on each host. If a remote host has more than one CPU,
you can use all of the CPUs on the remote host. For performance reasons, we recommend leaving
one CPU for system tasks and using only N-1 CPUs for parallel processing (if N is the number of CPUs
on the remote host).
Directories in which you want the child processes to be executed. It is recommend that you use
temporary directories like "/scratch" or "/tmp" on UNIX or "C:\TEMP" on PCs. These remote directories
are cleaned up automatically after the parallel processes finish running. Make sure that there is enough
disk space in the directory for the files created by the analysis.
An example jobname.hosts file looks like this:
# This file is used to specify those hosts that the ANSYS Parent may
# run children processes on.
#
# Each host entry is to be on its own line. The host entry consists of
# several fields which are space delimited.
#
# Field 1 - host IP address or name
# Field 2 - username to use for a remote shell to the host
# Field 3 - execution key:
# 0 - Use a remote shell to start the child process
# >1024 - Use a running ANSYS thin server on the host
# which is communicating on this port number
# Field 4 - the product code sequence to attempt to start ANSYS jobs
# on the host with. This is a colon delimited list.
# Field 5 - the number of jobs to run on this host
# Field 6 - The time in minutes to check again if the host is available.
# If this is zero then the host will not be checked again.
# Field 7 - If field 3 is nonzero then this is the local port number
# or range to start the communication with the ANSYS Thin
# Server on. This is tied to the "authentication" on the
# ANSYS Thin Server.
# If field 3 is zero then this should be set to zero.
# Field 8 - directory to create the children subdirectories in
#
# Example:
#
# UNIX box that has five processors and will first attempt to
# run with ANSYS Mechanical Batch Child and then ANSYS Mechanical
# zeus wjc 0 MEBACH:ANSYS 5 30 2000 /scratch/wjc
# XP box running the thin server on port 2010
# wjcpc wjc 2010 ANSYS 1 0 2000 C:\TEMP
alpha1 epc 0 MEBA:MEBACH 1 15 2000 /scratch/epc
athena epc 0 MEBA:MEBACH 1 15 2000 /scratch/epc
rs43p pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
rs260 pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
snoopy pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
hp160 pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
alpha24 pdstest 0 MEBA:MEBACH 2 15 2000 /home/pdstest
hp770 pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
us60 pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
ss60 pdstest 0 MEBA:MEBACH 1 15 2000 /home/pdstest
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
37
Probabilistic Design
UP - The Thin Server is up.
RUN - The ANSYS simulation is running on that host.
If you want to check diagnostics for a status code in parallel processing mode, select Diagnostics in
the Server Process Control dialog. The diagnostics system shows you the status of the slave machine(s),
and details to help you address any issues that arise. Suggested actions are presented.
Note
In batch mode, ANSYS will always attempt to start remote processing on all machines listed
in the Jobname.hosts file.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
mreht
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
39
Probabilistic Design
A full regression model uses the available sampling points to determine values for all regression coefficients c0 to c3. Of course the values for c0 to c2 will be zero or very close to zero; taking more coefficients
into account than really necessary reduces the degrees of freedom of the algebraic equation to be
solved to evaluate the coefficients. This in turn reduces the accuracy of the coefficients that are important
for the regression fit. The forward-stepwise regression analysis takes this into account and automatically
eliminates terms that are not needed.
The ANSYS PDS offers a variety of transformation functions that can be used to make the random response parameter to be more appropriately described by a quadratic function after the transformation
has been applied. These transformation functions can be found in Transformation of Random Output
Parameter Values for Regression Fitting in the Mechanical APDL Theory Reference.
i
B
is the
Here, yi is the value of a random output parameter obtained in the i-th sampling loop and
corresponding transformed value. The physical nature of the problem should indicate which transformation to use; for example, lifetime parameters (such as the number of cycles until low cycle fatigue occurs)
are usually transformed with a logarithmic transformation. If you do not have this kind of information,
then you should start with the Box-Cox transformation. The PDS automatically searches for an optimum
value for the Box-Cox parameter within the interval (-2,2). As guidelines:
If is close to -1.0 then the data is best transformed by a reciprocal transformation, which is a power
transformation with an exponent of -1.0.
If is close to zero then the data is best transformed by a logarithmic transformation.
If is close to 0.5 then use the square root transformation.
If is close to 1.0, then no transformation should be used.
If is not close to any of these specific values then the Box-Cox transformation is appropriate.
To fit a response surface:
Command(s): RSFIT
GUI: Main Menu> Prob Design> Response Surf> Fit Resp Surf
40
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
41
Probabilistic Design
Cumulative Distribution Function
Probabilities
Inverse Probabilities
Trends
Scatter Plot
Sensitivities
Correlation Matrix
Report
Print HTML Report
Plot Histogram
Use the PDHIST command to request a histogram plot of a design variable.
Command(s): PDHIST
GUI: Main Menu> Prob Design> Prob Results> Statistics> Histogram
You must specify the results set you want to use, the design variable you want to review, the number
of classes/points to use, and the type of histogram.
CumulativeDF
Use the PDCDF command to request a histogram plot of a design variable.
Command(s): PDCDF
GUI: Main Menu> Prob Design> Prob Results> Statistics> CumulativeDF
You must specify the results set you want to use, the design variable you want to review, and the
confidence level.
The confidence level is a probability expressing the confidence that the values for the cumulative distribution function are in fact between the confidence bounds. The larger the confidence level, the wider
the confidence bounds. Plotting of the confidence bounds only makes sense for the postprocessing of
Monte Carlo simulation results. Here, the confidence bounds represent the accuracy of the results and
with increasing sample size the width of the confidence bounds gets smaller for the same confidence
level. For response surface methods the number of simulations done on the response surface is usually
42
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Probabilities
Use the PDPROB command to request the value of a design variable at a specific point on the cumulative distribution curve.
Command(s): PDPROB
GUI: Main Menu> Prob Design> Prob Results> Statistics> Probabilities
You must specify the results set you want to use, the design variable you want to review, the relation
(greater than, less than), the limit value, and the confidence level.
Inverse Probabilities
Use the PDPINV command to request the value of a design variable at a specific point on the cumulative
distribution curve.
Command(s): PDPINV
GUI: Main Menu> Prob Design> Prob Results> Statistics> Inverse Prob
You must specify the results set you want to use, the design variable you want to review, the relation
(greater than, less than), the limit value, and the confidence level.
Scatter Plot
Use the PDSCAT command to request a scatter plot showing the correlation between two design
variables.
Command(s): PDSCAT
GUI: Main Menu> Prob Design> Prob Results> Trends> Scatter Plot
You must select the results set that you want to use, the design variables that you want to review, the
type of trendline curve to use (and if plotted, the polynomial order), and the maximum number of point
to include in the scatter plot.
Sensitivities
Use the PDSENS command to request the sensitivities of an output parameter to the input variables.
Command(s): PDSENS
GUI: Main Menu> Prob Design> Prob Results> Trends> Sensitivities
You must specify the results set and output parameter you want to use, the type of chart to plot, the
type of correlation coefficient, and the sensitivity level.
Correlation Matrix
Use the PDCMAT command to calculate the correlation coefficient matrix.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
43
Probabilistic Design
Command(s): PDCMAT
GUI: Main Menu> Prob Design> Prob Results> Statistics> Probabilities
You must specify the results set you want to use, which type of design variables you are looking at, the
specific design variable names, the type of correlation, the significance level, and whether you want to
see the probabilities with the correlation coefficients.
Report Options
Use the PDROPT command to request an HTML report.
Command(s): PDROPT
GUI: Main Menu> Prob Design> Prob Results> Report> Report Options
You must specify which statistics and trends to show in the report and in what form you want to see
them. See the PDROPT command for details, and see the other probabilistic results options for further
details.
Generate Report
Use the PDWRITE command to request the sensitivities of an output parameter to the input variables.
Command(s): PDWRITE
GUI: Main Menu> Prob Design> Prob Results> Report> Generate Report
You must enter a name for the report file, your first and last name, and whether links should be generated
between your report and the analysis file, each analysis loop, and the response surface output parameter
details (if the response surface method was used).
44
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The values and hints given below are simply guidelines; none are absolute. Always
verify this information with an expert in your organization and adapt it as needed to
suit your analysis.
45
Probabilistic Design
and maximum limit for a uniform distribution. In a mass production field, you probably have a lot of
data, in which case you could use a commercial statistical package that will allow you to actually fit a
statistical distribution function that best describes the scatter of the data.
PG
99.5%
2.5758
99.0%
2.3263
97.5%
1.9600
95.0%
1.6449
90.0%
1.2816
1.4.1.3.3. No Data
In situations where no information is available, there is never just one right answer. Below are hints
about which physical quantities are usually described in terms of which distribution functions. This
might help you with the particular physical quantity you have in mind. Also below is a list of which
distribution functions are usually used for which kind of phenomena. Keep in mind that you might need
to choose from multiple options.
46
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Geometric Tolerances
If you are designing a prototype, you could assume that the actual dimensions of the manufactured
parts would be somewhere within the manufacturing tolerances. In this case it is reasonable to use
a uniform distribution, where the tolerance bounds provide the lower and upper limits of the distribution function.
Sometimes the manufacturing process generates a skewed distribution; for example, one half of the
tolerance band is more likely to be hit than the other half. This is often the case if missing half of the
tolerance band means that rework is necessary, while falling outside the tolerance band on the other
side would lead to the part being scrapped. In this case a Beta distribution is more appropriate.
Often a Gaussian distribution is used. The fact that the normal distribution has no bounds (it spans
minus infinity to infinity), is theoretically a severe violation of the fact that geometrical extensions
are described by finite positive numbers only. However, in practice this is irrelevant if the standard
deviation is very small compared to the value of the geometric extension, as is typically true for
geometric tolerances.
Material Data
Very often the scatter of material data is described by a Gaussian distribution.
In some cases the material strength of a part is governed by the "weakest-link-theory". The "weakestlink-theory" assumes that the entire part would fail whenever its weakest spot would fail. for material
properties where the "weakest-link" assumptions are valid, then the Weibull distribution might be
applicable.
For some cases, it is acceptable to use the scatter information from a similar material type. Let's assume
that you know that a material type very similar to the one you are using has a certain material
property with a Gaussian distribution and a standard deviation of 5% around the measured mean
value; then let's assume that for the material type you are using, you only know its mean value. In
this case, you could consider using a Gaussian distribution with a standard deviation of 5% around
the given mean value.
For temperature-dependent materials it is prudent to describe the randomness by separating the
temperature dependency from the scatter effect. In this case you need the mean values of your material property as a function of temperature in the same way that you need this information to perform
a deterministic analysis. If M(T) denotes an arbitrary temperature dependent material property then
the following approaches are commonly used:
Multiplication equation:
M(T)rand = Crand
(T)
Additive equation:
M(T)rand =
(T) + Mrand
Linear equation:
M(T)rand = Crand
(T) + Mrand
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
47
Probabilistic Design
Here, (T) denotes the mean value of the material property as a function of temperature. In
the "multiplication equation" the mean value function is scaled with a coefficient Crand and this
coefficient is a random variable describing the scatter of the material property. In the "additive
equation" a random variable Mrand is added on top of the mean value function (T). The
"linear equation" combines both approaches and here both Crand and Mrand are random variables.
However, you should take into account that in general for the "linear equation" approach Crand
and Mrand are, correlated.
Deciding which of these approaches is most suitable to describing the scatter of the temperature dependent material property requires that you have some raw data about this material property. Only by
reviewing the raw data and plotting it versus temperature you can tell which approach is the better
one.
Load Data
For loads, you usually only have a nominal or average value. You could ask the person who provided
the nominal value the following questions: If we have 1000 components that are operated under real
life conditions, what would the lowest load value be that only one of these 1000 components is
subjected to and all others have a higher load? What would the most likely load value be, i.e. the
value that most of these 1000 components have (or are very close to)? What would the highest load
value be that only one of the 1000 components is subjected to and all others have a lower load? To
be safe you should ask these questions not only of the person who provided the nominal value, but
also to one or more experts who are familiar with how your products are operated under real-life
conditions. From all the answers you get, you can then consolidate what the minimum, the most
likely, and the maximum value probably is. As verification you can compare this picture with the
nominal value that you would use for a deterministic analysis. If the nominal value does not have a
conservative bias to it then it should be close to the most likely value. If the nominal value includes
a conservative assumption (is biased), then its value is probably close to the maximum value. Finally,
you can use a triangular distribution using the minimum, most likely, and maximum values obtained.
If the load parameter is generated by a computer program then the more accurate procedure is to
consider a probabilistic analysis using this computer program as the solver mechanism. Use a probabilistic design technique on that computer program to assess what the scatter of the output parameters are, and apply that data as input to a subsequent analysis. In other words, first run a probabilistic analysis to generate an output range, and then use that output range as input for a subsequent
probabilistic analysis.
Here, you have to distinguish if the program that generates the loads is ANSYS itself or your
own in-house program. If you have used ANSYS to generate the loads (for example, FLOTRAN
analysis calculating fluid loads on a structure or a thermal analysis calculating the thermal loads
of a structure) then we highly recommend that you include these load calculation steps in the
analysis file (and therefore in the probabilistic analysis). In this case you also need to model the
input parameters of these load calculation steps as random input variables. If you have used
your own in-house program to generate the loads, you can still integrate the load calculation
program in the analysis file (see the /SYS command for details), but you must have an interface
between that program and ANSYS that allows the programs to communicate with each other
and thus automatically transfer data.
You also have to distinguish if the load values are random fields or single random variables. If the load
is different from node to node (element to element) then it is most appropriate to include the program
calculating the load in the analysis file. If the load is described by one or very few constant values then
48
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Exponential Distribution
The exponential distribution is useful in cases where there is a physical reason that the probability
density function is strictly decreasing as the random input variable value increases. The distribution is
mostly used to describe time-related effects; for example, it describes the time between independent
events occurring at a constant rate. It is therefore very popular in the area of systems reliability and
lifetime-related systems reliability, and it can be used for the life distribution of non-redundant systems.
Typically, it is used if the lifetime is not subjected to wear-out and the failure rate is constant with time.
Wear-out is usually a dominant life-limiting factor for mechanical components, which would preclude
the use of the exponential distribution for mechanical parts. However in cases where preventive maintenance exchanges parts before wear-out can occur, then the exponential distribution is still useful to
describe the distribution of the time until exchanging the part is necessary.
Gamma Distribution
The Gamma distribution is again a more time-related distribution function. For example it describes
the distribution of the time required for exactly k events to occur under the assumption that the events
take place at a constant rate. It is also used to describe the time to failure for a system with standby
components.
Lognormal Distribution
The lognormal distribution is a basic and commonly used distribution. It is typically used to describe
the scatter of the measurement data of physical phenomena, where the logarithm of the data would
follow a normal distribution. The lognormal distribution is very suitable for phenomena that arise from
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
49
Probabilistic Design
the multiplication of a large number of error effects. It is also correct to use the lognormal distribution
for a random variable that is the result of multiplying two or more random effects (if the effects that
get multiplied are also lognormally distributed). It is often used for lifetime distributions; for example,
the scatter of the strain amplitude of a cyclic loading that a material can endure until low-cycle-fatigue
occurs is very often described by a lognormal distribution.
Uniform Distribution
The uniform distribution is a very fundamental distribution for cases where no other information apart
from a lower and an upper limit exists. It is very useful to describe geometric tolerances. It can also be
used in cases where there is no evidence that any value of the random variable is more likely than any
other within a certain interval. In this sense it can be used for cases where "lack of engineering knowledge" plays a role.
Triangular Distribution
The triangular distribution is most helpful to model a random variable when actual data is not available.
It is very often used to cast the results of expert-opinion into a mathematical form, and is often used
to describe the scatter of load parameters. However, regardless of the physical nature of the random
variable you want to model, you can always ask some experts questions like "What is the one-in-athousand minimum and maximum case for this random variable? and other similar questions. You
should also include an estimate for the random variable value derived from a computer program, as
described earlier. This is also described in more detail above for load parameters in Choosing a Distribution for a Random Variable (p. 45).
Weibull Distribution
In engineering, the Weibull distribution is most often used for strength or strength-related lifetime
parameters, and it is the standard distribution for material strength and lifetime parameters for very
brittle materials (for these very brittle material the "weakest-link-theory" is applicable). For more details
see Choosing a Distribution for a Random Variable (p. 45).
50
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
51
Probabilistic Design
they are built. One simulation loop represents one component that is subjected to a particular set of
loads and boundary conditions.
The Direct Monte Carlo Sampling technique is not the most efficient technique, but it is still widely
used and accepted, especially for benchmarking and validating probabilistic results. However, benchmarking and validating requires many simulation loops, which is not always feasible. This sampling
method is also inefficient due to the fact that the sampling process has no "memory."
For example, if we have two random input variables X1 and X2 both having a uniform distribution
ranging from 0.0 to 1.0, and we generate 15 samples, we could get a cluster of two (or even more)
sampling points that occur close to each other if we graphed the two variables (see figure below). While
in the space of all random input variables, it can happen that one sample has input values close to
another sample, this does not provide new information and insight into the behavior of a component
in a computer simulation if the same (or almost the same) samples are repeated.
Figure 1.7: Graph of X1 and X2 Showing Two Samples with Close Values
52
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
53
Probabilistic Design
in the PDEXE command for postprocessing. The PDS system does not check if the solution label in
the user-specified file and the one given in the PDEXE command match.
The second line of the file must contain the headers of the data columns. The first three column
headers must be ITER CYCL and LOOP respectively; then subsequent columns should contain
,
,
,
the names of the random variables. You must use one of the allowed separators as described above
between the column headers. No additional data is allowed on this line, and if found, will prompt an
error message. An error message is also issued if any of the required column headers are missing.
The random variable names in your file must match the names of the defined random variables. The
variable names that you specify must consist of all uppercase characters (regardless of the case used
in the defined variable names).
Columns four to n can be in arbitrary order. The ANSYS PDS tool determines the order for the random
variable data based on the order of the random variable names in the second line.
The third and subsequent lines must contain the order number for the iteration, cycle, and simulation
loop, then the random variable values for that loop. The iteration, cycle, and simulation loop numbers
must be in the first, second, and third columns respectively, followed by the random variable values.
The iteration and cycle numbers are used by the ANSYS PDS (internally) and for a user-defined
sampling method you will typically use a value of "1" for all simulation loops. The loop number is an
ascending number from 1 to the total number of loops requested. Additional data is not allowed,
and if found, will lead to an error message. An error message is also issued if any of the data columns
are missing.
You must be sure that the order of the random variable values in each line is identical to the order
of the random variable names in the second line.
The user-specified sampling file must contain a minimum of one data line for the random variable
values.
When the PDUSER command is issued, the PDS checks that the specified file exists, then verifies it for
completeness and consistency. Consistency is checked according to the possible minimum and maximum
boundaries of the distribution type of the individual random variables. An error message is issued if a
random variable value is found in the file that is below the minimum boundary or above the maximum
boundary of the distribution of that random variable. This also means that any value will be accepted
for a random variable if its distribution has no minimum or maximum boundary; for example, this is
the case for the Gaussian (normal) distribution. Apart from this check, it is your responsibility to provide
values for the random variables that are consistent with their random distribution.
Note
It is your responsibility to ensure that parameters defined as random input variables are actually input parameters for the analysis defined with the PDANL command. Likewise, you
must ensure that parameters defined as random output parameter are in fact results generated
in the analysis file.
Example
An excerpt of the content of a user-specified sampling file is given below. This example is based on
three random variables named X1, X2, and X3. A total of 100 simulation loops are requested.
USERSAMP
ITER CYCL
54
LOOP
X1
X2
X3
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
1
1
1
..
..
1
1
1
1
2
3
..
..
98
99
100
1.619379209e+000
2.237676559e-001
7.931615474e+000
...
...
1.797221666e+000
1.290815540e+001
4.699281922e+000
2.364528435e-001
5.788049712e-001
8.278689033e-001
...
...
3.029471373e-001
9.271606216e-001
6.526505821e-001
1.470789050e+000
1.821263115e+000
2.170793522e+000
...
...
1.877701081e+000
2.091047328e+000
1.901013985e+000
ji
is described by
x x
iaj
Ia i
x
i
Ia i
where c0 is the coefficient of the constant term, ci, i = 1,...NRV are the coefficients of the linear terms
and cij, i = 1,...NRV and j = i, ...,NRV are the coefficients of the quadratic terms. To evaluate these coefficients a regression analysis is used and the coefficients are usually evaluated such that the sum of
squared differences between the true simulation results and the values of the approximation function
is minimized.
Hence, a response surface analysis consists of two steps:
1. Performing the simulation loops to calculate the values of the random output parameters that correspond to the sample points in the space of random input variables.
2. Performing a regression analysis to derive the terms and the coefficients of the approximation
function.
The fundamental idea of Response Surface Methods is that once the coefficients of a suitable approximation function are found, then we can directly use the approximation function instead of looping
through the finite element model. To perform a finite element analysis might require minutes to hours
of computation time; in contrast, evaluating a quadratic function requires only a fraction of a second.
Hence, if using the approximation function, we can afford to evaluate the approximated response
parameter thousands of times.
A quadratic polynomial is sufficient in many cases of engineering analysis (for example, the evaluation
of the thermal stress mentioned above). For that evaluation, the Young's modulus and the thermal expansion coefficient both have a linear effect on the thermal stresses, which is taken into account in a
quadratic approximation by the mixed quadratic terms. However, there are cases where a quadratic
approximation is not sufficient; for example, if the finite element results are used to calculate the lifetime
of a component. For this evaluation, the lifetime typically shows an exponential behavior with respect
to the input parameters; thus the lifetime results cannot be directly or sufficiently described by a
quadratic polynomial. But often, if you apply a logarithmic transformation to the lifetime results, then
these transformed values can be approximated by a quadratic polynomial. The ANSYS PDS offers a
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
55
Probabilistic Design
variety of transformation functions that you can apply to the response parameters, and the logarithmic
transformation function is one of them.
Assuming the approximation function is suitable for your problem, the advantages of the Response
Surface Method are:
It often requires fewer simulation loops than the Monte Carlo Simulation method.
It can evaluate very low probability levels. This is something the Monte Carlo Simulation method
cannot do unless you perform a great number of simulation loops.
The goodness-of-fit parameters tell you how good the approximation function is (in other words,
how accurate the approximation function is that describes your "true" response parameter values).
The goodness-of-fit parameters can warn you if the approximation function is insufficient.
The individual simulation loops are inherently independent (the individual simulation loops do not
depend on the results of any other simulation loops). This makes Response Surface Method an ideal
candidate for parallel processing.
The disadvantages of the Response Surface Method are:
The number of required simulation loops depends on the number of random input variables. If you
have a very large number of random input variables (hundreds or even thousands), then a probabilistic analysis using Response Surface Methods would be impractical.
This method is not usually suitable for cases where a random output parameter is a non-smooth
function of the random input variables. For example, a non-smooth behavior is given if you observe
a sudden jump of the output parameter value even if the values for the random input variables vary
only slightly. This typically occurs if you have instability in your model (such as bulking). The same
might happen if the model includes a sharp nonlinearity such as a linear-elastic-ideal-plastic material
behavior. Or, if you are analyzing a contact problem, where only a slight variation in your random
input variables can change the contact situation from contact to non-contact or vice versa, then you
also might have problems using the Response Surface Method.
Note
To use Response Surface Methods, the random output parameters must be smooth and
continuous functions of the involved random input variables. Do not use Response
Surface Methods if this condition is not satisfied.
56
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The number of sample points (simulation loops) required for a central composite design as a function
of the number of random input variables is given in the table below:
Number of random input variables
Number of coefficients in
a quadratic function (with
cross-terms)
Factorial number f
Number sample
points (simulation
loops)
N/A
N/A
10
15
15
25
21
27
28
45
36
79
45
81
55
147
10
66
149
11
78
151
12
91
281
13
105
283
14
120
285
15
136
287
16
153
289
17
171
291
18
190
549
19
210
10
551
20
231
11
553
To use the Response Surface Method with a Central Composite Design, do the following:
Command(s): PDMETH,RSM,CCD PDDOEL,Name,CCD,...
GUI: Main Menu> Prob Design> Prob Method> Response Surface
PDDOEL allows you to specify design of experiment options.
See the Mechanical APDL Theory Reference for more details.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
57
Probabilistic Design
The number of sample points (simulation loops) required for a Box-Behnken design as a function of the
number of random input variables is given in the table below:
Number of random input variables
Number of coefficients in
a quadratic function (with
cross-terms)
Number sample
points (simulation
loops)
N/A
N/A
10
12
15
25
21
41
28
49
36
57
45
65
55
121
10
66
161
11
78
177
12
91
193
To use Response Surface Analysis with the Box-Behnken design option, do the following:
Command(s): PDMETH,RSM,BBM PDDOEL,Name,BBM,...
GUI: Main Menu> Prob Design> Prob Method> Response Surface
PDDOEL allows you to specify design of experiment options.
See the Mechanical APDL Theory Reference for more details.
58
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
59
Probabilistic Design
To review the simulation loop results:
Command(s): PDSHIS
GUI: Main Menu> Prob Design> Prob Results> Statistics> Sampl History
You need to select the results set label and the design variable, select a plot type, and set the confidence
level.
1.6.1.2. Histogram
A histogram plot is most commonly used to visualize the scatter of a probabilistic design variable. A
histogram is derived by dividing the range between the minimum value and the maximum value into
intervals of equal size. Then the PDS determines how many samples fall within each interval, that is,
how many "hits" landed in the intervals.
Most likely, you will use histograms to visualize the scatter of your random output parameters. The
ANSYS PDS also allows you to plot histograms of your random input variables so you can double check
that the sampling process generated the samples according to the distribution function you specified.
For random input variables, the PDS not only plots the histogram bars, but also a curve for values derived
from the distribution function you specified. Visualizing histograms of the random input variables is
another way to make sure that enough simulation loops have been performed. If the number of simulation loops is sufficient, the histogram bars will:
Be close to the curve that is derived from the distribution function
Be "smooth" (without large steps)
Not have major gaps
A major gap is given if you have no hits in an interval where neighboring intervals have many hits.
However, if the probability density function is flattening out at the far ends of a distribution (for example,
the exponential distribution flattens out for large values of the random input variable) then there might
logically be gaps. Hits are counted only as positive integer numbers and as these numbers gradually
get smaller, a zero hit can happen in an interval.
To plot histograms:
Command(s): PDHIST
GUI: Main Menu> Prob Design> Prob Results> Statistics> Histogram
60
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Ax@p
7HHI
7HW
7HV
7HU
7HT
7HS
7HR
7HQ
7HP
ix
FFF
7HI
7H
The value of the cumulative distribution function at the location x0 is the probability that the values of
X stay below x0. Whether this probability represents the failure probability or the reliability of your
component depends on how you define failure; for example, if you design a component such that a
certain deflection should not exceed a certain admissible limit then a failure event occurs if the critical
deflection exceeds this limit. Thus for this example, the cumulative distribution function is interpreted
as the reliability curve of the component. On the other hand, if you design a component such that the
eigenfrequencies are beyond a certain admissible limit then a failure event occurs if an eigenfrequency
drops below this limit. Thus for this example, the cumulative distribution function is interpreted as the
failure probability curve of the component.
The cumulative distribution function also lets you visualize what the reliability or failure probability
would be if you chose to change the admissible limits of your design.
Often you are interested in visualizing low probabilities and you want to assess the more extreme ends
of the distribution curve. In this case plotting the cumulative distribution function in one of the following
ways is more appropriate:
As a Gauss plot (also called a "normal plot"). If the probabilistic design variable follows a Gaussian
distribution then the cumulative distribution function is displayed as a straight line in this type of
plot.
As a lognormal plot. If the probabilistic design variable follows a lognormal distribution then the cumulative distribution function is displayed as a straight line in this type of plot
As a Weibull plot. If the probabilistic design variable follows a Weibull distribution then the cumulative
distribution function is displayed as a straight line in this type of plot.
The advantage of these plots is that the probability axis is scaled in a nonlinear fashion such that the
extreme ends of the distribution function are emphasized and more visible.
To plot the cumulative distribution function:
Command(s): PDCDF
GUI: Main Menu> Prob Design> Statistics> Prob Results> CumulativeDF
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
61
Probabilistic Design
1.6.2.1. Sensitivities
Probabilistic sensitivities are important in allowing you to improve your design toward a more reliable
and better quality product, or to save money while maintaining the reliability or quality of your product.
You can request a sensitivity plot for any random output parameter in your model.
There is a difference between probabilistic sensitivities and deterministic sensitivities. Deterministic
sensitivities are mostly only local gradient information. For example, to evaluate deterministic sensitivities you can vary each input parameters by 10% (one at a time) while keeping all other input parameters constant, then see how the output parameters react to these variations. As illustrated in the
figure below, an output parameter would be considered very sensitive with respect to a certain input
parameter if you observe a large change of the output parameter value.
62
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Y1
X
Y
Y2
These purely deterministic considerations have various disadvantages that are taken into consideration
for probabilistic sensitivities, namely:
A deterministic variation of an input parameter that is used to determine the gradient usually does
not take the physical range of variability into account. An input parameter varied by 10% is not
meaningful for the analysis if 10% is too much or too little compared with the actual range of
physical variability and randomness. In a probabilistic approach the physical range of variability is
inherently considered because of the distribution functions for input parameters. Probabilistic sensitivities measure how much the range of scatter of an output parameter is influenced by the scatter
of the random input variables. Hence, two effects have an influence on probabilistic sensitivities: the
slope of the gradient, plus the width of the scatter range of the random input variables. This is illustrated in the figures below. If a random input variable has a certain given range of scatter, then the
scatter of the corresponding random output parameter is larger, and the larger the slope of the
output parameter curve is (first illustration). But remember that an output parameter with a moderate
slope can have a significant scatter if the random input variables have a wider range of scatter (second
illustration).
Figure 1.13: Range of Scatter
range of
scatter X
Y1
range of
scatter X
range of
scatter Y
scatter
range of
Y2
scatter Y
Y2
range
Y2
X
Gradient information is local information only. It does not take into account that the output parameter
may react more or less with respect to variation of input parameters at other locations in the input
parameter space. However, the probabilistic approach not only takes the slope at a particular location
into account, but also all the values the random output parameter can have within the space of the
random input variables.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
63
Probabilistic Design
Deterministic sensitivities are typically evaluated using a finite-differencing scheme (varying one
parameter at a time while keeping all others fixed). This neglects the effect of interactions between
input parameters. An interaction between input parameters exists if the variation of a certain parameter has a greater or lesser effect if at the same time one or more other input parameters change
their values as well. In some cases interactions play an important or even dominant role. This is the
case if an input parameter is not significant on its own but only in connection with at least one other
input parameter. Generally, interactions play an important role in 10% to 15% of typical engineering
analysis cases (this figure is problem dependent). If interactions are important, then a deterministic
sensitivity analysis can give you completely incorrect results. However, in a probabilistic approach,
the results are always based on Monte Carlo simulations, either directly performed using you analysis
model or using response surface equations. Inherently, Monte Carlo simulations always vary all random
input variables at the same time; thus if interactions exist then they will always be correctly reflected
in the probabilistic sensitivities.
To display sensitivities, the PDS first groups the random input variables into two groups: those having
a significant influence on a particular random output parameter and those that are rather insignificant,
based on a statistical significance test. This tests the hypothesis that the sensitivity of a particular random
input variable is identical to zero and then calculates the probability that this hypothesis is true. If the
probability exceeds a certain significance level (determining that the hypothesis is likely to be true),
then the sensitivity of that random input variable is negligible. The PDS will plot only the sensitivities
of the random input variables that are found to be significant. However, insignificant sensitivities are
printed in the output window. You can also review the significance probabilities used by the hypothesis
test to decide which group a particular random input variable belonged to.
The PDS allows you to visualize sensitivities either as a bar chart, a pie chart, or both. Sensitivities are
ranked so the random input variable having the highest sensitivity appears first.
In a bar chart the most important random input variable (with the highest sensitivity) appears in the
leftmost position and the others follow to the right in the order of their importance. A bar chart describes
the sensitivities in an absolute fashion (taking the signs into account); a positive sensitivity indicates
that increasing the value of the random input variable increases the value of the random output parameter for which the sensitivities are plotted. Likewise, a negative sensitivity indicates that increasing
the random input variable value reduces the random output parameter value. In a pie chart, sensitivities
are relative to each other.
In a pie chart the most important random input variable (with the highest sensitivity) will appear first
after the 12 o'clock position, and the others follow in clockwise direction in the order of their importance.
Using a sensitivity plot, you can answer the following important questions.
How can I make the component more reliable or improve its quality?
If the results for the reliability or failure probability of the component do not reach the expected levels,
or if the scatter of an output parameter is too wide and therefore not robust enough for a quality
product, then you should make changes to the important input variables first. Modifying an input
variable that is insignificant would be waste of time.
Of course you are not in control of all random input parameters. A typical example where you have
very limited means of control are material properties. For example, if it turns out that the environmental
temperature (outdoor) is the most important input parameter then there is probably nothing you can
do. However, even if you find out that the reliability or quality of your product is driven by parameters
that you cannot control, this has importance it is likely that you have a fundamental flaw in your
product design! You should watch for influential parameters like these.
64
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
How can I save money without sacrificing the reliability or the quality of the product?
If the results for the reliability or failure probability of the component are acceptable or if the scatter
of an output parameter is small and therefore robust enough for a quality product then there is usually
the question of how to save money without reducing the reliability or quality. In this case, you should
first make changes to the input variables that turned out to be insignificant, because they do not effect
the reliability or quality of your product. If it is the geometrical properties or tolerances that are insignificant, you can consider applying a less expensive manufacturing process. If a material property turns
out to be insignificant, then this is not typically a good way to save money, because you are usually
not in control of individual material properties. However, the loads or boundary conditions can be a
potential for saving money, but in which sense this can be exploited is highly problem dependent.
Command(s): PDSENS
GUI: Main Menu> Prob Design> Prob Results> Trends> Sensitivities
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
65
Probabilistic Design
The effect of reducing and shifting the scatter of a random input variable is illustrated in the figures
below. "Input range before" denotes the scatter range of the random input variable before the reduction
or shifting, and "input range after" illustrates how the scatter range of the random input variable has
been modified. In both cases, the trendline tells how much the scatter of the output parameter is affected
and in which way if the range of scatter of the random input variable is modified.
Output range
after
Input range
after
Output
range
after
Output range
before
Input range
before
It depends on your particular problem if either reducing or shifting the range of scatter of a random
input variable is preferable. In general, reducing the range of scatter of a random input variable leads
to higher costs. A reduction of the scatter range requires a more accurate process in manufacturing or
operating the product the more accurate, the more expensive it is. This might lead you to conclude
that shifting the scatter range is a better idea, because it preserves the width of the scatter (which
means you can still use the manufacturing or operation process that you have). Below are some considerations if you want to do that:
Shifting the scatter range of a random input variable can only lead to a reduction of the scatter of a
random output parameter if the trendline shows a clear nonlinearity. If the trendline indicates a linear
trend (if it is a straight line), then shifting the range of the input variables anywhere along this straight
line doesn't make any difference. For this, reducing the scatter range of the random input variable
remains your only option.
It is obvious from the second illustration above that shifting the range of scatter of the random input
variable involves an extrapolation beyond the range where you have data. Extrapolation is always
difficult and even dangerous if done without care. The more sampling points the trendline is based
on the better you can extrapolate. Generally, you should not go more than 30-40% outside of the
range of your data. But the advantage of focusing on the important random input variables is that
a slight and careful modification can make a difference.
Note
ANSYS strongly recommends that you redo the entire probabilistic analysis using the new
and modified random input variables if you have reduced of shifted the scatter range of any
input variables using the procedure and recommendations above. To redo the probabilistic
analysis, save the probabilistic model using the PDSAVE command and clear the current
probabilistic analysis results using the PDCLR,POST command. Then you can modify the
random input variable definitions and redo the probabilistic analysis.
66
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
67
Probabilistic Design
If you perform all executions within the same session (or within the same batch input stream), the
procedure is very straightforward. After an execution, redefine the probabilistic design input as desired
and initiate the next execution. To initiate the execution:
Command(s): PDEXE
GUI: Main Menu> Prob Design> Run> Exec Serial> Run Serial
Main Menu> Prob Design> Run> Exec Parallel> Run Parallel
For more information on these execution options, see Execute Probabilistic Analysis Simulation
Loops (p. 25). If you left the program after performing an analysis and would like to continue your
probabilistic design analysis at some later time, you can do a save and restart as described next.
Caution
Currently, you must use the Save and Resume commands on the Prob Design menu to save
your work if interrupted or if you simply want to resume the project later.
68
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
In addition to probabilistic design data, the jobname is saved to the probabilistic design
database file (Jobname.PDS). Therefore, when a probabilistic design data file is resumed
(PDRESU), the jobname saved in that file will overwrite the current jobname (/FILNAME).
You can use the PDRESU command (Main Menu> Prob Design> Prob Database> Resume) in an interactive session to resume probabilistic design data that was created through a batch run, thus allowing
convenient interactive viewing of batch probabilistic design results.
69
Probabilistic Design
Young's modulus (E) = 30E6 psi
The following geometric properties are used for this analysis. These properties are the initial cross-sectional areas for each truss member:
ARE1 = 5.0
ARE2 = 5.0
ARE3 = 5.0
70
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
!
!
!
!
!
After the analysis file has been created, you can proceed with the probabilistic design analysis. You can
do this though commands or though the GUI. If you prefer using commands, the following input sets
up the probabilistic analysis for the three-bar truss example.
/inp,pds3bar,pdan
/com
/com, Enter PDS and specify the analysis file
/com
/PDS
! enter probabilistic design system
pdanl,pds3bar,pdan
/com
/com, Declare random input variables
/com
PDVAR,ARE1,GAUS,5,0.5
! define area1 with Gaussian distribution
! having mean of 5 and std. dev of 0.5
PDVAR,ARE2,tria,10,11,12
! define area2 with triangular distribution
! having low bound of 10, most likely point of 11
! and high bound of 12
PDVAR,ARE3,unif,5,6
! define area3 with uniform distribution
! with low bound of 5 and high bound of 6
/com
/com, Specify any correlations between the random variables
/com
PDCOR,ARE1,ARE3,0.25
! define a correlation coef of 0.25 between ARE1 and ARE3
/com
/com, Declare random output variables
/com
PDVAR,SIG1,resp
! define SIG1 a response parameter
PDVAR,SIG2,resp
! define SIG2 a response parameter
PDVAR,SIG3,resp
! define SIG3 a response parameter
PDVAR,TVOL,resp
! define TVOL a response parameter
/com
/com, Select the probabilistic design tool or method
/com
PDMETH,MCS,DIR
! specify direct Monte Carlo simulation
PDDMCS,100,NONE,ALL,,,,123457 ! use all 100 samples, initial seed of 123457
/com
/com, Execute the loops required for the probabilistic design analysis
/com
PDEXE,mcs3bar
! run analysis and define solution label 3bar_mcs
/com
/com, Review the results of the probabilistic analysis
/com
PDSENS,MCS3BAR,TVOL,BOTH,RANK,0.025 !Create Sensitivity plot
fini
/exit
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
71
Probabilistic Design
Step 2: Enter the Probabilistic Design Module and Identify Analysis File
First, enter the optimizer and identify the analysis file.
1. Select menu path Main Menu> Prob Design> Analysis File> Assign. The Assign Analysis File dialog
box appears.
2. In the Files list, click once on pds3bar.pdan and then clickOK.
72
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
73
Probabilistic Design
1. Select menu path Main Menu> Prob Design> Prob Results> Trends> Sensitivities. The Sensitivity of
a Response Parameter dialog box appears.
2. Select MCS3BAR in the Select Results Set field. Select SIG1 in the Select Response Parameter field.
Click OK.
3. A Rank-Order Correlation Sensitivity bar and pie chart is shown in the Graphics window.
Step 6: Exit
1. Click on Quit in the Toolbar.
2. Select an option to save, then click OK.
74
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
75
Variational Technology
Note
If the loading, material properties, or real constants defined with a data table (BF ,TBFIELD
command, or COMBIN14) are frequency-dependent and too few data points exist in the
table, the polynomial interpolation of the values will be poor. In case of abrupt variations
for a given frequency f2, ANSYS advises that you split the frequency range [f1 f3] into two
ranges (e.g., [f1 f2] ; [f2 f3] ) and perform separate analyses on each range.
High-frequency electromagnetic problems use ANSYS high-frequency elements HF119 or HF120 only.
The following Variational Technology topics are available for harmonic sweep:
2.1.1. Structural Elements Supporting Frequency-Dependent Properties
2.1.2. Harmonic Sweep for High-Frequency Electromagnetic Problems
2.1.3. Harmonic Sweep for Structural Analysis with Frequency-Dependent Material Properties
Element
PLANE182
EX, SDAMP
PLANE183
EX, SDAMP
SOLID185
EX, SDAMP
SOLID186
EX, SDAMP
SOLID187
EX, SDAMP
SOLID272
EX, SDAMP
SOLID273
EX, SDAMP
SOLID285
EX, SDAMP
76
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Define the port regions (flags) and boundary conditions. Use the
BFA, BFL, or BF commands for interior waveguide.
SF or SFA commands for exterior waveguide ports.
2.
Define the waveguide port type and excitation using the HFPORT command.
3.
4.
2.
Define a port excitation using BF, BFL, or BFA and a current density source. Offset this excitation from
the transmission line port plane.
3.
Define a voltage path from the central conductor to the ground on the transmission line port plane
using the PATH command. Save the path using the PASAVE command.
4.
5.
6.
Note
A large number of frequency steps will produce a very large database file.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
77
Variational Technology
78
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
|S21|
<|S21|
0.945 -52.364
0.935 -107.633
0.932 197.701
0.921 145.595
0.905
95.776
|S12|
<|S12|
0.947 -53.791
0.948 -106.050
0.959 200.509
0.966 146.452
0.964
93.750
|S22|
<|S22|
0.326 -69.987
0.337 -141.146
0.323 141.352
0.332
66.997
0.368
7.834
/batch,list
/title, Dielectric Post in a Rectangular Waveguide
/com, Waveguide Dimension: 22.86 x 10.16 mm^2 (Cutoff Frequency: 6.56 GHz)
/com, Dielectric Post: 12 x 10.16 x 6 mm^3 at the center of waveguide
/com,
epsr = 8.2
/com, Frequency Range: 8 - 12 GHz for TE10 mode
/com, PML Parameter: 5 Layers with -50 dB
/com, Mesh Size: Free Space Wavelength/15 at 12 GHz
/com, Numerical Model: IMPD Driven Port; PML Output Port
/com, Solution Target: S11 over frequency range using series expansion method
/prep7
ch=10.16e-3
! waveguide height
cw=22.86e-3
! waveguide width
c=12e-3
! post width
d=6e-3
! post length
epsr=8.2
! Dielectric constant
freq=10.e9
! Analysis frequency
fmesh=12e9
! Mesh frequency
lamda=3.e8/fmesh ! wavelength
h1=lamda/8
! mesh parameter
cl=5*d
! waveguide length
et,1,HF119,1
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
79
Variational Technology
et,2,HF119,1,,,1 ! PML option
mp,murx,1,1.
mp,perx,1,1.
mp,murx,2,1.
mp,perx,2,epsr ! Dielectric
block,-cw/2,cw/2,-ch/2,ch/2,-cl/2,cl/2
block,-c/2,c/2,-ch/2,ch/2,-d/2,d/2
vsbv,1,2,,delete,keep
block,-cw/2,cw/2,-ch/2,ch/2,cl/2,cl/2+lamda/5
vglue,all
esize,h1
type,1
mat,2
vmesh,2
mat,1
vmesh,4
type,2
vmesh,1
! Dielectric region
! Air region
! PML region
80
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
! STARTING FREQUENCY
! ENDING FREQUENCY
! NUMBER OF SUBSTEP
* ==============================================
*
*
Frequency dependent material properties
*
* ==============================================
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
81
Variational Technology
TB,ELASTIC,1,,2
MP,DENS,1,.10
BLOCK,0,10,0,2,0,2
LSEL,S,LOC,X,-.5,0.5
LESIZE,ALL,,,2
LSEL,S,LOC,X,9.5,10.5
LESIZE,ALL,,,2
LSEL,S,LOC,X,2,8
LESIZE,ALL,,,5
VMESH,ALL
FINISH
/SOL
LSEL,S,LOC,X,-.5,0.5
DL,all,,all
KSEL,S,LOC,X,8,12
KSEL,R,LOC,Y,-.5,.5
FK,ALL,FY,1000
ASEL,S,LOC,Y,1.8,2.2
SFA,ALL,,PRES,1000,
allsel
FINI
/com, * ==============================================
/com, *
/com, *
VT Harmonic Analysis
/com, *
/com, * ==============================================
/VT
VTCLR,ALL
VTRFIL
VTFREQ,frq,LF,UF,N
VTRSLT,disp,NODE,U,ALL,0.01,ALL
FINISH
/SOLU
ANTY,HARM
HROUT,OFF
! Print complex displacements as amplitude and phase angle
KBC,1
HROPT,VT,,,,HYST
Solve
FINISH
/show,post
/post26
nsol,10,57,u,y,d1
prvar,10
The graph below shows the difference in the Y displacement at the end of the beam, taking constant
material properties for 0Hz, constant properties for 500 Hz, and variable properties. Note that the first
peak is higher on the variable curve than the other two. This is because of the frequency step used to
create the frequency response.
82
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
83
84
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
85
Adaptive Meshing
2. Model your system using solid modeling procedures, creating meshable areas or volumes describing the
geometry of your system. It is not necessary to specify element sizes, nor do you need to mesh these
areas and volumes; the ADAPT macro initiates meshing automatically. (If you need to mesh your model
with both area and volume elements, create an ADAPTMSH.MAC subroutine.
3. You can either proceed to SOLUTION (/SOLU or menu path Main Menu> Solution) or remain in PREP7
to specify analysis type, analysis options, loads, and load step options. Apply only solid model loads and
inertia loads (accelerations, rotational accelerations, and rotational velocities) in a single load step. (Finite
element loads, coupling, and constraint equations can be introduced through the ADAPTBC.MAC subroutine. Multiple load steps can be introduced through the ADAPTSOL.MAC subroutine.)
4. If in /PREP7, exit the preprocessor (FINISH). (You can invoke the ADAPT macro from either SOLUTION
or the Begin level.)
5. Invoke the adaptive solution. To do so, use one of these methods:
Command(s): ADAPT
GUI: Main Menu> Solution> Solve> Adaptive Mesh
You can use the ADAPT macro in either a thermal or a structural analysis, but you cannot combine
the two disciplines in one adaptive solution. As the adaptive meshing iterations proceed, element
sizes are adjusted (within the limits set by FACMN and FACMX) to decrease and increase the elemental
error energies until the error in energy norm matches the target value (or until the specified maximum
number of solutions has been used).
After you have invoked the adaptive solution, this macro controls all program operations until the
solution is completed. The ADAPT macro will define element sizes, generate the mesh, solve, evaluate errors, and iterate as necessary till the target value of error in energy norm is met. All these
steps are performed automatically, with no further input required from you.
6. When adaptive meshing has converged, the program automatically turns element shape checking on
(SHPP,ON). It then returns to the SOLUTION phase or to the Begin level, depending on which phase you
were in when you invoked ADAPT. You can then enter /POST1 and postprocess as desired, using
standard techniques.
86
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Selective adaptivity can improve the performance of models having concentrated loads.
If you select a set of keypoints, the ADAPT macro will still include all your keypoints (that is, the ADAPT
macro will modify the mesh at both selected and non-selected keypoints), unless you also set KYKPS
= 1 in the ADAPT command (Main Menu> Solution> Solve> Adaptive Mesh).
If you select a set of areas or volumes, the ADAPT macro will adjust element sizes only in those regions
that are in the selected set. You will have to mesh your entire model in PREP7 before executing ADAPT.
87
Adaptive Meshing
mesh only areas. If you desire to mesh both volumes and areas, you can create a user subroutine, ADAPTMSH.MAC, to perform all the desired operations. You will need to clear any specially-meshed entities
before remeshing. Such a subroutine might look like this:
C*** Subroutine ADAPTMSH.MAC - Your name - Job Name - Date Created
TYPE,1
! Set element TYPE attribute for area meshing
ACLEAR,3,5,2
! Clear areas and volumes to be meshed by this subroutine
VCLEAR,ALL
AMESH,3,5,2
! Mesh areas 3 and 5 (no other areas will be meshed by ADAPT)
TYPE,2
! Change element type for volume mesh
VMESH,ALL
! Mesh all volumes
Please see the TYPE, ACLEAR, VCLEAR, AMESH, and VMESH command descriptions for more information.
We strongly recommend that you include a comment line (C***) to identify your macro uniquely. This
comment line will be echoed in the job printout, and will provide assurance that the ADAPT macro
has used the correct user subroutine.
This default sequence will solve only a single load step. You might be able to implement other solution
sequences by incorporating them into the user subroutine ADAPTSOL.MAC.
88
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
89
Adaptive Meshing
Mapped meshing (all hexahedral bricks) is available for 3-D solid analyses. Mapped volume meshing,
if possible for a given model, will probably give superior performance, compared to tetrahedral
meshing.
In general, midside-node elements will perform better than linear elements in an adaptive meshing
operation.
Do not use concentrated loads or other singularities (such as sharp re-entrant corners), because the
ADAPT procedure cannot show energy value convergence with mesh refinement in the vicinity of
these singularities. If a concentrated loading condition is present in your model, replace it with an
equivalent pressure load, applied over a small area. (Or, exclude the region of the concentrated load
from adaptive modification using the select options discussed previously.)
For many problems, it could be preferable to use a number of relatively small, simple regions in place
of fewer large, complicated regions, for best meshing performance.
If the location of maximum response is known or suspected beforehand, then a keypoint should be
placed near that location.
If the ADAPT procedure is used in an interactive run and the ANSYS program aborts abruptly without
issuing the proper error message, then the output file for the adaptive meshing portion of your run
(Jobname.ADPT) should be reviewed to determine the cause.
Similarly, if the ADAPT procedure is used in a batch run, then Jobname.ADPT should be saved and
examined in case of an abrupt abort.
If a model has an excessive curvature in some region, then your model might experience mesh failure.
In this case, use the SIZE field on the KESIZE command (Main Menu> Preprocessor> Meshing>
Size Cntrls> ManualSize> Keypoints> Picked KPs) to define the maximum element edge length
settings at keypoints near the excessive curvature region. FACMX should also be set equal to 1 (in
the ADAPT command) so that element size will not be permitted to grow in the vicinity of the excessive curvature.
You should save the results file (Jobname.RST or Jobname.RTH). In case of a program abort in
the middle of the ADAPT procedure, the results file will still have the entire analysis data from the
previous solution completed by the ADAPT procedure.
You should issue a SAVE (Utility Menu> File> Save as Jobname.db) before starting the adaptive
meshing operation. In case of system abort, Jobname.DB can then be used to restart the analysis.
90
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
91
92
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 4: Rezoning
In a finite large-deformation analysis, mesh distortion reduces simulation accuracy, causes convergence
difficulties, and can eventually terminate an analysis. Rezoning allows you to repair the distorted mesh
and continue the simulation.
The program's manual rezoning capability allows you to decide when to use rezoning, what region(s)
to rezone, and what remeshing method to use on the selected region(s).
During the rezoning process, the program updates the database as necessary, generates contact elements
if needed, transfers boundary conditions and loads from the original mesh, and maps all solved variables
(node and element solutions) to the new mesh automatically. Afterwards, with equilibrium achieved
based on the mapped variables, you can continue solving using the new mesh.
The following rezoning topics are available:
4.1. Benefits and Limitations of Rezoning
4.2. Rezoning Requirements
4.3. Understanding the Rezoning Process
4.4. Step 1: Determine the Substep to Initiate Rezoning
4.5. Step 2. Initiate Rezoning
4.6. Step 3: Select a Region to Remesh
4.7. Step 4: Perform the Remeshing Operation
4.8. Step 5: Verify Applied Contact Boundaries, Loads, and Boundary Conditions
4.9. Step 6: Automatically Map Variables and Balance Residuals
4.10. Step 7: Perform a Multiframe Restart
4.11. Repeating the Rezoning Process if Necessary
4.12. Postprocessing Rezoning Results
4.13. Rezoning Restrictions
4.14. Rezoning Examples
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
93
Rezoning
The simulation terminates at TIME = 0.44. Rezoning begins on the deformed mesh at substep 7 (TIME
= 0.40):
94
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Based on the new mesh, the simulation concludes successfully at TIME = 1.0:
95
Rezoning
difficult to determine when material is unstable, you can check the strain values, stress
states, and convergence patterns. A sudden convergence difficulty could mean that
material is no longer stable. The program also issues a warning at the beginning of the
solution indicating when hyperelastic material could be unstable, although such a
warning is very preliminary and applies only to cases involving simple stress states.
Unstable Structures
For some geometries and loads, a deformation may cause a "snap-through," or local
buckling. Such behavior can also manifest itself as a mesh distortion, but one that
rezoning cannot repair. The effect is usually easy to detect by closely checking the deformed region or the load-versus- time (displacement) curve.
Numerical Instabilities
A condition of numerical instability can occur when a problem is nearly overconstrained.
The constraints can include kinematic constraints such as applied displacements, CP,
and CE, and volumetric constraints introduced by fully incompressible material in mixed
u-P elements. In many cases, numerical instability is apparent even in the early stages
of an analysis.
For a successful rezoning, the new mesh must be of a higher quality than the old mesh. If the new
mesh is not better than the original mesh, rezoning cannot improve convergence, and can even worsen
convergence problems.
Requirements
PLANE182 -- B-bar method only (KEYOPT(1) = 0).
PLANE183
For both PLANE182 and PLANE183, all stress states (KEYOPT(3)) are allowed: plane strain, plane stress, axisymmetric, and generalized plane
strain. Pure displacement formulation (KEYOPT(6) = 0) or mixed u-P formulation (KEYOPT(6) = 1) is allowed.
SOLID185 -- B-bar method only (KEYOPT(1) = 0). Pure displacement formulation (KEYOPT(6) = 0) or mixed u-p formulation (KEYOPT(6) = 1) is
allowed in a nonlayered configuration (KEYOPT(3) = 0).
SOLID285
Contact elements
TARGE169
CONTA171 and CONTA172 with any of the following valid KEYOPT settings:
KEYOPT (1) = 0
96
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Requirements
Support Category
Requirements
KEYOPT (2) = 0, 1, 3, 4
KEYOPT (3) = 0
KEYOPT (4) = 0, 1, 2
KEYOPT (5) = 0, 1, 2, 3, 4
KEYOPT (7) = 0, 1, 2, 3
KEYOPT (8) = 0
KEYOPT (9) = 0, 1, 2, 3, 4
KEYOPT (10) = 0, 2
KEYOPT (11) = 0
KEYOPT (12) = 0, 1, 2, 3, 4, 5, 6
KEYOPT(14) = 0
TARGE170
CONTA173 and CONTA174 with any of the following valid KEYOPT settings:
KEYOPT (1) = 0
KEYOPT (2) = 0, 1
KEYOPT (4) = 0, 2
KEYOPT (5) = 0, 1, 2, 3, 4
KEYOPT (7) = 0, 1, 2, 3
KEYOPT (8) = 0
KEYOPT (9) = 0, 1, 2, 3, 4
KEYOPT (10) = 0, 2
KEYOPT (11) = 0
KEYOPT (12) = 0, 1, 2, 3, 4, 5, 6
KEYOPT(14) = 0
Contact pair
behavior
Materials
All structural materials are supported except for CAST (cast iron), CONCR
(concrete), MPLANE (microplane), SMA (shape memory alloy), and SWELL
(swelling).
For materials associated with user-defined state variables (TB,STATE), it
may be necessary to issue the MAPVAR command.
Analysis types
Loads and
boundary conditions (BCs)
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
97
Rezoning
Support Category
Requirements
Tabular displacements and forces are supported.
Tabular pressures that are functions of time only are supported.
Region to be
remeshed
The selected nodes inside the region must have the same nodal coordinate system.
Boundary nodes can have different nodal coordinate systems.
Elements must be of the same element type, material, element coordinate
system, and real constant.
If two regions with different attributes require remeshing, you must
remesh the regions separately. For more information, see Hints for
Remeshing Multiple Regions (p. 106).
Files
The conditions specified apply only to the region to be remeshed. No limitations exist for other regions,
although the analysis type itself must support rezoning.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Figure 4.2: Rezoning Using a Generic New Mesh Generated by Another Application
This flowchart shows a process for both 2-D and 3-D rezoning using a generic (.cdb
format) new mesh generated by a third-party application:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
99
Rezoning
100
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
Rezoning Comments
/CLEAR
Always clear the database first, before reentering the solution processor (/SOLU) and starting the rezoning process.
REZONE
Initiates rezoning
REMESH
Generates or
obtains the new
mesh required
for rezoning
You can remesh more than one region at the same specified
substep during the rezoning process.
After remeshing (REMESH,FINISH), the program generates contact elements (if any) and transfers loads and
boundary conditions automatically.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
101
Rezoning
Command
Description
Rezoning Comments
For more information, see Step 4: Perform the
Remeshing Operation (p. 104).
MAPSOLVE
102
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
103
Rezoning
A selected region that is too large may require more processing time and more subsequent remeshings.
If the selected region is too small to contain all of the highly distorted mesh areas, the new model after
rezoning may not be of sufficient quality to achieve convergence.
Using the GUI to Select a Region to Remesh
Select a region to remesh using either of the following methods available via the ANSYS Main Menu:
In the SOLU solution processor: Solution>Manual Rezoning>Select Rezone Elements
In the POST1 postprocessor: General Postproc>Manual Rezoning>Create Rezone Component
Select elements and group them into a component. During rezoning, import the component after
remeshing starts (REMESH,START).
This component can be imported into the program or exported to a third-party tool to create a new
mesh for the selected region.
If more than one region requires rezoning, see Remeshing Multiple Regions at the Same Substep (p. 116).
After you have selected the region(s) to remesh, proceed to Step 4: Perform the Remeshing Operation (p. 104).
104
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
105
Rezoning
106
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
107
Rezoning
Figure 4.4: Boundary Geometry of a Generic (CDB) New Mesh
If the rezoned part has contact/target elements, the program generates those elements automatically
for the new mesh, depending on whether the underlying old mesh had the same type of contact/target
elements. Isolated rigid target elements in the model remain the same throughout the analysis and
cannot be remeshed; however, all contact and target elements associated with solid elements are candidates for remeshing. While it is possible to read in the new contact/target elements of the new mesh
from the .cdb file, it is faster and more reliable to read in only the remeshed solid elements and allow
the program to generate the new contact/target elements.
The .cdb file of the new mesh must have no line breaks in the NBLOCK and EBLOCK records. Also,
while writing the mesh .cdb file, a block file format is necessary (CDWRITE,,,,,,,Fmat, where Fmat =
BLOCKED).
108
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The same problem appears in this example. However, the .cdb file for the new mesh
has node numbers defined from 28 through 43 and element numbers defined from 17
through 25. In this case, remeshing occurs using the REMESH command's KEEP option,
so the node and element numbers are not offset.
For more information, see the REMESH command documentation as it applies to remeshing using a
generic mesh created by another application.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
109
Rezoning
110
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
111
Rezoning
Figure 4.7: Splitting of Quadrilateral, Degenerate and Triangular Quadratic Elements ( LANE183)
P
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
113
Rezoning
Elements selected for splitting (ESEL) are marked as shown.
114
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Applicable Commands
Element size
Element shape
Element internal
angles
Refining
Other
The mesh-control commands are available at any point after issuing a REMESH,START command and
before issuing a REMESH,FINISH command.
Nodes at which force or isolated displacements are applied must remain during remeshing, as do nodes
at the boundaries of distributed displacements, pressures, or contact/target surfaces. The new elements
must have the same attributes as the old elements, such as element type, material, real constant, and
element coordinate system. The program rotates the new nodes to the same angles as the old nodes
in the region or on the boundary, as the case may be.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
115
Rezoning
After creating the new mesh, verify that the new elements cover the entire selected region and are
compatible with neighboring regions. (If another region requires remeshing, you can do so, but do not
issue another REMESH,START command.
116
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
dol ni egnhg
C C
C
stnemelpsid
deilppe
C C
C
C
C
hsem dlo
hsem wen
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
117
Rezoning
118
O P T I O N S
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
.
.
.
.
.
.
.
.
.
. .
UX
. .
. .
. .
. .
. .
. .
. .
. .
UY
. .
. .
. .
. .
. .
. .
. .
S T E P
. . . .
ROTZ
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
.2-D
.STATIC (STEADY-STATE)
.YES
.ON
.ON
.SPARSE
.PROGRAM CHOSEN
.SYMMETRIC
O P T I O N S
.
.
.
.
.
.
.
.
.
.
.
.
.
.
1
. 1.0000
.
ON
.
10
.
100
.
5
.
50
.
15
.
NO
.
ON
.YES (EXIT)
.USE DEFAULTS
.NO PRINTOUT
Mapping Substeps
The number of substeps used for mapping appears as follows:
*** LOAD STEP
1
*** TIME = 0.399110
SUBSTEP
10 COMPLETED.
REBALANCE FACTOR =
CUM ITER =
1.00000
101
Monitor File
The 6th and 7th columns of the monitor file indicate the rebalance factor (rather than the time) for
each mapping substep, as shown:
SOLUTION HISTORY INFORMATION FOR JOB: RznExample.mntr
13:06:16
LOAD
STEP
1
1
1
SUBSTEP
8
9
10
NO.
NO. TOTL
ATTMP ITER ITER
2
1
1
7
4
4
93
97
101
INCREMENT
REBALANCE
FACTOR
0.35000
0.35000
0.30000
09/20/2004
TOTAL
REBALANCE
FACTOR
0.35000
0.70000
1.0000
VARIAB 1
MONITOR
CPU
5.4300
8.0600
11.260
VARIAB 2
MONITOR
MxDs
-.64291
-.81383
-.93615
VARIAB 3
MONITOR
MxPl
0.78886E-30
0.78886E-30
0.78886E-30
Although multiple substeps may be necessary to balance all residuals, the program generates the restart
file and saves the result in a results file for only the last converged mapping substep. (Only the last
substep gives the balanced solution.) The program ignores any preexisting output specifications (set
via OUTRES or RESCONTROL commands, for example).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
119
Rezoning
120
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezone 0
(standard
run)
Rezone 1
Rezone 2
... Rezone 11
... Rezone nn
.rdb
.rdb
.rd01
.rd02
...
.rd11
...
.rdnn
.rxxx
.rxxx
.rxxx
.rxxx
...
.rxxx
...
.rxxx
...
.rs11
...
.rsnn
.ldhi
.ldhi
.rst
.rst
.rs01
.rs02
The maximum number of .rdb files and .rst files is 99. The .ldhi load history file contains information for all models created as a result of multiple rezonings.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
121
Rezoning
Displacement output for the rezoned model reports values as of the most recent rezoning, so the displacement will not seem to be continuous over multiple rezonings.
FILE SUFFIX
rst
rst
rst
rst
rst
rs01
rezone
rs01
rs01
The output from a SET,LIST command is also useful for simply determining what information is available.
For non-rezoning runs, for example, you can select substeps for further study by load step and substep
number, time, or set number.
Although the program creates multiple postprocessing data files for rezoning, you can consider them
a single file for POST1 processing (assuming that the appropriate jobname is already specified) because
there is no need to access individual files directly. The program automatically detects the specific file
needed according to the SET command issued, after which most of the usual postprocessing command
capability (PLESOL, PRNSOL, etc.) is available.
122
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Restrictions
During the macro's initial scan, the program stores the view location and size of every saved substep.
The program then combines the information to provide one fixed view (by internally setting the /FOCUS
and /DIST commands to fixed values).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
123
Rezoning
If multiple regions are read in (REMESH,READ) via a single .cdb file, the nodes of these regions must be
encapsulated in a single NBLOCK environment, and the elements must be encapsulated in a single EBLOCK
environment.
For a 2-D analysis, you can issue the following remeshing commands in the same session on multiple
regions:
AREMESH,AMESH (to create a program-generated new mesh)
REMESH,READ (to import a generic new mesh)
REMESH,SPLIT (to split an existing mesh)
For a 3-D analysis, remeshing is possible only by importing a new mesh created by another application
(REMESH,READ).
The POST1 load-case combination commands (such as LCASE, LCDEF, LCOPER, and LCSEL) are not
available.
Initial state loading (INISTATE) is not available for rezoning.
Rezoning is not available in Distributed ANSYS.
124
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Examples
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
125
Rezoning
k,2,xc,yc,yc
k,3,xc-rc,yc
k,4,0.0,0.0
k,5,2*rc,0.0
rect,0,b,0,h
circle,1,rc,2,3,360,1
/pnum,line,1
lplot
l,4,5
lplot
aplot
/com create solid elements
esize,el
mat,1
type,1
real,1
amesh,1
/com generate the 1st contact pair
mat,2
real,3
type,2
esize,h
lmesh,5,7
*get,PilotID,node,,num,max
PilotID=PilotID+1
nkpt,PilotID,1
tshap, pilo
e,PilotID
type,3
lsel,,,,2,3
nsll,,1
esln,,0
esurf
alls
/com generate the 2nd contact pair
real,4
type,4
lmesh,8
lsel,,,,8
esll
esurf,,reverse
alls
type,5
lsel,,,,1,2
nsll,,1
esln,,0
esurf
alls
/com apply boundary conditions and loads
d,PilotID,ux,0.0
d,PilotID,uy,PilotMove
d,PilotID,rotz,0.0
lsel,,,,4
nsll,,1
d,all,ux,0.0
alls
lsel,,,,8
nsll,,1
d,all,uy,0.0
alls
nlist
126
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Examples
elist
dlist
/com check the contact definition
cncheck
finish
/solution
rescontrol,,all,1
pred,off
nlgeom,on
time,1
NSUBST,10,100,5
outres,all,all
solve
finish
/post1
set,1,6
prns,u,comp
prns,s,comp
prns, cont
finish
remesh,start
esel,,,,65,128
aremesh
lesize,10,,,8
lesize,31,,,10
lesize,18,,,6
amesh,2
esel,all
nsel,all
remesh,fini
elist
dlist
mapsolve,50
! map solutions
finish
/solution
antype,,restart
solve
finish
! restart
127
Rezoning
run stops at t = 0.7875. Rezoning is applied at this time to achieve complete loading. The entire deformed
model at substep 4 is imported into ANSYS ICEM CFD, which generates a new mesh. After reading the
new mesh back in, the program creates the contact automatically when you issue the REMESH,FINISH
command.
The solid element used in the model is PLANE182 (using the B-Bar method with mixed u-P formulation).
CONTA171 and TARGE169 elements are also used. The material model used is a hyperelastic material
with a three-parameter OGDEN option.
This input results in a deformed mesh and causes the analysis to terminate at t = 0.7875 seconds:
/batch, list
/file,RznExample2
/prep7
h=4.6295
b=1.5
el=b/4
xc=0
yc=2.6295
rc=2.5
PilotMove= -yc
! ogden parameters
TB,HYPE,1,1,3,OGDE
TBTEMP,0
TBDATA,1,3.2084E-009,7.281,0.035198,3.0149,6.3712,2.0493
et,1,182
keyopt,1,3,0
keyopt,1,6,1
et,2,169
et,3,171
keyopt,3,9,0
keyopt,3,10,1
et,4,169
et,5,171
keyopt,5,9,0
keyopt,5,10,1
mp,mu,2,0.0
r,3
r,4
k,1,xc,yc
k,2,xc,yc,yc
128
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Examples
k,3,rc,yc
k,4,0.0,0.0
k,5,rc+1,0.0
rect,0,b,0,h
circle,1,rc,2,3,90,1
/pnum,line,1
lplot
l,4,5
lplot
aplot
esize,el
mat,1
type,1
real,1
amesh,1
/pnum,elem,1
/pnum,node,1
/com the 1st contact pair
mat,2
real,3
type,2
esize,h
lmesh,5
lsel,,,,5
esll
esurf,,reverse
alls
*get,PilotID,node,,num,max
PilotID=PilotID+1
nkpt,PilotID,1
tshap, pilo
e,PilotID
type,3
lsel,,,,2,3
nsll,,1
esln,,0
esurf
alls
/com the 2nd contact pair
real,4
type,4
lmesh,6
lsel,,,,6
esll
esurf,,reverse
alls
type,5
lsel,,,,1,2
nsll,,1
esln,,0
esurf
alls
d,PilotID,ux,0.0
d,PilotID,uy,PilotMove
d,PilotID,rotz,0.0
lsel,,,,4
nsll,,1
d,all,ux,0.0
alls
lsel,,,,6
nsll,,1
d,all,uy,0.0
alls
/solution
pred,off
rescontrol,,all,1,
eresx,no
nlgeom,on
time,1
NSUBST,10,100,5
outres,all,all
solve
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
129
Rezoning
finish
Following is the total elastic strain along the Y axis at t = 0.7875. The element distortion is apparent
and causes the problem to diverge.
!
!
!
!
!
!
!
With the deformed mesh corresponding to load step 1, substep 4 is shown next. ANSYS ICEM CFD uses
the boundary segments of this mesh next to generate the new mesh. The nodal discretization at the
boundary remains same for both the old and the new mesh.
130
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Examples
Next, the total elastic strain along the Y axis for this mesh is shown. This is one of the state variables
which is transferred to the new mesh when mapping solved node and element solutions from the original mesh to the new mesh (MAPSOLVE).
4.14.2.3. Importing the File into ANSYS ICEM CFD and Generating a New Mesh
At this stage of the rezoning process, start ANSYS ICEM CFD and read in the .cdb file. (Reminder: As
indicated in Exporting the Distorted Mesh as a CDB File (p. 130), only solid elements can be read in.)
Generate the new .cdb file as follows:
1. Import the .cdb file in ANSYS ICEM CFD as mesh (File Menu> Import Mesh> From Ansys).
2. Extract triangulated (STL) geometry from the mesh (Edit Menu> Mesh to Facets)
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
131
Rezoning
3. Set the global maximum element size in the order of the element size that you require (Mesh> Global
Mesh Setup> Global Mesh Size> Max Element).
4. Build the topology (Geometry> Repair geometry> Build topology)
5. Select the Respect line elements and Protect given line elements options (Mesh> Global Mesh Setup>
Shell Meshing Params).
6. Compute the new mesh (Mesh> Compute mesh> Surface mesh only> Mesh type: All Quad > Compute).
7. Select the Solve Options tab and write the input file. Do not include the bar elements.
8. Rename the new input file as a .cdb file.
The new mesh obtained from ANSYS ICEM CFD is shown here. Notice that the boundary discretization
remains the same as that of the old mesh.
!
!
!
!
!
!
!
After the MAPSOLVE command has executed (mapping the solved node and element solutions from
the original mesh to the new mesh), the total elastic strains along Y for the new mesh appears. Notice
that some expected nodal realignment has occurred in the new mesh.
132
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Rezoning Examples
!
!
!
!
Allow the analysis to complete. Following is a plot of the total elastic strain along the Y direction:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
133
134
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Requirements
PLANE182 -- Plane strain and axisymmetric states of stress only
PLANE183 -- Plane strain and axisymmetric states of stress only
Enhanced strain formulation is not supported.
Contact elements
TARGE169
CONTA171 and CONTA172 with any of the following valid KEYOPT settings:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
135
Requirements
KEYOPT (1) = 0
KEYOPT (2) = 0, 1, 3, 4
KEYOPT (3) = 0
KEYOPT (4) = 0, 1, 2
KEYOPT (5) = 0, 1, 2, 3, 4
KEYOPT (7) = 0, 1, 2, 3
KEYOPT (8) = 0
KEYOPT (9) = 0, 1, 2, 3, 4
KEYOPT (10) = 0, 2
KEYOPT (11) = 0
KEYOPT (12) = 0, 1, 2, 3, 4, 5, 6
KEYOPT(14) = 0
TARGE170
CONTA173 and CONTA174 with any of the following valid KEYOPT settings:
KEYOPT (1) = 0
KEYOPT (2) = 0, 1
KEYOPT (4) = 0, 2
KEYOPT (5) = 0, 1, 2, 3, 4
KEYOPT (7) = 0, 1, 2, 3
KEYOPT (8) = 0
KEYOPT (9) = 0, 1, 2, 3, 4
KEYOPT (10) = 0, 2
KEYOPT (11) = 0
KEYOPT (12) = 0, 1, 2, 3, 4, 5, 6
KEYOPT(14) = 0
Contact pair
behavior
Materials
All structural materials are supported except for CAST (cast iron), CONCR
(concrete), MPLANE (microplane), SMA (shape memory alloy), and SWELL
(swelling).
For materials associated with user-defined state variables (TB,STATE), it
may be necessary to issue the MAPVAR command. The defined size of
state variables must be such that full 3-D descriptions of stress and strain
tensors can be accommodated (that is, the defined stress and strain
tensors must be able to possess six independent components).
136
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Requirements
Displacements, forces, pressures, and nodal temperatures (BF,TEMP).
Loads applied on surface-effect elements (such as SURF153) are not
supported.
Tabular displacements and forces are supported.
Tabular pressures that are functions of time only are supported.
Files
The EEXTRUDE command documentation shows the topological mapping of 2-D elements (solid, contact,
and target) into specific 3-D elements.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
137
For information about the commands shown in the flowchart, see Key Commands Used in 2-D to 3-D
Analysis (p. 138)or the documentation for a given command in the Command Reference.
Description
/CLEAR
Always clear the database first, before reentering the solution processor (/SOLU) and starting the 2-D to 3-D analysis
process.
MAP2DTO3D,START
Initiates the
2-D to 3-D
analysis
138
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
EEXTRUDE
Extrudes
the 2-D
mesh to a
3-D mesh.
Deletes the
2-D rigid
targets and
regenerates
the 3-D rigid targets.
The entire 2-D mesh is selected automatically for extrusion. Partial mesh extrusion is not supported.
In the extrusion phase of a 2-D to 3-D analysis, the
command generates the 3-D contact elements. The
command also regenerates the loads, boundary conditions, and temperatures for the extruded 3-D mesh.
Maps the
loads, temperatures,
and boundary conditions from
the 2-D
mesh to
the corresponding 3D mesh.
Maps
MAP2DTO3D,FIN- node and
element soluISH
tions
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
139
5.4.3. Step 3: Extrude the 2-D Mesh to the New 3-D Mesh
Extrude the 2-D mesh to the 3-D mesh (EEXTRUDE).
The following figures illustrate how extrusion generates 3-D meshes from a plane strain and an
axisymmetric 2-D mesh, respectively:
Figure 5.2: 2-D Plane Strain to 3-D Solid Extrusion
Extrusion preserves the 2-D mesh topology on the generating plane and replicates it through the depth
of the extruded body.
140
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
141
142
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
143
144
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Cyclic Modeling
6.2.2. Edge Component Pairs
6.2.3. Model Verification (Preprocessing)
Define a basic sector model that is cyclically symmetric in any global or user-defined cylindrical coordinate
system. (For information about creating a model, see the Modeling and Meshing Guide.)
The angle (in degrees) spanned by the basic sector should be such that N = 360, where N is an integer.
The basic sector can consist of meshed or unmeshed geometry. The program allows user-defined
coupling and constraint equations on nodes that are not on the low or high edges of the cyclic sector.
(For more information about the cyclic sector's low and high edges, see Edge Component Pairs (p. 146).)
If meshed, the basic sector may have matching (as shown in Figure 6.3: Basic Sector Definition (p. 145))
or unmatched lower and higher angle edges. Matching means that corresponding nodes exist on each
edge, offset geometrically by the sector angle . The edges may be of any shape and need not be "flat"
in cylindrical coordinate space. For more information, see Identical vs. Dissimilar Edge Node Patterns (p. 149).
Figure 6.3: Basic Sector Definition
Z
Y
Sector
angle
X
CSYS = 1
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
145
146
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Cyclic Modeling
Figure 6.4: Full Cyclic Model
Figure 6.5: Cyclic Sector (p. 147) shows a cyclic sector from the model above. As you can see, this sector
model is leaning heavily in the direction of theta. The low and high boundaries in this figure have been
auto detected correctly using the default FACETOL value of 15 deg. This successful auto detection of
these boundaries can be seen in Figure 6.6: Successful Auto Detection with Default FACETOL = 15
Deg (p. 148).
Figure 6.5: Cyclic Sector
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
147
If you reset the FACETOL values in the figure above from 15 degrees to 30 degrees, auto detection
fails because the logic treats the left and bottom sides of the sector as a single face, and the top and
right sides as a second large face. This is because the bottom-left edge and top-right edge of this
model both contain a dihedral angle greater than 150 degrees (180-30), but less than 165 degrees (18015). These large angles add to the possibility of auto-detection errors.
If you start with a model tilted more than the sector in Figure 6.6: Successful Auto Detection with Default
FACETOL = 15 Deg (p. 148), it may fail at the default tolerance of 15 degrees. This may require you to
reduce FACETOL to 10 or even 5 degrees to get a successful result.
A FACETOL value that is set too low can also result in failure. A FACETOL value that is too low can
cause edges not to be detected along element boundaries on smooth surfaces. The valid range of FACETOL is model and mesh dependent, and may be dictated by a single edge shared between 2 elements.
Whenever auto detect fails for an element model, save the node groups for each element face that you
were working with as node components NPF_001, NPF_002, etc. for diagnostic purposes. Each node
component should represent exactly one face of the cyclic sector (cyclic boundary or not).
In Figure 6.7: Auto Detection Failure Due to Large Face Tolerance (p. 149) , NPF_001 is clearly too large.
As you can see from the nodes, the group has leaked across the lower left edge. This indicates that
FACETOL is too large for the given dihedral angle.
148
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Cyclic Modeling
Figure 6.7: Auto Detection Failure Due to Large Face Tolerance
149
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
Internal CEs cannot be listed or viewed after a cyclic symmetry solution.
During the solution stage of a cyclic symmetry analysis, the program automatically rotates the nodal
coordinate systems of all nodes on the low and high sector edges to be parallel with the cyclic coordinate
system.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
151
k
k
k
k
wo vf
o
wo v e
h gir f
h g ir e
(6.1)
where,
k = Harmonic index -- (0,1,2,,N / 2) when N is even, (0,1,2,,(N-1) / 2) when N is odd. (N is an integer
representing the number of sectors in 360.)
= Sector angle (2 / N)
U = Vector of displacement and rotational degrees of freedom
UALow represents the basic sector low side edge
UAHigh represents the basic sector high side edge
UBLow represents the duplicate sector low side edge
UBHigh represents the duplicate sector high side edge
The equation is a function of harmonic index k generating different sets of constraint equations for
each harmonic index. Therefore, for each harmonic index solution requested, the program creates the
appropriate constraint equations automatically, connects the edge-component nodes on basic sector
A and duplicate sector B, and solves.
Constraint equations that tie together the low and high edges of your model are generated from the
low- and high-edge components, and nowhere else. You should verify that automatically detected
components are in the correct locations and that you are able to account for all components; to do so,
you can list (CMLIST) or plot (CMPLOT) the components.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Commands
Loads
Comments
Nodal Loads
F, K
Surface
Loads
SF,
SFA,
SFE,
SFL
Inertia Loads
ACEL,
DOMEGA,
CMDOMEGA,
CMOMEGA,
OMEGA
153
SFL,ALL,PRES,%S1PRES%
When combined with other independent variables, SECTOR can be in positions 1, 2, or
3 only. Other independent variables operate as they do for non-cyclic data. (Think of X,
Y, and Z as ghost coordinates, behaving as though all sectors have been modeled with
actual nodes and elements.)
Example 6.3: Deleting a Sector Load
CYCOPT,LDSECT,3
154
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
! Apply a load (
value) on node 10 at sector
...
FDELE,10,UX
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
155
Only a harmonic index zero solution is valid for a static solution with cyclic loading.
For non-cyclically symmetric loading, the program supports linear static analysis only. The following
flowchart illustrates the process involved in a static cyclic analysis with non-cyclic loading.
Figure 6.10: Process Flow for a Static Cyclic Symmetry Analysis (Non-Cyclic Loading)
156
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
For a complicated structure exhibiting cyclic symmetry (for example, a turbine wheel), lines of zero
displacement may not be observable in a mode shape.
The harmonic index is an integer that determines the variation in the value of a single DOF at points
spaced at a circumferential angle equal to the sector angle. For a harmonic index equal to nodal diameter d, the function cos(d*) describes the variation. This definition allows a varying number of waves
to exist around the circumference for a given harmonic index, provided that the DOF at points separated
by the sector angle vary by cos(d*). For example, a harmonic index of 0 and a 60 sector produce
modes with 0, 6, 12, ... , 6N waves around the circumference.
The nodal diameter is the same as the harmonic index in only some cases. The solution of a given
harmonic index may contain modes of more than one nodal diameter.
The following equation represents the relationship between the harmonic index k and nodal diameter
d for a model consisting of N sectors:
(6.2)
k xm d
where m = 0, 1, 2, 3, ...,
For example, if a model has seven sectors (N = 7) and the specified harmonic index k = 2, the program
solves for nodal diameters 2, 5, 9, 12, 16, 19, 23, ....
The following table illustrates Equation 6.2 (p. 157), showing how the harmonic index, nodal diameter
and number of sectors relate to one another:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
157
2N
2N
...
N-1
N+1
2N - 1
2N + 1
...
N-2
N+2
2N - 2
2N + 2
...
N-3
N+3
2N - 3
2N + 3
...
N-4
N+4
2N - 4
2N + 4
...
...
...
...
...
...
...
...
N/2
N/2
3N / 2
3N / 2
5N / 2
...
(N - 1) /
2
(N + 1) /
2
(3N - 1)
/2
(3N + 1)
/2
(5N - 1)
/2
...
N/2
(N
is
even)
(N - 1) /
2
(N
is
odd)
Note
To avoid confusion, be aware that in some references mode refers to harmonic index as
defined here and nodal diameter describes the actual number of observable waves around
the structure.
Harmonic Index in an Electromagnetic Analysis For electromagnetic analyses, only the EVEN and
ODD harmonic index settings (see the CYCOPT command) are valid (for symmetric and antisymmetric
solutions, respectively).
Using VT Accelerator You can use the Variational Technology Accelerator (VT Accelerator) to speed
up the solve time needed to sweep over the range of values of the harmonic index. To activate VT Accelerator, issue CYCOPT, VTSOL prior to solving. You can use VT Accelerator only with matched node
pattern sectors in a modal cyclic symmetry analysis. Linear perturbation modal analysis does not support
VT Accelerator. You will see the most significant speed up for models with a large number of sectors
and/or a large number of eigenvalues. The benefit of using VT Accelerator is realized only for five or
more Harmonic Indices (HI). In addition, the level of performance improvement realized with VT Accelerator may also be dependent up on the problem. Using less than five sectors prevents a solution and
displays an error message.
The Variational Technology acceleration method is available for the perturbed modal cyclic symmetric
analysis.
158
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
A modal cyclic symmetry analysis allows only cyclically symmetric applied boundary constraints and
applied loads. Applied loads are not used in a modal analysis (which is a free-vibration problem). Eigensolutions are performed, looping on the number of harmonic indices specified (via the CYCOPT
command) at each load step.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
159
The modal cyclic symmetry solution occurs after the static cyclic symmetry solution. The modal solution
uses the same low- and high-edge components defined in the static cyclic analysis stage (via the CYCLIC
command). The analysis yields the eigenvectors of the structure in the prestressed state.
160
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Prestressed Modal
SOLUTION
Perform a static restart with perturbation:
ANTYPE,STATIC,RESTART,,PERTURB
PREP7
Define the basic sector
Perform a large-deflection:
NLGEOM,ON
POST1
Specify the result file for the linear
perturbation: FILE,,RSTP
For detailed information about the linear perturbation analysis procedure, see Linear Perturbation
Analysis in the Structural Analysis Guide. In addition, Example 9.4: Contact Status Control in a Linear
Perturbation Modal Analysis in that same document provides a complete example input listing for a
linear perturbation cyclic symmetry modal analysis.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
161
You can specify the SET command's ORDER option to sort the harmonic index results in ascending order
of buckling load multipliers.
162
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
163
164
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Meaning: 124 constraint equations are created, used, and then deleted to enforce cyclic symmetry
conditions between the low- and high-edge nodes. Every node on the low edge is precisely matched
to a corresponding node on the high edge, representing the best possible situation.
NUMBER OF CONSTRAINT EQUATIONS GENERATED =
124
(USING THE MATCHED NODES ALGORITHM -MAX NODE LOCATION ERROR = 0.73906E-02)
Meaning: 124 constraint equations are created, used, and then deleted to enforce cyclic symmetry
conditions between the low- and high-edge nodes. Every node on the low edge is matched to a corresponding node on the high edge within the current tolerance setting, but not all matches are precise.
The largest position mismatch is 0.0073906.
NUMBER OF CONSTRAINT EQUATIONS GENERATED =
(USING THE UNMATCHED NODES ALGORITHM)
504
Meaning: 504 constraint equations are created, used, and then deleted to enforce cyclic symmetry
conditions between the low- and high-edge nodes. At least one node on the low edge does not match
any node on the high edge within the current tolerance setting, so the program uses the unmatched
nodes algorithm.
165
Prestressed Modal
SOLUTION
Perform a static restart with perturbation:
ANTYPE,STATIC,RESTART,,PERTURB
PREP7
Define the basic sector
Perform a large-deflection:
NLGEOM,ON
SOLUTION
Define static analysis:
ANTYPE, STATIC
POST1
Specify the result file for the linear
perturbation: FILE,,RSTP
166
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The CYCPHASE command uses full model graphics (/GRAPHICS,FULL) to compute peak
values. Because of this, there may be slight differences between max/min values obtained
with CYCPHASE, and those obtained via /CYCEXPAND (/GRAPHICS,POWER).
For information about /CYCEXPAND and EXPAND command usage, see Expanding the Cyclic Symmetry
Solution (p. 167). For information about CYCPHASE command usage, see Phase Sweep of Repeated Eigenvector Shapes (p. 169).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
167
(6.3)
where,
Uj = Response of the full structure or assembly (displacement, stress, or strain) for sector number j
UA = Basic sector solution
UB = Duplicate sector solution
j = Sector number for response expansion -- j = 1,2,3,,N
k = Harmonic index -- (0,1,2,,N / 2) when N is even, (0,1,2,,(N-1) / 2) when N is odd. (N is an integer
representing the number of sectors in 360.)
= Sector angle (2/N)
The equation applies when expanding the static cyclic solution as well as the modal cyclic eigenvector
solution.
If cyclic expansion via the /CYCEXPAND command is active for a static or harmonic cyclic symmetry
analysis, the PLNSOL and PRNSOL commands have summation of all required harmonic index solutions
by default. (You can override the default behavior if necessary.)
In a static or harmonic cyclic symmetry analysis with non-cyclic loading, all applicable harmonic index
solutions are computed and saved in the results file as load step results.
A SET,LIST command lists the range of load step numbers in the group containing each solution. Each
load step post data header contains the first, last, and count of load steps from the given SOLVE command, as shown:
SET
1
2
3
4
5
6
...
GROUP
1-3
1-3
1-3
4-6
4-6
4-6
The SET command establishes which SOLVE load step group should display. Summation via /CYCEXPAND
is automatic (although you can override the default behavior). Plots and printed output show the
summation status.
With /CYCEXPAND turned on, the results are expanded at each load step and then combined to plot
the full solution as a complete sum. For example, in a four sector model where harmonic index results
0 through 3 are available in the results file, the plot command PLNSOL will display the results as STEP=1
THRU=3 COMPLETE SUM.
Accumulation occurs at the first applicable PLNSOL or PRNSOL command. After accumulation, the last
load step number of the current group becomes the new current load step number.
168
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Caution
Do not confuse the EXPAND command with /EXPAND.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
169
(6.4)
where,
c1 and c2 = Arbitrary constants
Ui and Ui+1 = Eigenvectors corresponding to fi and fi+1, respectively
The orientation of the combined mode shape U will be along a nodal diametral line that is neither
along that of Ui nor Ui+1. Because the full structure may have stress-raising features (such as bolt holes),
determining the eigenvector orientation that causes the most severe stresses, strains, or displacements
on the structure or assembly is critical.
To determine the peak value of stress, strain or displacement in the full structure or assembly, it is necessary to calculate U at all possible angular orientations in the range of 0 through 360. In the general postprocessor, the CYCPHASE command performs the computational task.
Because c1 and c2 are arbitrary constants, the CYCPHASE calculation rewrites Equation 6.4 (p. 170) as
follows:
(6.5)
Using the cyclic symmetry expansion of Equation 6.3 (p. 168) in Equation 6.5 (p. 170), the simplified phasesweep equation that operates on the cyclic sector solution (rather on the computation-intensive fullstructure expression in Equation 6.5 (p. 170)) is:
(6.6)
A phase sweep using the CYCPHASE command provides information about the peak values of stress,
strain and/or displacement components and the corresponding phase angle values. Using the phase
angle value further, you can expand the mode shape at that phase angle to construct the eigenvector
shape that produces the peak stress, strain and/or displacement. The expansion expression with the
phase angle used by the /CYCEXPAND command is:
(6.7)
where,
j = 1,2,3,...,N
Example:
To determine the eigenvector orientation that causes the highest equivalent stress, perform a phase
sweep on the stress via the CYCPHASE,STRESS command. Obtain a summary of the phase sweep via
the CYCPHASE,STAT command to determine the value of at which maximum equivalent stress occurred.
170
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The CYCPHASE command uses full model graphics (/GRAPHICS,FULL) to compute peak
values. Because of this, there may be slight differences between max/min values obtained
with CYCPHASE, and those obtained via /CYCEXPAND (/GRAPHICS,POWER).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
171
172
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
173
174
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
175
Description
Command
1.
---
2.
/INPUT,CYC_SYMM.INP
3.
CYCLIC
4.
AMESH,ALL
5.
/CYCEXPAND,,ON
6.
EPLOT
7.
CYCLIC,STATUS
8.
CYCOPT,STATUS
9.
/SOLU
SOLVE
10.
/POST1
RSYS,1
11.
SET,2,6
12.
PLNSOL,U,Y
SET,13,1
14.
PLNSOL,U,Y
Executing this step causes the struts of the assembly to bend out
of phase.
15.
SET,2,5
16.
PLNSOL,U,Y
176
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
Mode shape values may vary slightly depending on your computer system.
To view a traveling wave animation of your model, issue the ANCYC,24,,0.1 command. For more information, see Applying a Traveling Wave Animation to the Cyclic Model (p. 169).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
177
178
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
179
180
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
Command
1.
---
2.
/INPUT,buck_cyc_sym.inp
3.
CYCLIC
4.
/CYCEXPAND,,ON
5.
EPLOT
6.
CYCLIC,STATUS
7.
CYCOPT,STATUS
8.
9.
10.
11.
BUCOPT, LANB, 3
12.
MXPAND, 3, , , YES
13.
SOLVE
14.
Read the results from the smallest load factor. (This should correspond to the smallest frequency.)
15.
/SOLU
ANTYPE,STATIC
PSTRES,
ON
SOLVE
FINISH
/SOLU
ANTYPE,BUCKLE
FINISH
/POST1
SET, FIRST ,
, , , , , , ORDER
PLNSOL, U, SUM
---This step completes the example buckling cyclic symmetry analysis. Your results should match those shown in Figure 6.20: Example Buckling Cyclic Symmetry Analysis Results (p. 182).
The results of your analysis should match those shown here:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
181
Use the /PREP7 portion of the previous input file (buck_cyc_sym.inp) to construct your model.
After defining the model parameters--but before activating cyclic symmetry--define the arrays and the
programming operations, as follows:
*dim,Tstrut,array,10
*dim,Tfact,array,10
*do,I,1,10
/prep7
182
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
After you have defined the iterative parameters, proceed with the cyclic symmetry portion of the analysis:
cyclic
/cycexpand,,on
eplot
cyclic,status
cycopt,status
/solution
antype,static
pstres,on
solve
finish
/solution
antype,buckle
bucopt,lanb,3
mxpand,3,,,yes
solve
finish
/post1
set,first,,,,,,,order
plnsol,u,sum
*get,loadmult,active,,set,freq
Tfact(I)=loadmult
*enddo
The program then plots the data to determine the critical strut temperature:
*dim,data,table,10,2
data(0,1)=1
data(0,2)=2
*do,I,1,10
data(I,0)=I
data(I,1)=Tstrut(I)
data(I,2)=Tfact(I)
*enddo
/AXLAB,X,Strut Temperature
/AXLAB,Y,Load Factor
/GROPT,DIVX,5
/GROPT,DIVY,5
/XRANGE,100,200
/YRANGE,0.9,1.4
/GTHK,CURVE,1
/GMARKER,1,3
*VPLOT,data(1,1),data(1,2)
The eigenvalues (frequencies) calculated for the buckling analysis represent the buckling load factors.
The eigenvalues represent load factors for all applied loads.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
183
T (Strut)
Load Factor
100.00
1.3039
130.39
1.1845
154.44
1.0972
169.45
1.0461
177.27
1.0206
180.91
1.0089
182.52
1.0038
183.21
1.0016
183.50
1.0007
10
183.62
1.0003
184
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
185
186
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
set,,3
plns,u,sum,1
fini
! STEP #12
! Go to time/freq. history postprocessor
/post26
! STEP #13
! Store X and Y disp. of node #34 on sector #3"
nsol,2,34,U,x,UX_34,3
nsol,3,34,U,y,UY_34,3
! STEP #14
! Plot UX and UY disp. of node #34 as a function of freq.
plvar,2,3
Step Description
Description
1.
---
2.
/INPUT,HARM,CYC_SYMM.INP
3.
4.
AMESH,ALL
5.
/CYCEXPAND,,ON
6.
EPLOT
7.
CYCLIC,STATUS
8.
CYCOPT,STATUS
9.
/SOLU
SOLVE
10.
/POST1
11.
SET,1,3
12.
Plot the displacement sum contour. Figure 6.24: Contour Plot of Displacement Sum at
Frequency of 866 HZ (p. 189) shows a contour
plot of displacement sum at frequency = 866
HZ.
PLNSOL,U,SUM
13.
/POST26
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
187
NSOL,2,34,u,x,UX_34,3
NSOL,3,34,u,y,UY_34,3
15.
188
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
189
190
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Cf eshp liog
roto
gf litnegnt xulp
Ee eshp liog
rott
f
C
x
C
f
Ce eshp liog
f
pg rie
Ef eshp liog
d
d w
d
191
Apply Cylic
/title,Cyclic Symmetry Model for EMAG Analysis (Dual Coils with Iron Yoke)
/com
/com ***** Quarter Symmetry Model Expanded to Half Then to Full *****
/com
/com
/com
/nopr
/out,scratch
!!! Setup Model Parameters
_geomgen=1
p=1
!
!
!
!
!
alpha=22.5/p
beta=alpha+(45/p)
gamma=beta+(22.5/p)
r1=3
r2=4.5
r3=5
r4=7
r5=11
ncoil=(4*p)
i1=1
i2=2
*dim,alpha1,,ncoil
*dim,alpha2,,ncoil
*dim,current,,ncoil
*dim,coilname,string,ncoil
coilname(1)
coilname(2)
coilname(3)
coilname(4)
=
=
=
=
'coil1'
'coil2'
'coil3'
'coil4'
*do,i,1,ncoil
alpha1(i) = -alpha + (i-1)*(90/p)
alpha2(i) = alpha + (i-1)*(90/p)
*enddo
ii=0
*do,i,1,p
ii = ii + 1
current(ii) = i2
ii = ii + 1
current(ii) = i1
192
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
/prep7
ET,1,13,4
! Al Core
ASEL,S,LOC,X,r1,r2
CM,Outer_AL,AREA
AATT,2,,1,
! Air Gap
ASEL,S,LOC,X,r2,r3
CM,AIR,AREA
AATT,4,,1
! Coil 1
ASEL,S,LOC,X,r3,r4
ASEL,R,LOC,Y,0,alpha
CM,COIL1,AREA
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
193
/prep7
save,magtest,db
arsym,y,all
nummrg,all
csys,1
nsel,s,loc,x,r5
CM,extnode,NODE
! Apply BFE Current loads to each coil
*do,i,1,ncoil
asel,s,loc,x,r3,r4
asel,r,loc,y,alpha1(i),alpha2(i)
esla,s
cm,coilname(i),element
bfe,all,js,,,,current(i)
*enddo
csys,0
allsel
cmsel,s,extnode
d,all,az,0
194
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
!! Delete Bottom half of model and all loading attatched to bottom nodes
/prep7
allsel
nummrg,all
csys,1
nsel,s,loc,x,r5
D,all,az,0 ! AZ = 0 on outside nodes of arc
D,all,temp,25
!! Define Coils on Half Model
! Coil 1
ASEL,S,LOC,X,r3,r4
ASEL,R,LOC,Y,0,alpha
esla,s
CM,COIL1,ELEMENT
! Coil 2
ASEL,S,LOC,X,r3,r4
ASEL,R,LOC,Y,beta,(180-beta)
esla,s
CM,COIL2,ELEMENT
! Coil 3
ASEL,S,LOC,X,r3,r4
ASEL,R,LOC,Y,(180-alpha),180
esla,s
CM,COIL3,ELEMENT
!! Apply bfe loads to Half Model coils
cmsel,s,COIL1
bfe,all,js,,,,i2
cmsel,s,COIL2
bfe,all,js,,,,i1
cmsel,s,COIL3
bfe,all,js,,,,(-i2)
!!!
allsel
csys,0
cyclic,2
/solution
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
195
cycopt,hindex,odd
solve
FINISH
/post1
/out
Figure 6.28: Vector Plot of Cyclic Flux Density (B) - Half Model
196
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where [M], [C] and [K] are the structural mass, damping, and stiffness matrices, respectively.
[Kc] is the spin softening matrix due to the rotation of the structure. It changes the apparent stiffness
of the structure in a rotating reference frame (described in Rotating Structures in the Mechanical APDL
Theory Reference).
[G] is a damping matrix contribution due to the rotation of the structure. It is usually called Coriolis
matrix in a rotating reference frame, and gyroscopic matrix in a stationary reference frame (described
in Rotating Structures in the Mechanical APDL Theory Reference).
{F} is the external force vector in the stationary reference frame. In a rotating reference frame, it is the
sum of the external force and the effect of the angular rotational velocity force (as described in Acceleration Effect in the Mechanical APDL Theory Reference).
Without the inertia effect applied via the CORIOLIS command, ANSYS does not generate the [G] matrix,
and the usual effect of the angular rotation velocity specified by the OMEGA or CMOMEGA command
applies (as described in Acceleration Effect in the Mechanical APDL Theory Reference). An exception exists,
however, involving a nonlinear transient analysis using element MASS21; in this case, the inertia effect
due to rotation applied via an IC command (or a D command over an incremental time) is included
without having to issue the CORIOLIS and OMEGA or CMOMEGA commands.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
197
198
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
For a prestressed structure, set the Campbell key (CAMPBELL,ON) in the first solution pass.
Doing so allows a Campbell diagram analysis.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
199
Note
For a rotating structure meshed in shell elements lying in a plane perpendicular to the rotational velocity axis - such as a thin disk - the whirl effects are not plotted or printed by the
PRCAMP or PLCAMP commands. However, they can be visualized using the ANHARM
command.
To retrieve and store frequencies and whirls as parameters, use the *GET command.
Prestressed Structure
For a prestressed structure, set the Campbell key (CAMPBELL,ON) in the static solution portion of the
analysis. Doing so modifies the result file so that it can accommodate a subsequent Campbell diagram
analysis. In this case, static and modal solutions are calculated alternately and only the modal solutions
are retained.
Generating a Successful Campbell Diagram
To help you obtain a good Campbell diagram plot or printout, the sorting option is active by default
(PLCAMP,ON or PRCAMP,ON). ANSYS compares complex mode shapes and pairs similar mode shapes.
(Because eigenmodes at zero velocity are real modes, ANSYS does not pair them with complex modes.)
If the plot is unsatisfactory even with sorting enabled, try the following:
Start the Campbell analysis with a non-zero rotational velocity.
Modes at zero rotational velocity are real modes and may be difficult to pair with complex modes
obtained at non-zero rotational velocity.
Increase the number of load steps.
It helps if the mode shapes change significantly as the spin velocity increases.
Change the frequency window.
To do so, use the shift option (PLCAMP,,,FREQB or PRCAMP,,,FREQB). It helps if some modes fall
outside the default frequency window.
Overcoming Memory Problems
To run the Campbell analysis (PRCAMP or PLCAMP), the scratch memory needed may be important as
complex mode shapes are read from the result file for two consecutive load steps. If your computer
has insufficient scratch memory, try the following:
Decrease the number of extracted modes (MODOPT,,NMODE)
Generate the result file for a reduced set of selected nodes (for example, nodes on the axis of rotation).
Issue OUTRES,ALL,NONE and then OUTRES,Item,Freq,Cname where Item=NSOL, Freq=ALL and Cname
is the name of a node-based component.
200
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Real (VALUE)
Imaginary (VALUE2)
FY
F0cos
-F0sin
FZ
-F0sin
-F0cos
where:
F0 is the amplitude of the force. For unbalance, the amplitude is equal to the mass times
the distance of the unbalance mass to the spin axis.
is the phase of the force, needed only when several such forces, each with a different
relative phase, are defined.
If the forces are caused by an unbalance mass, multiplication of the amplitude of the static forces (F)
by the square of the spin velocity is unnecessary. ANSYS performs the calculation automatically at each
frequency step.
Because the rotational velocity commands (OMEGA and CMOMEGA) define only the orientation of the
spin axis, a harmonic analysis using the SYNCHRO command requires that you define the frequency
of excitation (HARFRQ) instead. For example, if the frequency of excitation is f, then:
= 2f/RATIO
where:
is the new magnitude of the rotational velocity vector used to calculate the gyroscopic
matrices.
RATIO is the ratio between the frequency of excitation and the frequency of the rotational velocity of the structure, as specified via the SYNCHRO command. If no RATIO
value is specified, an unbalance force is assumed; in all other cases, a general rotating
force is assumed.
Example Analysis
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
201
7.2.3. Orbits
When a structure is rotating about an axis and undergoes vibration motion, the trajectory of a node
executed around the axis is generally an ellipse designated as a whirl orbit.
In a local coordinate system xyz where x is the spin axis, the ellipse at node I is defined by semi-major
axis A, semi-minor axis B, and phase (PSI), as shown:
e
s
Angle (PHI) defines the initial position of the node (at t = 0). To compare the phases
of two nodes of the structure, you can examine the sum + .
Values YMAX and ZMAX are the maximum displacements along y and z axes, respectively.
You can print out the A, B, PSI, PHI, YMAX, and ZMAX values via a PRORB (print orbits)
command. Angles are in degrees and within the range of -180 through +180. The position
vector of local axis y in the global coordinate system is printed out along with the elliptical orbit characteristics. You can also animate the orbit (ANHARM) for further examination. For a typical usage example of these commands, see Example Unbalance Harmonic
Analysis (p. 210).
To retrieve and store orbits characteristics as parameters, use the *GET command after issuing the
PRORB command.
Note
The gyroscopic effect is not included in the dynamics equations expressed in a rotating reference frame. Therefore, if the structure contains a part with large inertia - such as a large
disk - the results obtained in the rotating reference frame may not compare well with stationary reference frame results.
202
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
203
Rotating structure can be part of a stationary structure in an analysis model (such as a gas turbine engine rotor-stator assembly). The stationary structure
and supports (such as bearings) need not be
axisymmetric.
204
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
See the CORIOLIS command for the list of elements supported in the Rotating Reference Frame.
Natural Frequencies
Natural frequencies differ according to the reference frame type. In most cases, natural frequencies are
known in a stationary reference frame through analytical expressions or experiment, for example. ANSYS
therefore recommends using the stationary reference frame for modal analyses.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
205
Description
ANSYS Command(s)
Set parameters.
lx=8
dia=0.2
2.
ET,
R,
MP,...
N,
E,
3.
D,...
4.
Set the analysis type (modal in this case). Use the QRDAMP eigensolver, request 8 modes, and specify complex eigensolutions.
ANTYPE,MODAL
MODOPT,QRDAMP,8,,,ON
MXPAND,...
5.
206
CORIOLIS,ON,,,ON
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
6.
ANSYS Command(s)
For each load step, define the rotation velocity, and then solve. OMEGA,
SOLVE
7.
PLCAMP
PRCAMP
Spin(rd/s)
1
2
3
4
5
6
7
8
BW
FW
BW
FW
BW
FW
BW
FW
0.000
30000.000
6.207
6.207
24.750
24.750
55.461
55.461
98.248
98.248
4.639
8.305
18.547
33.027
41.735
73.701
74.337
129.852
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
207
208
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
ANSYS Command(s)
ET,
MP,
2.
Set parameters.
thick1 = 0.03810
thick = thick1
hthick = thick/2.0
radius = 0.09525
inradius = radius
- hthick
outradius = radius + hthick
length = 0.254
3.
CYLIND,...
NUMMRG,...
ESIZE,...
VMESH,...
NSEL,...
4.
Set the analysis type (modal in this case) and use the QRDAMP ANTYPE,MODAL
eigensolver. Get 20 complex eigensolutions.
MODOPT,QRDAMP,20
6.
Define the angular rotation velocity and the number of modes OMEGA,
expanded.
MXPAND,
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
209
Description
ANSYS Command(s)
5.
CORIOLIS,ON
7.
Solve.
SOLVE
Disks are not visible in the plot because they are MASS21 elements.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
211
212
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
ANSYS Command(s)
MP,EX,1,2.1e+11
MP,DENS,1,7800
MP,PRXY,1,0.3
2.
ET,
SECTYPE,
SECDATA,
R,
N,
3.
TYPE,1
SECNUM,1
E,1,2
EGEN,7,1,1
TYPE,2
REAL,3
E,2
REAL,6
E,7
CM,inSpool,ELEM
4.
TYPE,1
SECNUM,2
E,9,10
EGEN,4,1,10
TYPE,2
REAL,4
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
213
Description
ANSYS Command(s)
E,10
REAL,5
E,12
ESEL,u,,,inSpool
CM,outSpool,ELEM
ALLSEL
6.
TYPE,3
REAL,7
E,1,101
REAL,8
E,9,109
REAL,9
E,6,13
REAL,10
E,8,108
5.
D,...
7.
f0 = 70e-6
F,7,FY,f0
F,7,FZ,,-f0
8.
ANTYPE,HARMIC
Harmonic analysis
SYNCHRO,,inSpool
NSUBST,500
500 substeps
HARFRQ,,14000/60
KBC,1
Step loading
DMPRAT,0.01
Damping ratio is 1%
CMOMEGA,inSpool,100.
CMOMEGA,outSpool,150.
214
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
ANSYS Command(s)
SOLVE
Note
The rotational velocities (CMOMEGA) are
not applied in the usual way. Rather, ANSYS considers only their direction cosines
and the velocity ratio between spools. For
more information, see the documentation
for the SYNCHRO command.
9.
/POST26
NSOL,2,7,U,Y,UY
NSOL,3,7,U,Z,UZ
REALVAR,4,2,,,UYR
REALVAR,5,3,,,UZR
PROD,6,4,4,,UYR_2
PROD,7,5,5,,UZR_2
ADD,8,6,7,,UYR_2+UZR_2
SQRT,9,8,,,AMPL7
!
NSOL,2,12,U,Y,UY
NSOL,3,12,U,Z,UZ
REALVAR,4,2,,,UYR
REALVAR,5,3,,,UZR
PROD,6,4,4,,UYR_2
PROD,7,5,5,,UZR_2
ADD,8,6,7,,UYR_2+UZR_2
SQRT,10,8,,,AMPL12
!
/GROPT,LOGY,1
/YRANGE,1.e-7,1.e-3
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
215
Description
ANSYS Command(s)
PLVAR,9,10
10.
/POST1
SET,1,262
/VIEW,,1,1,1
PLORB
Third output.
/ESHAPE,1
/GLINE,,-1
No element outline.
PLNSOL,U,SUM
11.
ANHARM
You can obtain the two critical frequencies (at which the amplitudes are largest) via
PRCAMP with SLOPE = 1.0.
216
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Orbits are represented in different colors. Orbits from the inner spool appear in sky blue,
and from the outer spool in purple. Spool lines appear in dark blue.
The following demo is presented as an animated GIF. Please view online if you are reading the PDF version
of the help. Interface names and other components shown in the demo may differ from those in the released
product.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
217
The animation of the whirls shown here is the third output resulting from the example
harmonic analysis with unbalance.
218
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 8: Submodeling
Submodeling is a finite element technique that you can use to obtain more accurate results in a particular region of a model. A finite element mesh may be too coarse to produce satisfactory results in a
given region of interest. The results away from this region, however, may be satisfactory.
Reanalyzing the entire model using a greater mesh refinement in order to obtain more accurate results
in one particular region is time-consuming and costly. Instead, you can use submodeling to generate
an independent, more finely meshed model of only the region (submodel) of interest and then analyze
it.
The following submodeling topics are available:
8.1. Understanding Submodeling
8.2. Using Submodeling
8.3. Example Submodeling Analysis Input
8.4. Shell-to-Solid Submodels
8.5. Where to Find Examples
Submodeling of a pulley hub and spokes: (a) coarsely meshed model, and (b) finely
meshed submodel (shown superimposed over coarse model)
Submodeling is also known as the cut-boundary displacement method (also known as the specified
boundary displacement method). The cut boundary is the boundary of the submodel which represents
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
219
Submodeling
a cut through the coarse model. Displacements calculated on the cut boundary of the coarse model
are specified as boundary conditions for the submodel.
Submodeling is based on St. Venant's principle, which states that if an actual distribution of forces is
replaced by a statically equivalent system, the distribution of stress and strain is altered only near the
regions of load application. The principle implies that stress concentration effects are localized around
the concentration; therefore, if the boundaries of the submodel are far enough away from the stress
concentration, reasonably accurate results can be calculated in the submodel.
The program does not restrict submodeling to structural (stress) analyses only. Submodeling is effective
in other disciplines as well. For example, in a magnetic field analysis, you can use submodeling to calculate more accurate magnetic forces in a region of interest.
Aside from the obvious benefit of yielding more accurate results in a region of your model, the submodeling technique has other advantages:
It reduces, or even eliminates, the need for complicated transition regions in solid finite element
models.
It enables you to experiment with different designs for the region of interest (different fillet radii, for
example).
It helps you in demonstrating the adequacy of mesh refinements.
The following restrictions apply to submodeling:
It is valid only for solid elements and shell elements.
The principle behind submodeling assumes that the cut boundaries are far enough away from the
stress concentration region. You must verify that this assumption is adequately satisfied.
220
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Submodeling
5. Verify that the distance between the cut boundaries and the stress concentration is adequate.
The steps are explained in detail next.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
221
Submodeling
database after solving the coarse model, but before going on to create the submodel. To save the
database, use either of these methods:
Command(s): SAVE
GUI: Utility Menu> File> Save as
Utility Menu> File> Save as Jobname.db
Specify appropriate node rotations. Node rotation angles on cut boundary nodes should not be
changed after they have been written to the node file in interpolation step 1 (see Perform CutBoundary Interpolation (p. 223)). To specify node rotations, use one of these methods:
Command(s): NROTAT
GUI: Main Menu> Preprocessor> Modeling> Create> Nodes> Rotate Node CS> To Active CS
Main Menu> Preprocessor> Modeling> Move/Modify> Rotate Node CS> To Active CS
222
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Submodeling
Be aware that node rotation angles might be changed by application of nodal constraints (DSYM), by
transfer of line constraints (SFL), or by transfer of area constraints (SFA), as well as by more obvious
methods (NROTAT and NMODIF).
The presence or absence of node rotation angles in the coarse model has no effect upon the submodel.
Loads and boundary conditions for the submodel are covered in subsequent steps.
At this point, it is worthwhile to discuss temperature interpolation. In an analysis with temperaturedependent material properties, or in a thermal-stress analysis, the temperature distribution
must be the same between the coarse model and the submodel. For such cases, you must also
interpolate the temperatures from the coarse model to all nodes in the submodel. To do so,
select all submodel nodes and write them to a different file using NWRITE,Filename,Ext.
Specify a file name; otherwise, your file of cut-boundary nodes will be overwritten. Step 7 shows
the command to do temperature interpolation.
2. Restore the full set of nodes, write the database to Jobname.DB, and leave PREP7. You must write
the database to Jobname.DB because you need to continue with the submodel later.
To restore the full set of nodes, use one of these methods:
Command(s): ALLSEL
GUI: Utility Menu> Select> Everything
To write the database to Jobname.DB, use one of these methods:
Command(s): SAVE
GUI: Utility Menu> File> Save as Jobname.db
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
223
Submodeling
3. To do the cut-boundary interpolation (and the temperature interpolation), the database must contain
the geometry for the coarse model. Therefore, you must resume the database using one of the
methods shown below, making sure to identify the name of the coarse-model database file:
Command(s): RESUME
GUI: Utility Menu> File> Resume from
For example, if the job name for the coarse-model analysis was COARSE, issue the command
RESUME,COARSE,DB.
4. Enter POST1, which is the general postprocessor (/POST1 or menu path Main Menu> General
Postproc). Interpolation can only be performed in POST1.
5. Point to the coarse results file (FILE or menu path Main Menu> General Postproc> Data & File
Opts).
6. Read in the desired set of data from the results file (SET or menu path Main Menu> General
Postproc> Read Results> option).
7. Initiate cut-boundary interpolation. To do so, use one of these methods:
Command(s): CBDOF
GUI: Main Menu> General Postproc> Submodeling> Interpolate DOF
By default, the CBDOF command assumes that the cut-boundary nodes are on file Jobname.NODE. The program then calculates the cut-boundary degree-of-freedom values and
write them in the form of D commands to the file Jobname.CBDO.
To do temperature interpolation, use one of these methods, being sure to identify the name
of the file containing all submodel nodes:
Command(s): BFINT
GUI: Main Menu> General Postproc> Submodeling> Interp Body Forc
Interpolated temperatures are written in the form of BF commands to the file Jobname.BFIN.
If real and imaginary data are involved, steps 6 and 7 must be performed twice. Issue the SET
command to get the real data, followed by the interpolation step (CBDOF and/or BFINT). Issue
the SET command with the field set to 1 to get the imaginary data, and repeat the interpolation
step, this time writing the interpolated imaginary data to a different file.
For load-history-dependent problems, steps 6 and 7 can be performed multiple times at different
data sets or substeps. The cut-boundary degree-of-freedom values of each data set will reside
in the same file but in different data blocks, separated by a :CBnn label (where nn is the cumulative iteration number of the results set used).
8. All interpolation work is now done, so leave POST1 (FINISH) and restore the submodel database
(RESUME or menu path Utility Menu> File> Resume from). (Be sure to use the submodel database
job name.)
224
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Submodeling
Then define the appropriate analysis type (usually static) and analysis options.
To apply the cut-boundary degree-of-freedom constraints, simply read in the file of D commands (created
by CBDOF) using one of these methods (for example, /INPUT,,CBDO):
Command(s): /INPUT
GUI: Utility Menu> File> Read Input from
Similarly, to apply the interpolated temperatures, read in the file of BF commands (created by BFINT)
using one of these methods (for example, /INPUT,,BFIN):
Command(s): /INPUT
GUI: Utility Menu> File> Read Input from
If real and imaginary data are involved, first read in the file(s) containing the real data, specify whether
degree-of-freedom constraint values and/or nodal body force loads are to be accumulated, and then
read in the file containing the imaginary data.
Specify that degree-of-freedom constraint values are to be accumulated:
Command(s): DCUM,ADD
GUI: Main Menu> Preprocessor> Loads> Define Loads> Settings> Replace vs Add> Constraints
Main Menu> Solution> Define Loads> Settings> Replace vs Add> Constraints
Specify that nodal body force loads are to be accumulated:
Command(s): BFCUM,,ADD
GUI: Main Menu> Preprocessor> Loads> Define Loads> Settings> Replace vs Add> Nodal Body
Ld
Main Menu> Solution> Define Loads> Settings> Replace vs Add> Nodal Body Ld
Be sure to reset the DCUM and BFCUM commands to their default status before proceeding.
It is important that you duplicate on the submodel any other loads and boundary conditions that existed
on the coarse model. Examples are symmetry boundary conditions, surface loads, inertia forces (such
as gravity), concentrated force loads, etc. (see Figure 8.5: Loads on the Submodel (p. 226)).
Then specify load step options (such as output controls) and initiate solution calculations using one of
these methods:
Command(s): SOLVE
GUI: Main Menu> Solution> Solve> Current LS
Main Menu> Solution> Run FLOTRAN
For most problems, exit the solution processor (FINISH) after the solution has completed.
For load-history-dependent problems, apply the cut-boundary degree-of-freedom constraints in the
next data block in the file created by cut-boundary interpolation (CBDOF), and solve (SOLVE) as a new
load step. Because it is a new load step, for any rate-dependent materials and transient analysis, the
time value must match the time when cut-boundary interpolation occurred. This step can be repeated
until all cut-boundary degree-of-freedom constraints are applied. Then, exit the solution processor.
The overall data flow for submodeling (without temperature interpolation) is shown in Figure 8.6: Data
Flow Diagram for Submodeling (Without Temperature Interpolation) (p. 226).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
225
Submodeling
Figure 8.5: Loads on the Submodel
Figure 8.6: Data Flow Diagram for Submodeling (Without Temperature Interpolation)
8.2.5.Verify the Distance Between the Cut Boundaries and the Stress Concentration
The final step is to verify that the cut boundaries of the submodel are far enough away from the concentration. You can do this by comparing results (stresses, magnetic flux density, etc.) along the cut
boundaries with those along the corresponding locations of the coarse model. If the results are in good
agreement, it indicates that proper cut boundaries have been chosen; otherwise, you will need to recreate
and reanalyze the submodel with different cut boundaries further away from the region of interest.
An effective way to compare results is to obtain contour displays and path plots, as shown in Figure 8.7: Contour Plots to Compare Results (p. 227) and Figure 8.8: Path Plots to Compare Results (p. 227).
226
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
227
Submodeling
/SOLU
ANTYPE,...
...
D,...
DSYMM,...
ACEL,...
...
SOLVE
SAVE
FINISH
! Create submodel:
/CLEAR
/FILNAME,submod
/PREP7
...
...
...
! Enter SOLUTION
! Analysis type and analysis options
! Loads and load step options
! Perform cut-boundary
NSEL,...
!
NWRITE
!
ALLSEL
!
NWRITE,temps,node
!
!
SAVE
!
FINISH
interpolation:
Select nodes on cut boundaries
Write those nodes to submod.node
Restore full sets of all entities
Write all nodes to temps.node (for
temperature interpolation)
Submodel database file submod.db
RESUME,coarse,db
/POST1
FILE,coarse,rst
SET,...
CBDOF
BFINT,temps,node
FINISH
!
!
!
!
!
!
!
!
!
!
RESUME
/SOLU
ANTYPE,...
...
/INPUT,submod,cbdo
/INPUT,submod,bfin
DSYMM,...
ACEL,...
...
SOLVE
FINISH
/POST1
...
...
...
FINISH
! Enter POST1
! Submodel solution
228
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
SAVE
FINISH
! Create submodel:
/CLEAR
!
/FILNAME,submod
!
/PREP7
!
...
...
!
...
! Perform cut-boundary
NSEL,...
!
NWRITE
!
ALLSEL
!
NWRITE,temps,node
!
!
SAVE
!
FINISH
RESUME,coarse,db
!
/POST1
!
FILE,coarse,rst
!
SET,...
!
CBDOF
!
!
BFINT,temps,node
!
!
!
SET,...
!
CBDOF
!
!
BFINT,temps,node
!
!
!
...
FINISH
!
RESUME
!
/SOLU
!
ANTYPE,...
!
...
/INPUT,submod,cbdo
!
/INPUT,submod,bfin
!
DSYMM,...
!
ACEL,...
...
TIME
SOLVE
!
/INPUT,submod,cbdo
!
/INPUT,submod,bfin
!
TIME
SOLVE
..
FINISH
/POST1
!
...
...
!
...
FINISH
Following is a simple model of a plate with a hole in its middle, with applied displacements in the x
direction. To simulate load-history dependency, cut boundary conditions are done at different substeps
of the coarse-mesh analysis, then applied as different load steps in the submodeling analysis. The material is assumed to be rate-dependent anisotropic plastic with isotropic hardening.
The problem is solved via the coarse-mesh model, the submodel, and the fine-mesh model, as shown
in this figure:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
229
Submodeling
Figure 8.9: Coarse-Mesh Model, Submodel, and Fine-Mesh Model
Both the coarse-mesh model and the fine-mesh model are solved in 20 substeps. The fine-mesh model
yields a benchmark solution for comparison purposes.
The submodel analysis is performed in four solutions, as follows:
1. With cut-boundary conditions at only the last substep of the coarse-mesh analysis, solved in one
load step with 20 substeps.
2. With cut-boundary conditions at five substeps (that is, at every four substeps) of the coarse-mesh
analysis, solved in five load steps each with about four substeps.
3. With cut-boundary conditions at 10 substeps (that is, at every two substeps) of the coarse-mesh
analysis, solved in 10 load steps, each with about two substeps.
4. With cut-boundary conditions at all 20 substeps of the coarse-mesh analysis, solved in 20 load steps,
each with about two substeps.
The specified load steps and substeps make the size of the substep in submodel solutions close to
those in whole-model solutions.
The following figure shows the equivalent plastic strains at point P in each solution:
230
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
When boundary conditions are cut at more substeps, the results are closer to the fine mesh (the correct
solution).
The following figure shows the contour plots of equivalent plastic strain at the end of loading for each
of the four solutions:
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
231
Submodeling
Figure 8.11: Equivalent Plastic Strain Distributions in a Submodeling Analysis with Load-History
Dependency
232
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Shell-to-Solid Submodels
The results with 20 boundary-condition cuttings most resemble the fine-mesh model.
As the example indicates, boundary conditions should be cut from more substeps of the coarse solution
to reflect load-history dependency more accurately.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
233
Submodeling
Figure 8.12: 3-D Solid Submodel Superimposed on Coarse Shell Model
The procedure for shell-to-solid submodeling is essentially the same as that for solid-to-solid submodeling, with these exceptions:
Shell-to-solid submodeling is activated by setting KSHS to 1 on the CBDOF command (Main Menu>
General Postproc> Submodeling> Interpolate DOF) and the BFINT command (Main Menu>
General Postproc> Submodeling> Interp Body Forc). This feature is not applicable to offsets used
with SHELL181 (SECOFFSET), or SHELL281 (SECOFFSET).
Cut boundaries on the submodel are the end planes that are normal to the shell plane (see Figure 8.13: Node Rotations (p. 235)). Nodes on these cut boundaries are written to the node file (NWRITE)
(Main Menu> Preprocessor> Modeling> Create> Nodes> Write Node File).
To determine the degree-of-freedom values at a cut-boundary node (CBDOF), the program first projects
the node onto the nearest element in the shell plane. The degree-of-freedom values of this projected
point are then calculated by interpolation and assigned to the corresponding node. Interpolated
temperatures (BFINT) are calculated based on the average temperature at the midplane of the nearest
shell element.
Note
The nodes on the cut boundary must lie within a distance of 0.75 times the average
thickness of the nearest shell element, as shown in Figure 8.13: Node Rotations (p. 235).
That is, the submodel should be approximately centered on the coarse model.
In a structural analysis, only translational displacements are calculated for the cut-boundary nodes,
but their values are based on both the translations and rotations of the projected point. Also, the
node is rotated such that the nodal UY direction is always perpendicular to the shell plane, as shown
234
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Node rotations: (a) before CBDOF command, (b) after CBDOF command
The .CBDO file written by the CBDOF command consists of two blocks:
a block of NMODIF commands (indicating node rotation angles) and DDELE commands (to delete
UY constraints)
a block of D commands (to apply the interpolated degree-of-freedom values).
The two blocks are separated by a /EOF command and a :CBnn label (where nn is the cumulative iteration number of the results set used).
You must read in the .CBDO file in PREP7, because the NMODIF command is only valid in PREP7. To
do so, enter the preprocessor, then use one of these methods:
Command(s): /INPUT
GUI: Utility Menu> File> Read Input from
Because the two blocks of commands are separated by a /EOF command, it is necessary to read in the
.CBDO file twice. The second time you read in the file, use the LINE field on /INPUT ("Optional line
number or label" in the GUI) to instruct the program to read the file starting with the :CBnn label, as
shown below:
/PREP7! The .CBDO file must be read in PREP7
/INPUT,,cbdo
! Reads Jobname.cbdo up to the /EOF command
/INPUT,,cbdo,,:cb1 ! Reads same file from the label :cb1
235
Submodeling
The Mechanical APDL Verification Manual contains the following submodeling test case:
VM142 - Stress Concentration at a Hole in a Plate
236
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 9: Substructuring
Substructuring is a procedure that condenses a group of finite elements into one element represented
as a matrix. The single-matrix element is called a superelement. You can use a superelement in an analysis as you would any other element type. The only difference is that you first create the superelement
by performing a substructure generation analysis.
Substructuring is available in the ANSYS Multiphysics, the ANSYS Mechanical, and the ANSYS Structural
products.
The following substructuring topics are available:
9.1. Benefits of Substructuring
9.2. Using Substructuring
9.3. Sample Analysis Input
9.4.Top-Down Substructuring
9.5. Automatically Generating Superelements
9.6. Nested Superelements
9.7. Prestressed Substructures
9.8. Where to Find Examples
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
237
Substructuring
The following figure shows the data flow for a complete substructure analysis and some of the files
involved.
Figure 9.1: Applicable Solvers in a Typical Substructuring Analysis
!F
f x
!hF
f %%
ph #$#"
gp Desrps
XAdeteles otu@ srevlos ellive
dohtem noitutitsuskf
e xysxei
p "
hF
PPvF hvF
F wF
f
238
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
9.2.1.2. Applying Loads and Creating the Superelement Matrices
gives the jobname GEN to all files produced by the generation pass. The default jobname is file or
whatever name was specified while entering the program.
Element Types - Most element types can be used to generate a substructure. In general, the only restriction
is that elements within the superelement are assumed to be linear and cannot use Lagrange multipliers.
If you include bilinear elements, they will be treated as linear elements (in their initial state).
Caution
Coupled-field elements used in a direct method coupled-field analysis with load vector
coupling are not valid in a substructure analysis. Other elements in the same shape family
should be used instead. See the Coupled-Field Analysis Guide for details. In addition, elements
with Lagrange multipliers cannot be used in substructuring. These type of elements include
MPC184, CONTA171, CONTA172, CONTA173, CONTA174, CONTA175, CONTA176, CONTA177,
and CONTA178 with appropriate KEYOPT(2) setting, and elements PLANE182, PLANE183,
SOLID185, SOLID186, SOLID187, SOLID272, SOLID273, and SOLID285 when using KEYOPT(6)
> 0.
Material Properties - Define all necessary material properties. For example, if the mass matrix is to be
generated, density (DENS) (or mass in some form) must be defined; if the specific heat matrix is to be
generated, the specific heat (C) must be defined; and so on. Again, because a superelement is linear,
any nonlinear material properties will be ignored.
Model Generation - In the generation pass, you are mainly concerned with creating the superelement
portion of the model. The nonsuperelement portion, if any, is defined later in the use pass. However,
you should plan the modeling approach for both portions before you start building the model. In particular, decide on how you want to connect the superelement to the other elements. To ensure the
connection, use the same node numbers at the interface. (Other methods requiring less effort on your
part are discussed in the use pass section later in this chapter.)
Edge Outline - Adjust the edge outline used to plot the superelement in the use pass [/EDGE]. A smaller
angle will produce more edges.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
239
Substructuring
One approach might be to develop the entire model, save it on a named database file, and select only
the portion to be substructured for the generation pass. In the use pass then, you can RESUME (Utility
Menu> File> Resume from) from the named database file, unselect the portion that was substructured
and replace it with the superelement matrix.
Note
Restarting a substructuring analysis is valid only if the backsubstitution method is chosen.
You cannot restart a run if the full resolve option is selected using the SEOPT command.
Name of the superelement matrix file - Specify the name (Sename) to be assigned to the superelement
matrix file. The program will automatically assign the extension SUB, so the complete file name will be
Sename.SUB. The default is to use the jobname [/FILNAME]. To specify the name of the superelement
matrix file:
Command(s): SEOPT
GUI: Main Menu> Solution> Analysis Type> Analysis Options
Equation Solver - The SPARSE solver is the only solver available for the generation pass of the substructure
analysis. To specify an equation solver:
Command(s): EQSLV
GUI: Main Menu> Solution> Analysis Type> Analysis Options
Matrices to be generated - You can request generation of just the stiffness matrix (or conductivity matrix,
magnetic coefficient matrix, etc.); stiffness and mass matrices (or specific heat, etc.); or stiffness, mass,
and damping matrices. The mass matrix is required if the use pass is a structural dynamic analysis or if
you need to apply inertia loads in the use pass. For the thermal case, the specific heat matrix is required
only if the use pass is a transient thermal analysis. Similar considerations apply to other disciplines and
to the damping matrix. To make your request, use the SEOPT command as described above.
240
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
Matrices to be printed - This option allows you to print out the superelement matrices. You can request
listing of both matrices and load vectors, or just the load vectors. The default is not to print any matrices.
To print out the matrices, use the SEOPT command:
Expansion Pass Method - Allows you to select the expansion pass method you plan to use during subsequent expansion passes with this superelement. The backsubstitution method (default) saves the
factorized matrix files needed to perform a backsubstitution of the master DOF solution during the expansion pass. The full resolve method does not save any factorized matrix files. The factorized matrix
files are named Sename.LNxx for the sparse solver.
Note
Factorized matrix files can become very large as the problem size increases, but are not
needed if the full resolve method is chosen during the expansion pass.
During the expansion pass, the full resolve method reforms the elements used to create the superelement,
reassembles the global stiffness matrix, and applies the master DOF solution as displacement boundary
conditions internally. These displacement boundary conditions are deleted upon finishing the expansion
pass solution.
Note
You cannot restart a substructure analysis with the full resolve expansion pass method
chosen.
Mass matrix formulation - Applicable only if you want the mass matrix to be generated. You can choose
between the default formulation (which depends on the element type) and a lumped mass approximation.
We recommend the default formulation for most applications. However, for dynamic analyses involving
"skinny" structures, such as slender beams or very thin shells, the lumped mass approximation has been
shown to give better results. To specify a lumped mass approximation, use one of these methods:
Command(s): LUMPM
GUI: Main Menu> Solution> Analysis Type> Analysis Options
Modes to be used - For superelements being used in a subsequent dynamic analysis [ANTYPE,MODAL,
HARMONIC, or TRANSIENT], you may include mode shapes as extra degrees of freedom to obtain better
accuracy [CMSOPT]. See the chapter on Component Mode Synthesis for more information.
Define master degrees of freedom using one of these methods
Command(s): M
GUI: Main Menu> Solution> Master DOFs> User Selected> Define
In a substructure, master DOFs serve three purposes:
They serve as the interface between the superelement and other elements. Be sure to define master DOFs
at all nodes that connect to nonsuperelements, as shown in Figure 9.2: Example of a Substructuring Application (p. 238). All degrees of freedom at these nodes should be defined as master DOFs (Lab = ALL on
the M command). Master DOFs must be defined even if you plan to have no elements in the model other
than a superelement.
If the superelement is to be used in a dynamic analysis, master DOFs characterize the dynamic behavior
of the structure if the Component Mode Synthesis method [CMSOPT] is not used. See Modal Analysis in
the Structural Analysis Guide for guidelines.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
241
Substructuring
If constraints [D] or force loads [F] are to be applied in the use pass, master DOFs must be defined at
those locations with the M command.
If this superelement is to be transformed [SETRAN] later in the use pass or used in a large deflection
analysis [NLGEOM,ON], then all nodes that have master DOFs must have all six DOFs (UX, UY, UZ, ROTX,
ROTY, ROTZ) defined and active.
For large deflections, master DOFs are typically defined at the joints of the flexible body and are at the
nodes connected to a joint element (MPC184), another rigid or flexible body node, or ground. At least
two master DOFs must be defined for each substructure, as the average rotation of the superelement
is computed from the average rotation of its master DOF. If only one node is a joint node, then another
must be chosen at the free end. See the Multibody Analysis Guide for more details.
Apply loads on the model You can apply all types of loads in a substructure generation pass (see
Table 9.1: Loads Applicable in a Substructure Analysis (p. 243)), but keep in mind the following points:
The program will generate a load vector that includes the effect of all applied loads. One load vector per
load step is written to the superelement matrix file. This load vector is the equivalent load for the combined
loads in the load step. A maximum of 31 load vectors are allowed.
Nonzero DOF constraints can be used in the generation pass and will become part of the load vector. (In
the expansion pass, if the load step being expanded contains nonzero DOF constraints, the database must
have matching DOF values. If it does not, the DOF constraints must be respecified [D] in the expansion
pass.)
Application of constraints [D] or force loads [F] can be postponed until the use pass, but master DOF must
be defined at those locations with the M command or corresponding GUI path.
Note
If a mass matrix is generated, apply the degree of freedom constraints in the use pass at
the master DOF (defined in the generation pass) to ensure that all mass is accounted for
in the substructure. For analyses with acceleration loadings, the load should be applied
in the generation pass and used in the use pass for greater accuracy, rather than apply
the acceleration load on the reduced mass matrix.
Similarly, application of linear and angular accelerations can be postponed until the use pass, but only if
a mass matrix is generated. This is recommended if you plan to rotate the superelement in the use pass,
because load vector directions are "frozen" and rotated with the superelement.
The Maxwell force flag (MXWF label on the SF family of commands) is normally used in a magnetic analysis to flag element surfaces on which the magnetic force distribution is to be calculated. This flag has
no effect (and therefore should not be used) for a superelement in a magnetic analysis.
If you intend to create an imaginary force vector, you should generate it as a real load vector, then use
it as an imaginary load vector in the use pass (SFE,,,,KVAL = 2) and expansion pass (SEEXP,,,ImagKy =
ON).
242
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
For large rotation analyses - Do not apply constraints to the model in this pass. You will apply constraints
for large rotation analyses in the use pass.
Table 9.1: Loads Applicable in a Substructure Analysis
Load Name
Load Category
Commands[1]
Solid Model Loads
Forces
Surface
Loads
Body Loads
ACEL, DOMEGA
1. The menu path used to access each command in the GUI will vary depending on the engineering discipline
of the analysis (structural, magnetic, etc.). For a list of menu paths, see the description of individual
commands in the Command Reference.
Specify load step options
options (damping).
The only options valid for the substructure generation pass are dynamics
243
Substructuring
Command(s): SAVE
GUI: Utility Menu> File> Save as Jobname.db
Start solution calculations using one of these methods:
Command(s): SOLVE
GUI: Main Menu> Solution> Solve> Current LS Output from the solution consists of the superelement
matrix file, Sename.SUB, where Sename is the name assigned as an analysis option [SEOPT] or the
jobname [/FILNAME]. The matrix file includes a load vector calculated based on the applied loads. (The
load vector will be zero if no loads are defined.)
Repeat for additional load steps (that is, to generate additional load vectors) The load vectors
are numbered sequentially (starting from 1) and appended to the same superelement matrix file. See
Loading in the Basic Analysis Guide for other methods for multiple load steps.
Leave SOLUTION using one of these methods
Command(s): FINISH
GUI: Main Menu> Finish
Caution
If you are using the command input method to clear the database, additional commands
may not be stacked on the same line (using the $ delimiter) as the /CLEAR command.
Be sure to define a jobname that is different from the one used for the generation pass. This way, you
can ensure that no generation pass files will be overwritten. To define a jobname, use one of these
methods:
Command(s): /FILNAME
GUI: Utility Menu >File> Change Jobname
244
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
Interface nodes between superelement and nonsuperelement must exactly match the master
node locations.
5. Define the superelement by pointing to the proper element type reference number and reading in
the superelement matrix. To point to the element type:
Command(s): TYPE
GUI: Main Menu> Preprocessor> Modeling> Create> Elements> Elem Attributes
Now read in the superelement matrix using one of these methods (you may first need to use
other commands to modify the matrix, as explained below):
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
245
Substructuring
Command(s): SE
GUI: Main Menu> Preprocessor> Modeling> Create> Elements> Superelements> From .SUB
File
a. If there are no nonsuperelements in the model, or if there are nonsuperelements and the interface
nodes have the exact same node numbers as the master nodes on the superelement, then simply
read in the superelement using the SE command:
TYPE,...! Element type reference number
SE,GEN! Reads in superelement from file GEN.SUB
The Sename field on the SE command shown above identifies the name of the superelement
matrix file. The extension .SUB is assumed, so the complete file name is Sename.SUB
(GEN.SUB in the above example). The superelement is given the next available element
number.
b. If there are nonsuperelements in the model and the interface nodes have a constant node
number offset from the master nodes, you must first create a new superelement matrix with
new node numbers and then read in the new matrix.
To create a new superelement matrix, use one of these methods:
Command(s): SETRAN
GUI: Main Menu> Preprocessor> Modeling> Create> Elements> Superelements> By
CS Transfer
Main Menu> Preprocessor> Modeling> Create> Elements> Superelements> By Geom
Offset
To read in the new matrix, use one of these methods:
Command(s): SE
GUI: Main Menu> Preprocessor> Modeling> Create> Elements> Superelements> From
.SUB File
For example, given an existing superelement matrix file GEN.SUB and a node number offset
of 2000, the commands would be:
SETRAN,GEN,,2000,GEN2,SUB
TYPE,...
SE,GEN2
c. If there are nonsuperelements in the model and the interface nodes have no relationship with
the master nodes (as would be the case with automatically meshed models), first observe the
following caution.
Caution
It is quite likely that the node numbers of the master nodes from the generation
pass overlap with node numbers in the use pass model. In such cases, reading
in the superelement [SE] will cause existing use pass nodes to be overwritten by
the superelement's master nodes. To avoid overwriting existing nodes, use a
node number offset [SETRAN] before reading in the superelement. In any case,
save the database [SAVE] before issuing the SE command.
Thus you should first save the database [SAVE], use the SETRAN command to create a new
superelement matrix with a node number offset, and then use the SE command to read in
246
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
the new matrix. The CPINTF command (Main Menu> Preprocessor> Coupling/Ceqn>
Coincident Nodes) can then be used to connect the pairs of nodes at the interface. For
example, given a superelement matrix file called GEN.SUB:
*GET,MAXNOD,NODE,,NUM,MAX
SETRAN,GEN,,MAXNOD,GEN2,SUB
SE,GEN2
NSEL,...
CPINTF,ALL
!
!
!
!
!
!
!
NSEL,ALL
d. If the superelement is to be transformed - moved or copied to a different position, or symmetrically reflected - you must first use the SETRAN command or SESYMM command (Main Menu>
Preprocessor> Modeling> Create> Elements> Superelements> By Reflection), with the appropriate node number offsets, to create new superelement matrix files and then use SE to read
in the new matrices. Connecting the superelements to the nonsuperelements is done the same
way as above - by using common node numbers, a constant node number offset, or the CPINTF
command.
Note
If you use SETRAN to transfer the superelement to a different coordinate system,
the superelement's master nodes are rotated with it by default. This is typically useful
if the original superelement's master nodes are rotated, into a cylindrical system for
example. (In this case, the transfer does not effect the superelement stiffness matrix.)
If the original superelement has no rotated nodes, it is likely that the transferred
superelement will not need rotated nodes either. You can prevent node rotation in
such cases by setting the NOROT field on SETRAN to 1. (The superelement stiffness
matrix and load vector are modified by the program for this type of transfer.)
6. Verify the location of the superelement using graphics displays and listings. Superelements are
represented by an edge outline display, the data for which are written to the matrix file in the
generation pass. To produce a graphics display:
Command(s): EPLOT
GUI: Utility Menu> Plot> Elements
To produce a listing:
Command(s): SELIST
GUI: Utility Menu> List> Other> Superelem Data
7. Save the complete model database:
Command(s): SAVE
GUI: Utility Menu> File> Save as Jobname.db
Leave PREP7 using one of these methods:
Command(s): FINISH
GUI: Main Menu> Finish
247
Substructuring
type of analysis. You should, of course, have the appropriate matrices generated during the generation
pass. For example, if you intend to do a structural dynamic analysis, the mass matrix must be available.
The procedure is as follows:
1. Enter SOLUTION using one of these methods:
Command(s): /SOLU
GUI: Main Menu> Solution
2. Define the analysis type and analysis options.
For large rotation analyses - turn large deformation effects on [NLGEOM,ON], and define the
proper number of substeps for the nonlinear analysis.
3. Apply loads on the nonsuperelements. These may consist of DOF constraints and symmetry conditions
[D family of commands], force loads [F family], surface loads [SF family], body loads [BF family], and
inertia loads [ACEL, etc.]. Remember that inertia loads will affect the superelement only if its mass
matrix was generated in the generation pass.
Note
For large rotation analyses, be sure to apply the proper constraints in this step.
4. Apply superelement load vectors (if any) using one of these methods:
Command(s): SFE
GUI: Main Menu> Solution> Define Loads> Apply> Load Vector> For Superelement
One load vector per load step (created during the generation pass) is available on the superelement matrix file, and is identified by its reference number:
SFE,63,1,SELV,0,0.75
applies, on element number 63, load vector number 1, with the load applied as a real load and
with a scale factor of 0.75. Thus the ELEM field represents the element number of the superelement, LKEY represents the load vector number (default = 1), Lab is SELV, KVAL is for a real
or imaginary load vector, and VAL1 represents the scale factor (default = 0.0). (See the SFE
command description for more information.)
Note
The load vector orientation is fixed (frozen) to the superelement, so if the superelement is used in a rotated position, the load vector rotates with it. The same applies
to the degree of freedom directions (UX, UY, ROTY, etc.). They too are fixed to the
superelement and will rotate with the superelement if it is rotated (unless NOROT =
1 on the SETRAN command, in which case the nodal coordinate systems will not
be rotated).
5. Specify load step options that are appropriate for the analysis type. Use the EQSLV command to
select an appropriate equation solver based on the chosen analysis type and the physics of the
problem.
6. Initiate the solution:
248
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
Command(s): SOLVE
GUI: Main Menu> Solution> Solve> Current LS
Results from the solution consist of the complete solution for nonsuperelements and the reduced
solution - DOF solution at masters - for the superelements. The complete solution for nonsuperelements is written to the results file (Jobname.RST, RTH, or RMG), which you can postprocess
using normal postprocessing procedures.
The reduced solution is written to the file Jobname.DSUB. You can review this file using one
of these methods:
Command(s): SEDLIST
GUI: Main Menu> General Postproc> List Results> Superelem DOF
Utility Menu> List> Results> Superelem DOF Solu
To expand the reduced solution to all elements within the superelement, you will need to
perform the expansion pass, explained next.
7. Leave SOLUTION:
Command(s): FINISH
GUI: Main Menu> Finish
Note
The displacement boundary conditions are automatically applied internally at the master
degrees of freedom during the expansion pass solution and are automatically deleted when
the solution completes.
The expansion pass logic for substructuring analyses first searches for the superelement .LN22 file
and, if found, chooses the sparse solver to perform a backsubstitution (the EQSLV command is ignored).
Otherwise, the program will stop the expansion pass and give a message suggesting an alternate expansion method.
If the .LN22 file is not detected for the specified superelement, the full resolve method is chosen. The
PCG solver is chosen by default for the full resolve method. You can select the sparse solver using the
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
249
Substructuring
EQSLV command to override the default. Other equation solvers cannot be used with the full resolve
method.
1. Clear the database:
Command(s): /CLEAR
GUI: Utility Menu> File> Clear & Start New
This has the same effect as leaving and re-entering the program.
2. Change the jobname to what it was during the generation pass. This way, the program can easily
identify the files required for the expansion pass:
Command(s): /FILNAME
GUI: Utility Menu> File> Change Jobname
3. Restore the generation pass database:
Command(s): RESUME
GUI: Utility Menu> File> Resume Jobname.db
4. Enter SOLUTION using one of these methods:
Command(s): /SOLU
GUI: Main Menu> Solution
5. Activate the expansion pass and its options:
Command(s): EXPASS
GUI: Main Menu> Solution> Load Step Opts> ExpansionPass
Expansion pass on or off - Choose "on."
Name of superelement to be expanded - Specify the name (Sename):
Command(s): SEEXP
GUI: Main Menu> Solution> Load Step Opts> ExpansionPass> Expand Superelem
(The complete name of the file is assumed to be Sename.SUB.)
Name of the reduced solution file from use pass - Specify the name (Usefil) using the SEEXP command (or the menu path shown above). The complete name of the file is assumed to be Usefil.DSUB.
Real or imaginary component of displacement - Applicable only if the use pass was a harmonic analysis. Use the Imagky key on the SEEXP command (or the menu path shown above). If all solutions
are to be expanded (NUMEXP,ALL), Imagky is ignored and both the real and imaginary solutions
are expanded.
6. Identify the use pass solution(s) to be expanded. You can either expand a single solution [EXPSOL] or a
range of solutions (including all) [NUMEXP]:
Single Solution - use either the load step and substep numbers or the time (or frequency) to identify
a solution:
Command(s): EXPSOL
GUI: Main Menu> Solution> Load Step Opts> ExpansionPass> Single Expand> By Load Step
250
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Using Substructuring
Main Menu> Solution> Load Step Opts> ExpansionPass> Single Expand> By Time/Freq
Note
If the load step being expanded contains nonzero DOF constraints, the database must
have matching DOF values. If it does not, the DOF constraints must be respecified [D] in
the expansion pass.
Range of Solutions - Identify the number of solution and time or frequency range:
Command(s): NUMEXP
GUI: Main Menu> Solution> Load Step Opts> ExpansionPass> Range of Solu's
7. Specify load step options. The only options valid for a substructure expansion pass are output controls:
Output Controls - These options control printed output, database and results file output, and extrapolation of results.
If you want to include any results data on the printed output file (Jobname.OUT):
Command(s): OUTPR
GUI: Main Menu> Solution> Load Step Opts> Output Ctrls> Solu Printout
If you want to control the data on the results file (Jobname.RST):
Command(s): OUTRES
GUI: Main Menu> Solution> Load Step Opts> Output Ctrls> DB/Results File
If you want to review element integration point results by copying them to the nodes instead of
extrapolating them (default):
Command(s): ERESX
GUI: Main Menu> Solution> Load Step Opts> Output Ctrls> Integration Pt
8. Start expansion pass calculations:
Command(s): SOLVE
GUI: Main Menu> Solution> Solve> Current LS
9. Repeat steps 6 to 8 for additional use pass solutions to be expanded. If you need to expand the solution
for a different superelement, you will need to leave and re-enter SOLUTION.
Note
If the superelement to be expanded contains contact elements and has multiple use
passes, use separate expansions (repeating steps 1 through 7) or issue the NUMEXP
command to expand multiple load steps together.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
251
Substructuring
11. Postprocess results in the superelement using standard techniques.
Note
An expansion pass is not valid if the use pass was a PSD analysis.
252
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Top-Down Substructuring
!
FINISH
! ... Review results in nonsuperelements
!
EXPANSION PASS
/CLEAR
! Clear the database
/FILNAME,GEN
! Change jobname back to generation pass jobname
RESUME
! Restore generation pass database
/SOLU
! Enter SOLUTION
EXPASS,ON
! Activate expansion pass
SEEXP,GEN,USE...
! Superelement name to be expanded (GEN, unless SETRAN used)
--! Load step options (mainly output controls)
--SOLVE
! Initiate expansion pass solution. Full
!
superelement solution written to GEN.RST (or
!
RTH or RMG).
FINISH
! ... Review results in superelements
For more information, see the ANTYPE, SEOPT, M, ET, SETRAN, SE, CPINTF, EXPASS, and SEEXP
command descriptions.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
253
Substructuring
connecting nodes between the superelements, because they were all generated from a single
model.
Enter SOLUTION and define the analysis type and analysis options. Apply loads on the nonsuperelements, read in load vectors (if any), specify load step options, and initiate the use pass
solution.
4. Perform the expansion pass. Start by restoring the full model database, with all elements and nodes
active. Then expand each superelement separately, using the appropriate jobnames and exiting and
re-entering SOLUTION each time. You can then review the results in each superelement using normal
postprocessing procedures. Use of the full database, FULL.DB, allows the reading in of multiple
superelement results:
RESUME,FULL,DB
/POST1
FILE,GEN1
SET,...
FILE,GEN2
SET,...!Will not clear previous superelement results
A sample input for top-down substructuring follows. This example assumes a model with one superelement and other nonsuperelements.
!
Sample
!
!
BUILD THE FULL
!
/FILNAME,FULL
/TITLE,...
/PREP7
--------SAVE
FINISH
!
GENERATION PASS
!
/FILNAME,GEN
/SOLU
ANTYPE,SUBSTR
SEOPT,GEN,...
ESEL,...
NSLE
M,...
D,...
------SOLVE
--SOLVE
----FINISH
!
USE PASS
!
/CLEAR
/FILNAME,USE
RESUME,FULL,DB
ESEL,...
NSLE
254
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
FINISH
!
EXPANSION PASS
!
/CLEAR
! Clear database for expansion pass
/FILNAME,GEN
! Change jobname back to generation pass jobname
RESUME,FULL,DB
! Restore full model database
/SOLU
ANTYPE,SUBSTR
EXPASS,ON
EXPSOL,...
SEEXP,GEN,USE,...
----SOLVE
! Enter SOLUTION
!
!
!
!
FINISH
! ... Review results in superelement
Please see the ANTYPE, SEOPT, M, ET, SE, EXPASS, and SEEXP command descriptions for more information.
If using the bottom-up substructuring method first create the part of the model that will become superelements. If using the top-down substructuring method, first create whole model, then select the part
of the model that will become superelements.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
255
Substructuring
2.
Perform the generation pass using SEOPT and any other /SOLU commands to define any necessary
options for the substructuring analysis.
3.
Use SEGEN to define the options for the automatic superelement generation process. If stopStage
= PREVIEW is selected, then the model is only broken into domains (superelements). No reduced matrices
are created and the superelements (.SUB files) are not actually created. You can then graphically
(visually) preview each domain by using /PNUM,DOMAIN. By default, master DOFs are automatically
defined at each of the following locations: all DOFs on the interfaces between each superelement, all
DOFs associated with contact elements (TARGE169 to CONTA177), and at all DOFs associated with nodes
having a point load defined. The option to manually define the master DOF only makes sense AFTER a
'preview pass' has been made, as the exact number of superelements and the superelement boundaries
for each superelement cannot be known until the process is completed at least once.
Note
Due to the heuristics in the automatic domain decomposer, which is used to create the
domains that will become superelements, the number of defined superelements may
exceed the number of requested superelements.
After completing a preview pass, you can then add master DOFs or remove master DOFs that were
automatically defined during the preview pass. At least one master DOF must be defined for each
superelement. Then set stopStage = GEN, and if any master DOFs were added or removed, set
mDof = YES, and solve the model.
4.
Use SOLVE to either preview or generate the automatically created superelements. Note that multiple
load steps are not supported with automatic superelement generation.
256
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
257
Substructuring
Mechanical APDL Verification Manual does not present them as step-by-step examples with lengthy
data-input instructions and printouts; however, most users with at least limited finite-element experience
should be able to fill in the missing details by reviewing each test case's finite element model and input
data with accompanying comments.
The Mechanical APDL Verification Manual contains the following substructuring cases:
VM125
VM141
VM147
VM153
258
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
259
Free (CMSOPT,FREE)
Interface nodes remain free durInterface nodes remain free during the CMS superelement gener- ing the CMS superelement generation pass.
ation pass.
For more information, see the discussion of component mode synthesis theory and methods in the
Mechanical APDL Theory Reference.
260
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
261
CMS Superelement
Matrix
Specifying the CMS method When specifying the CMS method, also specify the number of modes
and, optionally, the frequency range used to generate the superelement. The program supports the
fixed-interface (CMSOPT,FIX), free-interface (CMSOPT,FREE), and residual-flexible free-interface
(CMSOPT,RFFB) CMS methods. If using the free-interface method, also specify the rigid body modes
(CMSOPT,,,,FBDDEF). If using the residual-flexible free-interface method, specify pseudo-constraints
(D,,,SUPPORT).
Naming the superelement matrix file The program assigns the .sub extension to the superelement matrix file name that you specify (SEOPT,Sename); therefore, the complete file name is Sename.SUB. The default file name is the Jobname (/FILNAME).
Specifying the lumped mass matrix formulation Specify the lumped mass matrix formulation
(LUMPM) if necessary. For most applications, ANSYS, Inc. recommends the default formulation (depending upon the element type); however, for dynamic analyses involving "skinny" structures such
as slender beams or very thin shells, the lumped mass approximation typically yields better results.
Defining master DOFs In a substructure, master degrees of freedom (DOFs) serve as the interface
between the superelements or other elements. Define master DOFs (M) at all nodes that connect
to non-superelements (Lab1 = ALL), as shown in Example of a Substructuring Application. You must
define master DOFs even if you intend to have no elements in the model other than a superelement.
If this superelement is to be transformed (SETRAN) later in the use pass or used in a large deflection
analysis (NLGEOM,ON), all nodes that have master DOFs must have all six DOFs (UX, UY, UZ, ROTX,
ROTY, ROTZ) defined and active.
For large deflections, master DOFs are typically defined at the joints of the flexible body and are at
the nodes connected to a joint element (MPC184), another rigid or flexible body node, or ground.
At least two master DOFs must be defined for each substructure, as the average rotation of the su-
262
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
263
264
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
m IFH
mm S
htdiw
m SPHFH
suidr
m SQHFH
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
265
Q tr
P tr
I tr
efens
266
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
267
268
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
! g.
! All the active DOFs (that is, on the nodes which belong to "interface1")
! are set as masters
m,all,all
! h.
! Selects all the nodes attached to the selected elements
! (that is, elements which belong to "part1")
nsle
! i.
! solve the first CMS generation pass
solve
finish
! j.
! Save the generation pass 1 database
save
! Repeat the generation pass for "part2"
! Generation pass 2
/filnam,part2
/solu
antype,substr
seopt,part2,2
cmsopt,fix,10
cmsel,s,part2
cmsel,s,interface2
m,all,all
nsle
solve
finish
save
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
269
! c.
! A superelement element type is created
/prep7
et,1,matrix50
! d.
! Element type attribute pointer set to 1
type,1
! e.
! Brings in the three superelements created above
se,part1
se,part2
se,part3
finish
! f.
! A modal analysis is performed
/solu
antype,modal
! g.
! Specifies modal analysis options
modopt,lanb,10
! h.
! Expands 10 modes
mxpand,10
! i.
! Solve the modal analysis
solve
finish
! STEP #5 (a. through g.)
! Expansion pass
! Expansion pass 1
! a.
! Clears the database
/clear,nostart
! b.
! Changes the jobname to superelement 1 name
/filnam,part1
! c.
! resume the database
resume
! d.
! Specifies the expansion pass
/solu
expass,on
! e.
! Specifies superelement name and use pass name
seexp,part1,use
270
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
! f.
! Specifies the loadstep and substep to be expanded
expsol,1,4
! g.
! Solve the first expansion pass
solve
finish
! Repeat the expansion pass for "part2"
! Expansion pass 2
/clear,nostart
/filnam,part2
resume
/solu
expass,on
seexp,part2,use
expsol,1,4
solve
finish
! Repeat the expansion pass for "part3"
! Expansion pass 3
/clear,nostart
/filnam,part3
resume
/solu
expass,on
seexp,part3,use
expsol,1,4
solve
finish
/post1
cmsfile,add,part1,rst
cmsfile,add,part2,rst
cmsfile,add,part3,rst
! c.
! Reads the first data set
set,first
! d.
! Plots the displacement contour in the x direction
271
Description
Command(s)
1.
2.
3.
/INPUT,CMS_SAMPLE.INP
/FILNAME,PART1
/SOLU
ANTYPE,SUBSTR
SEOPT,PART1,2
CMSOPT,FIX,10
CMSEL
CMSEL,S,INTERFACE1
M,ALL,ALL
NSLE
SOLVE
SAVE
---As coded in the input file, generation passes for the remaining
parts PART2 and PART3 occur here. Steps a through j are repeated for PART2 and again for PART3. (The Jobname and
superelement matrix file name change accordingly.) Also, for
passes 2 and 3, the node component is INTERFACE2 and INTERFACE3, respectively.
4.
/CLEAR,NOSTART
/FILNAME,USE
/PREP7
ET,1,MATRIX50
TYPE,1
SE,PART1
SE,PART2
SE,PART3
FINISH
/SOLU
ANTYPE,MODAL
272
MODOPT,LANB,10
MXPAND,10
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
Command(s)
SOLVE
FINISH
5.
/CLEAR,NOSTART
/FILNAME,PART1
RESUME
/SOLU
EXPASS,ON
SEEXP,PART1,USE
EXPSOL,1,4
SOLVE
----
FINISH
/CLEAR,NOSTART
/POST1
CMSFILE,ADD,PART1,RST
CMSFILE,ADD,PART2,RST
CMSFILE,ADD,PART3,RST
SET,FIRST
PLNSOL,U,X
----
FINISH
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
273
274
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
275
Description
Command(s)
1.
...
2.
...
3.
/FILNAME,PART1
/SOLU
ANTYPE,SUBSTR
SEOPT,PART1,2
CMSOPT,FREE,10,,,FNUM,3
CMSEL
CMSEL,S,INTERFACE1
M,ALL,ALL
NSLE
SOLVE
SAVE
---As coded in the input file, generation passes for the remaining
parts PART2 and PART3 occur here. Steps a through j are repeated for PART2 and again for PART3. (The Jobname and
superelement matrix file name change accordingly.) Also, for
passes 2 and 3, the node component is INTERFACE2 and INTERFACE3, respectively.
4.
...
5.
...
6.
...
276
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
/solu
antype,substr
! c.
! Specifies the name to be assigned to superelement matrix file
! Strongly suggested to be the same as the active jobname
seopt,part1,2
! d.
! Specifies CMS options
cmsopt,rffb,10
! e.
! Selects element component named "part1"
cmsel,s,part1
! f.
! Selects node component named "interface1"
cmsel,s,interface1
! g.
! All the active DOFs (that is, on the nodes which belong to "interface1")
! are set as masters
m,all,all
! h.
! Selects all the nodes attached to the selected elements
! (that is, elements which belong to "part1")
nsle
! i.
! Specify only the minimum number of displacement constraints necessary
! to prevent rigid body motion: three constraints (or fewer, depending
! on the element type) for 2-D models and six (or fewer) for 3-D models.
d,430,all,support
d,440,ux,support
! j.
! solve the first CMS generation pass
solve
finish
! k.
! Save the generation pass 1 database
save
! Repeat the generation pass for "part2"
! Generation pass 2
/filnam,part2
/solu
antype,substr
seopt,part2,2
cmsopt,rffb,10
cmsel,s,part2
cmsel,s,interface2
m,all,all
nsle
d,705,all,support
d,715,ux,support
solve
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
277
Description
Command(s)
1.
...
2.
...
3.
/FILNAME,PART1
/SOLU
ANTYPE,SUBSTR
SEOPT,PART1,2
CMSOPT,RFFB,10
CMSEL
CMSEL,S,INTERFACE1
M,ALL,ALL
NSLE
i. Specify pseudo-constraints.
D,,,SUPPORT
278
SOLVE
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
Command(s)
SAVE
---As coded in the input file, generation passes for the remaining
parts PART2 and PART3 occur here. Steps a through j are repeated for PART2 and again for PART3. (The Jobname and
superelement matrix file name change accordingly.) Also, for
passes 2 and 3, the node component is INTERFACE2 and INTERFACE3, respectively.
4.
...
5.
...
6.
...
279
/clear,nostart
/filnam,use
/prep7
et,1,matrix50
type,1
! Define the three superelements to use in the model
se,part1
se,part2
*get,nmax,node,,num,max
sesymm,part2,x,nmax,part2sym,sub
se,part2sym
cpintf,all,0.001
finish
/solu
antype,modal
modopt,lanb,10
mxpand,10
solve
finish
! STEP #5
! Expansion pass
! Expansion pass 1
/clear,nostart
/filnam,part1
resume
/solu
expass,on
seexp,part1,use
expsol,1,4
solve
finish
! Expansion pass 2
/clear,nostart
/filnam,part2
resume
/solu
expass,on
seexp,part2,use
expsol,1,4
solve
finish
! Obtain the third part of the model from PART2
! Expand the solution in the transformed location
/assign,rst,part2sym,rst
/solu
expass,on
seexp,part2sym,use,,on
!
Offset node and element IDs in the new superelement
rstoff,node,10000
rstoff,elem,10000
expsol,1,4
solve
finish
! STEP #6
! Reads results for "load step 1 - substep 4"
280
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Description
Command(s)
1.
...
2.
...
3.
4.
SE,PART1
----
*GET,NMAX,NODE,,NUM,MAX
SESYMM,PART2,X,NMAX,PART2SYM,SUB
SE,PART2
SE,PART2SYM
CPINTF,ALL,0.001
FINISH
...
5.
/ASSIGN,RST,PART2SYM,RST
/SOLU
EXPASS,ON
SEEXP,PART2SYM,USE,,ON
RSTOFF,NODE,10000
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
281
Description
Command(s)
RSTOFF,ELEM,10000
EXPSOL,1,4
SOLVE
FINISH
6.
/CLEAR,NOSTART
/POST1
CMSFILE,ADD,PART1,RST
CMSFILE,ADD,PART2,RST
CMSFILE,ADD,PART2SYM,RST
...
...
282
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
283
284
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
285
286
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Exporting to ADAMS
The Select Interface Points dialog box appears first. From this dialog box, you must select two or more
interface points.
Do not select too many interface points, as one point gives rise to six degrees of freedom in ADAMS.
Too many interface points may lead to huge files and models.
After you confirm your selection, the Export to ADAMS dialog box appears.
Figure 11.2: Export to ADAMS Dialog Box
1 meter/millimeter
=
1000
Mass Factor =
1 kilogram/tonne
=
0.001
Force Factor
=
1 newton/newton
=1
Time Factor =
1 second/second
=1
2. Number of Modes to Extract: Input the number of normal modes to compute. Normal modes are the
eigenmodes of the component with all degrees of freedom of all interface points fixed. The number of
normal modes depends on the frequency range of the excitation you will apply in your ADAMS model.
You must select a sufficient number of modes to represent your structure in that frequency range. In
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
287
Note
Note that the algorithm used to compute the .MNF file adds constraints to the interface
points. If you create the .MNF file a second time using the same model in the same run, be
sure to delete all constraints on the interface points (or resume the database file Jobname.DBMNF) before you run it again.
!
!
!
!
!
See the ADAMS command description for more information. When you use command input to compute
the .MNF file, there is no option to change the file name. The default name of Jobname.MNF will be
used.
288
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
289
290
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
291
Note
If you use the import procedure a second time with the Add weak springs option, additional
weak springs will be added to the model. This will have only a small influence on the results.
!
!
!
!
!
!
!
!
Every subsequent call of the WSPRINGS command will apply weak springs. Therefore, this command
may be omitted when importing new loads.
292
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
The MSR toolkit features described here are not supported by the ANSYS-ADAMS Interface.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
293
294
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Link 3
U1
Input motion
Crank
U4
Link 2
U3
Link 1: Output
motion
U2
The following are dimensions and properties for the Link3 component.
Radius of holes (radh) = 6mm
Width of rectangular rod (width) = 25mm
Thickness of rectangular rod (thick) = 10mm
Length of rectangular rod (length) = 300mm + 4*Radius of holes = 324mm
Young's modulus for rod = 7.22 x 104 MPa
Poisson's ratio for rod = 0.34
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
295
296
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
At this point you may import the adamsout.mnf file into your ADAMS model and perform a rigidbody dynamics simulation. The ADAMS model should consist of the components shown in Figure 11.5: Linkage Assembly (p. 295). After the simulation is done, export the loads acting on the Link3
component at five arbitrary `time steps. Name the load file loads.lod.
Once you have exported the load file, you can perform a stress analysis for Link3 in Mechanical APDL
using the command input shown below.
RESUME,adamsout,db
! Resume model
/FILNAM,adamsin
! Change jobname
/TITLE,Import loads from ADAMS
! Change title
!
WSPRINGS
! Create weak springs
!
! Enter Solution and solve all load steps
/SOLU
/INPUT,loads,lod
! Read in 5 load steps written by ADAMS
*DO,i,1,5
! Use a do loop to solve each load step
LSREAD,i
! Read in load step
IRLF,1
! Activate inertia relief
SOLVE
! Solve current load step
*ENDDO
!
/POST1
! Enter the general postprocesser
! Write deformation and equivalent stress to graphics file
/VIEW,1,1,1,1
/AUTO,1
EPLOT
/TYPE,1,4
/SHOW,
EPLOT
*DO,i,1,5
SET,i
PLNSOL,u,sum
PLNSOL,s,eqv
*ENDDO
/SHOW,term
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
297
298
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
299
300
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
301
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
!
302
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Note
You can also use the ANSYS Parametric Design Language to write a macro to perform such
an operation. See the ANSYS Parametric Design Language Guide for more information.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
303
304
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Caution
By linking in your own Fortran routines, you are creating a custom, site-specific version of
the ANSYS program. When you use UPFs, you are using ANSYS in a nonstandard way, one
that ANSYS, Inc. verification testing does not cover. You are responsible for verifying that the
results produced are accurate and that the routines you link to ANSYS do not adversely affect
other, standard areas of the program.
Exercise care when using UPFs on parallel systems. Do not use the /CONFIG command or a
config145.ans file to activate parallelism on a system with UPFs.
The following topics concerning UPFs are available:
13.1.1. Understanding UPFs
13.1.2.Types of UPFs Available
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
305
PARAB now becomes a valid ANSYS command that simply calls the user routine USER01. You can call
up to ten such user routines as commands. By including /UCMD commands in your start-up file
(start145.ans), you can make the user routines available in all of your ANSYS sessions.
Allows you to define your own element type. You can add it to the element
library and use it as you would any other element. See User-Defined Elements
in the Element Reference, and Creating a New Element in the Guide to UserProgrammable Features.
User-defined materials
Allows you to define your own material model. See User-Defined Material
Model, and Subroutines for Customizing Material Behavior in the Guide to UserProgrammable Features.
306
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The contact elements (CONTA171 to CONTA178) allow you to define your own
interfacial behaviors between surfaces. See Subroutines for Customizing Contact
Interfacial Behavior in the Guide to User-Programmable Features.
User thickness
User stresses
Allows you to calculate plastic strains and form the tangent stress-strain matrix
at an integration point based on your own plasticity law.
If you need to account for swelling in an analysis (due to neutron bombardment, for example), you must write the appropriate swelling law as a user
routine. No built-in swelling laws are available.
Available for the layered elements (such as SOLID185 Layered Structural Solid,
SOLID186 Layered Structural Solid, and SHELL281). Up to nine user-defined
failure criteria can be supplied.
User viscosity
User loads
ANSYS as a subroutine
You can call the entire ANSYS program as a subroutine in your own program,
such as a user-written design optimization algorithm.
Allows you to evaluate results and perform any desired calculations during
solution.
USRSURF116
Allows you to modify SURF151 and SURF152 film coefficients and bulk temperatures based on information from FLUID116.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
307
308
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
309
! node 1022
! ux
spmwrite,modal,10,inputTab,inputLab,outputTab,outputLab,1
310
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
311
The state-space matrices file can be imported in Simplorer using: Simplorer Circuit > SubCircuit >
Add Mechanical Component
A transient analysis is performed in Simplorer for the 3 following reduced models:
All 32 modes are used.
First 3 modes are used.
First 3 modes plus residual vector are used (ANSYS input is listed above).
The schematic of the simulation is shown on the figure below:
312
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
A zoom on the graph showing the evolution of the spring force is given below:
Figure 14.2: Evolution of Spring Force
The reduced model based on the first 3 modes (red curve) is not accurate. When the residual vector is
included (purple curve), the spring force is similar to the reference given by the 32 modes reduced
model (dotted curve).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
313
314
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
315
Soil-Pile-Structure Analysis
American Petroleum Institute (API) Recommended Practice for Planning, Designing and Constructing Fixed
Offshore Platforms - Working Stress Design, RP2A-WSD, Twenty-first Edition, December 2000 (the Seventeenth
and Twentieth editions are also supported).
American Petroleum Institute (API) Recommended Practice for Planning, Designing and Constructing Fixed
Offshore Platforms - Load and Resistance Factor Design, RP2A-LRFD, First Edition, July 1, 1993.
International Standards Organization (ISO) 19902, Petroleum and Natural Gas Industries - Fixed Steel Offshore
Structures, First Edition, December 2007.
In order to calculate these unity checks, yield data must be supplied for the piles. For WSD checks, extreme environmental loads may be specified.
For a transient analysis, the pile calculation should be calculated at each time step, in order for the
updated loading (in the form of prescribed displacements) to be applied to the pile and a new stiffness
and loading to be calculated and applied to the pile cap. In this case FDELE should be used to ensure
that previously calculated nodal forces are not appended to.
Note
Multiple loadsteps and small time increments are necessary to allow for the nonlinear iterations
required for the soil-pile interaction equations, and to enable these calculations to converge.
316
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
The springs utilized for the soil model are characterized by a nonlinear force-deflection relationship of
the type shown in Figure 15.2: Soil/load Deflection Characteristics (p. 318). These are commonly known
as P-Y curves for lateral behavior and T-Z curves for axial behavior.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
317
Soil-Pile-Structure Analysis
Figure 15.2: Soil/load Deflection Characteristics
The nonlinear nature of the soil properties means an iterative solution technique is required. Each iteration takes an assumed, or calculated, value for the soil spring stiffness based upon the previous iteration.
Two solution methods exist for determining the representative soil stiffness from the defined P-Y and
T-Z curves, namely Tangent and Secant stiffness. These are shown diagramatically below:
For both types of soil representation, the curves are stored as a series of points as defined by the user
and depicted by crosses in the diagrams above. In order that the soil stiffness may be computed, the
soil properties are assumed to vary linearly between these points as shown above. For points beyond
the last y datum defined, linear extrapolation is utilized from the last two points supplied.
In general, the tangent stiffness method will converge faster than the secant stiffness approach, especially
as the pile approaches its limiting capacity. The tangent stiffness method is adopted by default. The
default may be overridden by the STIFF option in the PILEGEN command.
There is no limit to the number of points on any curve. For displacements that occur outside of the
range specified for a given curve, linear extrapolation of the last segment of the curve is assumed. If a
constant value is required, two points with similar P or T and different displacement must be input.
318
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Both symmetric and non-symmetric curves may be defined for soil layer data. If symmetric curves are
required, only that part of the curve for positive displacement values need be input. The program will
automatically generate the remaining information.
If any values are defined for a negative displacement, then the curve is assumed to be non-symmetric.
Care should be taken in providing non-symmetric P-Y data since the soil stiffnesses derived from the
curves are based upon local axes displacements, and these may vary from pile to pile and from iteration
to iteration. Shifted P-Y curves are not permitted, since this can result in undesired lateral deflections.
Use the SLID definition if shifted curves are required.
For T-Z and ENDB data, a positive local displacement for the purposes of soil stiffness formulation is
taken as being defined by the vector going from the pile cap to the pile tip. Thus, zero tensile stiffness
for end bearing forces may be modelled by supplying an ENDB definition with zero stiffness for negative
displacements.
The soil properties must be defined down to the full depth of the pile, or to a greater depth.
If both top and bottom layer depths are supplied on a P-Y or T-Z header, the data is taken as a constant
between these depths. If only one depth is supplied, the data is defining the properties at one depth
in the soil medium. P-Y and T-Z curves, either explicitly defined or generated from the soil properties,
are assumed to vary linearly between these depths.
If a sudden change in soil properties is required at a given depth, one of two options are available:
1. If one or both of the soil definitions represent a constant stratum, the given depth may be supplied for
both layers.
2. If both the soil definitions are single point, they must be separated by a finite distance so that the program
can identify which layer is uppermost. Provided the separation is less than the coordinate tolerance (0.1
x pile diameter) the program will utilize the data as though they represented coincident layers. The
higher of the two levels specified will be adopted as the point of the discontinuity. See also Automatic
Pile Subdivision (p. 316).
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
319
Soil-Pile-Structure Analysis
The procedures used for developing P-Y and T-Z curves from user defined soil properties are as indicated
in the list below. Detailed descriptions are not given but a typical curve for each procedure is shown.
Note
Overburden pressures used in the soil curve generation are computed from the soil densities
provided down the soil profile. If explicit P-Y and/or T-Z data is provided for any soil layer
within a soil profile, then overburden pressures will be computed based upon the soil
density local to the point of calculation.
The factor used for shaft friction is automatically calculated as specified in Clause 6.4.2 of
the API 21st edition. A limiting value of is set to 3.0. However no allowance is made for
pile length.
The resulting soil curves depend upon the pile to which it is to be applied since, in the
general form, pressures are generated. A different soil curve set will thus be produced for
each pile in an analysis. Where stepped piles are utilized, it is important that two soil definitions are provided at the step position(s) in order that the correct geometric data is utilized
for the curve generation.
Procedures for Developing P-Y and T-Z Curves
For sands, the method suggested by Reese, et. al. [1] is used.
Figure 15.3: P-Y Curve for Sand
For clays, the method suggested by Matlock [2], a modified Matlock procedure, and the basic recommendations as listed in API 15 [3] are available.
320
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
For T-Z and ENDB curves, the procedures recommended by Vijayvergiya [4] have been adopted.
For the purposes of the implementation in the ANSYS soil-pile analysis, the ultimate skin friction and
end bearing pressure that can be developed in cohesionless soils (sand) is limited to the values given
in the API RP2A code of practice. The user should also note that for plugged pile conditions it is assumed that any internal soil skin friction is sufficient to sustain the plug in position; there is no internal
check undertaken to check this requirement.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
321
Soil-Pile-Structure Analysis
Figure 15.6: T-Z Curve for Clay and Sand
Note
For P-Y, T-Z and ENDB curves, the curved section is divided into four divisions by default.
For T-Z curves, Zc is taken to be 0.01 x pile diameter.
For ENDB curves, Zc is taken to be 0.05 x pile diameter.
Whenever the soil is defined as banded, P-Y and T-Z curves are evaluated at the mid-point
of the band.
15.1.4. Mudslides
The soil is assumed to be initially in an undisturbed state and all pile displacements are referred to as
absolute soil origin. There are instances, however, where displacement of the soil itself may occur due
322
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
323
Soil-Pile-Structure Analysis
option in the PILEGEN command can be used to define extreme loadcases, which will affect the design
axial pile capacity computed because a different resistance factor is used.
The following utilization checks are undertaken when requested:
where
fa = axial stress
Fxc = inelastic local buckling stress
fby, fbz = bending stresses (including second order effects using the PDELTA option in the PILEGEN
command)
Fb = allowable bending stress
A one-third increase in the allowable stresses is included if an extreme loadcase is being processed.
where
!
" "
'%&
&$ #
where
324
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
=
#!
"! ! ! $ $
% %
=
where
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
325
Soil-Pile-Structure Analysis
It is assumed that the soil data file is available in the working directory. If this is not the case, the /COPY
command can be used to copy the data file into the working directory and rename it if necessary. For
example, the following command will copy the soil data file soil_1.txt from folder C:\soildata
to the working directory and rename it to be file_soil.txt.
/COPY,'C:\soildata\soil_1',txt,,file_soil
This data defines the inelastic and elastic properties and layering details for each soil profile required.
The data is divided into various categories, their use depending on the type of analysis being carried
out and the input requirements of the user. A complete soil profile is shown below followed by detailed
descriptions of each of the categories in turn.
326
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Parameters
MUDD
Keyword to identify this as mudline data.
muddat
Z coordinate of mudline datum relative to the global axes system. If omitted, the mudline is assumed
to pass through the global origin.
The MUDD command is used for defining a global mudline which does not coincide with the global
origin. This can be particularly useful if the structure has been defined with its origin at somewhere
other than the mudline and the user wishes to retain the same coordinate system.
If a particular soil profile requires a local scour to be modelled, both the P-Y and T-Z data should start
at a depth corresponding to the scour level. Dummy data does not have to be provided over this depth.
327
Soil-Pile-Structure Analysis
1. If the data supplied to the user are not in the correct form, (for example, pressure curves instead of force
per unit length or the wrong units such as N and mm instead of KN and m) then factors may be applied
to convert from one to the other.
2. If the data supplied is for a given diameter pile and a different diameter pile is to be analyzed, a factor
may be applied to P and T data equal to the ratio of the two diameters.
3. If a simplified group analysis is to be carried out utilizing Y and Z modifiers, the data can be factored
appropriately.
Parameters
PFAC
Keyword to signify the factor relates to P data.
YFAC
Keyword to signify the factor relates to Y data.
TFAC
Keyword to signify the factor relates to T data.
ZFAC
Keyword to signify the factor relates to Z data.
factor
Factor by which the data is to be multiplied (Real).
Note that both T-Z and ENDB data will be modified by the TFAC and ZFAC statements.
328
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Parameters
isoil
Identifying number of the soil profile (Integer). Only required on the first P-Y line of a given soil property.
P-Y
Keyword for P-Y data.
toplay
Depth to top of layer (Real).
botlay
Depth to bottom of layer (Real). If omitted, curve is defined at one level only. Depths specified are relative
to global mudline datum as defined using the MUDD data statement. A positive depth corresponds to
a point below the mudline datum.
P
Keyword defining P data ordinates from the P-Y curve.
funitl
Force per unit depth. Value taken from P-Y curve (Real). Only 80 characters are permitted per line. If
more data is required it can be entered on subsequent lines preceded with a colon.
Y
Keyword defining Y data ordinates from the P-Y curve.
disp
Displacement value taken from P-Y curve (Real). The order of the data must correspond with that used
for the P ordinates. Only 80 characters are permitted per line. If more data is required it can be entered
on subsequent lines preceded with a colon.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
329
Soil-Pile-Structure Analysis
Parameters
isoil
Identifying number of the soil profile (Integer). Only required on the first P-Y line of a given soil profile.
P-Y
Keyword for P-Y data.
toplay
Depth to top of layer (Real).
botlay
Depth to bottom of layer (Real). If omitted curve is defined at top layer level only. See Soil Representation (p. 317). Depths specified are relative to global mudline datum as defined using the MUDD data
statement. A positive depth corresponds to a point below the mudline datum.
REES, , MATL, , MMAT, , PY15
One of four keywords specifying the procedure to be used to generate the P-Y curve. The possible values
are shown below. See Soil Representation (p. 317) for more information.
REES use the Reese method for Sand soil types.
MATL the Matlock method can be used for Clay soil types.
MMAT the modified Matlock method can be used for Clay soil types.
PY15 the API method can be used for Clay soil types.
SAND, , CLAY
Keywords specifying the type of soil under consideration.
Soils may be of two general types: cohesionless (SAND) or cohesive (CLAY).
The following procedures and soil types are related as follows:
Soil type
Solution procedure
SAND
REES
CLAY
variables
Soil properties required to generate the P-Y curves. (Real)
For soil type SAND these are:
unit weight of soil
angle of internal friction
initial soil modulus
For soil type CLAY these are:
unit weight of soil
soil shear strength
330
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Parameters
isoil
Identifying number of the soil profile (Integer). Only required on the first T-Z line for a given soil profile.
T-Z
Keyword for T-Z data.
toplay
Depth to top of layer. (Real)
botlay
Depth to bottom of layer (Real). If omitted curve is defined at top layer level only. See Soil Representation (p. 317). Depths specified are relative to a global mudline datum as defined using the MUDD data
statement. A positive depth corresponds to a point below the mudline datum.
T
Keyword defining T data ordinates from the T-Z curve.
funitl
Force per unit depth. Value taken from T-Z curve (Real). Only 80 characters are permitted per line. If
more data is required it can be entered on subsequent lines preceded with a colon.
Z
Keyword defining Z data ordinates from the T-Z curve.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
331
Soil-Pile-Structure Analysis
disp
Displacement value taken from T-Z curve (Real). The order of the data must correspond with that used
for the T ordinates. Only 80 characters are permitted per line. If more data is required it can be entered
on subsequent lines preceded with a colon.
Parameters
isoil
Identifying number of the soil profile (Integer). Only required on the first P-Y data line of a given soil
profile.
T-Z
Keyword for T-Z data.
toplay
Depth to top of layer. (Real)
botlay
Depth to bottom of layer (Real). If omitted curve is defined at top layer level only. See Soil Representation (p. 317). Depths specified are relative to global mudline datum as defined using the MUDD data
statement. A positive depth corresponds to a point below the mudline datum.
V15O, , V15C
Keyword specifying the procedure to be used to generate the T-Z curve. See Soil Representation (p. 317).
The pile may be considered to be either plugged (V15C) or unplugged (V15O).
SAND, , CLAY
Keywords specifying the type of soil under consideration.
Soils may be of two general types: cohesionless (SAND) or cohesive (CLAY).
Either soil type may be used with each solution procedure.
variables
Soil properties required to generate the T-Z curves. (Real)
For soil type SAND these are:
unit weight of soil
angle of internal friction
For soil type CLAY these are:
332
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Parameters
isoil
Identifying number of the soil profile. (Integer)
ENDB
Keyword for end bearing data.
T
Keyword defining T data ordinates from ENDB curve.
force
Force value taken from ENDB curve (Real). Only 80 characters are permitted per line. If more data is required it can be entered on subsequent lines preceded with a colon.
Z
Keyword defining Z data ordinates from ENDB curve.
disp
Displacement value taken from ENDB curve (Real). The order of the data must correspond with that used
for the T ordinates. Only 80 characters are permitted per line. If more data is required it can be entered
on subsequent lines preceded with a colon.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
333
Soil-Pile-Structure Analysis
Parameters
isoil
Identifying number of the soil profile. (Integer)
ENDB
Keyword for ENDB data.
V15O, , V15C
Keyword specifying the solution procedure to be used to generate the ENDB curve. See Soil Representation (p. 317). The pile may be considered to be either plugged (V15C) or unplugged (V15O).
SAND, , CLAY
Keywords specifying the type of soil under consideration. Soils may be of two general types: cohesionless
(SAND) or cohesive (CLAY). Either soil type may be used with each solution procedure.
variables
Soil properties required to generate the ENDB curves (Real).
For soil type SAND these are:
unit weight of soil
angle of internal friction
For soil type CLAY these are:
unit weight of soil
soil shear strength
It is important that the properties are given in the order listed above.
SYMM, , ASYM
Keywords specifying the shape of the ENDB curves. These may be either symmetric (SYMM) or axisymmetric (ASYM). If omitted, a symmetric curve is assumed.
factor
Multiplication factor required for the tensile (negative) portion of the ENDB curve (Real). If omitted, a
value of 0.0 is assumed.
334
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Parameters
isoil
Identifying number of soil profile as used on P-Y and T-Z data. Only required on the first SLID command
line of a given soil profile.
SLID
Compulsory keyword to signify that soil displacement data is being defined.
depth
Depth from mudline to point where soil displacement data is being defined, positive downwards.
xdisp
Soil displacement in the global X axis.
ydisp
Soil displacement in the global Y axis.
zdisp
Soil displacement in the global Z axis.
Note
There are no limits to the number of soil movement definitions for a given soil.
The depths specified do not have to coincide with other soil data, such as P-Y and T-Z information. The program will interpolate values for unspecified depths. For pile depths outside the
values specified on the SLID definition, a constant value is adopted equal to the value defined
at an adjacent depth.
The vertical soil displacement (zdisp) may be omitted if zero. Similarly if both ydisp and zdisp
are zero then they may be omitted.
335
Soil-Pile-Structure Analysis
336
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
0.01 STAT
0.01 STAT
0.01 STAT
0.01 STAT
0.01 STAT
SYMM
SYMM
SYMM
SYMM
SYMM
SYMM
STAT
STAT
STAT
STAT
STAT
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
337
Soil-Pile-Structure Analysis
PILECALC
PILELOAD
PILEMASS
PILERUN
PILESEL
PILESTIF
338
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
!
!
!
!
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
339
Soil-Pile-Structure Analysis
PILEGEN,,ar10+1,ar10+1,1
FINISH
! Generate piles
/SOLU
PILECALC
Using the appropriate geometry, create a named selection that contains a single vertex locating the
pile cap; name the selection PileX, where X is the pile number (starting at 1).
Apply Forces and Standard Earth Gravity as required.
In the commands object, set input arguments as indicated below: the diameter, thickness, and length.
Length is the distance to the pile tip and will be applied to all piles; the length is vertically downwards
(Z direction) from the pile cap.
ARG1 Number of Piles (max 9 without Command Snippet customization)
340
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Soil-Pile-Structure Results
ARG2 Pile Diameter
ARG3 Pile Thickness
ARG4 Pile Length
In the analysis settings change the Solver Controls, Weak Springs to Off.
Equation Results
Available equation results types are:
DISPLACEMENT
REACTION
Valid associated result components are numbers 1 to 6, being X, Y, Z, RX, RY, and RZ respectively.
Element results
The available element results are listed below. Each surface is the subsection, and when the results
originate from a Splinter analysis initialized using the PILEGEN command, the element number is always
1 and the Structure used in forming the model name is the pile name entered via the StruPile argument. Components marked * are stored as character strings:
Type:
PILE DISPLACEMENT
PILE REACTION
RESULTS COMPONENT
DESCRIPTION
SEC.POSN
X value
Y value
Z value
RX
RX value
RY
RY value
RZ
RZ value
Type:
PILE FORCE/MOMENT
RESULTS COMPONENT
DESCRIPTION
SEC.POSN
FXX
X force (axial)
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
341
Soil-Pile-Structure Analysis
RESULTS COMPONENT
DESCRIPTION
FQY
Y shear force
FQZ
Z shear force
TXX
Torsional moment
MYY
Y bending moment
MZZ
Z bending moment
Type:
PILE STRESS
RESULTS COMPONENT
DESCRIPTION
SEC.POSN
SAX
Axial stress
SB
Bending stress
SXX.MAX
Type:
RESULTS COMPONENT
DESCRIPTION
SEC.POSN
UC.AX+BN
15.5. References
The following references are cited in this chapter:
1. Lymon C. Reese, William R. Cox and Francis D. Koop, Analysis of Laterally Loaded Piles in Sand, Paper
No 2080, Offshore Technology Conference, 1974.
2.
Hudson Mattock, Correlations for Design of Laterally Loaded Piles in Soft Clay, Paper No 1024 Offshore
Technology Conference, 1970.
3. American Petroleum Institute, Recommended Practice for Planning, Designing and Constructing Fixed
Offshore Platforms, 15th Edition ,1984.
4. V. N. Vijayvergiya, Load-Movement Characteristics of Piles, PORTS 77.
342
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
343
Note
The following points should be noted for the sequential solution approach:
The generalized matrices (stiffness, etc.) are computed based on the initial undeformed geometry
and assuming small displacement, linear behavior. It is thus implicitly assumed that the
foundation is linear with small deformations throughout the entire solution.
Likewise, when computing the generalized foundation external load time history in the substructure generation pass, it is assumed that the structural displacement and velocity are zero
since such information is not available when the loading is generated.
The hydrodynamic mass for the supporting structure is formed based on the water elevation
at the first time at which the solution is attempted.
If the supporting structure internal forces are recovered statically in step 4, the dynamic forces
(e.g. inertial force) in the foundation will be ignored. The dynamic effects can be accounted
for by running this step as a transient job. However, it should be noted that the following
points may affect the accuracy of the solution:
The generalized mass used in the aeroelastic analysis is only an approximation to the true
mass matrix (static reduction is exact but not dynamic).
Different time integration schemes may be adopted by aeroelastic code and Mechanical
APDL. Hence, while the displacement time histories are identical in both runs (for the prescribed displacement case), it may not be the case for the velocity and acceleration time
histories.
There should be little difference between applying forces or displacements to recover the
foundation forces. The two methods should yield identical results in a linear static analysis.
This approach should be much more efficient than the fully coupled approach as there is no
need to keep both the aeroelastic code and Mechanical APDL running simultaneously and keep
exchanging information every time step.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where ngenfr is the number of generalized freedoms, which is always 6 at present, and val is the
generalized matrix.
The row and column order in the generalized matrix corresponds to the order UX, UY, UZ, ROTX, ROTY,
ROTZ.
A time series of the generalized foundation external loading vector is written to another file. At each
solution time, the time (t) and the associated generalized load vector (f ) will be output to this file. The
force file has the following format:
ngenfr (1x, I6)
t1 (1x, E12.5)
(f1(i), i = 1,ngenfr) (6(1x, E12.5))
t2 (1x, E12.5)
(f2(i), i = 1,ngenfr) (6(1x, E12.5))
.
.
tn (1x, E12.5)
(fn(i), i = 1,ngenfr) (6(1x, E12.5))
The load values are ordered in the same way as the generalized matrices.
345
At this point the aeroelastic analysis can be run, using the output from the above analysis. Once complete,
a second Mechanical APDL analysis is run with a time series of forces and/or displacements at the interface node. These need to be converted from the aeroelastic output to Mechanical APDL compatible
output by the user (for example, using Excel) or the aeroelastic analysis program. Any ocean loading
346
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
347
348
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
349
where
*
text
where
GENE
350
water depth
density
water density
gravity
acceleration of gravity
symxy
symmetry integer
0: no symmetry
1: symmetry about x axis
2: symmetry about y axis
3:symmetry about both x and y axes
1.025000E+03
9.800000E+00
where
HYDR
npan
xk(j)
yk(j)
zk(j)
The nodes must be arranged in counterclockwise order when viewing the panel from the wetted side.
To define a triangular panel, specify the same coordinates for nodes 3 and 4.
Example 17.2: Hydrodynamic Surface Geometry
HYDR
12
15.0
15.0
15.0
15.0
7.5
7.5
7.5
7.5
15.0
15.0
15.0
15.0
-15.0
-15.0
-15.0
-7.5
-15.0
-15.0
-15.0
-7.5
-15.0
-15.0
-7.5
-7.5
0.0
-7.5
-15.0
-15.0
0.0
-7.5
-15.0
-15.0
0.0
-7.5
0.0
-7.5
7.5
7.5
7.5
7.5
0.0
0.0
0.0
0.0
15.0
15.0
15.0
15.0
-15.0
-15.0
-15.0
-7.5
-15.0
-15.0
-15.0
-7.5
-15.0
-15.0
-7.5
-7.5
0.0
-7.5
-15.0
-15.0
0.0
-7.5
-15.0
-15.0
-7.5
-15.0
-7.5
-15.0
7.5
7.5
7.5
7.5
0.0
0.0
0.0
0.0
15.0
15.0
15.0
15.0
-15.0
-15.0
-7.5
0.0
-15.0
-15.0
-7.5
0.0
-7.5
-7.5
0.0
0.0
-7.5
-15.0
-15.0
-15.0
-7.5
-15.0
-15.0
-15.0
-7.5
-15.0
-7.5
-15.0
15.0
15.0
15.0
15.0
7.5
7.5
7.5
7.5
15.0
15.0
15.0
15.0
-15.0
-15.0
-7.5
0.0
-15.0
-15.0
-7.5
0.0
-7.5
-7.5
0.0
0.0
-7.5
-15.0
-15.0
-15.0
-7.5
-15.0
-15.0
-15.0
0.0
-7.5
0.0
-7.5
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
351
where
PERD
nperd
period(i)
The data may be specified in one or more lines until all the periods are entered.
Example 17.3: Wave Periods
perd
1
18.0
where
DIRN
ndirn
heading(i)
the ith wave direction (in degrees, between -180 and +180)
The data may be specified in one or more lines until all the directions are entered.
Example 17.4: Wave Directions
dirn
2
0.0 90.0
where
PRES
iperiod
352
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
region
panel
magnitude
pressure amplitude
phase
associated phase angle in degrees (with respect to the wave at the center of gravity of the
vessel)
The associated phase angle (phase) determines the horizontal location of the loading. The magnitude
and phase together define the shape of the harmonic loading curve.
Pressure is positive when acting from the wetted surface towards the inner surface.
If no input pressure is given to a panel at a particular period and direction, the pressure on the panel
is assumed to be zero. If more than one set of pressures are defined, their effects are accumulative; that
is, the real and imaginary parts of each set are summed together.
Example 17.5: Panel Pressures
PRES
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
2
3
4
5
6
25.0
75.0
100.0
100.0
25.0
75.0
0.0
0.0
0.0
0.0
0.0
0.0
where
MORI
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
353
xk(j)
yk(j)
zk(j)
where
WKIN
iperiod
iheading
node
element
pmag
pphs
zsmag
zsphs
vxmag
x velocity amplitude
vxphs
vymag
y velocity amplitude
vyphs
vzmag
z velocity amplitude
vzphs
axmag
x acceleration amplitude
axphs
aymag
y acceleration amplitude
ayphs
azmag
z acceleration amplitude
azphs
Note that the data is defined for X, Y and Z in the global directions.
354
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
where
RAOS
xr yr zr
iperiod
iheading
xmag
x RAO amplitude
xphs
ymag
y RAO amplitude
yphs
zmag
z RAO amplitude
zphs
rxmag
rxphs
rymag
ryphs
rzmag
rzphs
1.042555E+00
2.027182E-07
1.042556E+00
1.605987E-07
4.346394E-02
0.000000E+00
4.340762E-02
0.000000E+00
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
355
where
MASS
mass(j,k)
0.000000E+00 -5.000000E+00
0.000000E+00 0.000000E+00
1.383750E+07 0.000000E+00
0.000000E+00 1.383750E+07
0.000000E+00 0.000000E+00
0.000000E+00 0.000000E+00
0.000000E+00 0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
7.780000E+08
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
7.780000E+08
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
7.780000E+08
356
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
An .ahd file named analysis_1, which contains the data from the hydrodynamic diffraction analysis,
was generated via the AQWA product's AQWA2NEUT utility, as follows:
aqwa2neut analysis 1
The diffracted wave loading data in the analysis_1.ahd file is imported (OCREAD) and applied to
the structural model for further analysis, as shown in the following example input:
/verify,difwaveeg
/FILNAME,difwaveg
/prep7
/TITLE,difwaveeg, EXAMPLE FOR AQWA LOADING
/com **************************************************************************
/com SHELL MODEL OF A HOLLOW BOX TUBE
/com CREATED 20/05/10
/com **************************************************************************
antype,static
nlgeom,off
! define shell section
sectype,11,shell
secdata,0.1
! define ocean
matwat=2
idwat=2
idcur=3
idwav=4
depth=250.0
! define geometry
/prep7
!wpoff,,,-15
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
357
358
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
! 1 SUBSTEPS
The following figure shows the analysis results using the pressures from the diffracted wave:
Figure 17.2: Analysis Results Using Pressure Data from a Hydrodynamic Diffraction Analysis
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
359
360
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
Index
Symbols
.mnf suffix, 293
/CYCEXPAND command
limitations, 172
2-D to 3-D analysis
benefits, 135
initial-state-based analysis, 141
introduction to, 135
key commands, 138
map solution variables, 141
performing, 139
requirements, 135
restrictions, 141
selecting a region for remeshing, 140
starting, 140
substep at which to initiate, 139
supported analyses, 135
understanding the process, 137
A
ADAMS (see ANSYS-ADAMS Interface)
ADAMS.MAC macro, 293
adaptive meshing
activating custom subroutines, 88
ADAPT macro: customizing (via UADAPT.MAC), 89
custom loads (via ADAPTBC.MAC), 88
custom meshing (via ADAPTMSH.MAC), 87
custom solution (via ADAPTSOL.MAC), 88
element types for, 85
example problems: finding, 90
hints, 89
modifying the basic process, 86
overview of, 85
prerequisites for, 85
process involved, 85
selective adaptivity, 86
user subroutines: employing, 87
aeroelastic analysis, 343
ANSYS-ADAMS Interface
ADAMS simulation: running, 289
ADAMS.MAC macro, 293
batch mode
exporting to ADAMS, 288
importing loads, 292
example analysis, 294
exporting to ADAMS, 286
interface points: modeling, 285
introduction, 283
loads
B
birth and death of elements
applying loads and solving, 301
building the model, 301
constraint equations, 300
controlled by analysis results, 303
example problems, 303
first load step: defining, 301
in linear models, 300
initial strain, 300-301
material properties: changing, 301
multiple load steps, 300
Newton-Raphson option, 300
nodal constraints and loads, 301
overview of, 299
reviewing results, 302
sample input, 302-303
status of an element: determining, 300
strain, 299
subsequent load steps: defining, 302
supported elements, 299
usage guidelines, 300
using in an analysis, 301
C
Campbell diagram, 199
CDB generic mesh file, 106
CMS (see component mode synthesis)
component mode synthesis
database: saving copy during generation pass, 261
definition of, 259
example analysis, 264
example: expanding a superelement in a transformed
location, 279
expanding superelements in transformed locations,
263
expansion pass, 263
generation pass, 261
introduction to, 259
lumped mass approximation, 261
master DOFs, 261
methods supported, 259
mode shapes: plotting or printing, 264
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
361
Index
offsets: applying to node and element IDs, 263
process (passes) involved, 261
solvers and files used, 260
specifying CMS results files for postprocessing, 264
superelement: creating, 261
understanding, 259
use pass, 263
constraint equations
in a cyclic symmetry analysis, 151
Coriolis effect, 197
cut boundary displacement method (see submodeling)
cyclic symmetry analysis
ANCYC command, 169
automated, 144
availability, 143
basic sector, 145
buckling, 162
commands affected by non-cyclic loading, 155
constraint equations for non-cyclically symmetric
loading, 152
coupling and constraint equations (CEs), 151
CYCLIC command, 146
database: saving after solution, 165
definition of cyclic symmetry, 143
deleting a sector load, 152
Duplicate Sector solution method, 151
edge component pairs, 146
edge node patterns: identical vs. dissimilar, 149
example: buckling, 177
example: modal, 172
EXPAND command, 169
expanding the solution, 167
geometric nonlinearity, 160
graphical picking of non-cyclic boundary conditions,
155
harmonic index, 157
inconsistent constraints, 152
introduction to, 143
large-deflection prestressed modal, 160
limitations, 171
linear buckling, 162
meshing the cyclic sector, 149
modal, 157
model only one sector, 144
model verification during preprocessing, 150
model verification during solution, 165
modeling, 144
most common commands, 144
naming convention for low- and high-edge component pairs, 146
nodal diameter, 157
non-cyclic loading, 155
362
D
damping
substructuring, 240
diffracted wave, 349
E
elements
birth and death of, 299
expansion pass
substructuring, 249
F
files
.mnf, 293
Jobname.LOOP, 3
Jobname.NODE, 223
Jobname.PDS, 3
probabilistic design analysis file, 9
sename.sub, 240
fixed-interface CMS method, 259
flexible body dynamics (see rotating structure analysis)
free-interface CMS method, 259
H
hydrodynamic analysis, 349
J
Jobname
in a submodeling analysis, 221
in a substructure analysis, 239
L
loads
substructuring, 240
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
M
macros
ADAMS.MAC, 293
material properties
changing to deactivate or reactivate elements, 301
meshing
adaptive, 85
modal neutral file, 293
model generation
adaptive meshing, 85
N
natural frequencies, 204
Newton-Raphson
birth and death analysis, 300
nonstandard uses of ANSYS, 307
O
ocean loading, 349
offsets
to node and element IDs in a CMS analysis, 263
P
PDS (see probabilistic design)
pile analysis, 315
probabilistic design
algorithms available, 24, 51
analysis file
creating, 9
specifying, 13
analysis flow, 3
analysis results: reviewing, 41
and traditional analysis methods: comparison of, 2
correlations between random variables, 20
cumulative distribution function: plotting, 42
definition of, 1
employing, 8
example analysis, 69
example problem, 10
histogram: plotting, 42
HTML report: generating, 67
inverse probabilities: requesting, 42
loops: executing, 25
Monte Carlo Simulation method, 51
parallel execution, 27
parameters: establishing, 12
PDS database
clearing, 69
description of, 3
saving, 68
postprocessing, 59
print correlation matrix, 43
R
residual-flexible free-interface CMS method, 259
rezoning
balancing residual forces, 118
benefits, 93
contract boundaries, 116
convergence difficulties during mapping operation,
120
example, 124, 127
file structures for repeated rezonings, 121
forces and isolated applied displacements, 116
horizontal rezoning, 106
introduction to, 93
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.
363
Index
key commands, 101
limitations, 95
mapping variables, 118
mesh control during, 114
multiframe restart after, 120
multiple regions remeshed at the same substep,106
nodal temperatures, 116
postprocessing results, 121
pressure and contiguous displacements, 116
remeshing, 104
repeating as needed, 121
requirements, 96
restrictions, 123
selecting a region for remeshing, 103
starting, 103
substep at which to initiate, 102
supported analyses, 96
understanding the process, 98
using a generic (third-party) mesh, 106
using manual mesh splitting, 110
variable mapping, 118
rigid-body dynamics, 283
rotating structure analysis
Campbell diagram, 199
Campbell diagram analysis example, 205
choosing the appropriate reference frame option,
204
Coriolis analysis example, 207
flexible body dynamics, 202
harmonic analysis for unbalance forces example,210
introduction to, 197
natural frequencies, 204
rotating reference frame, 202
rotordynamics, 198
stationary reference frame, 198
understanding, 197
rotordynamics (see rotating structure analysis)
S
soil data description, 325
soil-pile-structure analysis, 315
soil data description, 325
solid modeling
adaptive meshing, 85
mesh density, 85
State-Space Matrices
Export of, 309
submodeling, 219
analyzing the submodel, 224
coarse model: creating and analyzing, 221
creating the submodel, 222
cut-boundary interpolation: performing, 223
364
T
traveling wave animation, 169
U
UPFs (see user-programmable features)
user-programmable features
available types, 306
overview of, 305
process involved, 305
V
Variational Technology
example problem
transmission line, 78
waveguide, 79
forced-frequency simulations: solving, 76
overview, 75
W
wind turbine analysis, 343
sequentially coupled, 343
Release 14.5 - SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates.