| 查看: 725 | 回复: 0 | ||
[求助]
ANSYS预应力混凝土梁模拟节点耦合法
|
|
下面的两个受拉筋施加预应力，对它耦合，那么其他的箍筋还有受压筋采用切割实体的形式?模拟出来不收敛，模拟不出来，我把命令流放下面，希望各位大侠看看，有啥错误没，我在王新敏的命令流上改的 !预应力简支梁弹性分析(2007-武汉理工大学-预应力FRP筋受弯构件非线性有限元分析（CFRP+ANSYS）)!节点耦合法，初应变施加预应力!-------------------------------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.分别创建体和力筋线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分别划分单元网格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!体网格划分(此处手动操作)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.定义节点耦合自由度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 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
投票:
5