ANSYS中使用COMBIN39单元的实例

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,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,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!SPRINGDATARMORE,.3,.708,.4,1.056,.5,1.5RMORE,.6,2.064,.7,2.772,.8,3.648RMORE,.9,4.716,1.0,6.0R,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*2et,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,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,,23tbpt,,0.000005,hntra*0.0074627tbpt,,0.00001,hntra*0.00992537tbpt,,0.00005,hntra*0.0491293tbpt,,0.0001,hntra*0.0970149tbpt,,0.0002,hntra*0.19tbpt,,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,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.85mp,ex,2,1.94e5!纵向受拉钢筋mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,368.3,2000mp,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,LINEallsel,alllgen,2,13,,,,120lgen,2,14,,,,120lsel,s,line,,15,16lplotcm,zyj,lineallsel,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,lineWPAVE,0,0,0!CFRPwprota,,90!旋转工作平面yz，90°wpoffs,,,0.835!偏移工作平面zcfrp面与混凝土面之间距离为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,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,,,43LESIZE,42,50,,,,,,,0LESIZE,43,20,,,,,,,0m