下面的兩個受拉筋施加預應力,對它耦合,那么其他的箍筋還有受壓筋采用切割實體的形式?模擬出來不收斂,模擬不出來,我把命令流放下面,希望各位大俠看看,有啥錯誤沒,我在王新敏的命令流上改的
!預應力簡支梁彈性分析(2007-武漢理工大學-預應力FRP筋受彎構件非線性有限元分析(CFRP+ANSYS))!節(jié)點耦合法,初應變施加預應力!-------------------------------finish/clear/prep7!1.定義單元與材料性質--------------------et,1,solid65et,2,link8et,3,solid65mp,ex,1,3.25e4mp,prxy,1,0.2fc=19.1ft=1.71tb,concr,1tbdata,,0.5,1,1.71,-1tb,mkin,1tbtemp,,straintbdata,1,0.0002,0.0006,0.001,0.002,0.0033tbtemp,0tbdata,1,fc*0.34,fc*0.51,fc*0.75,fc,fc*0.85tbplot,mkin,1mp,ex,2,8.44e4 !預應力cfrp筋mp,prxy,2,0.23mp,ex,3,2.1e5mp,prxy,3,0.3 !受壓鋼筋tb,biso,3tbdata,,210,0mp,ex,4,2.06e5 !剛性墊塊mp,prxy,4,0.2pi=acos(-1)as=0.25*pi*9.5*9.5r,1,0r,2,as,70E3/(as*8.44e4) !預應力cfrp筋r,3,0.25*pi*6*6 r,4,0!2.分別創(chuàng)建體和力筋線blc4,,,160,240,2000blc4,,,160,-50,100 wpoffs,,,2000blc4,,,160,-50,-100wpcsys,-1wpoffs,,240,650blc4,,,160,50,100wpoffs,,,600blc4,,,160,50,100wpcsys,-1lsel,nonek,41,30,20,0k,42,30,20,2000l,41,42k,43,130,20,0k,44,130,20,2000l,43,44cm,sline,linewpcsys,-1*do,i,1,39$wpoff,,,50$vsbw,all$*enddowpcsys,-1wprota,,,90wpoffs,,,30vsbw,allwpoffs,,,100vsbw,allwpcsys,-1wprota,,90wpoffs,,,-20vsbw,allwpoffs,,,-200vsbw,allvglue,allnummrg,allnumcmp,all!3分別劃分單元網(wǎng)格elemsiz=50lsel,s,loc,x,30$lsel,r,loc,y,220$cm,JLJL,line$latt,3,3,2$lesize,all,elemsizlsel,s,loc,x,130$lsel,r,loc,y,220 $cm,JLJR,line$latt,3,3,2$lesize,all,elemsizlsel,s,tan1,zlsel,r,loc,x,30,130$lsel,r,loc,y,20,220 lsel,u,loc,z,0$lsel,u,loc,z,2000cm,gj,line$latt,3,3,2$lesize,all,elemsizcmsel,s,JLJLcmsel,a,JLJR$cmsel,a,gjlmesh,all$lsel,all!體網(wǎng)格劃分(此處手動操作)cmsel,s,sline !定義預應力CFRP筋單元尺寸,劃分單元latt,2,2,2lesize,all,50lmesh,alllsel,allesize,50vsel,s,loc,y,0,240vatt,1,4,1mshaps,0mshkey,1vmesh,allesize,50vatt,4,4,1lsel,s,loc,y,0,-50lsel,allmshaps,0mshkey,1vmesh,allesize,50vatt,4,4,1lsel,s,loc,y,240,290lsel,allmshaps,0mshkey,1vmesh,all!4.定義節(jié)點耦合自由度cmsel,s,slinensll,s,1cm,gsnode,node*get,senum,node,,count*dim,snode,,senum*dim,cnode,,senum*get,nodei,node,,num,minsnode(1)=nodei*do,i,2,senumsnode(i)=ndnext(snode(i-1))*enddoallsel,allcmsel,u,gsnode*do,i,1,senumnodei=snode(i)nodej=nnear(nodei)cnode(i)=nodejnsel,u,,,nodej*enddoallsel,all*do,i,1,senumcp,next,all,snode(i),cnode(i)*enddo!5.定義約束(手動操作)!6.施加荷載并求解(需要設置時間以及斜坡荷載,對于預應力的施加也需要查閱資料)/soluacel,,9800antype,0outres,all,allautots,onneqit,500nsubst,40cnvtol,u,,0.05cnvtol,f,,0.05solve
發(fā)自小木蟲Android客戶端 |