没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
索杆张力结构施工成形分析需要解决一个松弛态索杆系统的平衡形态求解问题,该问题可归结为一个给定构件原长的受荷索杆机构系统的找形问题。文中利用动力松弛法进行该类松弛索杆体系的找形分析。由于不建立刚度矩阵,避免了体系几何不稳定性所引起的刚度矩阵奇异问题。该方法采用悬链线索单元模型,可以考虑索大垂度效应。文中还推导了反映索原长和内力之间关系的悬链线索单元协调方程。通过一个正交松弛索网算例分析,考察了该找形方法的计算精度和收敛性。最后还模拟了一个索穹顶的施工成形过程,表明了该找形方法用于索杆张力结构施工成形分析的有
资源推荐
资源详情
资源评论
收稿日期 :2006‐05‐06 ;修改稿收到日期 :2007‐04‐03畅
基金项目 :国家自然科学基金(50578139)资助项目 畅
作者简介 :伍晓顺(1981‐) ,男 ,工学硕士 ;
邓 华
倡
(1971‐) ,男 ,工学博士 ,教授
(E‐mail :denghua@ zju .edu .cn)畅
第25卷第2期
2008 年 4 月
计 算 力 学 学 报
Chinese Journal of Computational Mechanics
Vol .25 ,No .2
April 2008
文章编号 :1007‐4708(2008)02‐0229‐08
基 于 动 力 松 弛 法 的 松 弛 索 杆 体 系 找 形 分 析
伍晓顺 , 邓 华
倡
(浙江大学 空间结构研究中心 ,杭州 310027)
摘 要 :索杆张力结构施工成形分析需要解决一个松弛态索杆系统的平衡形态求解问题 ,该问题可归结为一个
给定构件原长的受荷索杆机构系统的找形问题 。 文中利用动力松弛法进行该类松弛索杆体系的找形分析 。 由
于不建立刚度矩阵 ,避免了体系几何不稳定性所引起的刚度矩阵奇异问题 。 该方法采用悬链线索单元模型 ,可
以考虑索大垂度效应 。 文中还推导了反映索原长和内力之间关系的悬链线索单元协调方程 。 通过一个正交松
弛索网算例分析 ,考察了该找形方法的计算精度和收敛性 。 最后还模拟了一个索穹顶的施工成形过程 ,表明了
该找形方法用于索杆张力结构施工成形分析的有效性 。
关键词 :索杆张力结构 ;找形 ;施工模拟 ;悬链线索单元 ;动力松弛法
中图分类号 :TU393 .3 ;T U311 .3 文献标识码 :A
1 引 言
近些年来 ,以索穹顶结构
[1]
为代表的一类索杆
张力结构体系
[2]
受到结构工程界的关注并成为当
前空间结构的研究热点 。 索杆张力结构最主要的
特点是该类体系在几何上并不满足常规结构构成
准则(如 Maxwell 准则
[3]
) ,表观上属于“机构”系
统 。 但是索杆张力结构在合理预应力作用下 ,体系
可以稳定地承受外部荷载
[4 ,5]
,并且比常规结构具
备跨越更大空间的能力 。 索杆张力结构一个值得
关注的特点是其施工方法 ,无论是索穹顶还是一些
新型索杆张力结构 ,结构的施工通常借助于其自身
构件中的索单元 ,通过对这些索构件进行张拉以提
升整体结构 ,从而达到结构的最终成形 ,并且将预
应力施加到位
[1 ,6]
。
索杆张力结构的成形分析面临的是一个松弛
状态的强几何不稳定(可变)索杆系统 。 在提升过
程中 ,体系表现为一个大变形的性态 。 早期关于索
杆 张 力 结 构 施 工 成 形 的 研 究 基 本 针 对 于 索 穹
顶
[7 ,8]
,一般利用轴对称条件将空间问题简化为平
面问题 ,在结构分析时通过剔除直观的刚体位移以
克服刚度矩阵奇异 ,使得有限单元法得以实施 。 近
期的研究主要有两种思路 :一种是基于机构几何
运动学方法
[9]
建立索杆体系的运动解析 ,但通常忽
视荷载对形态的控制作用 ;一种是将运动驱动(张
拉索长度的缩短)独立于形态变量
[10]
,即将成形过
程看成由一系列离散的受荷机构平衡形态组成 。
从当前状况来看 ,基于第二种思路的研究已经取得
较为有效的进展 。
对于第二种思路 ,既然索杆张力结构的成形过
程可以描述为体系在张拉索长度变化下对应于一
系列施工步的平衡形态求解问题 ,那么对于每一个
特定状态(施工步) ,可以归结为一个已知构件原长
的松弛索杆体系在特定荷载(如自重)作用下的找
形问题 。 文献[10]中提出了一种基于“力密度”法
的求解策略 ,但当初始迭代形状与最终平衡形态偏
离较大时 ,基于力密度法的求解策略需要限制节点
位移增量的步长来保证计算的收敛性 ,而步长的减
小导致迭代次数的增加和计算效率的降低 。 同时
由于采用抛物线索单元 ,对于大垂度索单元将产生
较大误差 。 本文提出了一种基于动力松弛法的受
荷松弛索杆体系的找形策略 ,同时采用了精确的悬
链线单元以便考虑松弛态索段易出现大垂度的情
况 。 文中还引入了两个算例 ,从数值上考察该计算
方法的正确性和有效性 。
2 找形问题的特点和数学描述
从工程的角度来看 ,索杆张力结构的施工成形
分析需要求解一系列特定施工步骤或时刻下 ,体系
的形状以及相应的构件内力 。 施工过程体系形态
变化来自于运动的驱动 ,即张拉某些拉索使得其长
度缩短 ,因此运动驱动变量在数学上可以描述为张
拉索长度的变化 。 进而 ,可以将施工成形过程分析
落实在一系列特定驱动变量离散点(特定张拉索长
度)上的形态分析问题 。
应该指出的是 ,对于某个特定施工步骤或时
刻 ,构件的原长(即无伸长长度 ,或放样长度)是已
知的 ,即
s
0
=
s
1
-
Δ s
+
s
t
(1)
式中 s
1
为初始平衡态 (成形后的预应力状态) 的
几何长度 ,Δ s 为初始平衡态构件的弹性伸长量(对
于索) 或者缩短量(对于杆) ;s
t
为张拉索的牵引长
度 ,对于非张拉构件为零 。张拉索的 s
t
在张拉过程
中是变化的 ,其值需不断扣除千斤顶已拔出长度 。
对于某个特定时刻 ,体系形态分析实际上是一
些已知原长的构件根据特定的连接方式在其自重
作用下所达到的平衡形态求解问题(即找形问题) 。
因此对于具体的施工步 ,基本找形问题的数学描述
如下 :
(1) 已知条件 :该时刻各索段和杆的原长 s
0
,
构件的截面面积 A 和弹性模量 E ,构件之间的连接
关系(即拓扑关系) 以及边界约束情况 ;索上横向
荷载
q
(通常仅为索的自重) 和节点荷载
p
(包括实
际节点重量以及可能的外挂荷载) 。
(2) 求解内容 :在当前结构自重作用下的平衡
状态 ,具体包括结构形状(以节点坐标描述) 和构
件内力 。
3 单元分析
3 .1 基本假定
(1) 索为理想柔索 ,不受压且无弯曲刚度 。
(2) 杆为两节点单元 ,杆自重等效为两端节点
荷载 。
(3) 索上只有沿索长均布的竖向荷载(包括自
重)。
(4)构件材料符合虎克定律 ,且满足小应变假定 。
3 .2 悬链线索单元分析
图 1 所示空间索段 ,定义其局部坐标系 o
~
‐
x
~
z
~
,
其中 x
~
轴平行于索段在整体坐标系 o
‐
x
y
平面上的
投影线 ,z
~
轴和整体坐标系的 z 轴与索上均布荷载
同向 。根据悬索结构的基本理论
[11]
,索的平衡曲线
为悬链线 ,曲线方程为
z
~
(x) =
H
q
cosh
α
-
cosh
q
x
~
H
-
α
图 1 悬链线索单元
Fig .1 Catenary cable element
α
=
sinh
-
1
q
c
2 Hsinh
q
L
2 H
+
q
L
2 H
(2)
式中 L 是悬索两端节点的水平距离 ,c 是悬索两端
节点的竖向高差 ,
q
为沿索长的均布荷载 ,H 是悬
索张力的水平分量 ,为常数 。
将式(2)对 x
~
进行求导 ,可得悬链线的斜率方程
为
z
~
′
(x) = sinh
α
-
q
x
~
H
(3)
当 x
~
=
0 时 ,有 z
~
′
(0) = sinh
α
。索端张力 T
i
,T
j
与
z
~
′
(x
~
) 和 H 的关系式为
T
i
=
H 1
+
(z
~
′
(0))
2
, T
j
=
H 1
+
(z
~
′
(L ))
2
(4)
在局部坐标系 o
~
‐
x
~
z
~
下 ,索端张力在各坐标轴上的
分量为
F
i
x
~
= -
H , F
i
z
~
= -
H
·
z
~
′
(0)
F
j
x
~
=
H , F
j
z
~
= -
q
·
s
0
-
F
i
z
~
(5)
在整体坐标系下 ,索端张力对索两端 i 和
j
节
点产生的节点力分别为
F
i
x
=
H(x
j
-
x
i
)/L , F
i
y
=
H (
y
j
-
y
i
)/L
F
i
z
=
H z
~
′
(0) (6)
F
j
x
= -
H(x
j
-
x
i
)/L , F
j
y
= -
H (
y
j
- y
i
)/L
F
j
z
=
q
·
s
0
+
F
i
z
~
(7)
式中 (x
i
,
y
i
,z
i
) 和(x
j
,
y
j
,z
j
) 分别是单元两端节
点在整体坐标系下的坐标 。
根据几何条件 ,变形后索长 s 可按下式计算 :
s
=
∫
L
0
1
+
(z
~
′
(x))
2
dx
=
H
q
(sinh
α
-
sinh
β
)
β
=
sinh
-
1
q
c
2 Hsinh
q
L
2 H
-
q
L
2 H
(8)
引入小应变假定 ,索由于张力引起的变形伸长
量 Δs 可由下式求得
Δ s
=
H
EA
∫
L
0
ds
d x
~
2
dx
~
=
H
EA
∫
L
0
(1
+
z
~
′
2
(x
~
))dx
~
=
H
4
q
EA
(2 L
q
+
Hsinh(2
α
) - Hsinh(2
β
)) (9)
032
计 算 力 学 学 报
第 25 卷
剩余7页未读,继续阅读
资源评论
weixin_38678057
- 粉丝: 6
- 资源: 870
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功