求助，shell层合单元的问题

用ansys12，用关键点连接而成的曲线z轴旋转而成的层合壳，层厚度随总厚度均匀变化，壳体总厚度随着x，y的值而变化,使用shell63建模并划分网格后输入从网上找的下面代码即可实现变厚度，但是shell63不支持层合单元，所以打算改用shell99或shell91，但是下面代码就不适用了。。。。用section下面的lay-up也总是报错，求助大神们啊
MXNODE = NDINQR(0,14)
*DIM,THICK,,MXNODE
*DO,NODE,1,MXNODE
*IF,NDINQR(NODE,1),EQ,1,THEN
THICK(node) = 5000/1600/COS(200/(NX(NODE)**2+NY(NODE)**2)**(1/2))**2
*ELSE
THICK(node) = 0
*ENDIF
*ENDDO
RTHICK,THICK(1),1,2,3,4