/AUX15  
IGESIN,'Drawing10','iges','E:\desk\'  


/prep7
k,10000,0,0,0

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
et,1,solid45
mp,ex,1,2.06e5   !钢、不锈钢
mp,nuxy,1,0.3
mp,dens,1,7850e-9  


!GUI操作，划分10 的实体
vsel,all
vatt,1,1,1
ESIZE,10,10,
CM,_Y,VOLU  
VSEL, , , ,       1 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1 
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ET,2,MASS21
R,2,0,0,0,0,0,0,

!点划分网格
TYPE,2
REAL,2

!对10000这个点约束
KSEL,S,,,10000,,,
KMESH,ALL
ALLSEL
EPLOT
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!耦合,选取Z轴12到12的所有点与(0,0,12)点耦合
NSEL,S,LOC,y,0,1
NPLOT
CERIG,node(0,0,0),ALL,ALL,,,,
ALLSEL


/solu
allsel !全选

*do,i,161,191,15
d,i,ux,,,,,uy,uz
*enddo
*do,i,551,581,15
d,i,ux,,,,,uy,uz
*enddo
*do,i,941,971,15
d,i,ux,,,,,uy,uz
*enddo
*do,i,1331,1361,15
d,i,ux,,,,,uy,uz
*enddo

f,node(0,0,0),my,4842000
f,node(0,0,0),mx,3891000

f,node(0,0,0),Fz,16884
f,node(0,0,0),Fy,-4783
f,node(0,0,0),Fx,7449

allsel !全选
SOLVE