悬索桥集锦
一座悬索桥 本帖最后由 liushui 于 2013-6-3 23:40 编辑seuzhy 发表于 2007-2-2 08:52
我体会是:悬索桥是一种典型的柔性结构,具有很强的几何非线性,但一般不存在材料非线性问题。ANSYS分析的 ...
在下不才,愿意为之一解,有考虑不周之处,请各位多多海涵和指教。主缆的几何非线性比较明显,尤其是跨度大了之后。有限元进行动力分析的时候,单位分的细一点,可以更好的模拟该效应。在ansys中可以自己编写一个迭代的程序,直到模拟主缆的link单位的应变收敛位置,因为由弹性力学知道,应力等于应变乘以弹性模量。如果手上没有这个迭代程序,也不要紧,可以手动实现。具体作法为,先给link单位定义一个应变值,但是不要太离谱,这个值肯定不是最终值,先运行计算一遍。算完后,用ansys里面的etable这个表格提取应力的命令,用这个应力值代替命令流里初始的应变值。一般的桥,这样就差不多收敛了,因为ansys会帮你找到一个合理值。如果不放心,可以重复以上操作。当然以上作法是计算动力特性可以这么做,因为动力特性这样分析,误差很小。但是如果是施工设计,调主缆的内力,估计要考虑更多因素,还要模拟整个施工过程。本人就曾经这么手动调过主缆的初拉力,算出来的动力特性都非常精确。以上仅供参考,各位有更好的办法,也请赐教,刚好学习下,大家集思广益啊。
楼上谈论的问题好高深,表示膜拜! 转自关海土木
建模型
finish
/clear
/filename,xuansuoqiao
/title,xuansuodyna
/prep7
et,1,beam4
et,2,link8
mp,ex,1,2.1e11
mp,dens,1,14827
mp,nuxy,1,0.3
mp,ex,2,2e11
mp,dens,2,8923
mp,nuxy,2,0.3
mp,ex,3,1.5e11
mp,dens,3,7900
mp,nuxy,3,0.3
mp,ex,4,2.1e15
mp,dens,4,0
mp,nuxy,4,0.3
mp,ex,5,3.5e10
mp,dens,5,2600
mp,nuxy,5,0.167
mp,ex,6,2.1e11
mp,dens,6,7800
mp,nuxy,6,0.3
r,1,1.412,3.026,165.937,35.70,3.50
rmore,,10.571
r,2,0.3838
r,3,0.00622
r,4,66.80,553.00,911.00,11.00,11.5,,
rmore,,1190.00
r,5,62.70,507.00,838.00,11.00,11.50,,
rmore,,1090
r,6,58.60,461.00,765,11.00,11.5,,
rmore,,997
r,7,54.5,415,692,11.00,11.5,,
rmore,,899
r,8,50.4,369,619,11.00,11.5,,
rmore,,802
r,9,46,323,546,11.00,11.5,,
rmore,,704
r,10,41,277,474,11.00,11.5,,
rmore,,607
r,11,38,231,401,11.00,11.5,,
rmore,,510
r,12,34,185,328,11.00,11.5,,
rmore,,412
r,13,30,140,255,11.00,11.5,,
rmore,,316
r,14,27,116,219,6,8.5,,
rmore,,266
r,15,32,182,213,7.8,8.5,,
rmore,,318
r,16,31.5,170,212,6,8.5,,
rmore,,244
r,17,0.1288,0.0084,0.0084,2.00,8.50,,
rmore,,0.09911
!连续
n,1,
n,2,25,
n,70,1113
fill,2,70,67
n,71,1138
*do,i,1,71,1
n,100+i,nx(i),125.851-4*113.8*(1138-nx(i))*nx(i)/1138/1138,35.7/2
n,200+i,nx(i),125.851-4*113.8*(1138-nx(i))*nx(i)/1138/1138,-35.7/2
n,300+i,nx(i),,35.7/2
n,400+i,nx(i),,-35.7/2
*enddo
n,701,-330,-33,35.7/2
n,702,-330,-33,-35.7/2
n,703,1138+330,-33,35.7/2
n,704,1138+330,-33,-35.7/2
/view,1,1,1,1
gplot
type,1
mat,1
real,1
*do,i,1,70
e,i,i+1
*enddo
type,2
mat,2
real,2
*do,i,1,70
e,100+i,101+i
e,200+i,201+i
*enddo
e,701,101
e,702,201
e,703,171
e,704,271
type,2
mat,3
real,3
*do,i,2,70
e,100+i,300+i
e,200+i,400+i
*enddo
type,1
mat,4
real,1
*do,i,2,70
e,i,300+i
e,i,400+i
*enddo
gplot
n,501,0,-63.855,35.7/2
n,511,0,0,35.7/2
fill,501,511,9
n,528,0,125.851,35.7/2
fill,511,528,16
n,551,0,-63.855,-35.7/2
n,561,0,0,-35.7/2
fill,551,561,9
n,578,0,125.851,-35.7/2
fill,561,578,16
n,601,1138,-63.855,35.7/2
n,611,1138,0,35.7/2
fill,601,611,9
n,628,1138,125.851,35.7/2
fill,611,628,16
n,651,1138,-63.855,-35.7/2
n,661,1138,0,-35.7/2
fill,651,661,9
n,678,1138,125.851,-35.7/2
fill,661,678,16
*do,i,1,10 !qianhounanta
type,1
mat,5
real,3+i
e,500+i,501+i
*enddo
type,1
mat,5
real,14
*do,i,1,17
e,510+i,511+i
*enddo
*do,i,1,10
type,1
mat,5
real,3+i
e,550+i,551+i
*enddo
type,1
mat,5
real,14
*do,i,1,17
e,560+i,561+i
*enddo
*do,i,1,10 !qianhoubeita
type,1
mat,5
real,3+i
e,600+i,601+i
*enddo
type,1
mat,5
real,14
*do,i,1,17
e,610+i,611+i
*enddo
*do,i,1,10
type,1
mat,5
real,3+i
e,650+i,651+i
*enddo
type,1
mat,5
real,14
*do,i,1,17
e,660+i,661+i
*enddo
nummrg,node,,,,,low
n,500,
n,600,1138
type,1
mat,6
real,15
e,301,500
e,500,401
e,371,600
e,600,471
type,1
mat,6
real,17
e,577,524
e,574,521
e,574,527
e,571,524
e,677,624
e,674,627
e,674,621
e,671,624
type,1
mat,6
real,16
e,201,101
e,271,171
cp,1,uy,1,500
cp,2,uy,71,600
cp,3,uz,1,500
cp,4,uz,71,600
cp,5,rotx,1,500
cp,6,rotx,71,600
d,701,all
d,702,all
d,703,all
d,704,all
d,501,all
d,551,all
d,601,all
d,651,all
gplot
finish
[ 本帖最后由 水幽寒 于 2006-10-25 12:21 编辑 ] 动力分析部分
/solu
antype,2
modopt,lanb,20
mxpand,20,,,0
pstress,1
modopt,lanb,20,0,100,,off
solve
finish
/post1
fini
/config,nres,20000
*dim,aa,TABLE,3000,1 !输入三向天津地震波 (天津医院记录南北东西竖向)
*dim,bb,TABLE,3000,1
*dim,cc,TABLE,3000,1
*vread,aa(1,1),acex,txt,,1 !从文件中读取数值
(e16.6)
*vread,bb(1,1),acey,txt,,1 !从文件中读取数值
(e16.6)
*vread,cc(1,1),acez,txt,,1 !从文件中读取数值
(e16.6)
csys,0
/solu
outr,all,all
nsel,all
anty,trans
*do,tm,0.01,10,0.01
time,tm
acel,aa(tm)/100,bb(tm)/100,cc(tm)/100
lswrite,tm/0.01
*enddo
lssolve,1,1000
fini xuansuoqiao.JPG 哎哟,最后这个模型像是西侯门哟,:lol: 呵呵,这个并不是我做的,贴出来大家学习下
对啊,就是西堠门大桥
对啊,就是西堠门大桥 这个很好,楼主好这个很好,楼主好 楼主能否把西堠门的命令流上传一下 我体会是:悬索桥是一种典型的柔性结构,具有很强的几何非线性,但一般不存在材料非线性问题。ANSYS分析的最大难点就是成桥状态的找形,因为恒载作用下悬索桥的主缆和吊索就存在很大的张力,这种张力本身就是悬索桥的一种致关重要的刚度。在建模时,主缆和吊索的张力基本上要用Link单元的初应变来获得,可是,这个初应变值怎样准确估算,是比较困难的。本人正在研究这个问题,欢迎大家讨论! 这个问题确实很头疼,你在ansys里弄吗? 这个问题搞不懂,学习中~~~~~~~~ 有个3300m的悬索桥,实哪里的?只知道是个海峡大桥,是个小国家的。 我也很想知道这个问题,正为这个问题感到很苦恼 原帖由 yyxxgg000 于 2007-5-11 17:05 发表 http://www.bridgehead.com.cn/images/common/back.gif
有个3300m的悬索桥,实哪里的?只知道是个海峡大桥,是个小国家的。
呵呵,是Stretto di Messina Bridge,位于意大利半岛南端以西的地中海,给你一个效果图看看。
页:
[1]
2