完整word版ANSYS中使用COMBIN39单元的实例.docx
《完整word版ANSYS中使用COMBIN39单元的实例.docx》由会员分享,可在线阅读,更多相关《完整word版ANSYS中使用COMBIN39单元的实例.docx(12页珍藏版)》请在冰豆网上搜索。
完整word版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
TBDATA,,1。
97e6,,,,,
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,0,1,0。
4e6,—1,,
TBDATA,,,,,,,
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
mat,2
real,2
k,11,150,150
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
nsel,s,loc,z,1400
nsel,r,loc,x,0
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!
NONLINEARTRANSIENTDYNAMICANALYSIS
ET,1,COMBIN39,,,2!
ELEMENTWITHDISPLACEMENTALONGNODALY-AXIS
ET,2,MASS21,,,4!
MASSWITHOUTROTARYINERTIA
R,1,0。
0,0.0,。
1,.204,.2,.432!
SPRINGDATA
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!
MASSDATA
N,1
N,2
E,1,2
TYPE,2
REAL,2
E,2
D,1,ALL
D,2,UX
IC,2,UY,-1!
INITIALDISPLACEMENTANDVELOCITY
KBC,1!
STEPLOADING
SAVE
FINISH
/SOLU
SOLCONTROL,0
CNVTOL,F,1,1E—4!
FORCECONVERGENCECRITERIA
OUTRES,NSOL,1
NSUBST,5
OUTPR,BASIC,NONE
TIME,。
0002!
TIMETOALLOWINITIALCHANGEINACCELERATION
LSWRITE!
WRITELOADSTEPFILE1
NSUBST,40
OUTPR,BASIC,LAST
TIME,0.18!
TIMEARBITRARILYSELECTED
LSWRITE!
WRITELOADSTEPFILE2
LSSOLVE,1,2,1!
READIN2LOADSTEPSANDSOLVE
FINISH
/POST26
TIMERANGE,.003,。
18
NSOL,2,2,U,Y,2UY
PRVAR,2!
PRINTDISPLACEMENTS
*GET,PER,VARI,2,EXTREM,TMIN
*status,parm
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,1,0
keyopt,5,2,1
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 !
纵向受拉钢筋
mp,prxy,2,0。
3
tb,bkin,2,1,2,1
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
LPLOT
*do,ii,1,5
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
lplot
lATT,3,2,2
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
MSHAPE,0,3D
MSHKEY,1
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界面重合节点约束UXUY方向
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
time,1
nsubst,100
lnsrch,on
pred,on
outres,basic,all
ncnv,2
allsel,all
SOLVE
FINISH
有一些错误,请高手们帮我看看
***ERROR*** SUPPRESSEDMESSAGE CP= 7.375 TIME=21:
16:
15
Node154directionUYincludedincoupledsets696and275.