Model equations

mx'' + c x' + k u ( 1 - u x ) x = mg sin(angle)

u' = - dE(u,dp,T) / du

dT = T_surface / depth

T = T_surface/2

dp = f'(T)dT

x(0) = x₀,   x'(0) = v₀,   u(0) = u₀

Questions

• What are the initial conditions?
• What effects do parameters have?
• Does the model allow for (predict) fracture?
• What if several layers are coupled?