Materials/Failure

How to read 0001.out file? /FAIL/HASHIN + /PROP/TSH_COMP + Isolid=14

For Hashin criterion, a number of integration points present as FAILURE LAYER in Engine output file.

Example:

FAILURE LAYER # 11 ELEMENT # 1763 HASHIN MODE # 6 AT TIME #: 3.6405803218587E-0

The correspondence between number of integration point, ijk and the number of FAILURE LAYER in the output file depend on your property description. How many integration points each plane define and how many layers through thickness are defined in your property.

For example, if 2x2 integration points in the plane and 8 layers through the thickness are selected. Then you have 2x2x8=32 integration points, which will be present in the output file 32 FAILURE LAYER.
  • FAILURE LAYER # 1 ~ FAILURE LAYER # 8 means all integration points through layers in position of 1st integration point in plane
  • FAILURE LAYER # 9 ~ FAILURE LAYER # 16 means all integration points through layers in position of 2nd integration point in plane.
  • FAILURE LAYER # 17 ~ FAILURE LAYER # 24 means all integration points through layers in position of 3rd integration point in plane.
  • FAILURE LAYER # 25 ~ FAILURE LAYER # 32 means all integration points through layers in position of 4th integration point in plane.


    Figure 1.

How to get LAW24 (Concrete Material) data from real physical tests?

You get data for LAW24 from the following tests:
  • Get density by weighing a specimen
  • Get Young's modulus with a cylinder compression test
  • Poisson coefficient (usually assumed to be near 0.2)
  • Get ‘fc’ from Compression strength (using cube compression test)
With the following additional tests, material LAW24 will be more accurate (in fact, providing data fits the failure envelope):
  • Get ‘ft’ from tensile test:

    ‘ft’ is direct tensile strength. This test provides ‘ft/fc’ value for LAW24.

  • Get couple of (ft,Ht) from Splitting tensile test:

    ‘Fst’ is splitting tensile strength in Splitting tensile test. Assume that Fst=ft., then you need to model this test to fit limit strength (ft) by validating ‘Ht’ value. This test should be with very slow velocity, so use HA8 solid property for in Splitting tensile test.

  • Get ‘fb’ from Biaxial test and get (f2,s0) from confined test.

Which curve should be input in material LAW36 in order to characterize a material?

For elastic plastic laws in Radioss, one gives a true stress versus true strain curve, so you have to convert the experimental curve.

Refer to RD-E: 1100 Tensile Test on how to proceed.

For some other laws, other than the elastic plastic, such as /MAT/LAW38 (VISC_TAB), one gives directly an engineering stress versus engineering strain curve.

What do Warnings ID's 519 and 520 in the output file (Runname_0000.out) mean, relative to shell property TYPE11?

WARNING ID: 519
** WARNING IN SANDWICH SHELL INITIALIZATION
SHELL (ID=…) MASS (TYPE NUMBER …) 
SUM OF LAYER MASS DIFFER FROM TOTAL
WARNING ID: 520
** WARNING IN SANDWICH SHELL INITIALIZATION 
SHELL (ID=…) INERTIA (TYPE NUMBER …)
POTENTIAL INSTABILITY DUE TO LAYER INERTIA DISTRIBUTION

These messages are written as of Radioss V44. They concern a shell element using a /PROP/TYPE11 (SH_SANDW).

They mean that the mass (or inertia) of the shell, as it is computed from the characteristics of each layer (position, thickness, associated material), is not equal to the mass (or inertia) which is computed from the global material associated to the PART and the total thickness of the shell given in the property TYPE11.

Indeed, mass me and inertia I of the element which will be distributed to the nodes of the shell are correctly computed from the characteristics of each layer:(1)
m e = ( j = 1 , N p j t j ) S
Where,
S
Area of shell element
N
Number of layers
p j
Density of layer j
t j
Thickness of layer j
I = j = 1 , N p j t j S ( S + t j 2 12 + z j 2 )
z j
Distance to the middle plane of the shell element

The stability time step is computed from the global characteristics of the shell; it cannot be ensured if the mass and inertia computed for the shell are not close enough to the mass and inertia which is computed from the density of the global material and the total thickness given in the shell property. So these messages are written if the relative error with respect to the mass (or inertia) computed from the global characteristics is greater than 1%.

In order to ensure the stability, it is recommended to set a Young's modulus for the global material, at least equal to the maximum modulus of the materials associated to all layers.

In LAW25 there are different kinds of messages for failure in the output file (Runname_0001.out). What do they correspond to?

The possible messages in case of failure criteria for shell elements using LAW25 are:
FAILURE-1
Criteria ε m 1 for maximum tensile strain in first direction has been reached
FAILURE-2
Criteria ε m 2 for maximum tensile strain in second direction has been reached
FAILURE-P
Criteria W p max (maximum plastic work for failure) has been reached

In any case, the concerned element identifier and the number of the layer that has failed are written.

When the element is deleted (it depends on the failure of its different layers and on the flag Ioff - Total Element Failure Criteria - into the global material associated to the PART), the following message is written:
RUPTURE OF SHELL ELEMENT …

I use solid elements and ε p m a x in the corresponding material. After the criteria ε p m a x is reached, elements are not deleted; why?

Depending on the material law, the solid elements are not deleted after the criteria ε p m a x is reached (they do not appear as “deleted elements” in post-processors).

For Material Laws 2, 4 and 22, only the deviatoric part of the stress tensor is set to zero, the internal pressure of the solid is still computed.

On the other hand, for Material Laws 3, 23, 28 and 36 the solid elements are deleted when ε p m a x is reached.

LAW36 with beam /PROP/TYPE18: Why no element deletion with beam?

Beam element is compatible with LAW1 (only with /PROP/TYPE3), LAW2 and LAW36 (only with /PROP/TYPE18).

Besides, the max. strain failure criteria in LAW36 and LAW2 is not compatible with beam elements. Actually, the total strain in beam elements is not calculated, so any failure criteria based on max strain are not compatible with beams. That is why no element deletion with beam observed.

What is “global integration” in Radioss?

If number of integration points N=0 in /PROP/SHELL, then global integration approach is used.
  • Global integration (N=0) is only compatible with Material Laws 1, 2, 22, 36, 43 and 60.
  • Failure models are not available with global integration for shells.
  • Global integration is not compatible with saving stresses to state files using /STATE/DT which creates /INISHE and /INISH3N.

Up to V44, material LAW1 for shells does not use integration points but switches to the global formulation (corresponding to N=0); whatever the number of integration points N is asked for.

So up to V44, there is no way to use this material law with only one integration point and membrane only behavior.

Workaround: Use material LAW2 with integration point and a huge value for the yield stress. As of V51, material LAW1 for shells uses global formulation, except if one integration point is asked for in the property, then a membrane only behavior occurs.
Note: This can explain some differences in the results between versions up to V44 and V51.

Now, LAW1 is only available for global integration (N=0 or N > 1 in /PROP/SHELL) and membrane formulations (with NP=1 in /PROP/SHELL). Global integration option (N=0 in Radioss) directly computes the resultant stresses (F1, F2, F12, M1, M2, and M12) without considering the integration throughout thickness numerically (with classic strain, stress), so there is no value for SX_JJ, SY_JJ, SXY_JJ, SYZ_JJ, SZX_JJ with JJ=1,99 in /TH/SHELL. (Refer to How is the generalized stress tensor /ANIM/SHELL/TENS/MEMB and /ANIM/SHELL/TENS/BEND computed? in FAQs). The default value for Iplas in case of LAW2 and global integration (N=0 in shell property) is Iplas =2: radial return.

The default value for Iplas in case of LAW36 and global integration (N=0 in shell property) is Iplas =1: iterative projection.

Using Global Integration Approach

For elasto-plastic laws, this requires that yield criteria of plasticity be written with resultant stresses.

Advantage of this option is, of course, the computation cost. Which is the cost like N=3.

For drawback, it is the same behaviors which have been done throughout the thickness (especially for loading and unloading).

For simple load cases, like monotone in loading and unloading, you can still get good results with less computation times by using the global integration approach.

For complicated load cases, especially where loading and unloading happened differently throughout thickness, the precision of results might not be high enough by using the global integration approach.