combin39粘结滑移问题 ，急

tamxsy

我在做FRP加固混凝土梁剥离破坏模拟时，采用combin39单元模拟粘结滑移，单元keyopt(3)=UZ,模拟Z方向滑移，弹簧单元是在重合的节点上建立的，把x，y向进行耦合，为什么计算的时候 老是警告约束不够呢，而且计算很慢，但是当我也耦合节点UZ方向的时候，就没有那个警告了，速度也快多了，但是我这样做好像没有起到弹簧的作用，到底是怎么回事啊啊？

我的QQ：50343975


/config,nres,50000
/config,fsplit,750
...........................设置参数
/prep7
ET,1,SOLID65            混凝土
ET,2,LINK8  
ET,3,SHELL41            CFRP
ET,4,SOLID45            垫块
et,5,combin39        弹簧单元

                  keyopt,5,1,0
                  keyopt,5,2,1
                  KEYOPT,5,3,3
                  keyopt,5,4,0
keyopt,1,1,1
keyopt,1,7,1
keyopt,3,1,2

*SET,rd0,6  
*SET,rd1,12
*SET,hntra,20.1      钢筋参数
*SET,hntrl,2.01

r,1
r,2,acos(-1)*0.25*rd0*rd0   
r,3,acos(-1)*0.25*rd1*rd1  

                   r,5,0.0005,0.2934,0.001,0.41577,0.005,0.9297,         粘结滑移关系
                    rmore,0.01,1.31478,0.015,1.61028,0.02,1.85939,
                    rmore,0.025,2.07886,0.029214,2.24725
                    
                     MPTEMP,,,,,,,,  
                     MPTEMP,1,0  
                     MPDATA,EX,5,,3240   
                     MPDATA,PRXY,5,,0.3

mp,ex,1,3e4                    
mp,prxy,1,0.2    混凝土
tb,concr,1  
tbdata,,0.5,0.9,hntrl,-1   

tb,miso,1,,23   
tbpt,,0.000005,hntra*0.0074627
tbpt,,0.00001,hntra*0.00992537
tbpt,,0.00005,hntra*0.0491293
tbpt,,0.0001,hntra* 0.0970149

tbpt,,0.0002,hntra*0.19
tbpt,,0.0004,hntra*0.36
tbpt,,0.0006,hntra*0.51
tbpt,,0.0008,hntra*0.64
tbpt,,0.0010,hntra*0.75
tbpt,,0.0012,hntra*0.84
tbpt,,0.0014,hntra*0.91
tbpt,,0.0016,hntra*0.96
tbpt,,0.0018,hntra*0.99
tbpt,,0.0020,hntra*1.00
tbpt,,0.0022,hntra*0.98
tbpt,,0.0024,hntra*0.97
tbpt,,0.0026,hntra*0.95
tbpt,,0.0028,hntra*0.93
tbpt,,0.003,hntra*0.92  
tbpt,,0.0032,hntra*0.9  
tbpt,,0.0034,hntra*0.88
tbpt,,0.0036,hntra*0.87
tbpt,,0.0038,hntra*0.85

mp,ex,2,2e5    钢筋
mp,prxy,2,0.3
TB,BKIN,2,1,2,1
TBDATA,,335,2000,,,,

mp,ex,3,2.3e5      CFRP
mp,prxy,3,0.3   
R,4,0.167, , , , ,  

                       MPTEMP,,,,,,,,  
                     MPTEMP,1,0                       弹簧
                     MPDATA,EX,5,,3240   
                     MPDATA,PRXY,5,,0.3

mp,ex,4,2e5
mp,prxy,4,0.3   
tb,bkin,4,1,2          垫块
tbdata,,335,2000
                                 
建模

blc4,,,100,200,1500   混凝土

/VIEW,1,1,1,1   
/ANG,1  

k,9,25,25,50         钢筋
k,10,25,25,1450
l,9,10      
lgen,2,13,,,50,,
lgen,2,13,,,,150,
lgen,2,14,,,,150,

lsel,s,line,,13,16,1
LPLOT  
CM,zj,LINE  
allsel,all

k,17,25,25,120       箍筋
k,18,75,25,120
k,19,25,175,120
l,17,18
l,17,19
lgen,2,17,,,,150
lgen,2,18,,,50
lsel,s,line,,17,20,1
LPLOT
*do,ii,1,7
lgen,2,17,20,1,,,ii*180
ii=ii+1
*enddo
LPLOT
cm,gj,line

wpoffs,,,75           垫块     
blc4,,,100,-20,50
vgen,2,2,,,,,1300
vgen,2,2,,,,220,450
vgen,2,4,,,,,400

WPAVE,0,0,0       CFRP
wpoffs,,,150
wprota,,90
blc4,,,100,1200
                              
allsel,all
aplot
                                
     划分单元   
/PREP7  

VSEL,s,volu,,1     混凝土
VATT,1,1,1,0
LSEL,s,line,,12
LESIZE,12,50,,,,,,,0
LSEL,s,line,,7
LESIZE,7,25,,,,,,,0
LSEL,s,line,,8
LESIZE,8,25,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,1

VSEL,s,volu,,2,5,1   垫块
VATT,4,,4,0
LSEL,s,line,,55,91,12
LSEL,a,line,,59,95,12
LESIZE,all,25,,,,,,,0
LSEL,s,line,,54,90,12
LESIZE,all,10,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,all

LSEL,s,line,,13,14,1     纵筋
lATT,2,3,2,0
LSEL,s,line,,15,16,1
lATT,2,2,2,0
cmsel,s,zj,line
LESIZE,all,50,,,,,,,0
LMESH,all

cmsel,s,gj,line       箍筋
lplot
lATT,2,2,2,0
LESIZE,all,25,,,,,,,0
LMESH,all

asel,s,area,,31         CFRP
aplot
AATT,3,4,3,0,
LSEL,s,line,,98
LESIZE,98,50,,,,,,,0
LSEL,s,line,,99
LESIZE,99,25,,,,,,,0
mshape,1,2D                    
MSHKEY,1
aMESH,31
                        
                                          
                  
allsel,all
eplot
                      nsel,s,loc,y,0               弹簧单元
                      nsel,r,loc,z,150,1350      
                      nplot
                      TYPE,5   
                      MAT,5
                      REAL,5  
                      ESYS,0   
                      SECNUM,none,   
                      TSHAP,LINE  
                      eintf,0.0001,,low,,,,,
allsel,all
eplot



                      nsel,s,loc,y,0             混凝土和CFRP界面重合节点 约束UX UY方向
                      nsel,r,loc,z,150,1350
                 
                      nplot
                      CPINTF,UX,0.0001,   
                      CPINTF,Uy,0.0001,      
                                          
ALLSEL,ALL     其他节点耦合
NSEL,U,LOC,Y,0  
NPLOT                     
CPINTF,ALL,0.0001,  

nsel,s,loc,y,0
nsel,u,loc,z,150,1350
NPLOT
CPINTF,ALL,0.0001,        

VSEL,s,volu,,1
ESLV,S  
VSEL,s,volu,,2,5,1
NSLV,S,1
nsel,r,loc,y,0,200
nsel,u,loc,z,100
nsel,u,loc,z,1400
nsel,u,loc,z,550
nsel,u,loc,z,950
CEINTF,0.25,ux,uy,uz, , , ,0,   
eplot

VSEL,s,volu,,1
ESLV,S  
cmsel,s,zj,line
ESLL,a
cmsel,s,gj,line
NSLL,S,1
nsel,u,loc,z,300
nsel,u,loc,z,1200
CEINTF,0.25,ux,uy,uz, , , ,0,                     

                                   
allsel,all
eplot


求解设置
/solu
antype,0
nsel,s,loc,z,100
nsel,a,loc,z,1400
nsel,r,loc,y,-20
d,all,uy
d,all,ux
nsel,r,loc,z,1400
d,all,uz

*do,ii,0,60
allsel,all
nsel,s,loc,z,550
nsel,a,loc,z,950
nsel,r,loc,y,220
f,all,fy,-(500+ii*200)
allsel,all
autots,on
time,ii+1
nsubst,200
lnsrch,on
pred,on
neqit,3000
outres,basic,all
ncnv,2
lswrite,ii+1
ii=ii+1
*enddo

allsel,all
EPLOT

/solu
allsel,all
LSSOLVE,1,61,1