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 ! NONLINEAR TRANSIENT DYNAMIC ANALYSIS ET,1,COMBIN39,,,2 ! ELEMENT WITH DISPLACEMENT ALONG NODAL Y-AXIS

ET,2,MASS21,,,4 ! MASS WITHOUT ROTARY INERTIA R,1,0.0,0.0,.1,.204,.2,.432 ! SPRING DATA 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

N,2 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 SAVE 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 ARBITRARILY 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

PRVAR,2 ! PRINT DISPLACEMENTS *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界面重合节点 约束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 time,1 nsubst,100 lnsrch,on pred,on outres,basic,all ncnv,2 allsel,all SOLVE FINISH

有一些错误，请高手们帮我看看

*** ERROR *** SUPPRESSED MESSAGE CP = 7.375 Node 154 direction UY included in coupled sets 696 and 275.

TIME= 21:16:15