CN115130350B - 一种连续索的非线性滑移计算方法 - Google Patents
一种连续索的非线性滑移计算方法 Download PDFInfo
- Publication number
- CN115130350B CN115130350B CN202210790527.3A CN202210790527A CN115130350B CN 115130350 B CN115130350 B CN 115130350B CN 202210790527 A CN202210790527 A CN 202210790527A CN 115130350 B CN115130350 B CN 115130350B
- Authority
- CN
- China
- Prior art keywords
- cable
- pulley
- sides
- force difference
- rope
- 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.)
- Active
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供的一种连续索的非线性滑移计算方法,综合利用索的无应力总长度不变和索滑移刚度影响矩阵原理,通过几何非线性求解滑轮两侧索段的力差向量,并根据滑轮两侧的力差向量反复迭代模拟索滑移过程,以此建立连续非线性滑移索单元计算,根据力差向量、滑移刚度影响矩阵和变形向量的关系,求解各滑轮的变形向量,以此更新各滑轮两侧索段的无应力长度,再重新计算滑轮两侧力差向量,并结合滑移刚度影响矩阵,重新调整滑移点两侧索无应力长度情况,以此循环直至滑轮两侧力差向量收敛。本方法考虑了结构的几何非线性问题和拉索垂度问题,适应于多种工况,如多跨连续索和大垂度的结构计算,更贴近工程实际施工,且收敛更容易,计算效率更高。
Description
技术领域
本发明属于桥梁动力学技术领域,涉及一种连续索的非线性滑移计算方法。
背景技术
索作为空间结构重要的组成部分之一,其小自重、高效率、大跨度、省成本的优点,使其在张弦梁、缆索系统、弦支穹顶等轻型高效结构中占有重要的地位,广泛应用于大型体育馆屋面结构、斜拉桥、风景区游览索道以及拉索塔架等。连续索在节点处的连接有两种形式:一种是索和节点固定,节点两侧的索段在节点处不发生滑移;另一种是索在节点处采用穿管、滑轮等连接方式,索和节点之间能发生相对滑移。第二种连接形式称为接触连接,对应的节点称为接触点。第二种索连接形式无论是在张拉成型的施工阶段或施工完成后的使用阶段,始终无法避免索的滑移运动问题及在滑移过程中存在的多种滑移状态难题。
目前,针对滑移索单元的研究主要有以下几种方法:1)推导单个索段的滑移刚度来考虑索段滑移的影响,但需求解复杂的滑移刚度影响矩阵,且只对中间节点进行了求解,不便于扩展到多跨连续索结构中;2)利用线弹性理论推导单折点力筋单元,但是忽略了拉索的垂度,且将非线性问题线性化,使用范围受到限制,也不便于结构的受力全过程分析;3)滑移索系结构的冷冻升温法通过现有的有限元工具,通过虚加温度荷载的办法,来解决在加载过程中索的滑移问题,但是收敛速度不易控制;4)还有的索滑移的算法在求解各索段的原长时采用了不考虑弹性的索长计算式,对水平投影和竖直投影的修正也是非必要的;5)还有引入弯矩分配法概念,但采用的是二节点直线杆单元来代替索单元,不适用于大垂度索结构;6)还有用弹性悬链线单元模拟索段推导索在支撑点处的单侧滑移刚度,通过有限元迭代调整支撑点两侧的索原长,使索力相等,为了模拟索的滑移,采用通用有限元软件中的接触单元,但这种方法非常复杂且工作量很大。
发明内容
本发明的目的在于克服现有技术中在考虑索段滑移的影响的计算时,存在设计方法复杂、适用性窄、计算环节不易收敛等上述不足,提供一种连续索的非线性滑移计算方法。
为了实现上述发明目的,本发明提供了以下技术方案:
一种连续索的非线性滑移计算方法,包括如下步骤:
a、根据索结构的几何参数、材料参数和边界条件,将索与所有滑轮的连接均视为接触连接建立有限元模型,以索分别在1#,2#,...,i#滑轮滑移单位长度的情况下求解1#,2#,...,i#滑轮两侧索段的力差向量,获得滑移刚度影响矩阵:
b、结合1#,2#,...,i+1#索段的无应力长度,对索结构进行几何非线性分析,获取1#,2#,...,i+1#索段的力差向量得到第一力差向量其中,1#,2#,...,i+1#索段的初始无应力长度按照跨径比例分配,1#滑轮两侧分别为1#索段和2#索段,t1表示1#滑轮两侧1#索段和2#索段的力差向量,以此类推;
c、根据公式Δlslide=Kslide -1·ΔT求解1#,2#,...,i#滑轮的滑移长度向量Δlslide=[Δl1 Δl2 ... Δli]T,修正各索段的无应力长度并重复步骤a-b得到第二力差向量
d、重复步骤a-c,直至||(ΔT(n)-ΔT(n-1))/ΔT(n-1)||≤Δε,n表示迭代次数,提取对应有限元模型中的索的位移和各断面的内力,完成滑移计算。
优选的,Δε=0.0001-0.0005。
优选的,索结构的几何参数包含索的初始拟定截面尺寸、索的截面形状、各索段的跨径和各滑轮之间的高差,材料参数包含索的弹性模量、容重、抗拉强度和最小破断拉力。
优选的,几何非线性分析方法为牛顿-拉夫森迭代方法。
与现有技术相比,本发明的有益效果:综合利用索的无应力总长度不变和索滑移刚度影响矩阵原理,通过几何非线性求解滑轮两侧索段的力差向量,并根据滑轮两侧的力差向量反复迭代模拟索滑移过程,以此建立连续非线性滑移索单元计算,根据力差向量、滑移刚度影响矩阵和变形向量的关系,求解各滑轮的变形向量,以此更新各滑轮两侧索段的无应力长度,再重新计算滑轮两侧力差向量,并结合滑移刚度影响矩阵,重新调整滑移点两侧索无应力长度情况,以此循环直至滑轮两侧力差向量收敛。本方法在进行有限元建模时考虑了结构的几何非线性问题和拉索垂度问题,适应于多种工况和不同结构形式连续索计算,相较于现有技术仅适应于单跨拱桥,本方法可适应于多跨连续索和大垂度的结构计算,也可以适合穹顶、悬索桥和拱桥等的计算,而且可以同时考虑主索、牵引索、起重索各索组合作用下结构受力,更贴近工程实际施工,且对无应力长度的初始值依赖性小,无应力索长初始值只要不为负值就都能收敛,收敛更容易,计算效率更高。
附图说明
图1是本发明的一种连续索的非线性滑移计算方法的流程图;
图2是实施例1中的索鞍结构的结构示意图。
具体实施方式
下面结合实施例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
实施例1
一种连续索的非线性滑移计算方法,如图1所示,包括如下步骤:
a、根据索结构的几何参数、材料参数和边界条件,将索与所有滑轮的连接均视为接触连接建立有限元模型,以索分别在1#,2#,...,i#滑轮滑移单位长度的情况下求解1#,2#,...,i#滑轮两侧索段的力差向量,获得滑移刚度影响矩阵:
b、结合1#,2#,...,i+1#索段的无应力长度,对索结构进行几何非线性分析,获取1#,2#,...,i+1#索段的力差向量得到第一力差向量其中,1#,2#,...,i+1#索段的初始无应力长度按照跨径比例分配,1#滑轮两侧分别为1#索段和2#索段,t1表示1#滑轮两侧1#索段和2#索段的力差向量,以此类推;
c、根据公式Δlslide=Kslide -1·ΔT求解1#,2#,...,i#滑轮的滑移长度向量Δlslide=[Δl1 Δl2 ... Δli]T,修正各索段的无应力长度并重复步骤a-b得到第二力差向量
d、重复步骤a-c,直至||(ΔT(n)-ΔT(n-1))/ΔT(n-1)||≤Δε,n表示迭代次数,提取对应有限元模型中的索的位移和各断面的内力,完成滑移计算。
具体的,以图2的索鞍结构为例,图中共设有1#、2#、3#和4#滑轮,1#滑轮两侧分别为1#和2#索段,2#和3#滑轮之间为3#索段,4#滑轮两侧分别是3#和4#索段,根据索结构的几何参数包含索的初始拟定截面尺寸、索的截面形状、各索段的跨径和各滑轮之间的高差,材料参数包含索的弹性模量、容重、抗拉强度和最小破断拉力,将索与所有滑轮的连接均视为接触连接建立有限元模型,设定1#滑轮的滑移单位长度,获取1#滑轮、2#滑轮、...、i#滑轮对应两侧索段的索力差分别为k1,1、k2,1、...、ki,1,以各滑轮右侧索力减去左侧索力为索力差,比如1#滑轮在滑移单位长度时,1#滑轮右侧索力减去左侧索力的值为k1,1,2#滑轮右侧索力减去左侧索力的值为k2,1,以此类推得到相邻索段的力差向量{Ki1},即{Ki,1}=[k1,1 k2,1 ...ki,1]T,同理i#滑轮滑移单位长度时,各滑轮的相邻索段的力差向量为{Ki,i},即{Ki,i}=[k1,jk2,j ... ki,i]T,以此获得滑移刚度影响矩阵:
各索段初始无应力长度按照跨径比例分配,如图中的1#索段的无应力长度2#索段的无应力长度/>以此类推,S为索的总长,然后进行几何非线性分析,如采用牛顿-拉夫森迭代算法或弧长法。为方便计算,以滑轮右侧索力减去左侧索力大于零为正、小于零为负,如1#滑轮右侧索力减去左侧索力为下标1为1#滑轮,上标1为首次迭代,2#滑轮对应的索力差为/>...,i#滑轮对应的索力差为/>以此形成的第一力差向量/>
若||(ΔT(2)-ΔT(1))/ΔT(1)||≤Δε,取Δε=0.0001-0.0005,则提取当前有限元模型中的索的位移和各断面的内力,完成滑移计算。
若||(ΔT(2)-ΔT(1))/ΔT(1)||>Δε,则重复步骤a-c,获取并得到ΔT(3),以此类推,直至||(ΔT(n)-ΔT(n-1))/ΔT(n-1)||≤Δε达到收敛,进而获取收敛模型中的索的位移和各断面的内力。
以上仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以的权利要求的保护范围为准。
Claims (4)
1.一种连续索的非线性滑移计算方法,其特征在于,包括如下步骤:
a、根据索结构的几何参数、材料参数和边界条件,将索与所有滑轮的连接均视为接触连接建立有限元模型,以索分别在1#,2#,…,i#滑轮滑移单位长度的情况下求解1#,2#,…,i#滑轮两侧索段的力差向量,获得滑移刚度影响矩阵:
b、结合1#,2#,…,i+1#索段的无应力长度,对索结构进行几何非线性分析,获取1#,2#,…,i+1#索段的力差向量得到第一力差向量其中,1#,2#,…,i+1#索段的初始无应力长度按照跨径比例分配,1#滑轮两侧分别为1#索段和2#索段,t1表示1#滑轮两侧1#索段和2#索段的力差向量;
c、根据公式Δlslide=Kslide -1·ΔT求解1#,2#,…,i#滑轮的滑移长度向量Δlslide=[Δl1 Δl2 ... Δli]T,修正各索段的无应力长度并重复步骤a-b得到第二力差向量
d、重复步骤a-c,直至||(ΔT(n)-ΔT(n-1))/ΔT(n-1)||≤Δε,n表示迭代次数,提取对应有限元模型中的索的位移和各断面的内力,完成滑移计算。
2.如权利要求1所述的方法,其特征在于,Δε=0.0001-0.0005。
3.如权利要求1所述的方法,其特征在于,索结构的几何参数包含索的初始拟定截面尺寸、索的截面形状、各索段的跨径和各滑轮之间的高差,材料参数包含索的弹性模量、容重、抗拉强度和最小破断拉力。
4.如权利要求1-3任一所述的方法,其特征在于,几何非线性分析方法为牛顿-拉夫森迭代方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210790527.3A CN115130350B (zh) | 2022-07-06 | 2022-07-06 | 一种连续索的非线性滑移计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210790527.3A CN115130350B (zh) | 2022-07-06 | 2022-07-06 | 一种连续索的非线性滑移计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115130350A CN115130350A (zh) | 2022-09-30 |
CN115130350B true CN115130350B (zh) | 2023-06-20 |
Family
ID=83382914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210790527.3A Active CN115130350B (zh) | 2022-07-06 | 2022-07-06 | 一种连续索的非线性滑移计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115130350B (zh) |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4013142A (en) * | 1975-10-07 | 1977-03-22 | Westinghouse Electric Corporation | Elevator system having a drive sheave with rigid but circumferentially compliant cable grooves |
JP4310501B2 (ja) * | 2005-11-07 | 2009-08-12 | 学校法人慶應義塾 | 荷重伝達経路法に基づく数値構造解析装置 |
US7664354B2 (en) * | 2006-08-01 | 2010-02-16 | Nexans | System and method for loose tube tight buffer indoor/outdoor optical fiber cable |
CN101477340B (zh) * | 2009-01-20 | 2010-08-18 | 东南大学 | 基于索力监测的索结构中索系统的递进式健康监测方法 |
CN102733478B (zh) * | 2012-04-04 | 2013-07-10 | 中国航空规划建设发展有限公司 | 一种基于承载全过程研究的索穹顶结构设计指标确定方法 |
ES2550527B1 (es) * | 2014-05-09 | 2016-09-08 | Orona, S. Coop. | Elemento de suspensión y tracción, ascensor y procedimiento de control de la adherencia de dicho elemento a una polea |
CN111753435A (zh) * | 2020-07-04 | 2020-10-09 | 四川公路桥梁建设集团有限公司 | 基于分段悬链线和索力连续算法的缆索吊装系统计算方法 |
CN112035928B (zh) * | 2020-08-31 | 2022-04-15 | 中国石油天然气集团有限公司 | 大跨径悬索管道桥成桥主索线形计算方法 |
-
2022
- 2022-07-06 CN CN202210790527.3A patent/CN115130350B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN115130350A (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Analysis of cable-stayed bridges during construction by cantilever methods | |
CN111783201B (zh) | 一种三跨自锚式悬索桥动力特性的快速分析方法 | |
CN104899377B (zh) | 一种悬索桥索力优化方法 | |
CN111695188B (zh) | 一种地锚式悬索桥动力特性的快速精细分析方法 | |
CN105568864B (zh) | 确定斜拉桥合理施工索力的综合算法 | |
Greco et al. | On the force density method for slack cable nets | |
CN111523172A (zh) | 一种空间异形索面悬索桥主缆成桥线形分析方法 | |
CN107766670B (zh) | 周期性手征蜂窝结构材料等效弹性模量预测方法 | |
CN103850172A (zh) | 一种混合梁斜拉-悬索协作桥 | |
CN103758044A (zh) | 拱肋合拢方法与基于无应力状态法的拱肋合拢控制方法 | |
CN203755129U (zh) | 一种混合梁斜拉-悬索协作桥 | |
CN111666615A (zh) | 一种基于有限元的悬索桥空间缆索找形方法 | |
CN110807221A (zh) | 基于等效力位移法的索力计算方法 | |
Shi et al. | Axis optimisation of arch-shaped pylons for high-speed railway cable-stayed bridges | |
CN111222193A (zh) | 一种确定悬索桥主缆线形的方法 | |
CN115130350B (zh) | 一种连续索的非线性滑移计算方法 | |
Saka et al. | Optimum design of unbraced rigid frames | |
Zhang et al. | Reasonable completed state evaluation for hybrid cable-stayed suspension bridges: an analytical algorithm | |
CN113255029B (zh) | 一种活载作用下悬索桥结构变形及内力的确定方法 | |
CN111967085B (zh) | 大跨径悬索管道桥成桥风索线形计算方法 | |
CN115357965B (zh) | 自锚式悬索桥及其成桥线形确定方法 | |
CN109977565A (zh) | 一种基于fea的斜拉桥动力特性分析方法 | |
CN102561216B (zh) | 悬索桥吊拉组合体系加固结构的计算方法 | |
CN106960114A (zh) | 一种自锚式悬索桥体系转换过程模拟计算方法 | |
CN112035928B (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 |