如何用Shell99输入本构矩阵

yuanwucai

我在用Shell99自定义输入本构矩阵，计算出错如下：
*** WARNING ***                         CP =      43.602   TIME= 09:15:08
Input matrix of element 1 has terms on its main diagonal whose maximum  
ratio exceeds 1.0e9.  Accuracy may be poor.                             

*** ERROR ***                           CP =      43.899   TIME= 09:15:11
The value of UZ at node 11 is 1.343948114E+11.  It is greater than the  
current limit of 1000000.  This generally indicates rigid body motion   
as a result of an unconstrained model.  Verify that your model is      
properly constrained.                                                   

看起来是本构矩阵输入的时候有问题。但是我检查了一下，没发现什么问题，请各位大神帮我找找。
另外我发现，如果载荷是面内载荷，结果正确；但是面外载荷就会出错，如FZ。

完整代码如下：

FINISH  
/CLEAR,START

A11=1.0989
A22=1.0989
A12=0.3297
A66=0.3846       
D11=0.0916
D22=0.0916
D12=0.0275       
D66=0.0321                               

a = 2.7
b = 1.56

/PREP7
ET,1,SHELL99
KEYOPT,1,2,4
R,1,A11,A12,0,A22,0,A66,  
RMORE,0,0,0,0,0,0,  
RMORE,D11,D12,0,D22,0,D66,
RMORE,,,,,,,
RMORE, ,1, , , , ,  
!*  
BLC4,-a/2,-b/2,a,b
ESIZE,b/2,0,  
AMESH,all   
!*  
FINISH  
/SOLU   
ANTYPE,STATIC
nsel,s,loc,x,-a/2
D,all, , , , , ,UX,UY,UZ
nsel,all
*get,nnodez,node,0,count,
F,all,FZ,-1/nnodez
ALLSEL
SOLVE   
FINISH