Property and Elements

ERROR ID: 28 for /PROP/TYPE16.

In Radioss versions before 2017, this error indicates that the number of layers defined for /PROP/TYPE16 was set to N>1. To fix this issue, set N=1.

Strain formulations for /PROP/CONNECT

/PROP/CONNECT is used only for /MAT/LAW59 (CONNECT) and /MAT/LAW83. This element formulation is based on the relative element elongation in the normal and shear direction. This avoids the dependence of element time step on the element height. The element height can be even equal to 0.
Note: Material stiffness is defined in the normal and shear directions. The material has no stiffness in lateral direction. Therefore, it is always advised to attach the nodes of /PROP/CONNECT elements to some Lagrange components either directly or through a tied contact.
The parameter Ismstr in /PROP/TYPE43 (CONNECT) specifies how the nodal forces are calculated from stresses. When Ismstr=1, the initial middle area of the element is used in calculations. If Ismstr=4, the current middle area of the element is used. Ismstr=1 can be used to avoid a time step drop that happens when the connection element separates from the element it is connected to which can cause the connection element to expand laterally.

Layer thickness and position calculation.

In the properties /PROP/TYPE11 and /PROP/TYPE16, the properties, options Thick, ti and Ipos will affect the layer thickness and layer position.

Ipos = 0:

The layer positions Zi are automatically calculated with regard to layer thicknesses.

If the sum of layer thickness is different from the input value Thick, a warning message is displayed. Radioss then adjusts the layer thickness and layer position.

Example 1

Shell thickness Thick (1.8) is not equal to the sum of layer thickness (0.5+0.6+0.5=1.6) in the input.

During the computation, Radioss uses the sum of input layer thickness (1.6) to adjust each layer thickness using   ti = ti _ input T h i c k t i _ i n p u t and the layer thickness are automatically updated.

The position of each layer is also adjusted based on the new layer thickness.
Figure 1.




The Adjusted layer thickness is calculated:
Adjusted Layer Thickness Layer Position
ti = ti _ input T h i c k t i _ i n p u t
t 1 = 0.5 1.8 1.6 = 0.5625 Z 1 = ( Thick 2 t 1 2 ) = ( 1.8 2 0.5625 2 ) = 0.61875
t 2 = 0.6 1.8 1.6 = 0.675 Z 1 = 0
t 3 = 0.5 1.8 1.6 = 0.5625 Z 3 = Thick 2 t 3 2 = 1.8 2 0.5625 2 = 0.61875
Figure 2.


Note: The local Z-axis defines the location of the middle of the shell. Since the Zi values are real layer positions in the local Z axis, both negative and positive values of Zi are allowed.
Figure 3.


Ipos= 1:

All layer positions Zi in the element thickness MUST be user-defined.
  • If multiple layers have the same position and thus are overlapping, a warning message is displayed. In this case, the layer thickness will also be adjusted.
Example 2

If Ipos=1 and the layer position Zi is not defined or =0.0, then the layers are all in the same position and overlap each other.

During the computation, Radioss takes the thickness of the thickest layer (0.6) to adjust each layer thickness with ti = ti _ input T h i c k T h i c k n e s s   o f   t h i c k e s t   l a y e r .

The position of each layer will not be adjusted.
Figure 4.




The Adjusted layer thickness is calculated:
Adjusted Layer Thickness Layer Position
ti = ti _ input T h i c k T h i c k n e s s   o f   t h i c k e s t   l a y e r
t 1 = 0.5 1.8 0.6 = 1.5 Z 1 = 0
t 2 = 0.6 1.8 0.6 = 1.8 Z 1 = 0
t 3 = 0.5 1.8 0.6 = 1.5 Z 1 = 0
Figure 5.


  • If multiple layers may have a small overlap or a gap, a warning message is displayed. In this case, the layer thickness and layer position will be also adjusted.
Example 3

If Ipos=1 and you have defined Zi with a small overlap, layers thickness and layer position will be adjusted.

During the computation, Radioss calculates the distance between the upper and lower surfaces of layers:

Dist = | ( Z 1 + t 1 2 ) + ( Z 3 + t 3 2 ) | = | ( 0.5 + 0.6 2 ) + ( 0.5 + 0.6 2 ) | = 1.6

And takes this distance to adjust each layer thickness with ti = ti _ input T h i c k D i s t .
Figure 6.






The Adjusted layer thickness is calculated:
Adjusted Layer Thickness Layer Position
ti = ti _ input T h i c k D i s t
t 1 = 0.6 1.8 1.6 = 0.675 Z 1 = Thick 2 t 1 2 = 1.8 2 0.675 2 = 0.5625
t 2 = 0.6 1.8 1.6 = 0.675
t 3 = 0.6 1.8 1.6 = 0.675 Z 3 = ( Thick 2 t 3 2 ) = ( 1.8 2 0.675 2 ) = 0.5625
Note: The local Z-axis defines the location of the middle of the thickness shell. Since the Zi is user-defined, the layers are stacked according to the user definition.
Figure 7.


Zero or Negative Volume messages.

The Stater stops with the message:
** ERROR:  ZERO OR NEGATIVE 3D SOLID VOLUME
ZERO OR NEGATIVE VOLUME 3D-ELEMENT NB 1
Figure 8.


When using fully-integrated solid elements, the element is decomposed into sub-volumes associated to with each integration point. If the element is badly warped, one sub-volume could be negative even though the total volume of the element is positive. To fix this issue, remesh the part or switch to an element with 1 integration point.

During a simulation the Engine output message:
WARNING ZERO OR NEGATIVE VOLUME: 3D-ELEMENT ID: 1330, TOTAL (SMALL) STRAIN OPTION 
IS USED BY ELEMENT

Which indicates that the element has a zero or negative volume which would normally cause the analysis to top. By automatically switching to a small strain formulation, the analysis can continue without stopping. If the element originally used full geometric nonlinear strain (Ismstr = 2, 4), then the switch to small strain uses the element shape from cycle before the negative volume. If the element originally used Lagrange type total strain (Ismstr = 10, 12), then switch to total strain uses the element shape at time=0.0. This automatic switching can be disabled by using the Engine option /NEGVOL/STOP.

Sometimes a simulation will stop with the warning or error message:
Zero or Negative Volume

This happens when solid elements are very deformed and their characteristic length goes to 0 and; thus, the element has zero or negative volume. You may notice that before getting this error message, the time step of the element is reducing.

In the case of large strain formulation, the time step of an element goes to 0 when the element is compressed.
Figure 9.


To solve this problem, first check that the material properties are correct in the model. Often this issue is caused by a mistake in the material property input which causes the material to be too soft.

To prevent the drop in time step and negative volume issue, use Ismstr =2 or 12 in the solid property or in the option /DEF_SOLID. Then in the Radioss Engine file (Runname_0001.rad) use the option /DT/BRICK/CST which will set the time step value at which the solid elements will switch to small strain. The time step specified is recommended to be two to four times smaller than the constant time step is specified in /DT/NODA/CST.

This means that the solid elements using Ismstr =2 will use large strain formulation while their time step remains greater than Δ T min , and will then switch to small strain formulation. Solid elements using Ismstr =12 will use the Lagrange type total strain formulation while their time step remains greater than Δ T min , and will then switch to total small strain formulation.

Their volume will then remain constant and the element can even reverse its orientation. The drop in their time step normally stops except for some materials, especially viscous materials.

Integration points through the shell thickness

If only one integration point is used, a membrane only behavior will be obtained (except with LAW1, up to V44) which can be used to model fabric that has no bending stiffness.

When a material is elastic, three integration points give an exact solution and the bending moments are exactly integrated through the thickness of the shell.

For deformation that causes plasticity, the bending moments are not integrated exactly and more integration points are needed for an accurate solution. Therefore, it is recommended to use five integration points, /PROP/SHELL, N=5.
Figure 10.


Different results for shell properties

For shell elements that use the same material and through thickness integration points, the results can be different depending on which shell property is used. The through thickness integration used by /PROP/SHELL and /PROP/SH_ORTH defines the integration points and weights to integrate exactly the bending moments in the elastic case when three integration points are used.

The integration scheme which is used by /PROP/SH_COMP and /PROP/SH_SANDW is a step-by-step integration scheme and uses integration points at the center of each layer with weights that correspond to the relative thickness of each layer. An error occurs in the elastic field when there are a few layers or large differences in the thicknesses of the layers. One way to work around this problem is to subdivide the thicker layers. However, it is generally not well-suited when modeling the failure of the layers.
Figure 11.


  • Elastic case:
  • Stress is linear through the thickness; an integration of forces step-by-step is exact, but the integration of moments, step-by-step is not exact since σ ( z ) ; where, z is quadratic.

QEPH versus Belytschko shell

QEPH shells (Ishell=24) are more accurate for elastic or elasto-plastic loads, whatever the loading type - quasi-static or dynamic. The Belytschko shells (Ishell=1-4), will often hourglass.

QEPH shells (Ishell=24) will give better results if the mesh is fine enough. When using a coarse mesh, this formulation will be too stiff and some local buckling phenomena could be missed. In the case of a coarse mesh, the Belytschko (Ishell=4) shells often give better results.

Null Diameter SPH ERROR 174

MESSAGE ID :          174
** ERROR : NULL DIAMETER FOR SPH PARTICLE ID=52032255

This message means that an SPH particle is compressed so much that its diameter is zero which causes this error.

Most of the time, this is due to an input error, such as entering information in the wrong units. It is recommended to check the consistency between mass, density, diameter and particle pitch. Also, review that the material input is correct.

To see the particle diameter, review the “Diameter” animation contour is always available for the particles.