*parameter
transx = 0.0
transy = 0.0
phid = 0
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=z1,elset=one,file=libshell3d
10001,1001,1002,1003
*substructure property,elset=one
<transx>,<transy>,0

0,0,0,  0,1,0,  0,0,1
****************************
*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