CN111737849A - 一种拉索张拉施工过程中内力及变形跟踪方法 - Google Patents
一种拉索张拉施工过程中内力及变形跟踪方法 Download PDFInfo
- Publication number
- CN111737849A CN111737849A CN202010398077.4A CN202010398077A CN111737849A CN 111737849 A CN111737849 A CN 111737849A CN 202010398077 A CN202010398077 A CN 202010398077A CN 111737849 A CN111737849 A CN 111737849A
- Authority
- CN
- China
- Prior art keywords
- cable
- internal force
- value
- initial strain
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Load-Bearing And Curtain Walls (AREA)
Abstract
本发明属于拉索玻璃幕墙施工技术领域,具体提供了一种拉索张拉施工过程中内力及变形跟踪方法,S01:通过逆迭代法计算出满足预来力所需的初应变,再计算拉索初应变对应的实际内力和实际位移,所述逆迭代法包括改进张力补偿迭代法;S02:上一次迭代求解得到的初应变值之后,将所述初应变值做为下一次循环计算的初应变值;S03:当循环到k次时,若拉索实际内力与预来力的差值小于预设值时,循环结束,并记录上述每次计算得到应力值分别对应每个张拉阶段的拉索的实际内力,其中,k大于1。根据拉索在各阶段的目标索力来计算相对应的拉索内力和变形,为施工提供理论数据和参考依据,提高了拉索玻璃幕墙施工的安全可靠性。
Description
技术领域
本发明属于拉索玻璃幕墙施工技术领域,具体涉及一种拉索张拉施工过程中内力及变形跟踪方法。
背景技术
玻璃幕墙张拉控制施过程中,拉锁本身刚度的形成和大小,完全取决于施工成型时施加给各个索构件预应力的大小、张拉顺序、预应力松弛状况(张拉顺序引起的内力重分布)、补张拉等。拉索结构设计时,往往只给出设计状态下每根拉索的几何位形、应力、变形状态及承载要求,至于如何形成目标位形、应力和变形状态则是未知的,施工前需要进行复杂的计算机仿真仿真分析及迭代计算,以确定每根索的张拉应力和张拉顺序。索结构目标位形或者传力体系的些许改变、索夹具的不同(预应力摩擦损失)、张拉顺序的变化、预应力张拉时环境的不同、以及玻璃安装的顺序等,都会影响预应力索结构的预应力张拉顺序等施工工艺。
有关拉索的张拉控制,在工程实践中是不可回避的重要问题。预应力索的张拉施工,受各种因素的影响极大,其结果也有较大差别。考虑到拉索张拉工作量极大,施工条件所限,一般采取分批张拉的方法。预应力张拉的批次、预应力大小和张拉顺序对于预应力索的最终几何形状影响极大,因此,预应力分批、分级张拉是拉索施工成形的关键。
传统分析中,拉索玻璃幕墙应先进行施工找形分析,确定拉索的初始几何状态,然后再确定拉索初始状态的放样几何和内力分布。实际工程中大部分结构在进行找力分析后便进行张拉分析确定施工控制力进行施工。李咏梅针对逐环张拉成型的施工方法,利用索力为控制参数,编制了一种可以考虑施工进程的施工分析程序,该程序可以考虑几何非线性、后批预应力拉索对前批预应力拉索造成的损失;可以获得初始索力和位型。然而,该方法只是针对逐环张拉成型的无脚手架施工方法而提出的,且此种无脚手架施工方法本身只是用于中小跨度,对于大跨度结构不适用。除此之外,该程序只采用索力为控制参数,没有考虑位形控制参数。大跨度结构张拉施工过程中,位形的变化相对较大,采用双控进行施工控制是必要的。广义的施工找形分析应包括找力和找形,找力分析是确定对应预应力设计值的初应变值,找形分析是获得满足预应力设计态的放样态位形,二者是缺一不可,且密切联系,因此,在对结构进行施工前必须进行施工找形分析。
发明内容
本发明的目的是克服现有技术中拉锁张拉施工过程中的内力及变形判断误差大的问题。
为此,本发明提供了一种拉索张拉施工过程中内力及变形跟踪方法,包括以下步骤:
S01:通过逆迭代法计算出满足预来力所需的初应变,再计算拉索初应变对应的实际内力和实际位移,所述逆迭代法包括改进张力补偿迭代法;
S02:上一次迭代求解得到的初应变值之后,将所述初应变值做为下一次循环计算的初应变值;
S03:当循环到k次时,若拉索实际内力与预来力的差值小于预设值时,循环结束,并记录上述每次计算得到应力值分别对应每个张拉阶段的拉索的实际内力,其中,k大于1。
优选地,所述预来力的大小为拉索在实际工况下的目标索力。
优选地,通过获取施工最终状态下拉索的变形,然后通过变形计算得到所述目标索力。
优选地,所述改进张力补偿迭代法具体包括:
假设有n组索,张力设计值分别P1、P2、P3,…,Pn,进行k次循环。
各计算步骤为:
......
......
其中,k--循环序号;i--拉索组号;j--张拉批次号;Pi--第i组索的张力设计值;εi(k)--第k次循环计算中,由第i组索张力设计值计算出的初应应变;Pi(k)--第k次循环计算中,第i组索的张力;一第k次循环计算中,第i组索在第j批次张拉时的实际内力。
对于初应变εi(k),k=1时,由Pi计算得到值;当k≠1时,由Pi(k)计算得到。
由此,各组拉索均施加了初应变,各索内力变化如:
......
优选地,如果计算误差偏大,无法满足施工要求,则进行再一次循环计算。通过将内力变化值补偿给上一次循环时索的张力值,再换算成应变值,从而修改在第k+1次循环计算前各组索的张力控制值;
......
优选地,还包括找形分析,所述找形分析采用环向索的施工找形算法,具体地,根据设定的拉索张力值,计算拉索初始长度,获取拉索当前应力密度关系,采用逐渐迭代力密度法对拉索进行找形分析,在每次拉锁张拉迭代阶段过程中跟踪拉索实际变形,若拉索实际变形与拉索初始长度的差值小于预设值时,循环结束,并记录上述每次计算得到拉索实际变形。
优选地,计算分析用工具载体包括ANSYS或MATLAB。
本发明的有益效果:本发明提供的这种拉索张拉施工过程中内力及变形跟踪方法,S01:通过逆迭代法计算出满足预来力所需的初应变,再计算拉索初应变对应的实际内力和实际位移,所述逆迭代法包括改进张力补偿迭代法;S02:上一次迭代求解得到的初应变值之后,将所述初应变值做为下一次循环计算的初应变值;S03:当循环到k次时,若拉索实际内力与预来力的差值小于预设值时,循环结束,并记录上述每次计算得到应力值分别对应每个张拉阶段的拉索的实际内力,其中,k大于1。根据拉索在各阶段的目标索力来计算相对应的拉索内力和位移,找出结构在预应力施工过程中的不利状态,为施工提供理论数据和参考依据,提高了拉索玻璃幕墙施工的安全可靠性。
以下将结合附图对本发明做进一步详细说明。
附图说明
图1是本发明拉索张拉施工过程中内力及变形跟踪方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“中心”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征;在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
单向双索玻璃体系主要由三部分构成玻璃、连接件和索,也就是所谓的所玻璃面板、锚固结构和索结构。
1、索结构,是由拉索作为主要受力构件形成的预应力结构体系。由其承受玻璃幕墙的荷载,并将这些荷载传递给锚固结构。作为跨越幕墙支撑跨度的重要构建之一,索结构要悬挂于锚固结构上。
2、一锚固结构,就是指对锚杆、锚索等进行支撑的结构,它承受着索结构传来的荷载,由其传给上部主体结构和下部的地基基础。
3、玻璃面板,由安装在拉索上的连接件固定在索结构上,玻璃之间通过结构胶进行填充。
由此可见,以上三者是互相制约、相互依存、互相影响的,从而共同形成幕墙系统。这些结构能够承受相应承载和变形。
本发明实施例提供了一种拉索张拉施工过程中内力及变形跟踪方法,如图1所示,包括以下步骤:
S01:通过逆迭代法计算出满足预来力所需的初应变,再计算拉索初应变对应的实际内力和实际位移,所述逆迭代法包括改进张力补偿迭代法;
S02:上一次迭代求解得到的初应变值之后,将所述初应变值做为下一次循环计算的初应变值;
S03:当循环到k次时,若拉索实际内力与预来力的差值小于预设值时,循环结束,并记录上述每次计算得到应力值分别对应每个张拉阶段的拉索的实际内力,其中,k大于1。
在推导拉索的计算理论时,通常采用三个基本假设:(1)拉索是柔性结构,不能承受弯、压作用;(2)胡克定律适用于拉索材料;(3)外部载荷作用在节点上,索段呈直线状。公式推导如下:
方程式(1)和(2)就是单索问题的基本平衡微分方程。
由方程式(1)得
方程式(2)可以改写成
方程式(4)的物理意义就是:索曲线在某一点二阶导数与作用在该点的横向荷载集度成正比。
当qx沿跨度均布,,则qx=q(常量),公式(4)可改写成
当y=0时x=0
当y=L时x=C
由此可得
C2=0
拉索幕墙张力受力的特点:一是张拉过程中幕墙的受力状态、边界条件和几何形态是变化的;二是不同的张拉顺序和方法对张拉结构的受力状态很较大影响。因此,拉索幕墙的张拉模拟应该如实反映结构位移、边界条件等的变化,如实反映不同的施工顺序和方法的影响。
为此,数值模型应该确保幕墙构件的材料、几何尺寸与实际工程相同,预应力施加与实际张拉效果相同,在各阶段施工的荷载、边界条件等与施工施工状态相同。张拉时的数值模拟与使用阶段的数值模拟存在一些差别。
由于张拉力与模型的初应变不是一般的线性关系,应该通过逆迭代法计算出满足预来力所需的初应变,在计算拉索初应变对应的内力和位移。对于张拉方式采用拉索的张拉初应变求解可以采用改进张力补偿法进行迭代求解。接下来先找力,具体原理过程如下:
假设有n组索,张力设计值为君,P1、P2、P3,…,Pn,进行k次循环。
各计算步骤为:
......
……
其中,k--循环序号;i--拉索组号;j--张拉批次号;Pi--第i组索的张力设计值;εi(k)--第k次循环计算中,由第i组索张力设计值计算出的初应应变;Pi(k)--第k次循环计算中,第i组索的张力;一第k次循环计算中,第i组索在第j批次张拉时的实际内力。
对于初应变εi(k),k=1时,由Pi计算得到值;当k≠1时,由Pi(k)计算得到。
由此,各组拉索均施加了初应变,各索内力变化如:
......
在第k次循环完成后,如果 就是计算误差非常小时,循环计算即可结束。此时, 近似等于各索的设计张力,那么 就是ANYSY程序中第1,2,3......,n组索应施加的初应变。如果计算误差偏大,无法满足施工要求,则进行再一次循环计算。通过将内力变化值补偿给上一次循环时索的张力值,再换算成应变值,从而修改在第k+1次循环计算前各组索的张力控制值。
......
采用张力补偿法求解张拉控制力需要的初应变之后,分析,即跟踪张拉过程中每个张拉阶段的内力、变形。
接下来通过迭代算法进行找形,具体步骤如下:
在已知的初始几何形态和张拉力的条件下,计算结构初始状态下的内力和放样几何。它假设拉索幕墙的放样几何为{X}i,以一对方向相反、大小相同的力Tp代替张拉索,计算拉索预拉力Tp作用下结构的位形及其与设计几何{X}0的差值{ΔX}i。
{X}i+1={X}i-{ΔX}i (10)
并继续进行迭代过程。
上述的逆迭代法存在两个不足:一是该方法可行的前提是预拉力大小Tp必须事先确定,然而在找形分析之前,结构预拉力大小往往是未知的,需要在找形分析中确定;二是用一对大小相等、方向相反的拉力对结构施加预拉力,当结构变形较大时,这种模拟方法不够精确。
本文采用逐环张拉环索施工顺序的施工找形算法,并将其应用在拉索玻璃幕墙施工分析中。
该方法具有以下几个特点:
(1)索的力学计算建立在三个基本假设之上:①拉索是柔性结构,不能承受弯、压作用。②胡克定律适用于拉索材料。③外部荷载作用在节点上,索段呈直线状。
(2)根据拉索玻璃幕墙施工进程中的几何条件、荷载条件和边界约束条件的不同,分析了拉索玻璃幕墙施工中拉索张拉受力机理。
(3)对同一个拉索玻璃幕墙的初始状态施工找形分析,应包括找力和找形两个方面,找形主要的主要目的是计算放样几何。本文结合逆迭代法选用适用于本过程拉索的双控法施工找形方法。
以上例举仅仅是对本发明的举例说明,并不构成对本发明的保护范围的限制,凡是与本发明相同或相似的设计均属于本发明的保护范围之内。
Claims (7)
1.一种拉索张拉施工过程中内力及变形跟踪方法,其特征在于,包括以下步骤:
S01:通过逆迭代法计算出满足预来力所需的初应变,再计算拉索初应变对应的实际内力和实际位移,所述逆迭代法包括改进张力补偿迭代法;
S02:上一次迭代求解得到的初应变值之后,将所述初应变值做为下一次循环计算的初应变值;
S03:当循环到k次时,若拉索实际内力与预来力的差值小于预设值时,循环结束,并记录上述每次计算得到应力值分别对应每个张拉阶段的拉索的实际内力,其中,k大于1。
2.根据权利要求1所述的拉索张拉施工过程中内力及变形跟踪方法,其特征在于,所述预来力的大小为拉索在实际工况下的目标索力。
3.根据权利要求2所述的拉索张拉施工过程中内力及变形跟踪方法,其特征在于:通过获取施工最终状态下拉索的变形,然后通过变形计算得到所述目标索力。
4.根据权利要求1所述的拉索张拉施工过程中内力及变形跟踪方法,其特征在于,所述改进张力补偿迭代法具体包括:
假设有n组索,张力设计值分别P1、P2、P3,…,Pn,进行k次循环。
各计算步骤为:
......
......
其中,k--循环序号;i--拉索组号;j--张拉批次号;Pi--第i组索的张力设计值;εi(k)--第k次循环计算中,由第i组索张力设计值计算出的初应应变;Pi(k)--第k次循环计算中,第i组索的张力;一第k次循环计算中,第i组索在第j批次张拉时的实际内力。
对于初应变εi(k),k=1时,由Pi计算得到值;当k≠1时,由Pi(k)计算得到。
由此,各组拉索均施加了初应变,各索内力变化如:
......
6.根据权利要求1所述的拉索张拉施工过程中内力及变形跟踪方法,其特征在于,还包括找形分析,所述找形分析采用环向索的施工找形算法,具体地,根据设定的拉索张力值,计算拉索初始长度,获取拉索当前应力密度关系,采用逐渐迭代力密度法对拉索进行找形分析,在每次拉锁张拉迭代阶段过程中跟踪拉索实际变形,若拉索实际变形与拉索初始长度的差值小于预设值时,循环结束,并记录上述每次计算得到拉索实际变形。
7.根据权利要求1或6所述的拉索张拉施工过程中内力及变形跟踪方法,其特征在于:计算分析用工具载体包括ANSYS或MATLAB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010398077.4A CN111737849B (zh) | 2020-05-12 | 2020-05-12 | 一种拉索张拉施工过程中内力及变形跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010398077.4A CN111737849B (zh) | 2020-05-12 | 2020-05-12 | 一种拉索张拉施工过程中内力及变形跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111737849A true CN111737849A (zh) | 2020-10-02 |
CN111737849B CN111737849B (zh) | 2022-07-22 |
Family
ID=72647176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010398077.4A Active CN111737849B (zh) | 2020-05-12 | 2020-05-12 | 一种拉索张拉施工过程中内力及变形跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111737849B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120032879A1 (en) * | 2003-03-26 | 2012-02-09 | Mimic Technologies, Inc. | Method, Apparatus, and Article for Force Feedback Based on Tension Control and Tracking Through Cables |
CN106599509A (zh) * | 2016-12-26 | 2017-04-26 | 河北建筑工程学院 | 一种模拟弦支穹顶结构施加预应力的方法 |
CN108345729A (zh) * | 2018-01-26 | 2018-07-31 | 东南大学 | 一种基于群论的对称索杆结构找形方法 |
-
2020
- 2020-05-12 CN CN202010398077.4A patent/CN111737849B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120032879A1 (en) * | 2003-03-26 | 2012-02-09 | Mimic Technologies, Inc. | Method, Apparatus, and Article for Force Feedback Based on Tension Control and Tracking Through Cables |
CN106599509A (zh) * | 2016-12-26 | 2017-04-26 | 河北建筑工程学院 | 一种模拟弦支穹顶结构施加预应力的方法 |
CN108345729A (zh) * | 2018-01-26 | 2018-07-31 | 东南大学 | 一种基于群论的对称索杆结构找形方法 |
Non-Patent Citations (4)
Title |
---|
张海禄: "独塔斜拉桥钢绞线拉索施工过程中索力控制方法研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
王永泉等: "空间预应力钢结构拉索等效预张力确定方法研究", 《土木工程学报》 * |
王鹏 等: "自锚式悬索桥施工控制", 《结 构 工 程 师》 * |
赵眈崴等: "张弦梁结构张拉阶段受力分析及施工监测", 《建筑结构》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111737849B (zh) | 2022-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wakefield | Engineering analysis of tension structures: theory and practice | |
Symonds et al. | Recent progress in the plastic methods of structural analysis | |
CN106096257B (zh) | 一种非线性索单元分析方法及系统 | |
CN102622485B (zh) | 确定索杆梁系空间结构放样态的形态分析逆迭代法 | |
Xue et al. | Design optimization and experimental study on beam string structures | |
CN110700072B (zh) | 一种斜拉桥缩尺试验模型斜拉索的安装方法 | |
CN102789525B (zh) | 高层建筑结构非弹性时变分析方法 | |
CN107622174B (zh) | 基于影响矩阵修正的斜拉桥拉索张拉施工闭环控制方法 | |
Mesnil et al. | Stability of pseudo-funicular elastic grid shells | |
CN106599509A (zh) | 一种模拟弦支穹顶结构施加预应力的方法 | |
CN111737849B (zh) | 一种拉索张拉施工过程中内力及变形跟踪方法 | |
CN110387820B (zh) | 一种基于智能张拉系统的吊杆测控一体化施工方法 | |
Spiliopoulos et al. | An efficient mathematical programming method for the elastoplastic analysis of frames | |
CN102707623A (zh) | 一种预应力网格结构张拉全过程的反馈控制方法 | |
CN102733525B (zh) | 一种基于承载全过程开口式整体张拉结构多阶段设计方法 | |
Li et al. | Determination of reasonable finished state of self-anchored suspension bridges | |
Chen et al. | Practical second-order inelastic analysis for three-dimensional steel frames | |
Tian et al. | Nonlinear time-varying analysis algorithms for modeling the behavior of complex rigid long-span steel structures during construction processes | |
CN108319800A (zh) | 一种基于功能目标的张力结构通用体型设计方法 | |
CN111723418B (zh) | 预应力张拉控制方式的精确性判别方法 | |
Attiyah et al. | Finite element modelling of concrete shrinkage cracking in walls | |
Liu et al. | Cable tension preslack method construction simulation and engineering application for a prestressed suspended dome | |
CN112035928A (zh) | 大跨径悬索管道桥成桥主索线形计算方法 | |
Chan et al. | Nonlinear analysis of pre-tensioned glass wall facade by stability function with initial imperfection | |
CN112800538B (zh) | 一种任意薄壁梁结构的层合化分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |