boomsday 发表于 2006-5-4 22:11

求助!弦支穹顶的索预应力如何模拟?(附模型程序)

<P>我试了用降温法和初应变法作,但是总是不成。我把已经写好的程序发出来,那位大侠帮我看看,感激不尽。<BR>/prep7<BR>csys,2<BR>*do,i,1,56                  <BR>k,i,138,(i-1)*360/56,70.5<BR>*enddo<BR>*do,i,57,112<BR>k,i,138,(i-56)*360/56-180/56,72<BR>*enddo<BR>*do,i,113,168<BR>k,i,138,(i-113)*360/56,73.5<BR>*enddo<BR>*do,i,169,224<BR>k,i,138,(i-168)*360/56-180/56,75<BR>*enddo<BR>*do,i,225,280<BR>k,i,138,(i-225)*360/56,76.5<BR>*enddo<BR>*do,i,281,308<BR>k,i,138,(i-280)*360/28-360/56,78<BR>*enddo<BR>*do,i,309,336<BR>k,i,138,(i-309)*360/28,79.5<BR>*enddo<BR>*do,i,337,364<BR>k,i,138,(i-336)*360/28-360/56,81<BR>*enddo<BR>*do,i,365,392<BR>k,i,138,(i-365)*360/28,82.5<BR>*enddo<BR>*do,i,393,406<BR>k,i,138,(i-392)*360/14-360/28,84<BR>*enddo<BR>*do,i,407,420<BR>k,i,138,(i-407)*360/14,85.5<BR>*enddo<BR>*do,i,421,434<BR>k,i,138,(i-421)*360/14-360/28,87<BR>*enddo<BR>*do,i,435,441<BR>k,i,138,(i-435)*360/7+360/28,88.5<BR>*enddo<BR>k,442,138,0,90</P>
<P>et,1,beam188<BR>mp,ex,1,2.1e11<BR>mp,prxy,1,0.3<BR>mp,gxy,1,8e10<BR>sectype,1,beam,rect<BR>secdata,1.1,1.1<BR>sectype,2,beam,rect<BR>secdata,2.3,1.3<BR>sectype,3,beam,csolid</P>
<P><BR>*do,i,1,5<BR>*do,j,1,55,1<BR>   l,56*(i-1)+j,56*(i-1)+j+1<BR>*enddo<BR>l,i*56,i*56-55<BR>*enddo</P>
<P><BR>*do,i,1,4<BR>*do,j,1,27,1<BR>   l,280+28*(i-1)+j,280+28*(i-1)+j+1<BR>*enddo<BR>l,280+i*28,280+i*28-27<BR>*enddo</P>
<P><BR>*do,i,1,3<BR>*do,j,1,13,1<BR>   l,392+14*(i-1)+j,392+14*(i-1)+j+1<BR>*enddo<BR>l,392+i*14,392+i*14-13<BR>*enddo</P>
<P><BR>*do,j,1,6,1<BR>l,434+j,434+j+1<BR>*enddo<BR>l,441,435<BR></P>
<P>!建立径杆单元<BR>*do,i,1,4<BR>*do,j,1,56,1<BR>   l,56*(i-1)+j,56*i+j<BR>   *enddo<BR>*enddo</P>
<P>*do,j,2,56,1<BR>l,j,55+j<BR><BR>l,55+j,112+j<BR><BR>l,112+j,167+j<BR><BR>l,167+j,224+j<BR><BR>*enddo<BR>l,1,112<BR>l,112,113<BR>l,113,224<BR>l,224,225</P>
<P>*do,j,1,28,1<BR>l,223+2*j,280+j<BR>l,224+2*j,280+j<BR>*enddo<BR>*do,j,1,27,1<BR>l,225+2*j,280+j<BR><BR>*enddo<BR>l,225,308<BR>!建立第二道三层径干<BR>*do,i,1,3<BR>*do,j,1,28,1<BR>   l,280+28*(i-1)+j,280+28*i+j<BR>   *enddo<BR>*enddo<BR>*do,j,1,27,1<BR>l,280+j,309+j<BR>l,309+j,336+j<BR>l,336+j,365+j<BR>*enddo<BR>l,308,309<BR>l,309,364<BR>l,364,365<BR>!建立第三道径杆<BR>*do,j,1,14,1<BR>l,363+2*j,392+j<BR> l,364+2*j,392+j<BR>*enddo<BR>*do,j,1,13,1<BR>l,365+2*j,392+j<BR>*enddo<BR>l,365,406</P>
<P>!谢干<BR>*do,i,1,2<BR>*do,j,1,14,1<BR>   l,392+14*(i-1)+j,392+14*i+j<BR>   *enddo<BR>*enddo<BR>*do,j,1,13,1<BR>l,392+j,407+j<BR>*enddo<BR>l,406,407<BR>*do,j,1,13,1<BR>l,406+j,421+j<BR>*enddo<BR>l,420,421</P>
<P>!最后一圈<BR>*do,j,1,7,1<BR>l,419+2*j,434+j<BR>l,420+2*j,434+j<BR>*enddo</P>
<P>*do,j,1,6,1<BR>l,421+2*j,434+j<BR>*enddo<BR>l,421,441<BR>!屋顶节点杆<BR>*do,i,435,441,1<BR>l,i,442<BR>*enddo</P>
<P><BR>csys,0    !复制节点,要转换到卡式坐标系<BR>egjx=2e5 !Ey <BR>agjx=7600 !单根钢绞线面积 <BR>ehnt=4e4 !Eh <BR>xzxs=1.0e-5 !线胀系数 <BR>yjl1=3800000 !定义预加力<BR>yjl2=1700000<BR>yjl3=950000<BR>yjl4=650000<BR>yjl5=280000<BR>et,2,link8</P>
<P>r,2,agjx !定义link8单元的面积<BR>mp,ex,2,egjx !定义link8单元的弹性模量 <BR>mp,prxy,2,0.3 !定义link8单元的泊松系数 <BR>mp,alpx,2,xzxs !定义线膨胀系数 <BR>kgen,2,407,420,,,,-3.9,36!从网壳顶起第一道索<BR> !径向索</P>
<P><BR>*do,i,1,14,1      <BR>l,442+i,364+2*i<BR>*enddo<BR>*do,i,1,13,1      <BR>l,443+i,364+2*i<BR>*enddo<BR>l,443,392</P>
<P>!环向索<BR>*do,i,1,13,1<BR>l,442+i,443+i<BR>*enddo<BR>l,456,443   </P>
<P>type,1<BR>mat,1<BR>*do,i,1,14,1   !竖向撑杆<BR>l,442+i,406+i<BR>*enddo</P>
<P><BR>kgen,2,366,392,2,,,-3.9,91   !从网壳顶起第二道索<BR>!径向索</P>
<P><BR>*do,i,1,14,1      <BR>l,455+2*i,307+2*i<BR>*enddo<BR>*do,i,1,13,1      <BR>l,455+2*i,309+2*i<BR>*enddo<BR>l,483,309             !已经到了483号节点了<BR>!环向索<BR>*do,i,1,13,1<BR>l,455+2*i,457+2*i<BR>*enddo<BR>l,483,457</P>

<P><BR>mat,1<BR>*do,i,1,14,1   !竖向撑杆<BR>l,455+2*i,364+2*i<BR>*enddo</P>


<P>kgen,2,309,336,1,,,-3.9,175   !从网壳顶起第三道索</P>
<P><BR>*do,i,1,28,1                  !径向索<BR>   l,483+i,224+2*i<BR>*enddo<BR>*do,i,1,27,1<BR>   l,484+i,224+2*i<BR>*enddo<BR>l,484,280</P>
<P>*do,i,1,27,1                  !环向索<BR>l,483+i,484+i<BR>*enddo<BR>l,511,484            !倒511节点</P>

<P><BR>*do,i,1,28,1   !竖向撑杆<BR>l,483+i,308+i<BR>*enddo</P>

<P>kgen,2,226,280,1,,,-3.9,286   !从网壳顶起第四道索</P>
<P><BR>*do,i,1,28,1                  !径向索<BR>   l,510+2*i,111+2*i<BR>*enddo<BR>*do,i,1,27,1<BR>   l,510+2*i,113+2*i<BR>*enddo<BR>l,566,113</P>
<P>*do,i,1,27,1                  !环向索<BR>l,510+2*i,512+2*i<BR>*enddo<BR>l,566,512            !倒566节点</P>

<P><BR>*do,i,1,28,1   !竖向撑杆<BR>l,510+2*i,224+2*i<BR>*enddo</P>

<P>kgen,2,113,167,1,,,-3.9,454   !从网壳顶起第五道索</P>
<P><BR>*do,i,1,28,1                  !径向索<BR>   l,565+2*i,2*i<BR>*enddo<BR>*do,i,1,27,1<BR>   l,567+2*i,2*i<BR>*enddo<BR>l,567,56</P>
<P>*do,i,1,27,1                  !环向索<BR>l,565+2*i,567+2*i<BR>*enddo<BR>l,621,567            !倒621节点</P>

<P>*do,i,1,28,1   !竖向撑杆<BR>l,565+2*i,111+2*i<BR>*enddo</P>
<P><BR>lsel,s,line,,1,1267,1!划分网壳单元<BR>lesize,all,,,1<BR>LATT,1,,1, , , ,1<BR>lmesh,1,1267,1</P>
<P>lsel,none,line,,1,1267,1<BR>lsel,all<BR>lsel,s,line,,1268,1309,1    !划分索单元<BR>lesize,all,,,1<BR>LATT,1,2,2, , , ,3<BR>lmesh,1268,1309,1</P>
<P>lsel,none,line,,1268,1309,1<BR>lsel,all<BR>lsel,s,line,,1324,1365,1<BR>lesize,all,,,1<BR>LATT,1,2,2, , , ,3<BR>lmesh,1324,1365,1</P>
<P>lsel,none,line,,1324,1365,1<BR>lsel,all<BR>lsel,s,line,,1380,1463,1<BR>lesize,all,,,1<BR>LATT,1,2,2, , , ,3<BR>lmesh,1380,1463,1</P>
<P>lsel,none,line,,1380,1463,1<BR>lsel,all<BR>lsel,s,line,,1492,1575,1<BR>lesize,all,,,1<BR>LATT,1,2,2, , , ,3<BR>lmesh,1492,1575,1</P>
<P>lsel,none,line,,1492,1575,1<BR>lsel,all<BR>lsel,s,line,,1604,1687,1<BR>lesize,all,,,1<BR>LATT,1,2,2, , , ,3<BR>lmesh,1604,1687,1</P>
<P>lsel,none,line,,1604,1687,1<BR>lsel,all<BR>lsel,s,line,,1310,1323,1    !划分撑杆单元<BR>lesize,all,,,1<BR>LATT,1,,1, , , ,2<BR>lmesh,1310,1323,1</P>
<P>lsel,none,line,,1310,1323,1<BR>lsel,all<BR>lsel,s,line,,1366,1379,1<BR>lesize,all,,,1<BR>LATT,1,,1, , , ,2<BR>lmesh,1366,1379,1</P>
<P>lsel,none,line,,1366,1379,1<BR>lsel,all<BR>lsel,s,line,,1464,1491,1<BR>lesize,all,,,1<BR>LATT,1,,1, , , ,2<BR>lmesh,1464,1491,1</P>
<P>lsel,none,line,,1464,1491,1<BR>lsel,all<BR>lsel,s,line,,1576,1603,1<BR>lesize,all,,,1<BR>LATT,1,,1, , , ,2<BR>lmesh,1576,1603,1</P>
<P>lsel,none,line,,1576,1603,1<BR>lsel,all<BR>lsel,s,line,,1688,1715,1<BR>lesize,all,,,1<BR>LATT,1,,1, , , ,2<BR>lmesh,1688,1715,1<BR>finish<BR>lsel,none,line,,1688,1715,1<BR>lsel,all<BR>esel,all</P>
<P>/solu<BR>antype,0   !静力分析<BR>*do,i,1,56,1   !施加边界条件 <BR>   dk,i,ux<BR>   dk,i,uy<BR>   dk,i,uz<BR>*enddo<BR>*do,i,1,442,1    !施加外力<BR>   fk,i,fz,-1<BR>*enddo</P>
<P><BR>*do,i,1296,1309,1   <BR>bfl,i,temp,-yjl5/(xzxs*egjx*agjx) !对5钢绞线施加温度<BR>*enddo<BR>*do,i,1352,1365,1   <BR>bfl,i,temp,-yjl4/(xzxs*egjx*agjx) !对4钢绞线施加温度<BR>*enddo<BR>*do,i,1436,1463,1   <BR>bfl,i,temp,-yjl3/(xzxs*egjx*agjx) !对3钢绞线施加温度<BR>*enddo<BR>*do,i,1548,1575,1   <BR>bfl,i,temp,-yjl2/(xzxs*egjx*agjx) !对3钢绞线施加温度<BR>*enddo<BR>*do,i,1660,1687,1   <BR>bfl,i,temp,-yjl1/(xzxs*egjx*agjx) !对3钢绞线施加温度<BR>*enddo<BR>solve<BR>finish<BR><BR>不知道对不对。</P>

AaronSpark 发表于 2006-5-5 01:28

回复:(boomsday)求助!弦支穹顶的索预应力如何模拟...

从整个过程上看应该没什么大问题<BR>具体细节上你自己检查检查吧,别人没你对这个模型熟悉<BR>有什么具体不对的地方大家可以帮你看看
页: [1]
查看完整版本: 求助!弦支穹顶的索预应力如何模拟?(附模型程序)