Temperature field
disp u
coor x,y
shap %1 %2
gaus %3
mass %1
coef vn,un
mate ea eb gama 0.3 0.2 16

stif
dist=+[u/x;u/x]+[u/y;u/y]

load=+[u]*rc(un,vn,gama)*ea

end

fort
real*8 function rc(un,vn,gama)
implicit real*8 (a-h,o-z)
rc=(1-vn)*exp(gama*(1-1/un))
return
end

Concentration field of chemical constituents
disp v
coor x,y
shap %1 %2
gaus %3
mass %1
coef un,vn
mate ea eb gama 0.3 0.2 16

stif
dist=+[v/x;v/x]+[v/y;v/y]

load=+[v]*rc(un,vn,gama)*eb

end

fort
real*8 function rc(un,vn,gama)
implicit real*8 (a-h,o-z)
rc=(1-vn)*exp(gama*(1-1/un))
return
end

【Close】