完整word版ANSYS中使用COMBIN39单元的实例.docx
- 文档编号:14983422
- 上传时间:2023-06-29
- 格式:DOCX
- 页数:14
- 大小:18.23KB
完整word版ANSYS中使用COMBIN39单元的实例.docx
《完整word版ANSYS中使用COMBIN39单元的实例.docx》由会员分享,可在线阅读,更多相关《完整word版ANSYS中使用COMBIN39单元的实例.docx(14页珍藏版)》请在冰点文库上搜索。
完整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.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整 word ANSYS 使用 COMBIN39 单元 实例