15.11.2012 Views

Mechanical APDL Basic Analysis Guide - Ansys

Mechanical APDL Basic Analysis Guide - Ansys

Mechanical APDL Basic Analysis Guide - Ansys

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 5: Solution<br />

The PCG solver primarily solves for displacements/rotations (in structural analysis), temperatures (in thermal<br />

analysis), etc. The accuracy of other derived variables (such as strains, stresses, flux, etc.) is dependent upon<br />

accurate prediction of primary variables. Therefore, ANSYS uses a very conservative setting for PCG tolerance<br />

(defaults to 1.0E-8) The primary solution accuracy is controlled by the PCG. For most applications, setting<br />

the PCG tolerance to 1.0E-6 provides a very accurate displacement solution and may save considerable CPU<br />

time compared with the default setting. Use the EQSLV command to change the PCG solver tolerance.<br />

Direct solvers (such as the sparse direct solver) produce very accurate solutions. Iterative solvers, such as<br />

the PCG solver, require that a PCG convergence tolerance be specified. Therefore, a large relaxation of the<br />

default tolerance may significantly affect the accuracy, especially of derived quantities.<br />

The PCG solver does not support SOLID62 elements.<br />

With all iterative solvers you must verify that the model is appropriately constrained. No minimum pivot is<br />

calculated and the solver will continue to iterate if any rigid body motion exists.<br />

In a modal analysis using the PCG solver (MODOPT,LANPCG), the number of modes should be limited to<br />

100 or less for efficiency. PCG Lanczos modal solutions can solve for a few hundred modes, but with less<br />

efficiency than Block Lanczos (MODOPT,LANB).<br />

When the PCG solver encounters an indefinite matrix, the solver will invoke an algorithm that handles indefinite<br />

matrices. If the indefinite PCG algorithm also fails (this happens when the equation system is ill-conditioned;<br />

for example, losing contact at a substep or a plastic hinge development), the outer Newton-Raphson<br />

loop will be triggered to perform a bisection. Normally the stiffness matrix will be better conditioned after<br />

bisection and the PCG solver can eventually solve all the nonlinear steps.<br />

The solution time grows linearly with problems size for iterative methods so huge models can still be solved<br />

within very reasonable times. For modal analyses of large models (e.g., 10 million DOF or larger), MOD-<br />

OPT,LANPCG is a viable solution method if the number of modes is limited to approximately 100.<br />

Use MSAVE,ON (the default in most cases) for memory savings of up to 70 percent. The MSAVE command<br />

causes an element-by-element approach (rather than globally assembling the stiffness matrix) for the parts<br />

of the structure using SOLID185, SOLID186, SOLID187, SOLID272, SOLID273, and/or SOLID285 elements that<br />

have linear material properties. This feature applies only to static analyses or modal analyses using the PCG<br />

Lanczos method. (You specify these analysis types using the commands ANTYPE,STATIC, or ANTYPE,MODAL;<br />

MODOPT,LANPCG respectively.) When using SOLID186 and/or SOLID187, only small strain (NLGEOM,OFF)<br />

analyses are allowed. The solution time may be affected depending on the processor speed and manufacturer<br />

of your computer, as well as the chosen element options.<br />

5.2.3. The Jacobi Conjugate Gradient (JCG) Solver<br />

The JCG solver also starts with element matrix formulation. Instead of factoring the global matrix, the JCG<br />

solver assembles the full global stiffness matrix and calculates the DOF solution by iterating to convergence<br />

(starting with an initial guess solution for all DOFs). The JCG solver uses the diagonal of the stiffness matrix<br />

as a preconditioner. The JCG solver is typically used for thermal analyses and is best suited for 3-D scalar<br />

field analyses that involve large, sparse matrices.<br />

For some cases, the tolerance default value (set via the EQSLV,JCG command) of 1.0E-8 may be too restrictive,<br />

and may increase running time needlessly. The value 1.0E-5 may be acceptable in many situations.<br />

The JCG solver is available only for static analyses, full harmonic analyses, or full transient analyses. (You<br />

specify these analysis types using the commands ANTYPE,STATIC, HROPT,FULL, or TRNOPT,FULL respectively.)<br />

You cannot use this solver for coupled-field applications (SOLID5 or PLANE13).<br />

102<br />

Release 13.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information<br />

of ANSYS, Inc. and its subsidiaries and affiliates.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!