Rabu, 31 Oktober 2012

Running a Nonlinear Analysis in ANSYS

http://www.2doworld.com/ansys-and-etabs/running-a-nonlinear-analysis-ansys.html

Running a Nonlinear Analysis in ANSYS


PrinterPrinterEmailEmail
ANSYS employs an automatic solution control method that, based on the physics of your problem, sets various nonlinear analysis controls to the appropriate values. If you are not satisfied with the results obtained with these values, you can manually override the settings. The following commands are set to optimal defaults:
AUTOTS

PRED

MONITOR

DELTIM

NROPT

NEQIT

NSUBST

TINTP

SSTIF

CNVTOL

CUTCONTROL

KBC

LNSRCH

OPNCONTROL

EQSLV

ARCLEN

CDSWRITE

LSWRITE

These commands and the settings they control are discussed in later sections. You can also refer to the individual command descriptions in the ANSYS Commands Reference.
If you do choose to override the ANSYS-specified settings, or if you wish to use an input list from a previous release of ANSYS, issue SOLCONTROL,OFF in the/SOLUTION phase, or /CONFIG,NLCONTROL,OFF following the /BATCHcommand. See the SOLCONTROL command description for more details.
ANSYS' automatic solution control is active for the following analyses:

  • Single-field nonlinear or transient structural and solid mechanics analysis where the solution DOFs are combinations of UX, UY, UZ, ROTX, ROTY, ROTZ.
  • Single-field nonlinear or transient thermal analysis where solution DOF is TEMP.

1 Build the Model

This step is essentially the same for both linear and nonlinear analyses, although a nonlinear analysis might include special elements or nonlinear material properties. See Section 8.11, "Using Nonlinear (Changing-Status) Elements" and Section 8.3,"Modeling Material Nonlinearities" for more details. If your analysis includes large-strain effects, your stress-strain data must be expressed in terms of truestress and true (or logarithmic) strain. For more information on building models in ANSYS, see the ANSYS Modeling and Meshing Guide
.

2 Apply Loads and Obtain the Solution

After you've created a model in ANSYS, you define the analysis type and options, apply loads, specify load step options, and solve. A nonlinear solution will differ from a linear solution in that it often requires multiple load increments, and always requires equilibrium iterations. The general procedure follows. See Section 8.5 for a sample problem that walks you through a specific nonlinear analysis.
1. Specify a new or restart analysis and define the analysis type; in this case, static. The analysis type cannot be changed after the first load step (that is, after you issue your first SOLVE command).
Command(s):

    ANTYPE
GUI:
Main Menu>Solution>-Analysis Type-New Analysis/Restart>Static

    You will usually use New Analysis. Restarts are discussed in the ANSYS Basic Analysis Procedures Guide.
2. Specify if your problem has large deflections or large strain. This option cannot be changed after the first load step (that is, after you issue your first SOLVEcommand).
Command(s):

    NLGEOM
GUI:
Main Menu>Solution>Analysis Options

    Other analysis options are discussed in Section 8.4.2.1, "Advanced Analysis Options."
    Not all nonlinear analyses will produce large deformations. See Section 8.2 "Using Geometric Nonlinearities," for further discussion of large deformations.
3. Apply loads on the model. See Chapter 2, Table 2-3 of this guide and Chapter 2 of the ANSYS Basic Analysis Procedures Guide for load information. Remember that inertia and point loads will maintain constant direction, but surface loads will "follow" the structure in a large-deformation analysis. You can apply complex boundary conditions by defining a 1-dimensional table (TABLE type array parameter). See Section 2.5.2.3 "Applying Loads Using TABLE Type Array Parameters," for more information.
4. Specify time/frequency load step options: Time, number/size of time steps, and stepped/ramped loading.
Command(s):

    TIME
    NSUBST
    DELTIM
    KBC
GUI:
Main Menu>Solution>-Load Step Opts-Time/Frequenc>Time & Time Step/Time & Substep

    These options can be changed at any load step. See Chapter 2 of the ANSYS Basic Analysis Procedures Guide for more information on these options. Advanced time/frequency options, in addition to those listed here, are discussed in Section 8.4.2.2, "Advanced Load Step Options." Note that you only have to specify stepped or ramped loads (KBC) in a transient or dynamic analysis.
    A nonlinear analysis requires multiple substeps (or time steps; the two terms are equivalent) within each load step so that ANSYS can apply the specified loads gradually and obtain an accurate solution. The NSUBST and DELTIMcommands both achieve the same effect (establishing a load step's starting, minimum, and maximum step size), but by reciprocal means. NSUBST defines the number of substeps to be taken within a load step, whereas DELTIM defines the time step size explicitly. If automatic time stepping is off, then the starting substep size is used throughout the load step.
5. Specify nonlinear load step options as necessary: Solution termination criteria and arc-length option.
Command(s):

    ARCTRM
    ARCLEN
GUI:
Main Menu>Solution>-Load Step Opts-Nonlinear>Arc-Length Opts

    These options can be changed at any load step. See Chapter 2 of the ANSYS Basic Analysis Procedures Guide for more information on these options. Advanced nonlinear options, in addition to those listed here, are discussed in Section 8.4.2.2, "Advanced Load Step Options."
6. Specify birth and death options as necessary.
Command(s):

    EKILL
    EALIVE
    MPCHG
GUI:
Main Menu>Solution>-Load Step Opts-Other>Kill Elements
Main Menu>Solution>-Load Step Opts-Other>Activate Elements
Main Menu>Solution>-Load Step Opts-Other>Change Mat Props>Change Mat num

    The program "deactivates" an element by multiplying its stiffness by a very small number (which is set by the ESTIF command), and by removing its mass from the overall mass matrix. Element loads (pressure, heat flux, thermal strains, etc.) for inactive elements are also set to zero. You need to define all possible elements during preprocessing; you cannot create new elements in SOLUTION.
    Those elements to be "born" in later stages of your analysis should bedeactivated before the first load step, and then reactivated at the beginning of the appropriate load step. When elements are reactivated, they have a zero strain state, and (if NLGEOM,ON) their geometric configuration (length, area, etc.) is updated to match the current displaced positions of their nodes. See the ANSYS Advanced Analysis Techniques Guide for more information.
7. Specify output control options. See Chapter 2 of the ANSYS Basic Analysis Procedures Guide for more information on these options.
Command(s):

    OUTPR
    OUTRES
    ERESX
GUI:
Main Menu>Solution>-Load Step Opts-Output Ctrls>SoluPrintout
Main Menu>Solution>-Load Step Opts-Output Ctrls>DB/Results File
Main Menu>Solution>-Load Step Opts-Output Ctrls>Integration Pt

    Printed Output [OUTPR] includes any results data on the output file (Jobname.OUT).
    Database and Results File Output [OUTRES] controls the data on the results file (Jobname.RST). By default, only the last substep is written to the results file in a nonlinear analysis. To write all substeps, set the FREQ field on OUTRES to ALL.
    Only 1000 results sets (substeps) can be written to the results file, but you can use the command /CONFIG,NRES to increase the limit (see the ANSYS Basic Analysis Procedures Guide).
    Extrapolation of Results [ERESX] copies an element's integration point stress and elastic strain results to the nodes instead of extrapolating them, if nonlinear strains (plasticity, creep, swelling) are present in the element. The integration point nonlinear strains are always copied to the nodes.
    See the ANSYS Basic Analysis Procedures Guide for more information on how to use these commands properly.
8. Save a back-up copy of the database to a named file.
Command(s):

    SAVE
GUI:
Utility Menu>File>Save As
9. Solve.
Command(s):

    SOLVE
GUI:
Main Menu>Solution>-Solve-Current LS
10. If you need to define multiple load steps, repeat steps 4-8 for each additional load step. Other methods for multiple load steps-the load step file method and the array parameter method-are described in the ANSYS Basic Analysis Procedures Guide.
11. Leave the SOLUTION processor.
Command(s):

    FINISH
GUI:
Close the Solution menu.

2.1 Advanced Analysis Options

Stress Stiffness

To account for buckling, bifurcation behavior, ANSYS includes stress stiffness in all geometrically nonlinear analyses. If you are confident of ignoring such effects, you can turn stress stiffening off (SSTIF,OFF). This command has no effect with several ANSYS elements; see the ANSYS Elements Reference for the description of the specific elements you're using.

Newton-Raphson Option

ANSYS' automatic solution control will use the FULL Newton-Raphson option with adaptive descent off if there is a nonlinearity present. However, when node-to-node, node-to-surface contact elements are used for contact analysis with friction, then adaptive descent is automatically turned on (for example, PIPE20, BEAM23, BEAM24, and PIPE60). The underlying contact elements require adaptive descent for convergence.
Use this option only in a nonlinear analysis. This option specifies how often the tangent matrix is updated during solution. If you choose to override the default, you can specify one of these values:

  • Program-chosen (NROPT,AUTO): The program chooses which of the options to use, based on the kinds of nonlinearities present in your model. Adaptive descent will be automatically activated, when appropriate.
  • Full (NROPT,FULL): The program uses the full Newton-Raphson procedure, in which the stiffness matrix is updated at every equilibrium iteration.
  •  
      If adaptive descent is on (optional), the program will use the tangent stiffness matrix only as long as the iterations remain stable (that is, as long as the residual decreases, and no negative main diagonal pivot occurs). If divergent trends are detected on an iteration, the program discards the divergent iteration and restarts the solution, using a weighted combination of the secant and tangent stiffness matrices. When the iterations return to a convergent pattern, the program will resume using the tangent stiffness matrix. Activating adaptive descent will usually enhance the program's ability to obtain converged solutions for complicated nonlinear problems but is supported only for elements indicated under "Special Features" in Table 4.n.1 of the ANSYS Elements Reference.
  • Modified (NROPT,MODI): The program uses the modified Newton-Raphson technique, in which the tangent stiffness matrix is updated at each substep. The matrix is not changed during equilibrium iterations at a substep. This option is not applicable to large deformation analyses. Adaptive descent is not available.
  • Initial Stiffness (NROPT,INIT): The program uses the initial stiffness matrix in every equilibrium iteration. This option can be less likely to diverge than the full option, but it often requires more iterations to achieve convergence. It is not applicable to large deformation analyses. Adaptive descent is not available.
  • If a multistatus element is in the model, however, it would be updated at the iteration in which it changes status, irrespective of the Newton-Raphson option.

Equation Solver

ANSYS' automatic solution control activates the sparse direct solver (EQSLV, SPARSE) for most cases. Other choices include the frontal, sparse, or PCG solvers. For applications using solid elements (for example, SOLID92 or SOLID45), the PCG solver may be faster, especially when doing 3D modeling.
The sparse direct solver is in sharp contrast to the iterative solvers included in ANSYS, such as the PCG solvers. Although the PCG solver can solve indefinite matrix equations, when the PCG solver encounters an ill-conditioned matrix, the solver will iterate to the specified number of iterations and stop if it fails to converge. When this happens, it triggers bisection. After completing the bisection, the solver continues the solution if the resulting matrix is well-conditioned. Eventually, the entire nonlinear load step can be solved.
Use the following guidelines for selecting either the sparse or PCG solvers for nonlinear structure analysis:

  • If it is a beam/shell or beam/shell and solid structure, choose the sparse solver.
  • If it is a 3D solid structure and the number of DOF is relatively large, choose the PCG solver.
  • If the problem is ill-conditioned (triggered by poor element shapes), or has a big difference in material properties in different regions of the model, or has insufficient displacement boundary constraints, choose the sparse solver.

2.2 Advanced Load Step Options

Automatic Time Stepping

ANSYS' automatic solution control turns automatic time stepping on. An internal auto-time step scheme ensures that the time step variation is neither too aggressive (resulting in many bisection/cutbacks) nor too conservative (time step size is too small). At the end of a time step, the size of the next time step is predicted based on four factors:

  • Number of equilibrium iterations used in the last time step (more iterations cause the time step size to be reduced)
  • Predictions for nonlinear element status change (time step sizes are decreased when a status change is imminent)
  • Size of the plastic strain increment
  • Size of the creep strain increment

Convergence Criteria

The program will continue to do equilibrium iterations until the convergence criteria[CNVTOL] are satisfied (or until the maximum number of equilibrium equations is reached [NEQIT]). You can define custom criteria if the default settings are not suitable.
ANSYS' automatic solution control uses L2-norm of force (and moment) tolerance (TOLER) equal to 0.5%, a setting that is appropriate for most cases. In most cases, an L2-norm check on displacement with TOLER equal to 5% is also used in addition to the force norm check. The check that the displacements are loosely set serves as a double-check on convergence.
By default, the program will check for force (and, when rotational degrees of freedom are active, moment) convergence by comparing the square root sum of the squares (SRSS) of the force imbalances against the product of VALUE*TOLER. The default value of VALUE is the SRSS of the applied loads (or, for applied displacements, of the Newton-Raphson restoring forces), or MINREF (which defaults to 0.001), whichever is greater. The default value of TOLER is 0.005. IfSOLCONTROL,OFF, TOLER defaults to 0.001 and MINREF defaults to 1.0 for force convergence.
You should almost always use force convergence checking. You can also add displacement (and, when applicable, rotation) convergence checking. For displacements, the program bases convergence checking on the change in deflections (u) between the current (i) and the previous (i-1) iterations: u=ui-ui-1.
Note-If you explicitly define any custom convergence criteria [CNVTOL], the entire default criteria will be overwritten. Thus, if you define displacement convergence checking, you will have to redefine force convergence checking. (Use multipleCNVTOL commands to define multiple convergence criteria.)
Using tighter convergence criteria will improve the accuracy of your results, but at the cost of more equilibrium iterations. If you want to tighten (or loosen, which is not recommended) your criteria, you should change TOLER by one or two orders of magnitude. In general, you should continue to use the default value of VALUE; that is, change the convergence criteria by adjusting TOLER, not VALUE. You should make certain that the default value of MINREF=0.0011.0 makes sense in the context of your analysis. If your analysis uses certain sets of units or has very low load levels, you might want to specify a smaller value for MINREF.
Also, we do not recommend putting two or more disjointed structures into one model for the nonlinear analysis because the convergence check tries to relate these disjointed structures, often producing some unwanted residual force.
Checking Convergence in a Single and Multi-DOF System
To check convergence in a single degree of freedom (DOF) system, you compute the force (and moment) imbalance for the one DOF, and compare this value against the established convergence criteria (VALUE*TOLER). (You can also perform a similar check for displacement (and rotation) convergence for your single DOF.) However, in a multi-DOF system, you might want to use a different method of comparison.
The ANSYS program provides three different vector norms to use for convergence checking:

  • The infinite norm repeats the single-DOF check at each DOF in your model.
  • The L1 norm compares the convergence criterion against the sum of the absolute values of force (and moment) imbalance for all DOFs.
  • The L2 norm performs the convergence check using the square root sum of the squares of the force (and moment) imbalances for all DOFs. (Of course, additional L1 or L2 checking can be performed for a displacement convergence check.)
Example
For the following example, the substep will be considered to be converged if the out-of-balance force (checked at each DOF separately) is less than or equal to 5000*0.0005 (that is, 2.5), and if the change in displacements (checked as the square root sum of the squares) is less than or equal to 10*0.001 (that is, 0.01).

 CNVTOL,F,5000,0.0005,0
 CNVTOL,U,10,0.001,2

Maximum Number of Equilibrium Iterations [NEQIT]

ANSYS' automatic solution control sets this value to between 15 and 26 iterations, depending upon the physics of the problem. The idea is to employ a small time step with fewer quadratically converging iterations.
This option limits the maximum number of equilibrium iterations to be performed at each substep (default = 25 if solution control is off). If the convergence criteria have not been satisfied within this number of equilibrium iterations, and if auto time stepping is on [AUTOTS], the analysis will attempt to bisect. If bisection is not possible, then the analysis will either terminate or move on to the next load step, according to the instructions you issue in the NCNV command.

Predictor-Corrector Option [PRED]

ANSYS' automatic solution control will set PRED, ON if there are no beam or shell elements present. If the time step size is reduced greatly in the current substep,PRED is turned off. For transient analysis, the predictor is also turned off.
You can activate a predictor on the DOF solution for the first equilibrium iteration of each substep. This feature accelerates convergence and is particularly useful if nonlinear response is relatively smooth. It is not so helpful in analyses that incorporate large rotations or viscoelasticity. Using the predictor for large rotations can cause divergence and thus is not recommended for problems with large rotations.

Line Search Option [LNSRCH]

ANSYS' automatic solution control will toggle line search on and off as needed. For most contact problems, LNSRCH is toggled on. For most non-contact problems,LNSRCH is toggled off.
This convergence-enhancement tool multiplies the calculated displacement increment by a program-calculated scale factor (having a value between 0 and 1), whenever a stiffening response is detected. Because the line search algorithm is intended to be an alternative to the adaptive descent option [NROPT], adaptive descent is not automatically activated if the line search option is on. We do not recommend activating both line search and adaptive descent simultaneously.
When an imposed displacement exists, a run cannot converge until at least one of the iterations has a line search value of 1. ANSYS scales the entire U vector, including the imposed displacement value; otherwise, a "small" displacement would occur everywhere except at the imposed DOF). Until one of the iterations has a line search value of 1, ANSYS does not impose the full value of the displacement.

Creep Criteria [CRPLIM,CRCR]

If your structure exhibits creep behavior, you can specify a creep criterion for automatic time step adjustment. (If automatic time stepping [AUTOTS] is off, this creep criterion will have no effect.) The program will compute the ratio of creep strain increment (cr, the change in creep strain in the last time step) to the elastic strain (el), for all elements. If the maximum ratio is greater than the criterion CRCR, the program will then decrease the next time step size; if it is less, the program might increase the next time step size. (The program will also base automatic time stepping on the number of equilibrium iterations, impending element status change, and plastic strain increment. The time step size will be adjusted to the minimum size calculated for any of these items.) If the ratio cr/el is above the stability limit of 0.25, and if the time increment cannot be decreased, a divergent solution is possible and the analysis will be terminated with an error message. This problem can be avoided by making the minimum time step size sufficiently small [DELTIMNSUBST
].
Note-If you do not want to include the effects of creep in your analysis, set the time steps to be higher than the previous time step, but not more than 1.0e-6 higher.

Cutback Criteria [CUTCONTROL]

This command provides finer control over bisections and cutback in time step size. By default, the "plslimit" (maximum plastic strain increment limit) is set to 15%. This field is set to such a large value for avoiding unnecessary bisections caused by high plastic strain due to a local singularity which is not normally of interest to the user. The "crplimit" (creep increment limit) is set to 10%, as before; this is a reasonable limit for creep analysis. The number of points per cycle for second order dynamic equations is set to NPOINT=13 by default to gain efficiency at little cost to accuracy.

Time Step Open Control [OPNCONTROL]

This command is available for thermal analysis. It's primary use is in unsteady state thermal analysis where the final temperature stage reaches a steady state. In such cases, the time step can be opened quickly. The default is that if the TEMP increment is smaller than 0.1 in three (NUMSTEP=3) contiguous substeps, the time step size can be "opened-up" (value =0.1 by default). The time step size can then be opened continuously for greater solution efficiency.

Solution Monitoring [MONITOR]

This command provides a facility to monitor a solution value at a specified node in a specified DOF. The command also provides the user with a means to quickly review the solution convergence efficiency, rather than attempting to gather this information from a lengthy output file. For instance, if an excessive number of attempts were made for a substep, the information contained in the file provides hints to either reduce the initial time step size or increase the minimum number of substeps allowed through the NSUBST command to avoid an excessive number of bisections.

3 Review the Results

Results from a nonlinear static analysis consist mainly of displacements, stresses, strains, and reaction forces. You can review these results in POST1, the general postprocessor, or in POST26, the time-history postprocessor.
Remember that in POST1, only one substep can be read in at a time, and that the results from that substep should have been written to Jobname.RST. (The load step option command OUTRES controls which substep results are stored onJobname.RST.) A typical POST1 postprocessing sequence is described below.

3.1 Points to Remember

  • To review results in POST1, the database must contain the same model for which the solution was calculated.
  • The results file (Jobname.RST) must be available.

3.2 Reviewing Results in POST1

1. Verify from your output file (Jobname.OUT) whether or not the analysis converged at all load steps.

  • If not, you probably won't want to postprocess the results, other than to determine why convergence failed.
  • If your solution converged, then continue postprocessing.
2. Enter POST1. If your model is not currently in the database, issue RESUME
Command(s):

    /POST1
GUI:
Main Menu>General Postproc
3. Read in results for the desired load step and substep, which can be identified by load step and substep numbers or by time. (Note, however, that arc-length results should not be identified by time.
Command(s):

    SET
GUI:
Main Menu>General Postproc>-Read Results-load step

    You can also use the SUBSET or APPEND commands to read in or merge results data for selected portions of the model only. The LIST argument on any of these commands lists the available solutions on the results file. You can also limit the amount of data written from the results file to the database through theINRES command. Additionally, you can use the ETABLE command to store result items for selected elements. See the individual command descriptions in the ANSYS Commands Reference for more information.
Caution: If you specify a TIME value for which no results are available, the ANSYS program will perform a linear interpolation to calculate the results at that value of TIME. Realize that this interpolation will usually cause some loss of accuracy in a nonlinear analysis (see Figure 8-19). Thus, for a nonlinear analysis, you should usually postprocess at a TIME that corresponds exactly to the desired substep.
Figure 8-19 Linear interpolation of nonlinear results can introduce some error
4. Display the results using any of the following options:

3.3 Option: Display Deformed Shape

Command(s):

    PLDISP
GUI:
Main Menu>General Postproc>Plot Results>Deformed Shape

    In a large deformation analysis, you might prefer to use a true scale display [/DSCALE,,1].

3.4 Option: Contour Displays

Command(s):

    PLNSOL or PLESOL
GUI:
Main Menu>General Postproc>Plot Results>-Contour Plot-Nodal Solu orElement Solu

    Use these options to display contours of stresses, strains, or any other applicable item. If you have adjacent elements with different material behavior (such as can occur with plastic or multilinear elastic material properties, with different material types, or with adjacent deactivated and activated elements), you should take care to avoid nodal stress averaging errors in your results. Selecting logic (described in the ANSYS Basic Analysis Procedures Guide) provides a means of avoiding such errors.
    The KUND field on PLNSOL and PLESOL gives you the option of overlaying the undeformed shape on the display.
    You can also contour element table data and line element data:
Command(s):

    PLETABPLLS
GUI:
Main Menu>General Postproc>Element Table>Plot Element Table
Main Menu>General Postproc>Plot Results>-Contour Plot-Line Elem Res

    Use the PLETAB command (Main Menu>General Postproc>Element Table>Plot Element Table) to contour element table data and PLLS (Main Menu>General Postproc>Plot Results>Line elem Res) to contour line element data.

3.5 Option: Tabular Listings

Command(s):

    PRNSOL (nodal results)
    PRESOL (element-by-element results)
    PRRSOL (reaction data)
    PRETAB
    PRITER (substep summary data), etc.
    NSORT
    ESORT
GUI:
Main Menu>General Postproc>List Results>Nodal Solution
Main Menu>General Postproc>List Results>Element Solution
Main Menu>General Postproc>List Results>Reaction Solution

    Use the NSORT and ESORT commands to sort the data before listing them.

3.6 Other Capabilities

Many other postprocessing functions-mapping results onto a path, report quality listings, and so on-are available in POST1. See Chapter 5 of the ANSYS Basic Analysis Procedures Guide for details. Load case combinations usually are not valid for nonlinear analyses.

3.7 Reviewing Results in POST26

31.You can also review the load-history response of a nonlinear structure using POST26, the time-history postprocessor. Use POST26 to compare one ANSYS variable against another. For instance, you might graph the displacement at a node versus the corresponding level of applied load, or you might list the plastic strain at a node and the corresponding TIME value. A typical POST26 postprocessing sequence might follow these steps:
1. Verify from your output file (Jobname.OUT) whether or not the analysis converged at all desired load steps. You should not base design decisions on unconverged results.
2. If your solution converged, enter POST26. If your model is not currently in the database, issue RESUME.
Command(s):

    /POST26
GUI:
Main Menu>TimeHist Postpro
3. Define the variables to be used in your postprocessing session. The SOLUcommand will cause various iteration and convergence parameters to be read into the database, where you can incorporate them into your postprocessing.
Command(s):

    NSOL
    ESOL
    RFORCE
GUI:
Main Menu>TimeHist Postpro>Define Variables
4. Graph or list the variables.
Command(s):

    PLVAR (graph variables)
    PRVAR
    EXTREM (list variables)
GUI:
Main Menu>TimeHist Postpro>Graph Variables
Main Menu>TimeHist Postpro>List Variables
Main Menu>TimeHist Postpro>List Extremes

3.8 Other Capabilities

Many other postprocessing functions are available in POST26. See Chapter 6 of the ANSYS Basic Analysis Procedures Guide for details.
See the NLGEOMSSTIFNROPTTIMENSUBSTAUTOTSKBCCNVTOL,NEQITNCNVPREDOUTRES, and SOLU command descriptions for more information.

4 Terminating a Running Job; Restarting

You can stop a nonlinear analysis by creating an "abort" file (Jobname.ABT). See Chapter 3 of the ANSYS Basic Analysis Procedures Guide for details. The program will also stop upon successful completion of the solution, or if a convergence failure occurs.
You can often restart an analysis if it successfully completed one or more iterations before it terminated. Restart procedures are covered in Chapter 3 of the ANSYS Basic Analysis Procedures Guide.

Selasa, 21 Februari 2012

ANSYS 13.0 Enhanced Modal Analyses with Linear Perturbation

ANSYS 13.0 Enhanced Modal Analyses with Linear Perturbation

By Ted Harris24. March 2011 09:38
By Ted Harris

“Rest, rest, perturbed spirit!”  - William Shakespeare
If you have ever performed a large deflection prestressed modal analysis in ANSYS Mechanical APDL prior to version 13.0, your spirit might have been perturbed as well.  The procedure was not very user friendly, to sum it up.  For example, unless you were careful, the modal results would over-write the static preload results.  Thankfully, at 13.0 we have a smoother and more capable tool for handling large deflection prestressed modal analyses.  This new procedure is called Linear Perturbation.
We’ll focus on modal analyses in this article, but be aware that linear perturbation also applies to linear buckling analyses at 13.0, but only following a linear preload solution, and only in Workbench.  The capability for modal analyses is supported in both Workbench and Mechanical APDL.  Also, the preload, or ‘base’ analysis has to have multiframe restart capability turned on.  This will happen by default for a nonlinear analysis but needs to be manually activated in MAPDL for a linear analysis by issuing the command RESCONTROL,LINEAR.
In fairly simple terms, the prestress effects are included in a modal analysis via the change in the stiffness matrix that occurs during the prestress (typically nonlinear static) analysis.  This is the method that has been used in ANSYS for years.  What’s new at 13.0 is that the program keeps track of different components of the augmented tangent stiffness matrix.  The five possible contributing components are material property effects, stress stiffening effects, load stiffening effects, contact element effects, and spin softening effects.
While the material effects must remain linear, the contact stiffness can be altered, if desired, in the subsequent modal analysis.  More on that later.
The typical Mechanical APDL procedure to perform a nonlinear static structural prestress run followed by a modal analysis which utilizes those prestress effects is as follows:
! With static model prepared
/solu                            ! enter solution module
antype,0                         ! specify static analysis type
nlgeom,1                         ! turn on large deflection effects (nonlinear)
pstres,on                        ! turn on prestress effects for subsequent modal
nsub,10,10,10                    ! specify substep range
save                             ! save the database
solve                            ! solve the nonlinear static prestress case
finish                           ! leave the solution module
/solu                            ! re-enter solution so we can do a new analysis
antype,,restart,1,10,perturb     ! specify restart option for linear perturbation
                                 ! from last substep in this case
perturb,modal                    ! specify modal as next analysis
solve,elform                     ! calculate element formulation with solve command
modopt,lanb,12                   ! specify modal options for solution
mxpand,12                        ! specify number of modes for results calc
solve                            ! solve the prestress modal analysis
/post1                           ! enter general postprocessor
INRES,ALL                        ! make sure we read in all results from file
FILE,'nonlinear_static','rstp'   ! specify special results file for modal results
set,first                        ! read in results for first mode
plns,u,sum                       ! plot mode shape
set,next                         ! read in results for next mode
/repl                            ! plot mode shape, etc.

Note that the linear perturbation (modal) analysis has its own results file, the .rstp file.  Because of this, the preload results are still available in their own .rst file as it does not get overwritten by the modal step.
Here is a table of frequency results for a simple test case.  Three modal analyses were run:
1.  No prestress at all.
2.  With a linear static prestress state.
3.  With a nonlinear static prestress state.
image
Here is the model used for these runs in its initial configuration.  The block at the base was fixed in all DOF’s and the preload applied was a pressure load on one side of the vane.
undeformed
Here is the model with the deformed mesh due to the nonlinear prestress:
nonlin_prestress_deformed
Here is a mode shape plot for mode 12:
nonlin_prestress_mode12
The above example is all well and good but could have been done in prior versions of ANSYS using the old partial solve method.  What’s nice about the newer linear perturbation method is that it’s easy to get the mode shape plots relative to the deformed mesh from the prior prestress run, and you don’t need to worry about over-writing the prestress results with the modal results, since the corresponding results files are different.
Further, we can now perform modal analyses using different restart points in the static prestress run, assuming multiple restart points are available.
Finally, we can actually change some contact options between the static prestress solution and the modal solution.  For example, if the prestress analysis was run using frictional contact, the subsequent modal analysis can be run utilizing the prestressed state of the structure but with one of three contact states for the modal analysis: true status (that of the prior static analysis), force (to be) sticking, or  force (to be) bonded. The sticking option applies only to contact elements with a nonzero coefficient of friction. The bonded option will force contact pairs that are in contact in the static analysis to be bonded in the modal analysis.
The Mechanical APDL command sequence for this procedure would be something like this:

! first perform nonlinear static prestress run, then
/solu
antype,,restart,,,perturbation
perturb,modal,,BONDED,PARDELE                    ! pre-stress modal analysis, switch contact to bonded, delete
                                                                                              ! loads in case future MSUP
solve                                                                ! Generate matrices needed for perturbation analysis
! Next perform modal analysis
modopt,lanb,6
mxpand                    ! default expand in case of complex solution
solve
! modal results are now available for postprocessing
In Workbench Mechanical, the appropriate command sequence is sent to the solver when we link a modal analysis to a prior prestress analysis.  If the model involves contact, then in the modal analysis we’ll have choices for how the contact should be treated in the Pre-Stress branch under the Modal branch in the Outline Tree.  For frictional contact in the static prestress analysis, the choices in the Details view for the Pre-Stress branch in the modal analysis will be Use True Stress, Force Sticking, or Force Bonded as described above.
Here are some example plots for this scenario:
Two 3D plates subject to in plane bending, fixed at right ends, frictional contact between them.
image

Resulting contact status for static run (sliding is occurring)
image
Resulting static deformation:
def_c

Mode 6 result, “true” contact behavior:
mode6_nonlin_true_c.

Mode 6 result, “force bonded” contact behavior:
mode6_nonlin_force_bonded_c
Those last two images show a dramatic difference in modal results simply by changing the contact status behavior in the modal analysis.  In the first of those images, the contact status is set to ‘true’, meaning essentially the same as in the prestress analysis, subject to the linear nature of the modal analysis.  In this example, the frictional behavior in the static prestress run becomes no separation in the modal analysis, so the two plates can have mode shapes in which the plates slide relative to each other.  In the last plot, the contact status has been changed to ‘force bonded’ for the modal solution.  As the plot shows, mode shapes can only exist in which the two plates are bonded together.  Both modal analyses have the same prestress condition however.
Here is a frequency table comparing the first six modes of the two modal analyses.  Note that with the contact forced to be bonded we get a stiffening response as we might expect.
WB_freq_compare
So, although on the surface it might initially appear to be a black art, linear perturbation is a nice enhancement in ANSYS 13.0 that gives us a more robust and capable method for performing modal analyses with prestress effects included.  The prestress run is typically a linear or nonlinear static analysis, but it will also work with a full transient analysis to define the prestress state.  The ANSYS 13.0 Help has more information (see section 9.2 of the Mechanical APDL Structural Analysis Guide and section 17.8 of the Theory Reference).  We also recommend you try out the procedure on your own models.
Currently rated 5.0 by 5 people


ANSYS 13.0 Analisis Modal Ditingkatkan dengan kekacauan Linear

Oleh Ted Harris24.Maret 2011 09:38
Oleh Ted Harris

"Istirahat, istirahat, gelisah semangat!" - William Shakespeare
Jika Anda pernah melakukan defleksi besar pratekan analisis modal di APDL ANSYS Teknik sebelum versi 13,0, roh Anda mungkin telah terganggu juga. Prosedur ini tidak terlalu user friendly, untuk jumlah itu.Sebagai contoh, kecuali Anda berhati-hati, hasil modal akan over-menulis hasil preload statis. Untungnya, pada 13,0 kami memiliki suatu utilitas halus dan lebih mampu untuk menangani defleksi besar pratekan analisis modal. Prosedur baru ini disebut kekacauan Linear.
Kami akan fokus pada analisis modal dalam artikel ini, namun perlu diingat bahwa gangguan linier juga berlaku untuk analisis buckling linear pada 13,0, tapi hanya mengikuti solusi preload linier, dan hanya dalam Workbench. Kemampuan untuk analisis modal didukung di kedua Workbench dan APDL Mesin. Juga, preload, atau 'dasar' analisis harus memiliki multiframe restart kemampuan dihidupkan. Ini akan terjadi secara default untuk analisis nonlinier tetapi perlu diaktifkan secara manual di MAPDL untuk analisis linier dengan menerbitkan RESCONTROL perintah, LINEAR.
Dalam hal yang cukup sederhana, efek prategang termasuk dalam analisis modal melalui perubahan dalam matriks kekakuan yang terjadi selama prategang (biasanya nonlinier statis) analisis. Ini adalah metode yang telah digunakan dalam ANSYS selama bertahun-tahun. Apa yang baru di 13,0 adalah bahwa program melacak komponen yang berbeda dari matriks kekakuan tangen ditambah. Lima komponen memberikan kontribusi yang mungkin adalah efek properti material, efek stres kaku, efek beban kaku, efek kontak elemen, dan spin efek pelunakan.
Sementara dampak bahan harus tetap linear, kekakuan kontak dapat diubah, jika diinginkan, dalam analisis modal berikutnya. Lebih pada nanti.
Prosedur APDL khas Mesin untuk melakukan lari pratekan statis nonlinier struktural diikuti dengan analisis modal yang memanfaatkan efek-efek prategang adalah sebagai berikut:
! Dengan model statis disiapkan
/ Solu! memasukkan modul solusi
antype, 0! menentukan jenis analisis statis
nlgeom, 1! mengaktifkan efek defleksi besar (nonlinier)
pstres, pada! mengaktifkan efek prategang untuk modal berikutnya
nsub, 10,10,10! tentukan rentang substep
save! menyimpan database
memecahkan! menyelesaikan kasus pratekan statis nonlinier
selesai! meninggalkan modul solusi
/ Solu! masukkan kembali solusi sehingga kami dapat melakukan analisis baru
antype,, restart, 1,10, mengacaukan! tentukan pilihan untuk restart perturbasi linier
! dari substep terakhir dalam kasus ini
mengacaukan, modal! menentukan modal sebagai analisis berikutnya
memecahkan, elform! menghitung formulasi elemen dengan memecahkan perintah
modopt, lanb, 12! menentukan pilihan moda untuk solusi
mxpand, 12! menentukan sejumlah mode untuk hasil kalk
memecahkan! memecahkan analisis modal prategang
/ Post1! masuk postprocessor umum
INRES, SEMUA! pastikan kita baca dalam semua hasil dari file
FILE, 'nonlinear_static', 'RSTP'! menentukan file khusus hasil untuk hasil modal
ditetapkan, terlebih dahulu! dibaca dalam hasil untuk mode pertama
plns, u, jumlah! Plot modus bentuk
ditetapkan, berikutnya! dibaca dalam hasil untuk mode berikutnya
/ Repl! Plot modus bentuk, dll

Perhatikan bahwa perturbasi analisis (modal) linier memiliki file sendiri hasil, RSTP file.. Karena itu, hasil preload masih tersedia dalam sendiri. Terlebih dulu file mereka karena tidak mendapatkan ditimpa oleh langkah modal.
Berikut adalah tabel hasil frekuensi untuk kasus tes sederhana. Tiga analisis modal dijalankan:
1. Tidak ada pratekan sama sekali.
2. Dengan keadaan prategang linier statis.
3. Dengan keadaan nonlinier prategang statis.
image
Berikut adalah model yang digunakan untuk berjalan dalam konfigurasi awal. Blok di pangkalan itu tetap dalam semua itu DOF dan preload terapan adalah beban tekanan pada satu sisi baling-baling tersebut.
undeformed
Berikut adalah model dengan mesh cacat karena prategang nonlinier:
nonlin_prestress_deformed
Berikut ini adalah bentuk modus plot untuk modus 12:
nonlin_prestress_mode12
Contoh di atas adalah semua baik dan bagus tapi bisa dilakukan dalam versi sebelumnya dari ANSYS menggunakan metode parsial memecahkan tua. Apa yang baik tentang metode perturbasi lebih baru linier adalah bahwa hal itu mudah untuk mendapatkan plot modus getar relatif terhadap jaring cacat dari jangka prategang sebelumnya, dan Anda tidak perlu khawatir tentang over-menulis hasil prategang dengan hasil modal, karena hasil yang sesuai file yang berbeda.
Selanjutnya, kita sekarang dapat melakukan analisis moda menggunakan poin Restart berbeda dalam jangka pratekan statis, dengan asumsi poin Restart beberapa yang tersedia.
Akhirnya, kita sebenarnya dapat mengubah opsi kontak beberapa antara larutan pratekan statis dan solusi modal. Misalnya, jika analisis prategang tersebut diterjemahkan menggunakan kontak gesekan, analisis modal berikutnya dapat dijalankan memanfaatkan keadaan pratekan struktur tetapi dengan satu dari tiga negara bagian kontak untuk analisis modal: Status benar (bahwa dari analisis statis sebelumnya), gaya (menjadi) mencuat, atau kekuatan (menjadi) berikat. Pilihan menempel hanya berlaku untuk menghubungi elemen dengan koefisien nol gesekan. Pilihan terikat akan memaksa menghubungi pasangan yang berada dalam kontak dalam analisis statis untuk ditempelkan dalam analisis modal.
Urutan perintah Mesin APDL untuk prosedur ini akan menjadi sesuatu seperti ini:

! pertama kali melakukan lari prategang nonlinier statis, maka
/ Solu
antype,, restart,,, perturbasi
mengacaukan, modal,, terikat, PARDELE                     ! pra-stres analisis modal, beralih kontak ke berikat, menghapus
                                                                                              ! beban dalam MSUP kasus masa depan
memecahkan                                                                 ! Menghasilkan matriks yang diperlukan untuk analisis perturbasi
! Berikutnya melakukan analisis modal
modopt, lanb, 6
mxpand                     ! standar memperluas dalam hal solusi kompleks
memecahkan
! Hasil modal sekarang tersedia untuk postprocessing
Dalam Workbench Mesin, urutan perintah yang sesuai dikirim ke pemecah ketika kita menghubungkan analisis modal untuk analisis prategang sebelumnya. Jika model melibatkan kontak, maka dalam analisis modal kita akan memiliki pilihan untuk bagaimana kontak harus dirawat di cabang Pra-Stress bawah cabang Modal di Pohon Outline. Untuk kontak gesekan dalam analisis pratekan statis, pilihan pada tampilan Rincian untuk cabang Pra-Stress dalam analisis modal akan Gunakan Stres Benar, rekat Force, atau Angkatan Berikat sebagaimana dijelaskan di atas.
Berikut adalah contoh beberapa plot untuk skenario ini:
Dua 3D pelat dikenakan dalam pesawat lentur, tetap pada ujung kanan, kontak gesekan antara mereka.
image
Hasil status kontak untuk jangka statis (sliding ini terjadi)
image
Hasil deformasi statis:
def_c

Mode 6 hasil, "benar" perilaku hubungi:
mode6_nonlin_true_c.

Mode 6 hasil, "kekuatan ikatan" hubungi perilaku:
mode6_nonlin_force_bonded_c
Kedua terakhir gambar menunjukkan perbedaan yang dramatis dalam hasil modal hanya dengan mengubah perilaku kontak status dalam analisis modal. Pada bagian pertama dari gambar tersebut, status kontak diatur ke 'benar', yang berarti pada dasarnya sama seperti pada analisis prategang, sesuai dengan sifat linier dari analisis modal. Dalam contoh ini, perilaku gesekan dalam jangka pratekan statis menjadi ada pemisahan dalam analisis modal, jadi dua piring dapat memiliki bentuk modus di mana pelat geser relatif terhadap satu sama lain. Dalam plot terakhir, status kontak telah berubah menjadi 'kekuatan berikat' untuk solusi modal. Sebagai petak menunjukkan, bentuk modus hanya bisa ada di mana dua lempeng terikat bersama-sama. Kedua analisis modal memiliki kondisi prategang yang sama namun.
Berikut ini adalah tabel frekuensi membandingkan enam mode pertama dari dua analisis modal. Perhatikan bahwa dengan kontak dipaksa untuk ditempelkan kita mendapatkan respon kaku seperti yang kita harapkan.
WB_freq_compare
Jadi, meskipun di permukaan itu mungkin awalnya tampak seperti ilmu hitam, gangguan linear adalah perangkat tambahan yang bagus di ANSYS 13,0 yang memberi kita sebuah metode yang lebih kuat dan mampu untuk melakukan analisis modal dengan efek prategang disertakan. Jangka prategang biasanya sebuah analisis statis linier atau nonlinier, tetapi juga akan bekerja dengan analisis transien penuh untuk menentukan negara prategang. ANSYS 13,0 Bantuan memiliki informasi lebih lanjut (lihat bagian 9.2 dari Petunjuk Teknik Analisis APDL Struktural dan bagian 17,8 dari Referensi Teori). Kami juga merekomendasikan Anda mencoba prosedur pada model Anda sendiri.
Currently rated NaN by 5 people
  • Currently 5.0 b4b74173-b8ba-434f-8785-0d605bca78e0 /5 Stars.
  • 1
  • 2
  • 3
  • 4
  • 5

Add New Comment

Showing 0 comments

    blog komentar powered by DISQUS

    Berlangganan

    Untuk Email Marketing dapat Anda percaya

    PADT yang ANSYS Webinar Seri

    Jan 12, 2012 - 12:00 MST
    Update pada Seleksi Dinamakan untuk ANSYS Teknik R14
    Jan 26, 2012 - 12:00 MST
    Manajemen Memori di ANSYS
    Feb 9, 2012 - 12:00 MST
    Bekerja langsung dengan Nodes dan Elemen di ANSYS Teknik
    Feb 23, 2012 - 12:00 MST
    Majelis Meshing di ANSYS R14MEMBATALKAN
    8 Maret 2012 - 12:00 MST
    Pengenalan Workbench Kerangka Scripting - Mengontrol proyek, bahan, dan eksekusi solusi dengan python
    22 Maret 2012 - 12:00 MST
    Menguasai Manajer Jauh Solver (RSM) di R14

    Lihat daftar lengkap beserta link ke rekaman Webinars masa lalu di:
    dan klik "PADT ANSYS Webinar Seri"

    Dibawa ke anda oleh PADT

    PADT dengan bangga menerbitkan Fokus. Kami harap Anda menemukan nilai di dalamnya. Jika Anda melakukannya, kami meminta Anda mempertimbangkan kami bila Anda membutuhkan jasa. Meskipun Anda mungkin tahu kami untuk kegiatan simulasi kami, kami juga merupakan pengembangan produk layanan penuh dan penyedia prototyping cepat.

    Simulasi
    outsourcing
    Pelatihan
    Kustomisasi

    Rapid Protoypting 
    - SLA, SLS, FDM
    PolyJet
    - Tooling Lunak
    - Injection Molding
    Pengembangan Produk
    - Desain
    - Pengujian
    - Manufaktur
    - Comercialization
    - Spesialisasi dalam
    - Alat Kesehatan
    - Rotating Mesin
    - Energi Alternatif
    - Peralatan Semikonduktor

    Pelajari lebih lanjut diwww.PADTINC.com

    BATU HVPC Sistem

    Langkah Keluar dari Kotak 
    Langkah Ke Cube
    PADT menawarkan garis penuh sistem komputer yang dirancang khusus untuk pengguna FEA dan CFD. Komputer ini CUBE Kinerja Tinggi Nilai (HVPC) dikonfigurasi oleh staf TI PADT untuk memukul sweet spot antara performa dan biaya untuk pengguna simulasi.
    BARU: 
    Kami sekarang menawarkan Sistem Berbasis Intel 
    Mari kita Konfigurasi dan Quote satu untuk Anda hari ini!
    mini-Cluster
    96 core / 256 GB RAM / 3,6 Disk TB
    Handphone Rack / UPS / Monitor / Keyboard
    $ 43.250
    Hitung Server
    32 core / 128 GB RAM / 3 Disk TB
    $ 12.300
    Simulasi Workstation
    12 core / 64 GB RAM / 1,5 Disk TB
    $ 5.800
    Simulasi fileserver
    10 TB Disk / eSATA Eksternal
    $ 5.800
    Lebih Sistem Tersedia
    Mengunjungi
    www.CUBE-HVPC.com
    Untuk Informasi lebih

    Pedal: Editor APDL

    Dapatkan Sosial dengan PADT

     
    PADT di Facebook PADT di Web
    PADT di Twitter PADT E-Mail Langganan
    PADT pada Linked In Blog Fokus

    Latihan

    Dapatkan Dilatih oleh Orang yang Menulis Fokus
    Cari tahu mengapa pengguna produk ANSYS, Inc dari seluruh dunia datang ke PADT untuk pelatihan mereka. Kami menawarkan kelas di hampir setiap produk dengan harga yang bersaing dan dapat membuat kelas kustom sesuai kebutuhan. Teknisi kami dapat datang ke fasilitas Anda atau Anda bisa datang ke kantor kami di Tempe AZ atau Littleton CO
    Silakan melihat kami jadwal kami padahalaman pelatihan , atau hanyamenghubungi kami dan kami dapat berbicara tentang memenuhi kebutuhan spesifik Anda.

    Bulan Daftar

    • 2002
    • 2003
    • 2004
    • 2005
    • 2006
    • 2007
    • 2008
    • 2009
    • 2010
    • 2011
    • 2012

    Halaman Daftar