TITLE 'Thaw Consolidation' COORDINATES cartesian1 SELECT {ngrid=10} !errlim=1e-20 {regrid=on} VARIABLES u(100) { pore water pressure } vx xm = move(x) DEFINITIONS alpha=0.15 Po=100 { Overburden pressure - kpa=kg/m*s^2} !gamma=1300 { Submerged unit weight - kg/m^3} cv=0.05 { Coefficien of consolidation - m^2/j} L=0.001 { Soil column depth (total) - m} endtime=50 lmove = L+ alpha*t^0.5 !THIS IS WHERE THE PROBLEM SEEMS TO BE !u0=Po*x/L INITIAL VALUES xm=L vx=0 EQUATIONS u: cv*dxx(u)=dt(u) vx: div(grad(vx))=0 xm: dt(xm) = vx ! CONSTRAINTS { Integral constraints } BOUNDARIES { The domain definition } REGION 1 { For each material region } START(0) point value(u)=0 point value(vx)= 0 point value(xm)= 0 LINE TO (L) point value(u)=-cv*dx(u)/dt(lmove)+Po point value(vx)=dt(lmove) point value(xm)= lmove TIME 0 BY 1e-10 TO endtime { if time dependent } !MONITORS { show progress } PLOTS { save result displays } FOR TIME 0 BY 0.01 TO endtime {ELEVATION(u) FROM (0) to (D)} ELEVATION(u) FROM (0) to (lmove) !ELEVATION(u/Po) FROM (0) to (lmove) fixed range (0,1) !ELEVATION(dx(u)) FROM (0) to (lmove) !HISTORY(u) at (0) at (3*lmove/4) at (lmove/2) at (lmove/4) at (lmove) !HISTORY(dx(u)) at (0) (3*lmove/4) (lmove/2) (lmove/4) (lmove) fixed range (0,-1000) HISTORY (lmove) END