% distort.prb P nu 0.0001 D 1 all D 2 all D 3 all D 4 all A stress e - D_j{nu}D_j U1_i - D_j{nu}D_iU1_j b 1 V301 = -X2 + X1*(1-X1) b 2 V301 = -X1 + 1 + X2*(1-X2) b 3 V301 = X2 - 1 + X1*(1-X1) b 4 V301 = X1 + X2*(1-X2) b all V200 = normal b all V400 = [grad] {V301} b all V302 = V400@V200 b all V500 = V302, V302 b all V500 = V500 * V200 b all V302 = V301 b all {V200*V300}_i+{V500 t V200}_ijU1_j < run mapmesh extent = max V101 V101 = V101/extent mapmesh step step step > < step stress show L1 V100 solve mapm 200 V200=V200+V100 mapm 200 clear prim popp >