Project

Rupture plate

Rupture Plate:

Aim

  1. Become familiar with format of RADIOSS Engine and Starter files
  2. Know how to use HyperWorks Solver Run Manager to execute a RADIOSS simulation
  3. Analyze starter and engine output files
  4. Be able to post-process key results from an explicit analysis with HyperView
  5. To learn about materials law 2,law 27 & law 36 with its various failure modes.

Case 1:

Open Failure Johnson model.Run the model as it is. Rename it as Law2Eps_maxfailure.

About the model

The material undergoes isotropic elasto-plastic behavior which can be reproduced by a Johnson-Cook model with or without damage (/MAT/LAW27 and /MAT/LAW2, respectively).The tabulated material law (/MAT/LAW36) is also studied.

The model is an aluminum plate with 3mm thick.

This law represents an isotropic elasto-plastic material using the Johnson-Cook material model. This model expresses material stress as a function of strain, strain rate and temperature. A built-in failure criterion based on the maximum plastic strain is available.

The values a,b & n are calculated.

The EPSpmax are given as 0.151.So the rupture will happen at the failure strain of 0.151.

/FAIL/XFEM_JOHNS (Obsolete) is an XFEM (eXtended Finite Element Method) failure model.A failure model is created.Ifail_sh(Shell failure flag) = 2. (For each integration point, the stress tensor is set to zero, if D=∑Δεpεf≥1\">D=∑Δεpεf≥1D=∑Δεpεf≥1 ; i.e the platic strain is greater than the failure strain so the shell is cracked (not deleted), for all integration points or layers.)The constant D1,D2 & D3 which are real constants are also provided in order to calculate the εf (failure strain).

The material also fails when plastic strain becomes equal to rupture strain


So when the material is loaded, the von mises increases upto the EPS limit, then rupture takes place.

The plate initially has a rupture and the damages the element, the plate at the centre breaks into small fragments

From the graph, we get to know about the effect of energy involved.

Initially, The rigid wall is moving towards the aluminium plate,the plate deforms,a punch force which is created due to the initial velocity,so the internal energy gradually increases, upto to the rupture and then it becomes constant.Here the internal energy is the same as the external work,no hourglass energy is present.

Once the plate is ruptured,the punch force decreases become it has no resistance.The kinetic energy decreases and there are some ripples because,the it has to pierce the plate and the elements are cracked.The internal energy constant becomes constant, because of no resistance in the plate.

Case 2:

In Fail/Johnson, change Ifail_sh = 1,Dadv = 1, Ixfem = 1.Run the model as Law2epsmaxcrack

Where Dadv is the criterion for crack propagation

Shell failure flag = 1, Shell is cracked (not deleted), if D=∑Δεpεf≥1\">D=∑Δεpεf≥1D=∑Δεpεf≥1 for one integration point or layer.

XFEM formulation (Ixfem=1) is only compatible with BT Q4 (Ishell=1, 2, 3 or 4), and QEPH (Ishell=24) shell elements. If XFEM is activated (Ixfem=1), the failure criteria will lead to element cracking instead of element or layer deletion.

  • If /PROP/SHELL (TYPE1) is used, then mono-layer XFEM will be applied.

In this case, the whole element thickness is considered as a single layer. The failure criterion is calculated in each integration point, but only one single crack can appear in this element. This approach is compatible with all values of the shell flag (Ifail_sh=1 or 2). The crack direction is determined by the principal constraints in the last failed integration point.

So when the material is loaded, the von mises stress increases upto the EPS limit, then rupture takes place.

The plate initially has a rupture and the damages the element, the plate at the centre makes a cut into 4 parts,here the plate is much more stiffer than case 1,the elements hold together and just makes way for the spherical ball to pass through

From the graph, we get to know about the effect of energy involved.

Here the energy graphs when compared to case 1 has a smooth failure, which is in the sense that the plate at the centre splits into 4 equal halves,instead to small fragments, so the the crack propagation is uniform,less resistance when compared to case 1.The rupture gets delayed when compared with case 1 & no hourglass energy

Case 3:

For case 1 & 2 we used Fail/Johnson model

Using the Johnson-Cook failure model /FAIL/JOHNSON along with the Johnson-Cook material model provides more user control of the failure
Johnson-Cook failure strain definition:

Failure of the element happens when 1 integration point reaches εpmax

The previous models uses cumulative damage to compute failure.

The energy dissipation is similar to the previous cases.

Even though the failure card is deleted but still the failure occurs due to EPS_p_max

1. EPS_p_max - Plastic Strain beyond this will cause element deletion

2. EPS_t1 - Strain beyond which the line connecting t1 and t2 will be used

3. EPS_t2 - Strain beyond which the element will be separated from the connected elements

Since the model is modelled with plastic strain Epsmax values not zero, it ruptures

Case 4:

Delete the EPS_p_max and run the model as Law2

It means that failure strain is not defined,the plate once it is loaded , is subjected to elastic and plastic strain,the plate has more elastic strain than plastic strain,so the internal energy gradually increases and goes on and kinetic energy increases and becomes stable, hourglass energy is zero & the punch gradually increases and never drops because of no failure.

This is is mainly because the EPS_p_max is zero, where the failure takes place, the plate behaves much more ductile, internal energy is stored in the plate

    Case 5:

    In this case, the material is changed to Law1.Only E,Poisson ratio & Density are specified.THe material will have only elastic strain.No failure occurs.The internal energy increases gradually and the kinetic energy also increases and never drops.

    Case 6:

    Open Law27 file and change the properties to recommended values.Here the material used is law27

    /MAT/LAW27 (PLAS_BRIT) Elasto-plastic Material with Brittle Failure
    This law combines an isotropic elasto-plastic Johnson-Cook material model with an orthotropic brittle failure model. Material damage is accounted for prior to failure. Failure and damage occur only in tension. This law is applicable only for shells.Yield surface definition same as Johnson-Cook (Law2)

    Damage and rupture defined with 4 parameters for each principal direction:
    et = Strain at the beginning of tensile failure
    em = Maximum tensile strain at which the stress in the element is set to a value dependent on dmax1
    dmax = Maximum damage factor
    ef = Maximum tensile strain for element deletion

    Law 27 is used to simulate material damage following a Johnson-Cook plasticity law. Thus, a damage model is incorporated into the material law to take into account the damage evolution with stress decreasing up to element rupture.

    The damage parameters are:

    Tensile rupture strain t1: damage starts if the highest principal strain reaches this tension value.
    Maximum strain m1: the element is damaged if the highest principal strain is above the tension value. The element is not deleted.
    Maximum damage factors max: this value should be kept at its default value (0.999).
    Failure strain f1: the element is deleted if the highest principal strain reaches the tension value.

    From the energy graph we find that the graphs are similar to case 1 & 2, kinetic energy has a sharp increase , this is the brittle nature,i.e the area around the corners of the plate are not deformed much, and the plate has failed as if it is drilled, when compared to case 1 & 2.No hourglass occurs.

    Case 7:

    The material is changed to Law36

    This is a tabulated law; therefore, the true stress versus plastic strain function can be directly used. The rupture phase can be simulated by adding points to this hardening function.The point-by-point definition of the hardening curve in law 36 enables to bypass the adaptation difficulties when using the Johnson-Cook model.

    This law models an isotropic elasto-plastic material using user-defined functions for the work-hardening portion of the stress-strain curve (for example, plastic strain vs. stress) for different strain rates.

    Note: The first point of yield stress functions (plastic strain vs stress) should have a plastic strain value of zero. If the last point of the first (static) function equals 0 in stress, default value of ep max is set to the corresponding value of ep

    Isotropic elasto-plastic material
    • User-defined function for the work-hardening portion of the stress-strain curve
    • Available for brick and shell elements
    • Elastic portion of material stress-strain curve defined by Young modulus and Poisson’s ratio
    • Material plasticity curves can be given for up an arbitrary number of strain rates
    • Linear interpolation of strain-stress curve
    • For a given strain rate
    • For a given plastic strain

    From this material, we can get the true stress, we can get the actual stress strain curve.The

    From the energy curves, we find that the curves are similar to the earlier curves , the failure happens earlier when compared to the earlier cases.

    The plate is damaged at the centre,and breaks into small fragments.

    The comparison of .out file values for the 7 cases

    Material type CYCLE TIME TIME-STEP ELEMENT ERROR I-ENERGY K-ENERGY T K-ENERGY R EXT-WORK MAS.ERR
    Law2epsmaxfailure 49383 5.01 1.01E-04 SHELL 159 0.80% 2.63E+04 94.22 63.87 2.62E+04 0
    Law2epsmaxcrack 49406 5.01 1.01E-04 SHELL 159 3.70% 3.06E+04 122.9 62.23 2.96E+04 0
    Law2epsmaxnofail 49406 5.01 1.01E-04 SHELL 159 3.70% 3.06E+04 122.9 62.23 2.96E+04 0
    Law2 48736 5.01 1.05E-04 SHELL 159 3.00% 6.49E+04 2.976 2.38E-02 6.30E+04 0
    Law1 47968 5.01 1.10E-04 SHELL 159 0.30% 8.52E+05 2647 4090 8.57E+05 0
    Law27 49355 5.01 1.01E-04 SHELL 159 1.10% 3.98E+04 45.51 14.73 3.94E+04 0
    Law36 52228 5.01 9.59E-05 SHELL 159 -1.20% 4.22E+04 287.6 235.6 4.32E+04 0

    From the above table, we find that there is no mass error

    Cycle time increases for Law 36 and decreases for Law1

    Error increases for Law2epsmaxnofail

    But overall the error % is comparatively small and the simulations are ok

    Law2epsmaxfailure has maximum and law1 has minimum internal energy

    Conclusions:

    The rupture of a circular plate, due to the impact of a rigid sphere was studied and several failure models with different simulation parameters were compared. The results obtained highlight the sensitivity of the numerical models to simulate the failure.

    When compared with the graphs,all cases are almost similar

    Law2 (case 1,2,3) - Plate ruptures

    Case 1 element damages and ruptures due to failure model and epsmax value

    Case 2 Creates a crack and cuts into 4 equal halves, rupture is smooth due to failure model and epsmax value (smooth due to Ixfem ,Dadv values)

    Case 3 Deformation is more, ruptures due to epsmax value (failure card deleted )

    Case 4 Rupture almost negligible (epsmax value 0)

    Case 5 No rupture (Law 1 has only elastic strain)

    Case 6 Ruptures (Due to law27 material\'s brittle nature)

    Case 7 Ruptures (Due to law 36 , close to true stress & strain since true stress and strain values are inputted but takes more cycles i.e more time to complete simulation

    Hourglass energy is zero in all the cases, this is due to the use of recommended parameters

    Your portfolio, next

    Ship work like this
    in your own course.

    Talk to a career counsellor

    Free 20-minute call. We'll map the right course to where you want to work.

    By submitting you agree to our Terms and Privacy Policy.

    Talk to a career counsellor

    Free 20-minute call. We'll map the right course to where you want to work.

    By submitting you agree to our Terms and Privacy Policy.