《ANSYS中使用COMBIN39单元的实例(11页).doc》由会员分享,可在线阅读,更多相关《ANSYS中使用COMBIN39单元的实例(11页).doc(11页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-ANSYS中使用COMBIN39单元的实例-第 11 页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
2、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 mp,ex,3,36534 r,3,0.012,340.98,0.023144013,676.822,0.041007667,10
3、09.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,
4、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,14
5、00 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,
6、all nsubst,10 nlgeom,on solve fdele,all *enddo ANSYS中使用COMBIN39单元的实例2/PREP7ANTYPE,TRANS ! NONLINEAR TRANSIENT DYNAMIC ANALYSISET,1,COMBIN39,2 ! ELEMENT WITH DISPLACEMENT ALONG NODAL Y-AXISET,2,MASS21,4 ! MASS WITHOUT ROTARY INERTIAR,1,0.0,0.0,.1,.204,.2,.432 ! SPRING DATAR,2,2588E-6 ! MASS DATAN,1N,
7、2E,1,2TYPE,2REAL,2E,2D,1,ALLD,2,UXIC,2,UY,-1 ! INITIAL DISPLACEMENT AND VELOCITYKBC,1 ! STEP LOADINGSAVEFINISH/SOLU SOLCONTROL,0CNVTOL,F,1,1E-4 ! FORCE CONVERGENCE CRITERIAOUTRES,NSOL,1NSUBST,5OUTPR,BASIC,NONETIME,.0002 ! TIME TO ALLOW INITIAL CHANGE IN ACCELERATIONLSWRITE ! WRITE LOAD STEP FILE 1NS
8、UBST,40OUTPR,BASIC,LASTTIME,0.18 ! TIME ARBITRARILY SELECTEDLSWRITE ! WRITE LOAD STEP FILE 2LSSOLVE,1,2,1 ! READ IN 2 LOAD STEPS AND SOLVEFINISH/POST26TIMERANGE,.003,.18NSOL,2,2,U,Y,2UY PRVAR,2 ! PRINT DISPLACEMENTS*GET,PER,VARI,2,EXTREM,TMIN*status,parmFINISHfrp钢筋混凝土加固梁的粘结滑移问题我最近在做frp钢筋混凝土加固梁的剥离问题,
9、采用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
10、,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,1 KEYOPT,5,3,3r,6,0,0,0.0005,29
11、340,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,possible tb,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.
12、0491293tbpt,0.0001,hntra* 0.0970149tbpt,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,
13、hntra*0.97 tbpt,0.0026,hntra*0.95 tbpt,0.0028,hntra*0.93 tbpt,0.003,hntra*0.92tbpt,0.0032,hntra*0.9tbpt,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.3tb,bkin,2,1,2,1tbdata,368.3,2000mp,ex,3,1.97e5 !横向箍筋,受压钢筋材料属性mp,prxy,3,0.3tb,bkin,3,1,2,1tbdata,4
14、10.6,2000!定义FRPmp,ex,4,ef mp,prxy,4,0.3mp,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,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,2
15、0,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,90wpoffs,0.835 !偏移工作平面zcfrp面与混凝土面之间距离为0.835mmblc4,100,-1000wpoffs,-0.835 !支座面与cfp下表面距离为0.835mmwpoffs,-1100blc4,100,100,20 !支座处混凝土垫
16、块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,0 MSHAPE,0,3DMSHKEY,1VMESH,1cmsel,s,zlj,line !主la筋划分为50mmlATT,2,1,2LESIZE,all,50,0
17、 LMESH,allcmsel,s,zyj,line !主ya筋划分为50mmlATT,3,2,2LESIZE,all,50,0 LMESH,allcmsel,s,gj,line !箍筋为20mmlplotlATT,3,2,2LESIZE,all,20,0 LMESH,allasel,s,area,7 !CFRP划分单元为四边形20*50aplotAATT,3,4,3lsel,s,42lsel,a,43LESIZE,42,50,0 LESIZE,43,20,0 mshape,0,2D MSHKEY,1AMESH,7allsel,all !划分支座垫块为25*50*20VSEL,s,2 !选择体
18、2VATT,5,3,4LESIZE,50,50,0LESIZE,51,20,0LESIZE,55,20,0LESIZE,47,50,0LESIZE,49,20,0LESIZE,53,20,0 MSHAPE,0,3DMSHKEY,1VMESH,all/VIEW,1,1,1,1 allsel,all !划分支座垫块为25*50*20VSEL,s,3 !选择体2VATT,5,3,4LESIZE,67,50,0LESIZE,62,20,0LESIZE,61,20,0LESIZE,66,50,0LESIZE,57,20,0LESIZE,58,20,0 MSHAPE,0,3DMSHKEY,1VMESH,a
19、llallsel,allnsel,s,loc,y,-0.835,0 !弹簧单元0nplotTYPE,5REAL,6ESYS,0eintf,0.84,low !结点号码从最低点到最高点编号allsel,allnsel,s,loc,y,-1.67,0!混凝土和CFRP界面重合节点 约束UX UY方向CPINTF,UX,1.67CPINTF,Uy,1.67nsel,s,loc,y,-1.67,-0.835!cfrp内部结点耦合CPINTF,Uz,0.835ALLSEL,ALL !其他节点耦合CPINTF,UXCPINTF,UyCPINTF,Uzallsel,all!求解设置/soluantype,0
20、allsel,all asel,s,9 !右边界约束da,all,uxda,all,uyallsel,allasel,s,loc,z,0 !对称面约束da,all,symmallsel,all !施加外部载荷asel,s,17 sfa,all,pres,2.456allsel,allautots,ontime,1nsubst,100lnsrch,onpred,onoutres,basic,allncnv,2allsel,allSOLVE FINISH 有一些错误,请高手们帮我看看* ERROR * SUPPRESSED MESSAGE CP = 7.375 TIME= 21:16:15Node 154 direction UY included in coupled sets 696 and 275.