求助：出现刚体移动，但觉得约束没问题啊

qyt913

我在一个非线性材料的分析，但是施加载荷后，出现errror： The value of UX at node 34 is 3.997319433E+18.  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.                                                   



但觉得约束时没问题的! 不知道有人遇到过没，并且解决掉的呢？能否帮忙下呢？
我的qq：303133639

qyt913

命令流文件，需要的话，留下邮箱，我发给你……

ferriswanglei

为什么不直接传上来呢？

qyt913

3# ferriswanglei


/prep7  
bosongbi=0.40
et,1,solid185
et,2,shell63
MP,EX,1,1.36e5  
MP,PRXY,1,bosongbi  
MP,DENS,1,1200  
TB,HYPE,1,1,3,OGDE  
TBDATA,,-7.594e4,4.941,1.138e4,6.425,6.572e4,4.712  
TB,PRONY,1,1,2,SHEAR
TBDATA,,0.148,2.123,0.252,9.371,,   
MP,EX,2,8e4
MP,PRXY,2,bosongbi
MP,DENS,2,1200  
TB,HYPE,2,1,3,OGDE  
TBDATA,,-7.594e4,4.941,1.138e4,6.425,6.572e4,4.712  
TB,PRONY,2,1,2,SHEAR
TBDATA,,0.148,2.123,0.252,9.371,,   
MP,EX,3,3.4e4   
MP,PRXY,3,bosongbi
MP,DENS,3,1200  
TB,HYPE,3,1,3,OGDE  
TBDATA,,-4.895e4,5.511,9.889e3,6.571,3.964e4,5.262   


R,1,1e-10,1e-10,1e-10,1e-10,

BLOCK,0,0.010,0,0.010,0,0.006,  
CYLIND,0.001, ,0,0.006,0,90
vptn,all
wpoff,0,0,0.001
vsbw,all
wpoff,0,0,0.001
vsbw,all
wpcs,-1
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL  


VSEL,S,LOC,Z,0,0.001
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,4   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       1,       1,   1,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  

FLST,5,4,4,ORDE,4   
FITEM,5,1   
FITEM,5,5   
FITEM,5,19  
FITEM,5,-20
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,10, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,14  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,4, , , , ,1  
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,4   
FITEM,5,15  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30, , , , ,1
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,5,4,4,ORDE,4   
FITEM,5,23  
FITEM,5,25  
FITEM,5,32  
FITEM,5,-33
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20,0.9, , , ,1   
!*  

FLST,5,4,4,ORDE,4   
FITEM,5,2   
FITEM,5,-3  
FITEM,5,29  
FITEM,5,31  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,15, , , , ,1

FLST,2,2,4,ORDE,2   
FITEM,2,29  
FITEM,2,31  
LCCAT,P51X  
FLST,2,2,4,ORDE,2   
FITEM,2,2   
FITEM,2,-3  
LCCAT,P51X  
CM,_Y,VOLU  
VSEL, , , ,       4
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
!!!!!!!!!!!!!!!!!!!!!!!!!! 10*10*30，四边形对边网格数分别为30和20，其中20的比例为0.9


VSEL,S,LOC,Z,0.001,0.002
FLST,5,2,6,ORDE,2   
FITEM,5,2   
FITEM,5,6   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       2,       1,   1,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  

FLST,5,1,4,ORDE,1   
FITEM,5,17  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,4, , , , ,1  
!*  
FLST,5,2,6,ORDE,2   
FITEM,5,2   
FITEM,5,6   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


allsel,all

FLST,5,2,6,ORDE,2   
FITEM,5,3   
FITEM,5,5   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
VATT,       3,       1,   1,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,34  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,10, , , , ,1
!*  
FLST,5,2,6,ORDE,2   
FITEM,5,3   
FITEM,5,5   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  

ALLSEL,ALL  
/REPLOT


/SOL
ANTYPE,0
NLGEOM,1

!!!!!!!!!求解设置

asel,s,loc,x,0  

FLST,2,6,5,ORDE,6   
FITEM,2,5   
FITEM,2,-6  
FITEM,2,12  
FITEM,2,18  
FITEM,2,24  
FITEM,2,26  
DA,P51X,SYMM


ALLSEL,ALL  
asel,s,loc,y,0  

FLST,2,6,5,ORDE,6   
FITEM,2,3   
FITEM,2,11  
FITEM,2,14  
FITEM,2,19  
FITEM,2,23  
FITEM,2,27  
DA,P51X,SYMM

ALLSEL,ALL  
GPLOT   

asel,s,loc,z,0.006  
APLOT   

FLST,2,2,5,ORDE,2   
FITEM,2,10  
FITEM,2,16  

DA,P51X,ALL,0  

ALLSEL,ALL  
SAVE



FLST,2,1,5,ORDE,1   
FITEM,2,9   
SFA,P51X,1,PRES,101325

ferriswanglei

可以求解
/prep7  
bosongbi=0.40
et,1,solid185
et,2,shell63
MP,EX,1,1.36e5  
MP,PRXY,1,bosongbi  
MP,DENS,1,1200  
TB,HYPE,1,1,3,OGDE  
TBDATA,,-7.594e4,4.941,1.138e4,6.425,6.572e4,4.712
TB,PRONY,1,1,2,SHEAR
TBDATA,,0.148,2.123,0.252,9.371,,   
MP,EX,2,8e4
MP,PRXY,2,bosongbi
MP,DENS,2,1200  
TB,HYPE,2,1,3,OGDE  
TBDATA,,-7.594e4,4.941,1.138e4,6.425,6.572e4,4.712


TB,PRONY,2,1,2,SHEAR
TBDATA,,0.148,2.123,0.252,9.371,,   
MP,EX,3,3.4e4   
MP,PRXY,3,bosongbi
MP,DENS,3,1200  
TB,HYPE,3,1,3,OGDE  
TBDATA,,-4.895e4,5.511,9.889e3,6.571,3.964e4,5.262

  
R,1,1e-10,1e-10,1e-10,1e-10,
BLOCK,0,0.010,0,0.010,0,0.006,  
CYLIND,0.001, ,0,0.006,0,90
vptn,all
wpoff,0,0,0.001
vsbw,all
wpoff,0,0,0.001
vsbw,all
wpcs,-1
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL  
VSEL,S,LOC,Z,0,0.001
VATT,1,1,1,0
ALLSEL
LSEL,S, , ,1,5,4
LSEL,A,,,19,20
LESIZE,ALL, , ,10, , , , ,1
ALLSEL
LSEL,S, , ,14
LESIZE,ALL, , ,4, , , , ,1  
ALLSEL
LSEL,S, , ,4,15,11
LESIZE,ALL, , ,30, , , , ,1
ALLSEL
VSWEEP,1  
LSEL,S, , ,23,25,2
LSEL,A,,,32,33
LESIZE,ALL, , ,20,0.9, , , ,1   
ALLSEL
LSEL,S, , ,2,3
LSEL,A,,,29,31,2
LESIZE,ALL, , ,15, , , , ,1
ALLSEL
LCCAT,29,31  
LCCAT,2,3
VSWEEP,4
!!!!!!!!!!!!!!!!!!!!!!!!!! 10*10*30，四边形对边网格

数分别为30和20，其中20的比例为0.9
VSEL,S,LOC,Z,0.001,0.002
VATT,2,1,1,0
ALLSEL
LESIZE,17, , ,4, , , , ,1  
VSEL,S, , ,2,6,4
VSWEEP,ALL  
allsel,all
VSEL, , , ,3,5,2
VATT,3,1,1,0
ALLSEL
LESIZE,34, , ,10, , , , ,1
VSEL, , , ,3,5,2
VSWEEP,ALL  
/REPLOT
/SOLU
ANTYPE,0
NLGEOM,1
NSUBST,50,100,20
TIME,1
!!!!!!!!!求解设置
asel,s,loc,x,0  
DA,ALL,SYMM
ALLSEL,ALL  
asel,s,loc,y,0  
DA,ALL,SYMM
ALLSEL,ALL  
asel,s,loc,z,0.006  
APLOT   
DA,ALL,ALL,0  
ALLSEL,ALL  
SFA,9,1,PRES,101325

qyt913

5# ferriswanglei

我的为什么还是不能计算啊？依然是
The value of UZ at node 2 is 1.499044699E+12.  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.

qyt913

5# ferriswanglei
你只是把命令流简化下，但我是我求解时，依然有问题……希望您能再指导指导

ferriswanglei

我这里可以计算

qyt913

8# ferriswanglei 你的ansys是什么版本呢？那可真的很奇怪，斑竹，可以留个更直接交流的方式吗？我的qq：303133639，可以吗？我的网格，什么的都觉得没错……

qyt913

8# ferriswanglei
那我可不可以看看你的结果，比如z方向的变形云图！可以吗？