《ANSYS中使用COMBIN39单元的实例.doc》由会员分享，可在线阅读，更多相关《ANSYS中使用COMBIN39单元的实例.doc（10页珍藏版）》请在金锄头文库上搜索。
1、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 lesize,i,50 *enddo *do,i,9,12
2、,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,1009.222 rmore,0.060508355,1329
3、.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,3.232973,4385.644667,3.244492
4、,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,1400 nsel,r,loc,x,0 d,all,all a
5、llsel 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,all nsubst,10 nlgeom,on solve
6、 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 DATARMORE,.3,.708,.4,1.056,.5,1.5RMORE,.6,2.064,.7,2.772,.8,3
7、.648RMORE,.9,4.716,1.0,6.0R,2,2588E-6 ! MASS DATAN,1N,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 CH
8、ANGE IN ACCELERATIONLSWRITE ! WRITE LOAD STEP FILE 1NSUBST,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*sta
9、tus,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,
10、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!法向弹簧
11、keyopt,5,1,0keyopt,5,2,1 KEYOPT,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,possible tb,conc,1 !混凝土自带的破坏准则tbdata,0.5,0.9,hntrl,-1!混凝土本构关系 tb,miso,1,23tbpt,0.000005,hntra*0.007
12、4627tbpt,0.00001,hntra*0.00992537tbpt,0.00005,hntra*0.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,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.0
