*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=z1,elset=one,file=libshell3d
10001,1001,1002,1003
*substructure property,elset=one
<transx>,<transy>,0
0,0,0,  0,0,1,  <phid>
****************************
*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,
*output,history,freq=0
*substructure path,enter=10001
*output,field
*node output
u,
*substructure path,leave
*end step