*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*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
*substructure generate, type=z1, overwrite, 
recovery matrix=yes, mass matrix=yes,libr=lib_misc_pert
*RETAINED NODAL DOFS,SORTED=NO
1001,1,6
1002,1,6
1003,1,6
*substructure matrix output, stiffness=yes,
mass=yes, recovery=yes, file=lib_misc_pert
*END STEP