RD-E: 0800 Hopkinson Bar

The high strain rate tensile behavior of the 7010 aluminum alloy is studied using the Hopkinson pressure bar technique (stress wave).

Precise data for high strain rate materials is necessary to enable the accurate modeling of high-speed impacts. The high strain rate characterization of materials is usually performed using the split Hopkinson Pressure Bar within the strain rate range 100-10000 s-1. Using the one-dimensional analysis of the Hopkinson bar experiment, it is assumed that the object deforms under uni-axial stress, the bar object interfaces remain planar at all times, and the stress equilibrium in the object is achieved using travel times. The Radioss explicit finite element code is used to investigate these assumptions.

rad_ex_8_Hopkinson_bar
Figure 1.

Options and Keywords Used

  • Axisymmetrical analysis (/ANALY) and quad elements
  • High strain rate and Split Hopkinson Pressure Bar (SHPB)
  • Wave propagation and stress pulse
  • Elastic model (/MAT/LAW1 (ELAST)) and Johnson-Cook elasto-plastic model (/MAT/LAW2 (PLAS_JOHNS))
  • Boundary conditions (/BCS)
  • Imposed velocities (/IMPVEL)

Low extremity nodes of the output bar are fixed in the Z direction. The axisymmetrical condition on the revolutionary symmetry axis requires the blocking of the Y translation and X rotation.

The projectile is modeled using a steel cylinder with a fixed velocity in the direction Z. The required strain rate is taken into account by applying two imposed velocities, 1.7 ms-1 and 5.8 ms-1 in order to produce strain rate ranges in the of 80 s-1 and 900 s-1 (low and high rates) object.

True Stress, True Strain and True Strain Rate Measurement from Time History

rad_ex_fig_8-6
Figure 2. Nodes and Quads Saved for Time History

In the experiment, the strain gauge is attached to the object. In simulation, the true strain will be determined from 9040 and 6 nodes' relative Z displacements ( l 0 = 3.83638 mm).

The true stress can be given using two data sources. The first methodology consists of using the equation previously presented, based on the assumption of the one-dimensional propagation of bar-object forces. The engineering strain ε t associated with the output stress wave is obtained from the Z displacement of nodes located on the output bar. The true plastic strain is extracted from the quads on the object, saved in the Time History file. True stress can also be measured directly from the Time History using the average of the Z stress quads 6243, 6244, 6224 and 6235. It should be noted that the section option is not an available option with the quad elements.

The strain rate can be calculated from either the true plastic strain of quads saved in /TH/QUAD or from the true strain ε true .
Table 1. Relations Used in the Analysis
  High Rate Testing
True stress σ t r u e ( t ) = S b a r E b a r S O b j e c t ε T ( t ) exp ( ε p l ( t ) ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaHdpWCdaWgaaWcbaGaamiDaiaadkhacaWG1bGaamyzaaqabaGc daqadaqaaiaadshaaiaawIcacaGLPaaacqGH9aqpdaWcaaqaaiaado fadaWgaaWcbaGaamOyaiaadggacaWGYbaabeaakiaadweadaWgaaWc baGaamOyaiaadggacaWGYbaabeaaaOqaaiaadofadaWgaaWcbaGaam 4taiaadkgacaWGQbGaamyzaiaadogacaWG0baabeaaaaGccqaH1oqz daWgaaWcbaGaamivaaqabaGcdaqadaqaaiaadshaaiaawIcacaGLPa aaciGGLbGaaiiEaiaacchadaqadaqaaiabew7aLnaaBaaaleaacaWG WbGaamiBaaqabaGcdaqadaqaaiaadshaaiaawIcacaGLPaaaaiaawI cacaGLPaaaaaa@6036@ Z stress average from quads saved in /TH
True strain ε t r u e ( t ) = ln ( l i ( t ) l 0 ) l i ( t ) = l 0 + Δ l = l 0 + ( u 9040 ( t ) u 6 ( t ) ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaH1oqzdaWgaaWcbaGaamiDaiaadkhacaWG1bGaamyzaaqabaGc daqadaqaaiaadshaaiaawIcacaGLPaaacqGH9aqpciGGSbGaaiOBam aabmaabaWaaSaaaeaacaWGSbWaaSbaaSqaaiaadMgaaeqaaOWaaeWa aeaacaWG0baacaGLOaGaayzkaaaabaGaamiBamaaBaaaleaacaaIWa aabeaaaaaakiaawIcacaGLPaaacaaMf8UaamiBamaaBaaaleaacaWG PbaabeaakmaabmaabaGaamiDaaGaayjkaiaawMcaaiabg2da9iaadY gadaWgaaWcbaGaaGimaaqabaGccqGHRaWkcaqGuoGaamiBaiabg2da 9iaadYgadaWgaaWcbaGaaGimaaqabaGccqGHRaWkdaqadaqaaiaadw hadaWgaaWcbaGaaGyoaiaaicdacaaI0aGaaGimaaqabaGcdaqadaqa aiaadshaaiaawIcacaGLPaaacqGHsislcaWG1bWaaSbaaSqaaiaaiA daaeqaaOWaaeWaaeaacaWG0baacaGLOaGaayzkaaaacaGLOaGaayzk aaaaaa@6956@
True strain rate ε ˙ = Δ ε p l ( t ) Δ t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacuaH1oqzgaGaaiabg2da9maalaaabaGaaeiLdiabew7aLnaaBaaa leaacaWGWbGaamiBaaqabaGcdaqadaqaaiaadshaaiaawIcacaGLPa aaaeaacaqGuoGaamiDaaaaaaa@450F@ ε ˙ = Δ ε t r u e ( t ) Δ t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacuaH1oqzgaGaaiabg2da9maalaaabaGaaeiLdiabew7aLnaaBaaa leaacaWG0bGaamOCaiaadwhacaWGLbaabeaakmaabmaabaGaamiDaa GaayjkaiaawMcaaaqaaiaabs5acaWG0baaaaaa@46FD@

Input Files

The input files used in this example include:
High_strain_rate
<install_directory>/hwsolvers/demos/radioss/example/08_Hopkinson_Bar/High_strain_rate/SHPB_H*
Low_strain_rate
<install_directory>/hwsolvers/demos/radioss/example/08_Hopkinson_Bar/Low_strain_rate/SHPB_L*

Model Description

In order to model and predict the behavior of material during impact, the responses at very high strain rates should be studied.

The Split Hopkinson Bar is an inexpensive device for performing high strain-rate experiments. 1

This equipment consists of four long pressure bars:
  • the striker bar
  • the incident bar
  • the transmission bar
  • the drop bar

The object is sandwiched between the transmission and the incident bar. Assuming that the wave propagation in the bar is non-dispersive, the force and displacement upon contact between the bar and the object can be obtained from the strains measured through experience. In this example, the dynamic tensile behavior, achieved through experience of the 7010 aluminum alloy with a Split Hopkinson Pressure Bar (SHPB) is compared to numerical simulations. Two cases are studied at the strain rates of 80 s-1 (low rate) and 900 s-1 (high rate), respectively. At high strain rates, experience shows that the stress flow significantly increases by more than 30% with the strain rate increasing; thus demonstrating strain rate dependence in aluminum alloys in general. For the strain rates' range applied here, an existing Johnson-Cook model is used to describe the stress flow as a strain and strain rate function. Failure is not taken into account.

The Split Hopkinson Pressure Bar technique corresponds to a high strain rate deformation of the aluminum alloy at high stress. Figure 3 shows a diagram of the basic Hopkinson bar setup. It consists of two cylindrical bars of the same diameter, respectively called Input and Output bars.

rad_ex_fig_8-1
Figure 3. Hopkinson Bar Device

The objects material undergoes an isotropic elasto-plastic behavior which can be reproduced using a Johnson-Cook model (/MAT/LAW2). The steel bars and the striker follow a linear elastic law (/MAT/LAW1).

The following system is used: mm, ms, g, N, MPa.

rad_ex_fig_8-2
Figure 4. Object Geometry and Cross-section (dimensions in mm)

Taking into account the geometry's revolution symmetry, the material and the kinematic conditions, an axisymmetrical model is used (N2D3D = 1 in /ANALY set up in the Starter file). Y is the radial direction and Z is the axis of revolution.

Johnson-Cook Model

The Johnson-Cook model describes the stress in relation to the plastic strain and the strain rate using the following equation:

rad_ex_8_Johnson-Cook
Figure 5.
Where,
ε ˙
Strain rate
ε ˙ 0
Reference strain rate
ε p
Plastic strain (true strain)
a
Yield stress
b
Hardening parameter
n
Hardening exponent
c
Strain rate coefficient

The two optional inputs, strain rate coefficient and reference strain rate, must be defined for each material in /MAT/LAW2 in order to take account of the strain rate effect on stress, that is the increase in stress when increasing the strain rate. The constants a, b and n define the shape of the strain-stress curve.

In CRAHVI, G4RD-CT-2000-00395, D.1.1.1, Material Tests - Tensile properties of Aluminum Alloys 7010T7651 and AU4G Over a Range of Strain Rates, the behavior of the 7010 aluminum alloy can be described according to the relations:
σ = ( 496 + 225 ε 0.35 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaHdpWCcqGH9aqpdaqadaqaaiaaisdacaaI5aGaaGOnaiabgUca RiaaikdacaaIYaGaaGynaiabew7aLnaaCaaaleqabaGaaGimaiaac6 cacaaIZaGaaGynaaaaaOGaayjkaiaawMcaaaaa@4749@
Strain rates below 80 s-1
σ = ( 496 + 225 ε 0.35 ) ( 1 + 0.16 ln ( ε ˙ 0.08 ) ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaHdpWCcqGH9aqpdaqadaqaaiaaisdacaaI5aGaaGOnaiabgUca RiaaikdacaaIYaGaaGynaiabew7aLnaaCaaaleqabaGaaGimaiaac6 cacaaIZaGaaGynaaaaaOGaayjkaiaawMcaamaabmaabaGaaGymaiab gUcaRiaaicdacaGGUaGaaGymaiaaiAdaciGGSbGaaiOBamaabmaaba WaaSaaaeaacuaH1oqzgaGaaaqaaiaaicdacaGGUaGaaGimaiaaiIda aaaacaGLOaGaayzkaaaacaGLOaGaayzkaaaaaa@556B@
Strain rates exceeding 80 s-1 up to 3000 s-1

rad_ex_fig_8-3
Figure 6. Yield Curve of the Johnson-Cook Model: σ = ( 496 + 225 ε 0.35 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqaHdpWCcqGH9aqpdaqadaqaaiaaisdacaaI5aGaaGOnaiabgUca RiaaikdacaaIYaGaaGynaiabew7aLnaaCaaaleqabaGaaGimaiaac6 cacaaIZaGaaGynaaaaaOGaayjkaiaawMcaaaaa@4749@
The material properties of the object are:
Material Properties
Young's modulus
73000 [ MPa ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaai Gac2eacaGGqbGaaiyyaaGaay5waiaaw2faaaaa@3BE6@
Poisson's ratio
0.33
Density
0.0028 [ g m m 3 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaadEgaaeaacaWGTbGaamyBamaaCaaaleqabaGaaG4maaaa aaaakiaawUfacaGLDbaaaaa@3BBC@
The material used for the bars and projectile is TYPE1 (linear elastic) with the following properties:
Material Properties
Young's modulus
210000 [ MPa ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaai Gac2eacaGGqbGaaiyyaaGaay5waiaaw2faaaaa@3BE6@
Poisson's ratio
0.33
Density
0.0078 [ g m m 3 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaadEgaaeaacaWGTbGaamyBamaaCaaaleqabaGaaG4maaaa aaaakiaawUfacaGLDbaaaaa@3BBC@
The geometrical characteristics of the bars and projectile are:
Bars
Length
4 m
Diameter
12 mm
Projectile
Radius
12 mm
Weight
170 g

High Strain Rate Test Method

The object is screwed in between the incident and transmission bars. A stress pulse is introduced into the input bar through impact from a steel projectile on the steel disc attached to one end of the input bar. The impact generates a tensile wave which propagates along the input bar. Part of the wave is reflected and a part is transmitted via the object's interface. The stress pulse continues through the object and into the transmitted bar. The wave reflections inside the sample enable the stress to be homogenized during the test. The strain associated with the output or transmitted stress wave is measured by the strain gauges on the output or transmitted bar. The strain gauges attached to the object gauge length, provide direct measuring of the true strain, and the true plastic strain in the object during the experiment. The transmitted elastic wave provides a direct force measurement to the bar object interfaces by way of the following relation:(1) F ( t ) = S b a r E b a r ε T ( t )
Where,
E b a r
Modulus of the output bar,
ε T
Strain associated with the output stress wave
S b a r
Cross-section of the output bar

If the two bars remain elastic and wave dispersion is ignored, then the measured stress pulses can be assumed to be the same as those acting on the object.

The engineering stress value in the object can be determined by the wave analysis, using the transmitted wave:(2) σ e n g i n e e r i n g ( t ) = F ( t ) S o b j e c t = S b a r E b a r S o b j e c t ε T ( t )
Engineering stress can also be found by averaging out the force applied by the incident that is the reflected and transmitted wave, as shown in the equation:(3) σ e n g i n e e r i n g ( t ) = S b a r E b a r 2 S o b j e c t [ ε I ( t ) + ε R ( t ) + ε T ( t ) ]
Where,
ε I and ε R
Strains associated with input stress wave
ε T
Strain associated with output stress wave
True stress in the object is computed using the following relation (refer to Example 11 - Tensile Test for further details):(4) σ t r u e = σ e n g i n e e r i n g exp ( ε t r u e )
The true strain rate is given by:(5) ε ˙ = Δ ε t r u e Δ t
True stress and true strain are evaluated up to the failure point.

rad_ex_fig_8-4
Figure 7. 1D Analysis
Interface 1
F 1 = S b a r ( σ I ( t ) + σ R ( t ) ) = S b a r E b a r ( ε I ( t ) + ε R ( t ) )
Interface 2
F 2 = S b a r σ T ( t ) = S b a r E b a r ε T ( t )
Balance in object
F 1 = F 2 ; ε I ( t ) + ε R ( t ) = ε T ( t )
Engineering stress in object
σ o b j e c t ( t ) = F 1 / S o b j e c t = F 2 / S o b j e c t

Model Method

The mesh is made of 12054 2D solid elements (quads). The quad dimension is about 2 mm.

rad_ex_fig_8-5
Figure 8. Axisymmetrical Model Mesh with Imposed Velocities on the Top of the Input Bar

Strain Rate Filtering

Because of the dynamic load, strain rates cause high frequency vibrations which are not physical. Thus, the stress-strain curve may appear noisy. The strain rate filtering option enables to dampen such oscillations by removing the high frequency vibrations in order to obtain smooth results. A cut-off frequency for strain rate filtering, Fcut = 30 kHz was used in this example. Good engineering judgment should be used to determine a good value for each simulation. Refer to Example 11 - Tensile Test for further details.

Results

The purpose of this test is to obtain the results observed in experiments with a Johnson-Cook model. The increase of stress is expected to equal approximately 30% compared to the low strain rate test.

Experimental Data

Experimental results show that the variation of the true tensile flow stress compared with the true strain is approximately equivalent to a strain rate between 80 s-1 and 100 s-1. The reference strain, ε in the Johnson-Cook model is set to 0.08 ms-1. At higher rates, the true flow stress increases significantly compared with the strain rate. The 7010 aluminum alloy exhibits an increase in the flow stress by a typical 30% at high strain rates (900 s-1 to 3000 s-1) compared to static values.

Results are given at the specific true strains of 0.02, 0.05 and 0.10. The influence of the strain rate on stress can be seen in Figure 9 1:


ex8_fig11
Figure 9. Variation of True Stress Compared with True Strain for 7010 Alloy using 2 Different Rates (experimental data)
For the test performed with a strain rate of 900 s-1, the flow stress reaches 850 [ MPa ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqGqFfpeea0xe9vq=Jb9 vqpeea0xd9q8qiYRWxGi6xij=hbba9q8aq0=yq=He9q8qiLsFr0=vr 0=vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaai Gac2eacaGGqbGaaiyyaaGaay5waiaaw2faaaaa@3BE6@ at a 0.25 strain.
Table 2. True Stress at Specific Strains using Both Strain Rates (experimental data)
  Strain rate: 80 s-1 Strain rate: 900 s-1
True strain 0.02 0.05 0.1 0.02 0.05 0.1 0.25
True stress (MPa) 550 600 610 625 775 800 850

Johnson-Cook Model

Figure 10 shows the variation of true stress in time in relation to the wave propagation along the bars. Stresses are evaluated on the input bar, the object and the output bar.

rad_ex_fig_8-8
Figure 10. Stress Measurement Localizations

rad_ex_fig_8-9
Figure 11. Stress Waves in the Input Bar, the Output Bar and the Object (imposed velocities = 5.8 ms-1 )
The stress-time curve shows the incident, reflected and transmitted signals.

rad_ex_fig_8-10
Figure 12. SHPB of the Motion in Time of the Tensile Pulse

rad_ex_fig_8-11
Figure 13. von Mises Stress Wave Propagation Along Bars (imposed velocities = 5.8 ms-1 )
The speed of wave, C MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qaaaa@36BE@ along the bars is calculated using the relation:(6) C = E ρ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qaiabg2 da9maakaaabaWaaSaaaeaacaWGfbaabaGaeqyWdihaaaWcbeaaaaa@3A79@
Where,
C MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qaaaa@36BE@
5189 ms-1
E MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qaaaa@36BE@
Young's modulus
ρ
Density of the bars

The time step element is controlled by the smallest element located in the object. It is set at 5x10-5 ms. The stress wave thus reaches the object in 0.77 ms and travels 0.26 mm along the bar for each time step. Obviously, it remains lower than the element length of the smallest dimension (0.88 mm).

An imposed velocity of 5.8 ms-1 produces a strain rate in the object of approximately 900 s-1, while a strain rate of approximately 80 s-1 is achieved using an imposed velocity of 1.7 ms-1. A simulation is performed for each velocity value. It should be noted that the study on low rates is more limited in time than on high rates due to the reflected wave generated on top of the output bar.

Figure 14 shows the true stress and true strain as a function of the strain rate.

rad_ex_fig_8-12
Figure 14. Variation of True Stress with True Strain for High and Medium Strain Rates
At a high strain rate (900/s), an increase in the flow stress is observed, being approximately 30% higher than the stress obtained for a low strain rate (80/s). The Johnson-Cook model used provides precise results compared with the experimental data.

ex_8_stressZ
Figure 15. Stress Z and Plastic Strain on Object at 0.6 ms
The true stresses determined from both methodologies are shown side-by-side. This validates the analysis based on a transmitted wave. Typical curves for a model having imposed velocities equal to 5.8 ms-1 (Figure 16).

rad_ex_8_truestress
Figure 16. True Stress Comparison in the Object

rad_ex_8_strainrate
Figure 17. True Strain Rate in the Object (using both computations)

Either data sources used to evaluate the strain rate give similar results.

The results show:
  • the strain rate effect on stress, with or without the cut-off frequency for smoothing (100 kHz);
  • the influence of the strain rate coefficient (comparison with experimental data).

rad_ex_8_truestress-mpa
Figure 18. Strain Rate Effect

rad_ex_8_truestress-mpa2
Figure 19. Influence of the Strain Rate Coefficient c

These studies are performed for the high strain rate model ( ε = 900 s-1).

Figure 20 compares the distribution of the von Mises stress on the object, with and without the strain rate filtering at time t=0.6 ms.

rad_ex_fig_8-14
Figure 20. Comparison of the Distribution of the von Mises Stress at Time t=0.6 ms

More physical flow stress distribution is obtained using filtering. Explicit is an element-by-element method, while the local treatment of temporal oscillations puts spatial oscillations into the mesh.

1 CRAHVI, G4RD-CT-2000-00395, D.1.1.1, Material Tests - Tensile properties of Aluminum Alloys 7010T7651 and AU4G Over a Range of Strain Rates.