*parameter
transx = 1.0
transy = 2.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*cosphi - yt1*sinphi
y1 = xt1*sinphi + yt1*cosphi
x2 = xt2*cosphi - yt2*sinphi
y2 = xt2*sinphi + yt2*cosphi
x3 = xt3*cosphi - yt3*sinphi
y3 = xt3*sinphi + yt3*cosphi
x4 = xt4*cosphi - yt4*sinphi
y4 = xt4*sinphi + yt4*cosphi
*node
1001,<x1>,<y1>
1002,<x2>,<y2>
1003,<x3>,<y3>
1004,<x4>,<y4>
*Element, type=S4R,elset=shell
1,  1001,  1002,  1003,  1004
*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