Displacement field
disp u v
coor x y
func ex ey exy
mate pe pv fx fy 2.0e5;0.3;0.0;0.0;
shap q 4
gaus q
load fx fy
$c6 fact = pe/(1.+pv)/(1.-pv)
vect ed ex ey
matr de 2 2
1 pv
pv 1
func
ex=+[u/x]
ey=+[v/y]
exy=+[u/y]+[v/x]
stif
dist =+[ed_i;ed_j]*de_i_j*fact+[exy;exy]*fact*((1.-pv)/2)
end
Stress field
coef u v
coor x y
disp sx,sy,sxy,
mate pe pv fx fy 2.0e5;0.3;0.0;0.0;
shap q 4
gaus q
mass q
load fsx fsy fsxy
$c6 fact = pe/(1.+pv)/(1.-pv)
matr de 2 2
1 pv
pv 1
vect ed ex ey
vect fs fsx fsy
stif
$cv ex=+{u/x}
$cv ey=+{v/y}
$cv fs_i=+de_i_j*ed_j*fact
$cv exy=+{u/y}+{v/x}
$c6 fsxy=exy*((1.-pv)/2)*fact
dist =+[sx;sx]*0.0
end