(完整word版)ANSYS中使用COMBIN39单元的实例

1、院整院整word版）ANSYS中使用COMBIN39单元的实例ANSYS中使用COMBIN39单元的实例1finish/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

2、,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,333!定义Z方向弹簧mp,ex,3,36534r,3,0.012,340.98,0.023144013,676.822,0.041007667,1009.222rmore,0.060508355,1329.6,0.078506,1

3、669.288667,0.13106,2006.533333rmore,0.200204,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

4、,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,

5、allallseltime,0*do,i,1,11,1time,if,737,fz,5000*icnvtol,f,0.1,2nropt,fulloutres,all,allnsubst,10nlgeom,on院整院整word版）ANSYS中使用COMBIN39单元的实例院整院整word版）ANSYS中使用COMBIN39单元的实例RMORE,.9,4.716,1.0,6.0RMORE,.9,4.716,1.0,6.0RMORE,.9,4.716,1.0,6.0RMORE,.9,4.716,1.0,6.0solvefdele,all*enddo*do,i,12,25,1f,737,fz,5500

6、0+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!SPRINGDAT

7、ARMORE,.3,.708,.4,1.056,.5,1.5RMORE,.6,2.064,.7,2.772,.8,3.648院整院整word版）ANSYS中使用COMBIN39单元的实例院整院整word版）ANSYS中使用COMBIN39单元的实例R,2,2588E-6!MASSDATAN,1N,2E,1,2TYPE,2REAL,2E,2D,1,ALLD,2,UX!INITIALDISPLACEMENTANDVELOCITY!INITIALDISPLACEMENTANDVELOCITY!STEPLOADINGKBC,1SAVEFINISH/SOLUSOLCONTROL,0CNVTOL,F,1,

8、1E-4OUTRES,NSOL,1NSUBST,5!FORCECONVERGENCECRITERIAOUTPR,BASIC,NONETIME,.0002TIME,.0002!TIMETOALLOWINITIALCHANGEINACCELERATIONLSWRITE!WRITELOADSTEPFILE1NSUBST,40OUTPR,BASIC,LASTTIME,0.18!TIMEARBITRARILYSELECTED院整院整word版冶在丫5中使用COMBIN39单元的实例院整院整word版冶在丫5中使用COMBIN39单元的实例院整院整word版）ANSYS中使用COMBIN39单元的实例*S

9、ET,hntrl,2.35*SET,hntrl,2.35LSWRITE!WRITELOADSTEPFILE2LSSOLVE,1,2,1LSSOLVE,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单元来模

10、拟frp与混凝土之间的粘结滑移，有一些错误，希望高手们有时间帮我看看。多谢了很多人说，采用ansys算粘结滑移很难算到收敛，但是不知道有没有算出来的，只要有可能，我就不放弃了。命令流如下：!第一组/prep7!定义参数和单元类型!混凝土参数*SET,hntra,20.5*SET,ec,28800*SET,possible,0.3!FRP参数SET,ff,1800!抗拉强度*SET,ef,220000!弹性模量!FRP纤维布的极限拉应变SET,ckf,ff/ef!FRP厚度SET,tf,0.121*2et,1,solid65!混凝土单元KEYOPT,1,1,1!不考虑大变形KEYOPT,1,7,

11、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!垫块院整院整word版）ANSYS中使用COMBIN39单元的实例院整院整word版）ANSYS中使用COMBIN39单元的实例院整院整word版）ANSYS中使用COMBIN39单元的实例allsel,allallsel,alltbpt,0.0002,hntra*0.19tbpt,0.000

12、2,hntra*0.19ET,5,COMBIN39!法向弹簧keyopt,5,1,0keyopt,5,2,1KEYOPT,5,3,3r,6,0,0,0.0005,29340,0.001,41577!粘结滑移关系rmore,0.005,92970,0.01,131478,0.015,161028rmore,0.02,185939,0.025,207886,0.029214,224725!定义材料属性mp,ex,1,ec!混凝土材料属性mp,prxy,1,possibletb,conc,1!混凝土自带的破坏准则tbdata,0.5,0.9,hntrl,-1!混凝土本构关系tb,miso,1,23t

13、bpt,0.000005,hntra*0.0074627tbpt,0.00001,hntra*0.00992537tbpt,0.00005,hntra*0.0491293tbpt,0.0001,hntra*0.0970149tbpt,0.0004,hntra*0.36tbpt,0.0004,hntra*0.36tbpt,0.0006,hntra*0.51tbpt,0.0008,hntra*0.64tbpt,0.0010,hntra*0.75tbpt,0.0012,hntra*0.84tbpt,0.0014,hntra*0.91tbpt,0.0016,hntra*0.96tbpt,0.0018,

14、hntra*0.99tbpt,0.0020,hntra*1.00tbpt,0.0022,hntra*0.98tbpt,0.0024,hntra*0.97tbpt,0.0026,hntra*0.95tbpt,0.0028,hntra*0.93tbpt,0.003,hntra*0.92tbpt,0.0032,hntra*0.9tbpt,0.0034,hntra*0.88tbpt,0.0036,hntra*0.87tbpt,0.0038,hntra*0.85!纵向受拉钢筋mp,ex,2,1.94e5!纵向受拉钢筋mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,368.3,200

15、0mp,ex,3,1.97e5!横向箍筋，受压钢筋材料属性mp,prxy,3,0.3tb,bkin,3,1,2,1tbdata,410.6,2000定义FRPmp,ex,4,efmp,prxy,4,0.3mp,ex,5,ec!垫块材料属性mp,prxy,5,possible!建模blc4,100,200,-1100k,9,20,40,0!钢筋k,10,20,40,-1100l,9,10lgen,2,13,60lsel,s,line,13,14LPLOTCM,zlj,LINElgen,2,13,120lgen,2,14,120lsel,s,line,15,16lplotcm,zyj,lineal

16、lsel,allk,17,20,40,-500!箍筋k,18,80,40,-500k,19,20,160,-500l,17,18l,17,19lgen,2,17,120lgen,2,18,60lsel,s,line,17,20,1LPLOT*do,ii,1,5lgen,2,17,20,1,-ii*100ii=ii+1*enddoLPLOTcm,gj,line院整院整word版）ANSYS中使用COMBIN39单元的实例院整院整word版）ANSYS中使用COMBIN39单元的实例WPAVE,0,0,0!CFRPwprota,90!旋转工作平面yz,90wpoffs,0.835!偏移工作平面zc

17、frp面与混凝土面之间距离为0.835mmblc4,100,-1000wpoffs,-0.835!支座面与cfp下表面距离为0.835mmwpoffs,-1100blc4,100,100,20!支座处混凝土垫块100*100*20wpoffs,600wpoffs,-200wprota,-90blc4,100,20,-100!荷载处混凝土垫块是100*100*20allsel,all!划分单元VSEL,s,volu,1!混凝土网格划分较细为25*50*25VATT,1,3,1LSEL,s,line,12LSEL,a,line,7LSEL,a,line,8LESIZE,12,50,0LESIZE,

18、7,20,0LESIZE,8,20,0MSHAPE,0,3DMSHKEY,1VMESH,1cmsel,s,zlj,line!主la筋划分为50mmlATT,2,1,2LESIZE,all,50,0LMESH,allcmsel,s,zyj,line!主ya筋划分为50mmlATT,3,2,2LESIZE,all,50,0LMESH,allcmsel,s,gj,line!箍筋为20mmlplotlATT,3,2,2LESIZE,all,20,0LMESH,allasel,s,area,7!CFRP划分单元为四边形20*50aplotAATT,3,4,3lsel,s,42lsel,a,43LESIZ

19、E,42,50,0LESIZE,43,20,0mshape,0,2DMSHKEY,1AMESH,7allsel,all!划分支座垫块为25*50*20VSEL,s,2!选择体2VATT,5,3,4LESIZE,50,50,0LESIZE,51,20,0LESIZE,55,20,0LESIZE,47,50,0LESIZE,49,20,0LESIZE,53,20,0MSHAPE,0,3DMSHKEY,1VMESH,all/VIEW,1,1,1,1allsel,all!划分支座垫块为25*50*20VSEL,s,3!选择体3VATT,5,3,4LESIZE,67,50,0LESIZE,62,20,0LESIZE,61,20,0LESIZE,66,50,0LESIZE,57,20,0LESIZE,58,20,0MSHAPE,0,3DMSHKEY,1