ANSYS中使用COMBIN39单元的实例11

第页ANSYS中使用COMBIN39单元的实例1转贴个例子吧finish/clear/units,si/prep7antype,staticet,1,solid65r,1,MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,1.97e9MPDATA,PRXY,1,,0.27TB,BKIN,1,1,2,1TBTEMP,0TBDATA,,1.97e6,,,,,TB,CONC,1,1,9,TBTEMP,0TBDATA,,0,1,0.4e6,-1,,TBDATA,,,,,,,blc4,,,300,300,1400*do,i,1,8,1lesize,i,50*enddo*do,i,9,12,1lesize,i,100*enddovmesh,1et,2,link8!是否应用PIPE单元r,2,0.000576mp,ex,2,2.1e11mp,prxy,2,0.3type,2mat,2real,2k,11,150,150k,12,150,150,1400l,11,12lesize,13,100lmesh,13et,3,combin39keyopt,3,3,3mp,ex,3,36534r,3,0.012,340.98,0.023144013,676.822,0.041007667,1009.222rmore,0.060508355,1329.6,0.078506,1669.288667,0.13106,2019.533333rmore,0.201904,2332.511333,0.297857,2678.866667,0.381036,3004.466667rmore,0.491887,3340.466667,0.634276,3667.711333,2.961887,3864.733333rmore,2.989726,3982.733333,3.05323,4143.8,3.141894,4204.355333rmore,3.16088,4294.266667,3.232973,4385.644667,3.244492,4473.044667rmore,3.333723,4047.6,type,3mat,3real,3e,37,736e,86,737*do,i,567,579,1e,i,i+171*enddosaveet,4,shell63r,4,0.015,0.015,0.015,0.015mp,ex,4,2.5e4mp,prxy,4,0.3type,4mat,4real,4asel,,,,3,6,1aatt,4,4,4,amesh,3,6,1allselfinish/solunsel,s,loc,z,1400nsel,r,loc,x,0d,all,allallselnsel,s,loc,z,1400nsel,r,loc,x,300d,all,allallseltime,0*do,i,1,11,1time,if,737,fz,5000*icnvtol,f,,0.1,2nropt,fulloutres,all,allnsubst,10nlgeom,onsolvefdele,all*enddo*do,i,12,25,1f,737,fz,55000+2500*(i-11)time,icnvtol,f,,0.05,2nropt,fulloutres,all,allnsubst,10nlgeom,onsolvefdele,all*enddoANSYS中使用COMBIN39单元的实例2/PREP7ANTYPE,TRANS!NONLINEARTRANSIENTDYNAMICANALYSISET,1,COMBIN39,,,2!ELEMENTWITHDISPLACEMENTALONGNODALY-AXISET,2,MASS21,,,4!MASSWITHOUTROTARYINERTIAR,1,0.0,0.0,.1,.204,.2,.432!SPRINGDATAR,2,2588E-6!MASSDATAN,1N,2E,1,2TYPE,2REAL,2E,2D,1,ALLD,2,UXIC,2,UY,-1!INITIALDISPLACEMENTANDVELOCITYKBC,1!STEPLOADINGSAVEFINISH/SOLUSOLCONTROL,0CNVTOL,F,1,1E-4!FORCECONVERGENCECRITERIAOUTRES,NSOL,1NSUBST,5OUTPR,BASIC,NONETIME,.0002!TIMETOALLOWINITIALCHANGEINACCELERATIONLSWRITE!WRITELOADSTEPFILE1NSUBST,40OUTPR,BASIC,LASTTIME,0.18!TIMEARBITRARILYSELECTEDLSWRITE!WRITELOADSTEPFILE2LSSOLVE,1,2,1!READIN2LOADSTEPSANDSOLVEFINISH/POST26TIMERANGE,.003,.18NSOL,2,2,U,Y,2UYPRVAR,2!PRINTDISPLACEMENTS*GET,PER,VARI,2,EXTREM,TMIN*status,parmFINISHfrp钢筋混凝土加固梁的粘结滑移问题我最近在做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

!选择体2

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,