求助，预应力碳纤维布加固混凝土梁不收敛

lfj_0129

预应力碳纤维布加固混凝土，怎么就不收敛呢，请教！！！命令流如下：

/FILNAME,PFRP ou,0     
/TITLE,PCFRP BEAM
/com,structural
/prep7
et,1,link8    !定义钢筋单元
et,2,solid65  !定义混凝土单元
et,3,plane42  !  ?
et,4,shell41  !定义FRP单元
et,5,solid45  !定义垫块单元
KEYOPT,4,1,2  !
KEYOPT,2,7,1  !
r,1,33.18   !定义箍筋的面积（直径6.5 mm）
r,2,153.94  !定义受拉(压)筋面积（直径14mm）
r,3         !定义第3种实常数：混凝土
r,4,0.167   !定义CFRP的厚度（mm）
r,5         !定义第5种实常数：刚支座
mp,ex,1,2e5       !定义箍筋（和受压钢筋）的弹性模量
mp,prxy,1,0.25    !定义箍筋（和受压钢筋）的泊松比
mp,dens,1,7800e-9 !定义箍筋（和受压钢筋）的质量密度
tb,bkin,1,1,2,1   !定义非线性材料特性表, Bkin双线性随动强化:定义箍筋（和受压钢筋）的应力—应变曲线（采用经典双线性随动强化模型BKIN）
tbdata,,210       !单位N/mm2（和tbmodif命令有什么区别？），给当前数据表定义数据（配合tbtemp,及tb使用）
/axlab,x,strain   !
/axlab,y,stress
tbplot,Bkin,1,,,yes  !绘制应力—应变曲线
mp,ex,2,2e5        !定义受拉筋的弹性模量
mp,prxy,2,0.3      !定义受拉筋的泊松比
mp,dens,2,7800e-9  !定义受拉筋的质量密度
tb,bkin,2,1,2,1    !定义受拉筋的应力—应变曲线（采用经典双线性随动强化模型BKIN）
tbdata,,300        !定义受拉筋的极限强度
/axlab,x,strain
/axlab,y,stress
tbplot,Bkin,2,,,yes
mp,ex,3,3.3e4        !定义混凝土的弹性模量N/mm2(C45是3.35e4)
mp,prxy,3,0.2        !定义混凝土的泊松比
mp,dens,3,2500e-9    !定义混凝土的质量密度（单位Kg／mm3）
tb,conc,3,1,9
tbdata,,0.7,1.0,2.77,-1   !定义混凝土强度准则,张口裂缝的剪切传递系数0.7(祝0.3),闭合裂缝的剪切传递系数1.0,单轴抗拉强度2.77,-1表示该混凝土材料不具有抗压能力(是关闭压碎选项?)
!混凝土第2个Data，定义其本构关系,多线性随动强化(MISO)
tb,miso,3,,20      !
tbpt,,0.0001,3.3   !定义混凝土的本构关系，应力单位是Mpa
tbpt,,0.0003,8.523
tbpt,,0.0005,13.45
tbpt,,0.0007,17.83
tbpt,,0.0009,21.545
tbpt,,0.001,23.138
tbpt,,0.0012,25.793
tbpt,,0.0014,27.751
tbpt,,0.0016,29.057
tbpt,,0.0018,29.779
tbpt,,0.002,30   !即fc 为30 N/mm2，fcu,k达到65 N/mm2，E应该是3.65以上？ tbpt,,0.0022,29.809
tbpt,,0.0024,29.294
tbpt,,0.0026,28.534
tbpt,,0.0028,27.601
tbpt,,0.003,26.551
tbpt,,0.0032,25.434
tbpt,,0.0034,24.285
tbpt,,0.0036,23.13
tbpt,,0.0038,22
/axlab,x,strain
/axlab,y,stress
tbplot,miso,3,,,yes     !混凝土材料的编号为3
tb,concr,3       !定义混凝土的特征值
tbmodif,2,3,0.4
tbmodif,3,3,0.6
tbmodif,4,3,2.2         !35号混凝土开裂应力？
tbmodif,5,3,-1
!CFRP材料参数定义：
MP,EX,4,2.37e5
      MP,EY,4,14816  !根据拉伸弹性模量2.37×105Mpa推得
       MP,EZ,4,14816
       MP,PRXY,4,0.17   !根据其它资料得到CFRP泊松比
       MP,PRYZ,4,0.17   
       MP,PRXZ,4,0.17   
       MP,GXY,4,7408  !根据拉伸弹性模量2.37×105Mpa推得
       MP,GYZ,4,7408
       MP,GXZ,4,7408
mp,dens,4,1880e-9    !定义CFRP的质量密度（单位Kg／mm3）
mp,alpx,4,-7.0e-7      !FRP线膨胀系数（为负，后面采用升温法施加预应力）,预应力施加后，Read Results，向上的位移大概0.031mm？
tbtemp,0             !为Data Table定义温度
!强度值在破坏准则中定义，具体为：
        tb,fail,4,            !4为材料号
       tbtemp,,crit
        tbdata,2,1         !1为最大应力破坏准则
        tbtemp,
        tbdata,10,4330,-2172,57.6,-316.8,57.6,- 316.8  !（纵向拉伸强度、纵向压缩强度、横向拉伸强度、横向压缩强度、Z向—厚度方向拉伸强度、Z向压缩强度）  !根据拉伸拉伸强度4.33×103Mpa推得
       tbdata,16,48,48,48   !（XY、YZ、XZ向剪切强度）   !根据拉伸拉伸强度4.33×103Mpa推得
/axlab,x,strain  
/axlab,y,stress  
!定义刚支座垫块材料
mp,ex,5,2e5
mp,prxy,5,0.3
/pnum,node,1
/pnum,elem,1

n,1                          !定义节点1
n,6,150                      !定义节点9,梁宽150
fill,1,6                     !利用节点填充命令生成中间7个节点
ngen,9,6,1,6,1,,31.25         !节点Y向复制,梁高250mm
ngen,29,1000,1,54,1,,,-100   !节点Z向复制,梁长2800mm
!另外产生左支座节点,后面没用？
NGEN,2,100000,1,9,1,,-25
NGEN,2,200000,1001,1009,1,,-25
NGEN,2,300000,2001,2009,1,,-25
!另外产生右支座节点
NGEN,2,100000,24001,24009,1,,-25
NGEN,2,100000,25001,25009,1,,-25
NGEN,2,100000,26001,26009,1,,-25
/view,1,1,1,1   
nplot         

type,1  !定义单元之前先声明即将定义单元的类型为Link8  
real,1  !箍筋（直径6.5 mm）
mat,1   !钢筋
*do,i,8,10,1 !在节点11和16间依次建立箍筋单元，
e,i,i+1
*enddo  
*do,i,44,46,1  !在节点83和88间依次建立单元，
e,i,i+1
*enddo  
*do,i,8,38,6   !在节点11和74间依次建立单元,
e,i,i+6
*enddo  
*do,i,11,41,6   !在节点17和80间依次建立单元,
e,i,i+6
*enddo  
egen,29,1000,1,28,1 !
!受压钢筋单元建立
type,1  !定义单元之前先声明即将定义单元的类型为Link8  
real,2  !受压钢筋（直径14mm）
mat,1   !钢筋
*do,i,44,27044,1000  !
e,i,i+1000  
*enddo  
*do,i,47,27047,1000
e,i,i+1000  
*enddo  
!受拉钢筋单元建立
type,1  !Link8
real,2  !主筋（直径14mm）
mat,2   !钢筋
*do,i,8,27008,1000
e,i,i+1000  
*enddo  
*do,i,11,27011,1000
e,i,i+1000  
*enddo  
type,4  !Shell41
real,4  !CFRP
mat,4   !CFRP
*do,i,2001,2005,1      !
*do,j,i,23000+i,1000   !
e,j,j+1,j+1000         !
e,j+1,j+1000,j+1000+1  !
*enddo  
*enddo  
/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1   
eplot  

k,1
k,2,150
k,3,150,250
k,4,,250
a,1,2,3,4            ！在上述关键点上生成一个面
lsel,s,loc,y,0       !定义新的选择集（选择沿y轴上y＝0的面上的线）
lsel,a,loc,y,250     !定义新的选择集（选择沿y轴上y＝250的面上的线）
lesize,all,,,5       !设定直线网格划分的等分数为
lsel,all
lsel,s,loc,x,0       !定义新的选择集（选择沿x轴上x＝0的面上的线）
lsel,a,loc,x,150     !定义新的选择集（选择沿x轴上x＝150的面上的线）
lesize,all,,,8      !设定直线网格划分的等分数为
type,3               !？
amesh,all            !对前面生成的面划分网格
type,2  !Solid65
real,3  !混凝土
mat,3   !混凝土
extopt,esize,28  !?
extopt,aclear,1  !?
vext,all,,,,,-2800

!左刚支座单元的建立
k,9
k,10,150
k,11,150,-31.25
k,12,,-31.25
a,9,10,11,12       !在上述关键点上生成一个面
lsel,s,loc,y,0       !定义新的选择集（选择沿y轴上y＝0的面上的线）
lsel,a,loc,y,-31.25     !定义新的选择集
lesize,all,,,5       !设定直线网格划分的等分数为
lsel,all
lsel,s,loc,x,0       !定义新的选择集（选择沿x轴上x＝0的面上的线）
lsel,a,loc,x,150     !定义新的选择集（选择沿x轴上x＝150的面上的线）
lesize,all,,,1       !设定直线网格划分的等分数为1
type,3               
amesh,7              !对7号面划分网格（前面拉伸后已经有6个面了，故从7开始编号）
type,5  !Solid45
real,5  !刚支座
mat,5   !刚支座
extopt,esize,2   
extopt,aclear,1  
vext,7,,,,,-200  !支座长是200  !定义拉伸的面,应该是前面划分网格的7号面
!右刚支座单元的建立
k,17,,,-2800
k,18,150,,-2800
k,19,150,-31.25,-2800
k,20,,-31.25,-2800
a,17,18,19,20    !在上述关键点上生成一个面
lsel,s,loc,y,0       !定义新的选择集（选择沿y轴上y＝0的面上的线）
lsel,a,loc,y,-31.25     !定义新的选择集（）
lesize,all,,,5       !设定直线网格划分的等分数为
lsel,all
lsel,s,loc,x,0       !定义新的选择集（选择沿x轴上x＝0的面上的线）
lsel,a,loc,x,150     !定义新的选择集（选择沿x轴上x＝150的面上的线）
lesize,all,,,1       !设定直线网格划分的等分数为1
type,3               
amesh,13          !对13号面划分网格（前面拉伸后已经有12个面了，故从13开始编号）
type,5  !Solid45
real,5  !刚支座
mat,5   !刚支座
extopt,esize,2   
extopt,aclear,1  
vext,13,,,,,200  !支座长是200  !定义拉伸的面,应该是前面划分网格的13号面

!左垫块单元的建立
k,25,0,250,-1000
k,26,150,250,-1000
k,27,150,281.25,-1000
k,28,0,281.25,-1000
a,25,26,27,28    !在上述关键点上生成面19
lsel,s,loc,y,250     !定义新的选择集（选择沿y轴上y＝250的面上的线）
lsel,a,loc,y,281.25     !定义新的选择集（选择沿y轴上y＝的面上的线）
lesize,all,,,5       !设定直线网格划分的等分数为
lsel,all
lsel,s,loc,x,0       !定义新的选择集（选择沿x轴上x＝0的面上的线）
lsel,a,loc,x,150     !定义新的选择集（选择沿x轴上x＝150的面上的线）
lesize,all,,,1       !设定直线网格划分的等分数为1
type,3               
amesh,19             !对19号面划分网格
type,5  !Solid45
real,5  !刚支座
mat,5   !刚支座
extopt,esize,2   
extopt,aclear,1  
vext,19,,,,,-200

!右垫块单元的建立
k,33,0,250,-1800
k,34,150,250,-1800
k,35,150,281.25,-1800
k,36,0,281.25,-1800
a,33,34,35,36    !在上述关键点上生成面25
lsel,s,loc,y,250     !定义新的选择集（选择沿y轴上y＝250的面上的线）
lsel,a,loc,y,281.25     !定义新的选择集（选择沿y轴上y＝的面上的线）
lesize,all,,,5       !设定直线网格划分的等分数为
lsel,all
lsel,s,loc,x,0       !定义新的选择集（选择沿x轴上x＝0的面上的线）
lsel,a,loc,x,150     !定义新的选择集（选择沿x轴上x＝150的面上的线）
lesize,all,,,1       !设定直线网格划分的等分数为1
type,3               
amesh,25             !对25号面划分网格
type,5  !Solid45
real,5  !刚支座
mat,5   !刚支座
extopt,esize,2   
extopt,aclear,1  
vext,25,,,,,200  !支座长是200  !定义拉伸的面25

/pnum,mat,1
/pnum,node,0
/pnum,elem,0
/view,1,1,1,1   
eplot   
allsel  
!nummrg,all             !用all或node均是相同的警告，对计算结果影响不大
nummrg,node, , , ,low   !合并FRP单元和混凝土单元节点，支座和CFRP的怎么合并\或者该处不设CFRP？
Nummrg,KP     !如果多次发生合并命令，一定要先合并节点，再合并关键点。合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题
!CPINTF,ALL,0.0001   !(耦合),但生成许多约束符号，后面有警告？
numcmp,all   !压缩相同编号的节点，复杂计算可以提高计算速度
eplot   
FINISH  

/SOL
nsel,s,loc,y,-31.25  
nsel,r,loc,z,-100  
d,all,uy
d,all,uz
d,all,ux                !
nsel,s,loc,y,-31.25  
nsel,r,loc,z,-2700  
d,all,uy
d,all,uz
d,all,ux
allsel  
gplot      !绘制所有的实体
nsel,s,node,,1669,1669,1   !
nsel,a,node,,1670,1670,1   !a-再选择某些节点，加入Active节点中
f,all,fy,-27000               !集中力，单位是N
nsel,s,node,,1687,1687,1
nsel,a,node,,1688,1688,1
f,all,fy,-27000
allsel   save
esel,s,type,,4     !用此命令为什么需要单独一行一行的输入或者执行两次bfe才能执行
eplot  
bfe,all,temp,,2610
bfe,all,temp,,2610
nlgeom,on            !打开大位移选项
nropt,full,          !指定完全的牛顿-拉普森选项求解,祝202（对结果没有影响），郝239
eqslv,spar, ,0,      !选用稀疏矩阵求解器SPARSE, 适用于实对称和非对称的矩阵， antype,0             !求解类型为静态分析
nsubst,200           !指定此荷载步的子步数200，
OUTRES,ERASE
outres,all,all       !输出全部子步的结果
neqit,50           !在非线性分析中指定平衡迭代的最大次数（越大越容易收敛）
pred,on              !在非线性分析中打开预测选项
ncnv,1               !终止分析选项(0 如果求解不收敛，也不终止分析, 1 如果求解不收敛，终止分析和程序（缺省）,2如果求解不收敛，终止分析,但不终止程序)
cnvtol,U,,0.5,0,,!位移收敛控制(则应再用力收敛)，容易收敛，0无穷范数控制收敛？MINREF值1.0）
cnvtol,f,,0.02,2,1       !设置迭代收敛准则（一般用力收敛控制即可，Value默认，Toler收敛标准为0.02, Specifies norm selection－NORM（残余力的第2范数控制收敛？），MINREF默认1）
allsel
acel,,9.8            !定义加速度
AUTOTS,1   !
LNSRCH,1   !
TIME,2     !
solve

hqx9912

我也是用的降温法也不收敛，把预应力去掉就可以收敛了，是不是降温法的问题？新手，见笑了。

sqh031566

楼主能直接把模型弄上来不？  一点都不直观……