*parameter
transx = 0.0
transy = 0.0
phid = 90
phi = pi*phid/180
sinphi = sin(phi)
cosphi = cos(phi)
dlen = 1.5
xt1 = transx
yt1 = transy
xt2 = xt1 + dlen 
yt2 = yt1
xt3 = xt1 + dlen
yt3 = yt1 + dlen
xt4 = xt1 
yt4 = yt1 + dlen
**
x1 = -xt1
y1 =  yt1
x2 = -xt2
y2 =  yt2
x3 = -xt3
y3 =  yt3
x4 = -xt4
y4 =  yt4
*node
1001,<x1>,<y1>
1002,<x2>,<y2>
1003,<x3>,<y3>
1004,<x4>,<y4>
*Element, type=S4R,elset=shell
1,  1001,  1004,  1003,  1002
*shell section, elset=shell,material=elas
0.5,5
*Material, name=Elas
*Density
40
*Elastic
1e+08,0.
*step,nlgeom
*static
0.1,1.0
*controls,param=field
,,,1.0
*boundary
1001,1,3
1001,4,5
1001,6,6,0.2
*output,field
*node output
u,
*end step