Appendix
Basic Relations
E,ν | E,G | E,B | G,ν | G,B | B,ν | λ,μ | |
---|---|---|---|---|---|---|---|
E | E | E | E | 2(1+ν)G | 9BG3B+G | 3(1−2ν)B | (3λ+2μ)μλ+μ |
G=μ | E2(1+ν) | G | 3EB9B−E | G | G | 3(1−2ν)B2(1+ν) | μ |
B=K | E3(1−2ν) | EG9G−3E | B | 2(1+ν)G3(1−2ν) | B | B | 3λ+2μ3 |
ν | ν | E−2G2G | 3B−E6B | ν | 3B−2G6B+2G | ν | λ2(λ+μ) |
D11 | E(1−ν)(1+ν)(1−2ν) | (4G−E)G3G−E | 3B(3B+E)9B−E | 2G(1−ν)1−2ν | 3B+4G3 | 3B(1−ν)1+ν | λ+2μ |
D12=λ | Eν(1+ν)(1−2ν) | (E−2G)G3G−E | 3B(3B−E)9B−E | 2Gν1−2ν | 3B−2G3 | 3Bν1+ν | λ |
C11 | E1+ν2 | 4GG4G−E | 36BE36−(3−EB)2 | 2G1−ν | 4G(3B+G)3B+4G | 3B(1−2ν)1−ν2 | |
C12 | Eν1+ν2 | (E−2G)2G4G−E | 6E(3−EB)36−(3−EB)2 | 2Gν1−ν | 2G(3B−2G)3B+4G | 3B(1−2ν)1−ν2 |
Hook Law 3D (principal stress and strain)
σ=Dε
σ1=D11ε1+D12ε2+D13ε3
σ1=(λ+2μ)ε1+λ(ε2+ε3)
σ1=λ(ε1+ε2+ε3)+2με1
σ1=Kεkk+2μe1 with εkk=ε1+ε2+ε3 and e1=ε1−13(ε1+ε2+ε3)
Hook Law 2D (plane stress)
σ=Cε
σ1=C11ε1+C12ε2
Unit Systems
Length | Time | Mass | Force | Pressure | Velocity | ρ | Energy | G |
---|---|---|---|---|---|---|---|---|
m | s | Kg | Kg m/s2 | N/m2 | m/s | Kg/m3 | Kmg2/s2 | 9.81 |
m | s | Kg | N | Pa | m/s | m Kg/l | J | 9.81 |
m | s | g | mN | mPa | m/s |
μ
Kg/l |
mJ | 9.81 |
m | s | Mg (ton) | KN | KPa | m/s | Kg/l | KJ | 9.81 |
m | ms | Kg | MN | MPa | Km/s | m Kg/l | MJ | 9.81e-6 |
m | ms | g | KN | KPa | Km/s |
μ
Kg/l |
KJ | 9.81e-6 |
m | ms | Mg (ton) |
GN | GPa | Km/s | Kg/l | GJ | 9.81e-6 |
mm | s | Kg | mN | KPa | mm/s | M Kg/l | mJ | 9.81e+3 |
mm | s | g | mN | Pa | mm/s | K Kg/l | nJ | 9.81e+3 |
mm | s | Mg (ton) |
N | MPa | mm/s | G Kg/l | mJ | 9.81e+3 |
mm | ms | Kg | KN | GPa | m/s | M Kg/l | J | 9.81e-3 |
mm | ms | g | N | MPa | m/s | K Kg/l | mJ | 9.81e-3 |
mm | ms | Mg (ton) |
MN | TPa | m/s | G Kg/l | KJ | 9.81e-3 |
cm | ms | g | daN | 105Pa bar |
dam/s | Kg/l | dJ | 9.81e-4 |
cm | ms | Kg | 104
N (KdaN) |
108Pa (Kbar) |
dam/s | K Kg/l | hJ | 9.81e-4 |
cm | ms | Mg (ton) |
107 N (MdaN) |
1011 Pa (Mbar) |
dam/s | M Kg/l | 105 J | 9.81e-10 |
cm | μs | g | 107
N (MdaN) |
1011
Pa (Mbar) |
104 m/s | Kg/l | 105 J | 9.81e-10 |
Filtering
Often it is useful to filter results in a material or failure law to remove numerical noise. The most common filter is an exponential moving average filter. This is especially important for material models that include strain rate effects.
In most materials, the flag Fsmooth = 1 must be defined to enable the filtering and the cutoff frequency entered using Fcut. For the case of filtering strain rates, use:
- ˙εfiltered(t)
- Filtered strain rate.
- ˙ε(t)
- Strain rate at the current timestep before filtering.
- α
- Degree of weighting decrease, a constant smoothing factor between 0 and 1. A higher value discounts previous values faster which results in less filtering.
- dt
- Timestep of the simulation.
- ˙εfiltered(t−dt)
- Filtered strain rate at the previous time step.
For materials laws where Fcut can be entered.
Where, Fcut is the cutoff frequency.
Thus,
The cutoff frequency is a function of the model timestep. Experience shows that the speed of the deformation is important also. For slower speeds, like a car crash, 1 – 10 kHz (1000 – 10,000 Hz) is a good value, but for high-speed events, like ballistic, less filtering should be used - so 1 – 10 GHz is appropriate. Good engineering judgment should be used to determine a reasonable value for each simulation. Refer to RD-E: 1102 Strain Rate Effect for an example of strain rate filtering.