P 1 nu 0.3 % Cu P 1 E 110000 P 1 Expan 1.65e-5 P 2 nu 0.33 % Fe P 2 E 209000 P 2 Expan 1.1e-5 P 1 mu E_1/(2+2*nu_1) P 1 lambda mu_1*nu_1*2/(1 -nu_1) P 2 mu E_2/(2+2*nu_2) P 2 lambda mu_2*nu_2*2/(1-nu_2) P dt .5 P 1 k 393 P 2 k 80 P 1 rhocp 396*8930 P 2 rhocp 450*7860 A displace e D_j{mu}D_j U1_i + D_j{mu}D_i U1_j +\ D_i{lambda} D_j U1_j = D_i {V201*E*Expan/(1-nu)} b 3 U1={0.,0.} A tempdiff e {rhocp}U1 - D_j {k*dt} D_j U1 = {rhocp*V202} b 3 U1=400. < rr V201=100 V202=0. > < run step step step step > < step heat stretch black show V400 V400=V400+V100 cont 202 show V200 > < stretch displace improve V200=V100 mapm V100=v100+v200 mapm popp show V100 > < heat tempdiff solve V301=V101-V302 V302=V101 V202=V202+v101 AvTemp= sum V202/Nodes show AvTemp popp > < improve V600=V100 solve V200=V100 show V100 redo redo > < redo mult V100=V600-V100 resolve show v100 V100=V100+V200 V200=V100 >