| 查看: 731 | 回复: 0 | ||
[求助]
ansys初学者自己编写的命令流，出现大变形结构不收敛，求大神帮我找找问题
|
|
附命令流，请大神解答 FINISH /CLEAR /PREP7 !定义单元类型 ET,1,beam188 !定义钢筋的单元类型 sectype,1,beam,csolid secdata,6,, MP,EX,1,2.06E5 MP,PRXY,1,0.25 MP,DENS,1,7850 ET,2,beam188 sectype,2,beam,csolid secdata,4,, MP,EX,2,2.06E5 MP,PRXY,2,0.25 MP,DENS,2,7850 TB,BISO,1,1,2,1 !双线性各向同性强化模型 TBDATA,,235,0 TB,BISO,2,1,2,1 TBDATA,,300,0 ET,3,SOLID65 !定义混凝土的单元类型 MP,EX,3,13585 MP,PRXY,3,0.2 fc=14.3$ft=1.43 tb,concr,3$tbdata,,0.5,0.95,ft,-1 tb,miso,3,,11 tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64$tbpt,,0.0010,fc*0.75$tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99 ET,4,solid45 !定义钢管的单元类型 MP,EX,4,2.06E5 MP,PRXY,4,0.25 MP,DENS,4,7850 TB,BISO,4,1,2,1 !双线性各向同性强化模型 TBDATA,,235,0 !建立几何模型(钢筋) k,1,20,20,820 k,2,20,20,1280 l,1,2 k,3,20,80,820 k,4,20,80,1280 l,3,4 k,5,80,20,820 k,6,80,20,1020 l,5,6 k,7,80,20,1080 k,8,80,20,1280 l,7,8 k,9,80,80,820 k,10,80,80,1020 l,9,10 k,11,80,80,1080 k,12,80,80,1280 l,11,12 k,13,280,20,1020 l,6,13 k,14,280,80,1020 l,10,14 k,15,280,20,1080 l,7,15 k,16,280,80,1080 l,11,16 k,17,20,20,860 k,18,20,80,860 k,19,80,80,860 k,20,80,20,860 l,17,18 l,18,19 l,19,20 l,20,17 lglue,11,12,13,14 lgen,3,11,14,1,,,60 k,29,20,20,1120 k,30,20,80,1120 k,31,80,80,1120 k,32,80,20,1120 l,29,30 l,30,31 l,31,32 l,32,29 lglue,23,24,25,26 lgen,3,23,26,1,,,60 k,41,120,20,1020 k,42,120,80,1020 k,43,120,80,1080 k,44,120,20,1080 l,41,42 l,42,43 l,43,44 l,44,41 lglue,35,36,37,38 lgen,3,35,38,1,60 !建立几何模型(混凝土) k,53,,100 k,54,100,100 k,55,100 k,56,,, a,53,54,55,56 vext,all,,,,,2100 k,100,100,100,1100 k,101,100,100,1000 k,102,100,,1000 k,103,100,,1100 a,100,101,102,103 vext,7,,,1000 vglue,all !钢筋划分 elemsiz=20 lsel,s,line,,1,10 latt,1,,1$lesize,all,elemsiz lsel,s,line,,11,46 latt,2,,1$lesize,all,elemsiz lsel,s,type,,1 lmesh,all !建立外包钢管 k,201,-2,102 k,202,102,102 k,203,102,-2 k,204,-2,-2 k,205,,100 k,206,100,100 k,207,100, k,208,, a,201,202,203,204 a,205,206,207,208 asba,3,4 vext,5,,,,,1000 vgen,2,1,,,,,1100 k,209,102,102,1102 k,210,102,102,998 k,211,102,-2,998 k,212,102,-2,1102 k,213,102,100,1100 k,214,102,100,1000 k,215,102,,1000 k,216,102,,1100 a,209,210,211,212 a,213,214,215,216 asba,34,35 vext,36,,,998 wpcsys,-1$wprota,,-90$vsbw,1$vsbw,4$vsbw,5 wpcsys,-1$wpoff,,100$wprota,,-90$vsbw,7$vsbw,8$vsbw,9 wpcsys,-1$wprota,,,-90$vsbw,11$vsbw,13$vsbw,1$vsbw,6 wpcsys,-1$wpoff,100$wprota,,,-90$vsbw,16$vsbw,17$vsbw,18$vsbw,19 wpcsys,-1$wpoff,,,1000$vsbw,4$vsbw,15 wpcsys,-1$wpoff,,,1100$vsbw,24$vsbw,25 wpcsys,-1$wpoff,1000$wprota,,,-90$asbw,all wpcsys,-1 esize,20 vsel,s,loc,x,,100$vsel,r,loc,y,,100$vsel,r,loc,z,0,2100 vatt,3,,3$vsweep,all vsel,all vsel,s,loc,x,100,1100$vsel,r,loc,y,,100$vsel,r,loc,z,1000,1100 vatt,3,,3$vsweep,all vsel,all vsel,s,loc,x,-2,0$vsel,r,loc,y,-2,102$vsel,r,loc,z,0,2100 vatt,4,,4$vsweep,all vsel,all vsel,s,loc,x,60,102$vsel,r,loc,y,-2,102$vsel,r,loc,z,0,2100 vatt,4,,4$vsweep,all vsel,all vsel,s,loc,x,-2,2100$vsel,r,loc,y,100,102$vsel,r,loc,z,0,2100 vatt,4,,4$vsweep,all vsel,all vsel,s,loc,x,-2,2100$vsel,r,loc,y,-2,0$vsel,r,loc,z,0,2100 vatt,4,,4$vsweep,all vsel,all vsel,s,loc,x,100,1100$vsel,r,loc,y,-2,102$vsel,r,loc,z,998,1000 vatt,4,,4$vsweep,all vsel,all vsel,s,loc,x,100,1100$vsel,r,loc,y,-2,102$vsel,r,loc,z,1100,1102 vatt,4,,4$vsweep,all nummrg,node /solu asel,s,loc,z,0 da,all,all asel,s,loc,z,2100 da,all,all allsel,all sfa,169,1,pres,10 allsel,all antype,0 nsubst,100 outres,all,all autos,on neqit,100 cnvtol,u,,0.05 solve |
回复此楼
