PERTURB, Type
, MatKey
, ContKey
, LoadControl
Sets linear perturbation analysis options.
Type
Type of linear perturbation analysis to be performed:
STATIC | — | Perform a linear perturbation static analysis. |
MODAL | — | Perform a linear perturbation modal analysis. |
BUCKLE | — | Perform a linear perturbation eigenvalue buckling analysis. |
HARMONIC | — | Perform a linear perturbation full harmonic analysis. |
SUBSTR | — | Perform a linear perturbation substructure generation pass. |
OFF | — | Do not perform a linear perturbation analysis (default). |
MatKey
Key for specifying how the linear perturbation analysis uses material properties, valid for all structural elements except contact elements. For more information, see Linear Perturbation Analysis in the Mechanical APDL Theory Reference.
AUTO | — | The program selects the material properties for the linear perturbation analysis automatically (default). The materials are handled in the following way:
|
TANGENT | — | Use the tangent (material Jacobian) on the material constitutive curve as the material property. The material property remains linear in the linear perturbation analysis and is obtained at the point of the base analysis where restart occurs. The materials are handled in the following way:
|
SPOFF | — | Provide the same values as AUTO, but set the spin softening matrix to zero (ignoring the spin softening effect). The spin softening effect is excluded in all the linear perturbation analysis types except for linear perturbation buckling. Note that although the spin softening effect is excluded in linear perturbation analysis, it is still included in the base static or full transient analysis if NLGEOM,ON is issued in the base analysis. |
ContKey
Key that controls contact status for the linear perturbation
analysis. This key controls all contact elements (TARGE169, TARGE170, and CONTA171 through CONTA178) globally for all contact
pairs. Alternatively, contact status can be controlled locally per
contact pair by using the CNKMOD command. Note
that the contact status from the base analysis solution is always
adjusted by the local contact controls specified by CNKMOD first and then modified by the global sticking or bonded control
(ContKey
= STICKING or BONDED). The tables
in the Notes section show how
the contact status is adjusted by CNKMOD and/or
the ContKey
setting.
CURRENT | — | Use the current contact status from the restart snapshot (default). If the previous run is nonlinear, then the nonlinear contact status at the point of restart is frozen and used throughout the linear perturbation analysis. |
STICKING | — | For frictional contact pairs (MU > 0), use sticking contact (e.g., MU*KN for tangential contact stiffness) everywhere the contact state is closed (i.e., status is sticking or sliding). This option only applies to contact pairs that are in contact and have a frictional coefficient MU greater than zero. Contact pairs without friction (MU = 0) and in a sliding state remain free to slide in the linear perturbation analysis. |
BONDED | — | Any contact pairs that are in the closed (sticking or sliding) state are moved to bonded (for example, KN for both normal and tangential contact stiffness). Contact pairs that have a status of far-field or near-field remain open. |
LoadControl
Key that controls how the load vector of {Fperturbed} is calculated. This control is provided for convenience of load generation for linear perturbation analysis. In general, a new set of loads is required for a linear perturbation analysis. This key controls all mechanical loads; it does not affect non-mechanical loads. Non-mechanical loads (including thermal loads) are always kept (i.e., not deleted).
ALLKEEP | — | Keep all the boundary conditions (loads and constraints) from the end of the load step of the current restart point. This option is convenient for further load application and is useful for a linear perturbation analysis restarted from a previous linear analysis. For this option, {Fend} is the total load vector at the end of the load step at the restart point. |
INERKEEP | — | Delete all loads and constraints from the restart step, except for displacement constraints and inertia loads (default). All displacement constraints and inertia loads are kept for convenience when performing the linear perturbation analysis. Note that nonzero and tabular displacement constraints can be considered as external loads; however, they are not deleted when using this option. |
PARKEEP | — | Delete all loads and constraints from the restart step, except for displacement constraints. All displacement constraints are kept for convenience when performing the linear perturbation analysis. Note that nonzero and tabular displacement constraints can be considered as external loads; however, they are not deleted when using this option. |
DZEROKEEP | — | Behaves the same as the PARKEEP option, except that all nonzero displacement constraints are set to zero upon the onset of linear perturbation. |
NOKEEP | — | Delete all the loads and constraints, including all displacement constraints. For this option, {Fend} is zero unless non-mechanical loads (e.g., thermal loads) are present. |
Linear perturbation analysis is disabled (Type
= OFF) by default. When the linear perturbation analysis is enabled,
linear material property behavior is assumed for stress calculations;
contact status for all contact pairs from the point of restart is
used by default; and all loads and constraints from the restart step
are deleted, except for displacement constraints and inertia loads,
by default.
This command controls options relating to linear perturbation analyses. It must be issued in the first phase of a linear perturbation analysis.
This command is also valid in PREP7.
A linear perturbation analysis consists of two phases (two SOLVE commands). The first phase is a restart from a base analysis. This base analysis must be a linear or nonlinear static analysis or full transient analysis. The first phase starts with the ANTYPE,,RESTART,,,PERTURB command and ends with the SOLVE,ELFORM command. The purpose of the first phase is to re-establish a snapshot of the stiffness matrices at the specified restart point. The second phase, ending with the second SOLVE command, is for the actual linear perturbation analysis.
The total perturbed loads are calculated as follows:
{Fperturbed} = {Fend} + {Fadd}
where:
{Fend} = total loads at the end of the load step of the current restart
point (load applications are read from the .LDHI file). By default, all of the loads in {Fend} are deleted except for displacement boundary conditions and inertia
loads (see the description of LoadControl above). |
{Fadd} = Additional (new) loads prescribed by the user in the second phase of the linear perturbation analysis (after the first SOLVE command is invoked). |
In the first phase of a linear perturbation analysis, the ANTYPE,,RESTART command resumes the Jobname.RDB database and reads in the .LDHI file to establish the {Fend} load. New load application (adding to {Fadd}) or load removal (changing {Fend}) can be done only in the second phase of the linear perturbation analysis (after the first SOLVE command), allowing flexibility in controlling the final {Fperturbed} to be used.
For Type
= STATIC, {Fperturbed} is the actual external load for the static analysis.
For Type
= MODAL, {Fperturbed} is calculated and stored in the .FULL and .MODE files for a subsequent mode-superposition, PSD,
or other type of modal-based linear dynamic analysis. Linear dynamic
options such as multiple load generations (MODCONT,ON), enforced motion (MODCONT, ,ON), and residual
vector methods (RESVEC,ON) can be used in a linear
perturbation analysis. For these methods, the MODCONT or RESVEC command must be invoked in the second
phase (after the first SOLVE) of the linear perturbation
procedure. For the enforced motion method, the base identification
number should be specified (D command) in the second
phase of the linear perturbation analysis. This base identification
number is used later in the downstream mode-superposition or other
mode-superposition based analysis.
For Type
= BUCKLE, {Fperturbed} is the actual linear buckling load which is used to generate the
linear stress stiffening matrix for the buckling analysis.
For Type
= HARMONIC, {Fperturbed} is the actual external load for the full harmonic
analysis. In this case, {Fperturbed} can be
frequency dependent and can use complex input.
For Type
= SUBSTR,
{Fperturbed} is used to generate the first
reduced external load vector of the substructure.
In most cases involving linear perturbation analysis, Matkey
= AUTO is the best option for controlling material
behavior. Matkey
= TANGENT is often the
better option, however, in special cases such as the following:
A linear perturbation buckling analysis, to introduce preferred buckling modes into a subsequent post-buckling nonlinear analysis.
A linear perturbation modal analysis, to introduce preferred modes into a subsequent bifurcation analysis.
If the TANGENT option is used in conjunction with a modal restart solution that follows a linear perturbation modal analysis, then the AUTO option is assumed and linear material properties are used for stress calculations in the modal restart solution. This occurs because the TANGENT material properties are not available during the modal restart phase due to a data architecture limitation. Furthermore, linear material properties are used for the stress calculation in any downstream analysis that uses the modal restart solution.
For more information about the automatic and tangent options, see Specifying Material Behavior in Linear Perturbation in the Element Reference.
You can control the contact status for the linear perturbation
analysis by using the ContKey
field on
this command and/or the CNKMOD command. The first
table shows the effects of using only the ContKey
setting on the PERTURB command. The second table
shows the effects of using both the CNKMOD command
and the ContKey
setting on PERTURB.
Table 207: Adjusted Contact Status when PERTURB Command Is Issued
Contact Status from the Base Analysis Solution at the Restart Point | ContKey Setting on PERTURB Command | |
ContKey Value | Adjusted Contact Status | |
0 - far-field | any | 0 - far-field |
1 - near-field | any | 1 - near-field |
2 - sliding | CURRENT or STICKING (mu=0) | 2 - sliding |
STICKING (mu>0) or BONDED | 3 - sticking | |
3 - sticking | any | 3 - sticking |
Table 208: Adjusted Contact Status when Both CNKMOD and PERTURB Are Issued
Contact Status from the Base Analysis Solution at the Restart Point | CNKMOD,ITYPE ,12,Value | ContKey Setting
on PERTURB Command | ||
KEYOPT(12) Value | Adjusted Contact Status | ContKey Value | Final Adjusted Contact Status | |
0 - far-field | any | 0 - far-field | any | 0 - far-field |
1 - near-field | 0, 1, 2, 3, 6 | 1 - near-field | any | 1 - near-field |
4 | 1 - near-field (if outside of the adjusted pinball region) | any | 1 - near-field | |
2 - sliding (if inside of the adjusted pinball region) | CURRENT or STICKING (mu=0) | 2 - sliding | ||
STICKING (mu>0) or BONDED | 3 - sticking | |||
5 | 1 - near-field (if outside of the adjusted pinball region) | any | 1 - near-field | |
3 - sticking (if inside of the adjusted pinball region) | any | 3 - sticking | ||
2 - sliding | 0, 2, 4 | 2 - sliding | CURRENT or STICKING (mu=0) | 2 - sliding |
STICKING (mu>0) or BONDED | 3 - sticking | |||
1, 3, 5, 6 | 3 - sticking | any | 3 - sticking | |
3 - sticking | any | 3 - sticking | any | 3 - sticking |
When ContKey
is set to CURRENT, all
contact related parameters (such as normal stiffness and tangential
stiffness) will remain unchanged throughout the linear perturbation
analysis. However when ContKey
is set to
STICKING or BONDED, the program will re-evaluate the contact normal
and tangential stiffness in order to perform the linear perturbation
analysis based on the actual sticking behavior regardless of the friction
coefficient value.
Note that the CNKMOD command allows you to take points in the base analysis that are near contact (within the pinball region) and modify them to be treated as "in contact" in the linear perturbation analysis; see the "1 - near-field" row in the above table with KEYOPT(12) values set to 4 or 5. CNKMOD also allows you to take points that are sliding in the base analysis and treat them as sticking in the linear perturbation analysis, irrespective of the MU value; see the "2 - sliding" row in the above table with KEYOPT(12) values set to 1,5, or 6.
If an open gap exists at the restart point of the base static/transient solution and the contact status is adjusted as sliding or sticking due to a “bonded” or “no separation” contact behavior definition, then the program will treat it as near-field contact when executing the CNKMOD command in a downstream linear perturbation analysis.