ANSYS中使用COMBIN39单元的实例

ANSYS中使⽤COMBIN39单元的实例ANSYS中使⽤COMBIN39单元的实例1
finish
/clear
/units,si
/prep7
antype,static
et,1,solid65
r,1,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.97e9
MPDATA,PRXY,1,,0.27
TB,BKIN,1,1,2,1
TBTEMP,0
TBDA TA,,1.97e6,,,,,
TB,CONC,1,1,9,
TBTEMP,0
TBDA TA,,0,1,0.4e6,-1,,
TBDA TA,,,,,,,
blc4,,,300,300,1400
*do,i,1,8,1
lesize,i,50
*enddo
*do,i,9,12,1
lesize,i,100
*enddo
vmesh,1
et,2,link8 !是否应⽤PIPE 单元
r,2,0.000576
mp,ex,2,2.1e11
mp,prxy,2,0.3
type,2
k,12,150,150,1400
l,11,12
lesize,13,100
lmesh,13
et,3,combin39
keyopt,3,3,3 !定义Z⽅向弹簧
mp,ex,3,36534
r,3,0.012,340.98,0.023144013,676.822,0.041007667,1009.222
rmore,0.060508355,1329.6,0.078506,1669.288667,0.13106,2006.533333
rmore,0.200204,2332.511333,0.297857,2678.866667,0.381036,3004.466667
rmore,0.491887,3340.466667,0.634276,3667.711333,2.961887,3864.733333
rmore,2.989726,3982.733333,3.05323,4143.8,3.141894,4204.355333
rmore,3.16088,4294.266667,3.232973,4385.644667,3.244492,4473.044667 rmore,3.333723,4047.6, type,3
mat,3
real,3
e,37,736
e,86,737
*do,i,567,579,1
e,i,i+171
*enddo
save
et,4,shell63
r,4,0.015,0.015,0.015,0.015
mp,ex,4,2.5e4
mp,prxy,4,0.3
type,4
mat,4
real,4
asel,,,,3,6,1
aatt,4,4,4,
amesh,3,6,1
allsel
finish
/solu
d,all,all
allsel
nsel,s,loc,z,1400
nsel,r,loc,x,300
d,all,all
allsel
time,0
*do,i,1,11,1
time,i
f,737,fz,5000*i
cnvtol,f,,0.1,2
nropt,full
outres,all,all
nsubst,10
nlgeom,on
solve
fdele,all
*enddo
*do,i,12,25,1
f,737,fz,55000+2500*(i-11)
time,i
cnvtol,f,,0.05,2
nropt,full
outres,all,all
nsubst,10
nlgeom,on
solve
fdele,all
*enddo
ANSYS中使⽤COMBIN39单元的实例2
/PREP7
ANTYPE,TRANS ! NONLINEAR TRANSIENT DYNAMIC ANAL YSIS ET,1,COMBIN39,,,2 ! ELEMENT WITH DISPLACEMENT ALONG NODAL Y-AXIS
RMORE,.3,.708,.4,1.056,.5,1.5
RMORE,.6,2.064,.7,2.772,.8,3.648
RMORE,.9,4.716,1.0,6.0
R,2,2588E-6 ! MASS DATA
N,1
E,1,2
TYPE,2
REAL,2
E,2
D,1,ALL
D,2,UX
IC,2,UY,-1 ! INITIAL DISPLACEMENT AND VELOCITY
KBC,1 ! STEP LOADING
SA VE
FINISH
/SOLU
SOLCONTROL,0
CNVTOL,F,1,1E-4 ! FORCE CONVERGENCE CRITERIA OUTRES,NSOL,1
NSUBST,5
OUTPR,BASIC,NONE
TIME,.0002 ! TIME TO ALLOW INITIAL CHANGE IN ACCELERATION LSWRITE ! WRITE LOAD STEP FILE 1
NSUBST,40
OUTPR,BASIC,LAST
TIME,0.18 ! TIME ARBITRARIL Y SELECTED
LSWRITE ! WRITE LOAD STEP FILE 2
LSSOLVE,1,2,1 ! READ IN 2 LOAD STEPS AND SOLVE
FINISH
/POST26
TIMERANGE,.003,.18
NSOL,2,2,U,Y,2UY
PRV AR,2 ! PRINT DISPLACEMENTS
*GET,PER,V ARI,2,EXTREM,TMIN
FINISH
frp钢筋混凝⼟加固梁的粘结滑移问题
我最近在做frp钢筋混凝⼟加固梁的剥离问题，采⽤ansys建⽴模型。

solid65模拟混凝⼟，combin39单元来模拟frp与混凝⼟之间的粘结滑移，有⼀些错误，希望⾼⼿们有时间帮我看看。

多谢了
很多⼈说，采⽤ansys算粘结滑移很难算到收敛，但是不知道有没有算出来的，只要有可能，我就不放弃了。

命令流如下：
!第⼀组
/prep7
!定义参数和单元类型
!混凝⼟参数
*SET,hntra,20.5
*SET,ec,28800
*SET,possible,0.3
*SET,hntrl,2.35
!FRP参数
*SET,ff,1800 !抗拉强度
*SET,ef,220000 !弹性模量
!FRP纤维布的极限拉应变
*SET,ckf,ff/ ef
!FRP厚度
*SET,tf,0.121*2
et,1,solid65 !混凝⼟单元
KEYOPT,1,1,1 !不考虑⼤变形
KEYOPT,1,7,1 !考虑应⼒松弛，有助于计算收敛
et,2,link8 !钢筋单元
et,3,shell43 !模拟frp布
et,4,solid45 !刚性垫块
!定义实常数
r,1,50.3 !纵向受拉钢筋实常数
r,2,28.3 !纵向受压钢筋，箍筋实常数
r,3 !混凝⼟
R,4,tf,tf,tf,tf,,, !frp的实常数
RMORE,,
r,5 !垫块
ET,5,COMBIN39 !法向弹簧
KEYOPT,5,3,3
r,6,0,0,0.0005,29340,0.001,41577 !粘结滑移关系rmore,0.005,92970,0.01,131478,0.015,161028 rmore,0.02,185939,0.025,207886,0.029214,224725
!定义材料属性
mp,ex,1,ec !混凝⼟材料属性
mp,prxy,1,possible
tb,conc,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,1.94e5 !纵向受拉钢筋
tbdata,,368.3,2000
mp,ex,3,1.97e5 !横向箍筋，受压钢筋材料属性mp,prxy,3,0.3 tb,bkin,3,1,2,1
tbdata,,410.6,2000
!定义FRP
mp,ex,4,ef
mp,prxy,4,0.3
mp,ex,5,ec !垫块材料属性
mp,prxy,5,possible
!建模
blc4,,,100,200,-1100
k,9,20,40,0 !钢筋
k,10,20,40,-1100
l,9,10
lgen,2,13,,,60
lsel,s,line,,13,14
LPLOT
CM,zlj,LINE
allsel,all
lgen,2,13,,,,120
lgen,2,14,,,,120
lsel,s,line,,15,16
lplot
cm,zyj,line
allsel,all
k,17,20,40,-500 !箍筋
k,18,80,40,-500
k,19,20,160,-500
l,17,18
l,17,19
lgen,2,17,,,,120
lgen,2,18,,,60
lsel,s,line,,17,20,1
lgen,2,17,20,1,,,-ii*100
ii=ii+1
*enddo
LPLOT
cm,gj,line
WPAVE,0,0,0 !CFRP
wprota,,90 !旋转⼯作平⾯yz，90°
wpoffs,,,0.835 !偏移⼯作平⾯z cfrp⾯与混凝⼟⾯之间距离为0.835mm blc4,,,100,-1000 wpoffs,,,-0.835 !⽀座⾯与cfp下表⾯距离为0.835mm wpoffs,,-1100
blc4,,,100,100,20 !⽀座处混凝⼟垫块100*100*20 wpoffs,,600
wpoffs,,,-200
wprota,,-90
blc4,,,100,20,-100 !荷载处混凝⼟垫块是100*100*20 allsel,all
!划分单元
VSEL,s,volu,,1 !混凝⼟⽹格划分较细为25*50*25 VATT,1,3,1
LSEL,s,line,,12
LSEL,a,line,,7
LSEL,a,line,,8
LESIZE,12,50,,,,,,,0
LESIZE,7,20,,,,,,,0
LESIZE,8,20,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,1
cmsel,s,zlj,line !主la筋划分为50mm
lATT,2,1,2
LESIZE,all,50,,,,,,,0
LMESH,all
cmsel,s,zyj,line !主ya筋划分为50mm
lATT,3,2,2
LESIZE,all,50,,,,,,,0
LMESH,all
cmsel,s,gj,line !箍筋为20mm
LESIZE,all,20,,,,,,,0
LMESH,all
asel,s,area,,7 !CFRP划分单元为四边形20*50 aplot AATT,3,4,3
lsel,s,,,42
lsel,a,,,43
LESIZE,42,50,,,,,,,0
LESIZE,43,20,,,,,,,0
mshape,0,2D
MSHKEY,1
AMESH,7
allsel,all !划分⽀座垫块为25*50*20
VSEL,s,,,2 !选择体2
VATT,5,3,4
LESIZE,50,50,,,,,,,0
LESIZE,51,20,,,,,,,0
LESIZE,55,20,,,,,,,0
LESIZE,47,50,,,,,,,0
LESIZE,49,20,,,,,,,0
LESIZE,53,20,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,all
/VIEW,1,1,1,1
allsel,all !划分⽀座垫块为25*50*20
VSEL,s,,,3 !选择体3
VATT,5,3,4
LESIZE,67,50,,,,,,,0
LESIZE,62,20,,,,,,,0
LESIZE,61,20,,,,,,,0
LESIZE,66,50,,,,,,,0
LESIZE,57,20,,,,,,,0
LESIZE,58,20,,,,,,,0
VMESH,all
allsel,all
nsel,s,loc,y,-0.835,0 !弹簧单元0
nplot
TYPE,5
REAL,6
ESYS,0
eintf,0.84,,low !结点号码从最低点到最⾼点编号
allsel,all
nsel,s,loc,y,-1.67,0 !混凝⼟和CFRP界⾯重合节点约束UX UY⽅向CPINTF,UX,1.67 CPINTF,Uy,1.67
nsel,s,loc,y,-1.67,-0.835 ！cfrp内部结点耦合
CPINTF,Uz,0.835
ALLSEL,ALL !其他节点耦合
CPINTF,UX
CPINTF,Uy
CPINTF,Uz
allsel,all
!求解设置
/solu
antype,0
allsel,all
asel,s,,,9 !右边界约束
da,all,ux
da,all,uy
allsel,all
asel,s,loc,z,0 !对称⾯约束
da,all,symm
allsel,all !施加外部载荷
asel,s,,,17
sfa,all,,pres,2.456
allsel,all
autots,on
nsubst,100
lnsrch,on
pred,on
outres,basic,all
ncnv,2
allsel,all
SOLVE
FINISH
有⼀些错误，请⾼⼿们帮我看看
*** ERROR *** SUPPRESSED MESSAGE CP = 7.375 TIME= 21:16:15 Node 154 direction UY included in coupled sets 696 and 275.。