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

!定義參數(shù)和單元類型

!混凝土參數(shù)

*SET,hntra,20.5

*SET,ec,28800

*SET,possible,0.3

*SET,hntrl,2.35

!FRP參數(shù)

*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

!剛性墊塊

!定義實常數(shù)

r,1,50.3

!縱向受拉鋼筋實常數(shù)

r,2,28.3

!縱向受壓鋼筋，箍筋實常數(shù)

r,3

!混凝土

R,4,tf,tf,tf,tf,,,

!frp的實常數(shù)

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

!混凝土網(wǎng)格劃分較細為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,