去评论
距米网-精简版

How to Define AUTODYN solution controls

JUMU
2019/11/17 19:05:19
AUTODYN  Controls

This panel lets you define solution controls for your model.

Wrapup Criteria

These controls are the only ones visible when you open this panel for the first time. This is because you must specify some of the parameters.
(all other control inputs have default values that are usually OK).

Cycle limit

Enter the maximum number of cycles you want your problem to run (enter a large integer if you don't want your problem to stop on a cycle limit).

Time limit

Enter the maximum time you want your problem to run (enter a large number - e.g. 1.0E20 - if you don't want your problem to stop on a time limit).

Energy fraction

Enter the energy fraction you want to use to stop your problem if the energy error gets too large. The default value of 0.05 causes you problem to stop if the energy error exceeds 5%.

Energy reference cycle

Enter the cycle # at which you want AUTODYN to start checking for an energy error.


Press  alongside any of the remaining options to reveal the input fields for that option :

Timestep Options

These options control the timestep AUTODYN uses to run your problem.

Start time

Enter the start time for your problem.

Minimum timestep

Enter the minimum timestep you want to use.

If the timestep drops below this value, your problem will be terminated.

If you leave this field zero, the minimum timestep will be set to 1/10 of the Initial timestep.

Maximum timestep

Enter the maximum timestep you want to use.

AUTODYN will use the minimum of this value or the computed stability timestep.

Initial timestep

Enter the initial timestep you want to use.

If you leave this field zero, the initial timestep will be set to 1/2 the computed stability timestep.

Safety factor

It is not wise to run problem at the stability limit, so a safety factor is applied to the computed stability timestep.

Enter the stability factor you want to use..

The default value (0.6666) is good for all problems, but you may want to increase this value to 0.9 in some cases to use a larger timestep (0.9 is a good value for most Lagrangian calculations).

Method of calculating characteristics zone

Three options are available for calculating the effective zone size used in the timestep calculation.  By default, all models will use a zone size calculation based on cell diagonals.  If the user wishes to obtain improved efficiency, especially in 3D simulations, the zone size calculation can be switched to one of two alternative options (based on either distances between opposing face centers or closest face centres).

Experience to date has shown that these options can significantly improve the efficiency of Lagrange simulations, especially in 3D.  However, in certain circumstances when cells become highly distorted instabilities have been observed causing the calculation to terminate with high energy errors.  The correct choice of erosion strain can reduce these problems.  It is therefore recommended that users only utilize either of the new zonesize options if efficiency is critical.

Automatic Mass Scaling

Mass scaling is an artificial  (numerical) mechanism for increasing the CFL (Courant-Friedrichs-Lewy) timestep of individual elements that govern the maximum allowed timestep of explicit transient dynamic solutions in AUTODYN.  Increasing the timestep has the obvious benefit of reducing the number of cycles required to run a simulation to a given point in time.  Educated use of this option can therefore result in significant improvements in efficiency.

Damping Options



These options control the magnitude of damping applied to your problem.

The default values are best for most problems.

Be very careful if you increase these values, they can significantly change your solution.

Quadratic viscosity

Enter the coefficient of quadratic viscosity you want to apply to shocks (compression waves).

Linear viscosity

Enter the coefficient of linear viscosity you want to apply to shocks (compression waves).

Linear viscosity in expansion

Choose whether you also want to apply the linear viscosity in expansion.

Hourglass damping

Enter the coefficient of houglass damping  you want to apply to Lagrangian elements.

Static damping

A static damping constant may be specified which changes the AUTODYN calculation from a dynamic solution to a relaxation iteration converging to a state of stress equilibrium.  For optimal convergence, the value chosen for the constant, R, may be defined by:

R=2*timstep/T

where timestep is the expected average value of the timestep and T is longest period of vibration for the system being analyzed.

Solver Options



Method for Lagrange/ALE density update

The element density can be updated each cycle using one of two methods; incrementally from the change in volume relative to the last cycle or recalculated based on the current mass and volume ofthe element.  If the incremental change in volume over a cycle is small, it can be more accurate to use the incremental approach rather than the total approach to updating the density.  This especially applies to the single precision version of AUTODYN.

Hence the options are:

Automatic      - AUTODYN decides whether it needs to do an incremental update or not based on the change in volume.

Incremental    - Force AUTODYN to always use the incremental update.

Total               - Force AUTODYN to always recalculate the density.

Method for Euler Strain rate calculation

For multi-material Euler there are two methods to compute strain rates.

Weighted.  For multi-material Euler analysis, involving materials with strength, strain rates need to be calculated to compute the increment in strain from cycle to cycle.  The strain rates are derived from the velocity gradients, which are obtained by an element surface integration of the velocity.  Using the Weighted option the face velocity that is used in the integration is a density weighted average of the element velocity on both sides of the face.

This is the default option and provides the most accurate strain rate calculation.

Averaged.  With this option the face velocity used in the integration is a simple average of the element velocity on both sides of the face and provides less accurate results than the Weighted option.  Until AUTODYN Version 6.1 the Averaged option was the only option available, and can be used to reproduce old results.  Also models that are generated before Version 6.1 will still run using the Averaged method.

Method for Euler pressure calculation

For multi-material Euler cells containing two or more materials, there are two methods for calculating the resultant cell pressure:

Equilibrium        -        An iterative procedure is used for each material in a cell to establish a single pressure consistent with the equation of state and conditions of each material in the cell.  This method is currently only applied if there are two materials in the cell and these materials are defined by ideal gas, or fully burnt JWL material.  If this criteria is not met in a cell the pressure averaging option is used.

Average        -        The resultant cell pressure is calculated from a weighted average of the pressure of the individual materials in the cell.  The weighting algorithm takes account of both the volume fraction of each material and the relative stiffness of each material.  This is the default method.

Global Cutoffs

These options let you set global cutoffs.

The default values are usually good if you are using the default units (mm,mg,ms). You may want to change one or more of these defaults if you are using some other set of units.

Minimum Velocity

Enter the minimum velocity you want to use.

If any velocity drops below this value it is set to zero.

Maximum Velocity

Enter the maximum velocity you want to use.

If any velocity rises above this value it is set to this value.

This option is sometimes useful if your problem develops high velocities, which result in a small timestep, in regions that are not important to your solution.

Radius cutoff

Enter a value for the radius cuttoff for 2D axial models.

If a Lagrangian node is inside this radius at the start of your calculation, it is placed on the axis of symmetry.

If a Lagrangian node is outside this radius at the start of your calculation, it will not be allowed to come closer than this radius to the axis of symmetry as your calculation proceeeds.

SPH Cutoffs

These options let you set SPH cutoffs.

The default values are usually good if you are using the default units (mm,mg,ms). You may want to change one or more of these defaults if you are using some other set of units.

Limit density / Delete node

If the density of an SPH node drops below the minimum density or above the maximum density, you have the choice of either limiting the density to the cutoff density or deleting the node. Make your choice here.

Minimum density

Enter the minimum cutoff density for SPH nodes.

Maximum density

Enter the maximum cutoff density for SPH nodes.

Gravity

This option lets you set gravitational components for your model.

X / Y / Z  acceleration

Enter the X Y and Z components of gravity in these fields.

Transport



These option relate to the transport of material in various Solvers.

Material velocity timestep safety factor

Enter the timestep safety factor for material velocities.

In the velocity used to calculate the stability timestep is c+u, where  c is the soundspeed of the material and u is its physical velocity.

If u is of the order of or greater that c, then using c + u/f can give more stable results, where f is the safety factor.

Euler-FCT

Select the type of transport algorithm that is used for Euler-FCT calculations from this pull-down menu.  The default setting is multi-dimensional where all three spatial dimensions are solved simultaneously which is more efficient and helps to maintain symmetric flow fields.  The alternative is Operator Split where each spatial dimension is solved in turn.

ALE/Euler Energy

During a computational Euler cycle mass, momentum and energy is being transported (fluxed) from element to element.  This option determines if total energy or internal energy is being transported.

Total. This option conserves the total energy (kinetic + internal energy) during transport.  This is the most robust and accurate option and is the default.

Internal. This option conserves the internal energy during transport.  In cases where the total energy transport doesn't provide an accurate answer, e.g. material overheating or when the kinetic energy is predominant, the Internal option might be chosen as an alternative.

Euler

This option determines which transport algorithm will be used for multi-material transport.

SLIC. When material is being transported out of an element that contains multiple materials or voids, it is necessary to use special rules to determine the amount of each material being transported, because the only information available for the transport algorithm is the fraction of each material in the individual elements.  In AUTODYN the Simple Lime Interface Calculation (SLIC) method is used which uses a small number of stereotypes of mixed material element topologies.  The SLIC algorithm is the default and works well for all multi-material Euler analysis.

Updated SLIC. Isolated small parts of material smaller than an element do not move through the mesh at the correct velocity.  They tend to stop moving while surrounding material rushes by.  The Updated SLIC option will transport these small material fractions with statistically the correct velocity.  This option should be used with care, because it can have some undesirable side effects.  It is possible for these small material fragments to "tunnel" through other multi-material Euler objects like walls or plates and emerge at the other side.