求助！弦支穹顶的索预应力如何模拟？（附模型程序）

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

 mat,1 *do,i,1,14,1 !竖向撑杆 l,455+2*i,364+2*i *enddo

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

 *do,i,1,28,1 !竖向撑杆 l,483+i,308+i *enddo

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

 *do,i,1,28,1 !竖向撑杆 l,510+2*i,224+2*i *enddo

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

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

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

从整个过程上看应该没什么大问题 具体细节上你自己检查检查吧，别人没你对这个模型熟悉 有什么具体不对的地方大家可以帮你看看

账号		自动登录	找回密码
密码			我要加入

求助！弦支穹顶的索预应力如何模拟？（附模型程序）

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（boomsday）求助！弦支穹顶的索预应力如何模拟...