Appendix

Basic Relations

E,ν E,G E,B G,ν G,B B,ν λ,μ
E E E E 2(1+ν)G 9BG3B+G 3(12ν)B (3λ+2μ)μλ+μ
G=μ E2(1+ν) G 3EB9BE G G 3(12ν)B2(1+ν) μ
B=K E3(12ν) EG9G3E B 2(1+ν)G3(12ν) B B 3λ+2μ3
ν ν E2G2G 3BE6B ν 3B2G6B+2G ν λ2(λ+μ)
D11 E(1ν)(1+ν)(12ν) (4GE)G3GE 3B(3B+E)9BE 2G(1ν)12ν 3B+4G3 3B(1ν)1+ν λ+2μ
D12=λ Eν(1+ν)(12ν) (E2G)G3GE 3B(3BE)9BE 2Gν12ν 3B2G3 3Bν1+ν λ
C11 E1+ν2 4GG4GE 36BE36(3EB)2 2G1ν 4G(3B+G)3B+4G 3B(12ν)1ν2
C12 Eν1+ν2 (E2G)2G4GE 6E(3EB)36(3EB)2 2Gν1ν 2G(3B2G)3B+4G 3B(12ν)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=ε113(ε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)=α˙ε(t)+(1α)˙εfiltered(tdt)

Where,
˙ε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(tdt)
Filtered strain rate at the previous time step.

For materials laws where Fcut can be entered.

α=2π dt Fcut

Where, Fcut is the cutoff frequency.

Thus,

Fcut=α2π dt 

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.