CN107194080B - 一种轮辐式张拉体系的结构确定方法 - Google Patents

一种轮辐式张拉体系的结构确定方法 Download PDF

Info

Publication number
CN107194080B
CN107194080B CN201710378618.5A CN201710378618A CN107194080B CN 107194080 B CN107194080 B CN 107194080B CN 201710378618 A CN201710378618 A CN 201710378618A CN 107194080 B CN107194080 B CN 107194080B
Authority
CN
China
Prior art keywords
spoke type
parameter
tensioning system
parameters
structural model
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
Application number
CN201710378618.5A
Other languages
English (en)
Other versions
CN107194080A (zh
Inventor
张月强
丁洁民
张峥
黄卓驹
李璐
周旋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Architecture Design and Research Institute of Tongji University Group Co Ltd
Original Assignee
Architecture Design and Research Institute of Tongji University Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Architecture Design and Research Institute of Tongji University Group Co Ltd filed Critical Architecture Design and Research Institute of Tongji University Group Co Ltd
Priority to CN201710378618.5A priority Critical patent/CN107194080B/zh
Publication of CN107194080A publication Critical patent/CN107194080A/zh
Application granted granted Critical
Publication of CN107194080B publication Critical patent/CN107194080B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling 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)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Complex Calculations (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种轮辐式张拉体系的结构确定方法,所述方法包括下列步骤:根据轮辐式张拉体系的结构参数,确定轮辐式张拉体系的结构模型;对结构模型进行受力分析,得到相关的线性方程组;判断线性方程组是否得到最优解,若是则得到轮辐式张拉体系的结构,若否则修改轮辐式张拉体系的结构参数并返回修改结构模型。与现有技术相比,本发明具有结构确定准确、实现简单以及便于推广等优点。

Description

一种轮辐式张拉体系的结构确定方法
技术领域
本发明涉及结构找形领域,尤其是涉及一种轮辐式张拉体系的结构确定方法。
背景技术
轮辐式张拉结构是在张拉结构的基础上根据自行车车轮辐条结构概念提出来的,这种结构体系具有张拉结构的受力特点,同时克服了一般张拉结构体系复杂、传力不直接和容易形成机构的缺点。这种结构轻巧美观,能够与建筑造型和建筑美学很好的融合。同时有研究表明,随着跨度的增加,其造价增加并不大,这在大跨度空间结构方面具有很大的优势。近年来该结构体系深受建筑师的喜爱,在大型体育场得到广泛的应用,但确定该结构体系的几何形态和相应的预应力分布比较困难。
发明内容
本发明的目的是针对上述问题提供一种轮辐式张拉体系的结构确定方法。
本发明的目的可以通过以下技术方案来实现:
一种轮辐式张拉体系的结构确定方法,所述方法包括下列步骤:
1)根据轮辐式张拉体系的结构参数,确定轮辐式张拉体系的结构模型;
2)对步骤1)中建立的结构模型进行受力分析,得到相关的线性方程组;
3)判断步骤2)得到的线性方程组是否得到最优解,若是则得到轮辐式张拉体系的结构,若否则修改轮辐式张拉体系的结构参数并返回步骤1)。
所述步骤1)具体为:
11)接收外界输入的轮辐式张拉体系的结构参数;
12)根据步骤11)接收的结构参数,通过参数化设计软件确定轮辐式张拉体系的结构。
所述轮辐式张拉体系的结构参数包括上压环参数、内拉环参数、下压环参数、索桁架榀数、承重索参数、吊索参数、稳定索参数、上压环环索参数和吊索数。
所述参数化设计软件通过Python语言实现。
所述步骤2)具体为:
21)接收外界输入的预应力分析初始参数;
22)根据步骤21)接收的预应力分析初始参数,计算步骤1)中建立的结构模型的预应力分布;
23)根据步骤22)计算得到的结构模型的预应力分布,计算结构模型在初始态和荷载态下的受力性能,得到相关的线性方程组。
所述预应力分析初始参数包括主动索初始拉力和被动索初始拉力。
所述步骤22)具体为:通过矩阵平衡法计算,结合步骤21)接收的预应力分析参数,得到步骤1)中建立的结构模型的预应力分布。
所述计算结构模型在初始态和荷载态下的受力性能具体为:通过牛顿-拉夫逊法将非线性的结构模型的预应力分布变换为线性方程组。
所述判断步骤2)得到的线性方程组是否得到最优解具体为:根据结构模型的受力结果参数最小化原则判断步骤2)中得到相关的线性方程组是否得到最优解,若是则表明受力分析结果达到最优,若否则表明受力分析结果未达到最优。
所述受力结果参数包括应变能或预应力。
与现有技术相比,本发明具有以下有益效果:
(1)通过确定轮辐式张拉体系的结构参数,对轮辐式张拉体系进行结构找形分析,继而再通过受力分析对轮辐式张拉体系进行结构找力分析,通过结构找力分析和结构找形分析的结合,可以充分确定最符合当前情况的轮辐式张拉体系结构,确保得到结果的准确性。
(2)对于轮辐式张拉体系进行结构找形分析时,通过输入参数得到其初始结构,并结合找形分析的结果对参数进行优化,从而对结构进行修正,保证了轮辐式张拉体系可以满足受力最优,从而使得得到的结果最好且具有实用性。
(3)轮辐式张拉体系的结构参数包括上压环参数、内拉环参数、下压环参数、索桁架榀数、承重索参数、吊索参数、稳定索参数、上压环环索参数和吊索数,充分考虑了轮辐式张拉体系的所有关键参数,从而提升了最后结果的准确性。
(4)参数化设计软件通过Python语言实现,便于实现的同时也保证了参数化设计软件的高运算速度。
(5)对于轮辐式张拉体系进行结构找力分析时,首先接收预应力分析的初始参数,充分考虑到了轮辐式张拉结构为全柔性结构,在未施加初始预应力之前,结构并不具有刚度和承载能力,且初始参数是根据结构所处地区的实际情况确定的,因而得到的结果也更具有实际意义和应用价值。
(6)通过矩阵平衡法计算预应力分布,该方法准确度高且便于计算,实用性能强。
(7)最后通过判断应变能或预应力是否达到最小来判断结构是否达到最优,可以根据结构所处环境选择判断方法,灵活性强,适应性好。
附图说明
图1为轮辐式张拉体系的结构确定方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本实施例提供的是一种轮辐式张拉体系的结构确定方法,包括下列步骤:
1)根据轮辐式张拉体系的结构参数,确定轮辐式张拉体系的结构模型:
11)接收外界输入的轮辐式张拉体系的结构参数;
12)根据步骤11)接收的结构参数,通过参数化设计软件确定轮辐式张拉体系的结构;
2)对步骤1)中建立的结构模型进行受力分析,得到相关的线性方程组:
21)接收外界输入的预应力分析初始参数;
22)根据步骤21)接收的预应力分析初始参数,计算步骤1)中建立的结构模型的预应力分布:
通过矩阵平衡法计算,结合步骤21)接收的预应力分析参数,得到步骤1)中建立的结构模型的预应力分布;
23)根据步骤22)计算得到的结构模型的预应力分布,计算结构模型在初始态和荷载态下的受力性能,得到相关的线性方程组,具体为:
通过牛顿-拉夫逊法将非线性的结构模型的预应力分布变换为线性方程组。;
3)判断步骤2)得到的线性方程组是否得到最优解,若是则得到轮辐式张拉体系的结构,若否则修改轮辐式张拉体系的结构参数并返回步骤1),具体为:
根据结构模型的受力结果参数最小化原则判断步骤2)中得到相关的线性方程组是否得到最优解,若是则表明受力分析结果达到最优,若否则表明受力分析结果未达到最优。
其中,轮辐式张拉体系的结构参数包括上压环参数、内拉环参数、下压环参数、索桁架榀数、承重索参数、吊索参数、稳定索参数、上压环环索参数和吊索数,参数化设计软件通过Python语言实现,预应力分析初始参数包括主动索初始拉力和被动索初始拉力,受力结果参数包括应变能或预应力。
根据上述步骤来具体确定轮辐式张拉体系的结构的过程为:参数化设计采用Python语言编制参数化的模块,最终通过调节若干个参数可以创建结构的整体模型;将得到的模型导入有限元软件,通过编制的找形分析有限元插件利用矩阵平衡法找出结构的预应力分布,并计算结构在初始态和荷载态下的受力性能,通过参数化设计软件和有限元软件的若干次循化设计,根据最小应变能原则或最小预应力方差原则找到线性方程组的最优解,最终找出结构的合理形态和预应力分布。实际操作过程即为:将参数化的模型转化为DXF格式,通过编制的导入插件一键导入到有限元软件中,然后通过编制的找形分析插件进行找力分析。找力分析前需要设置主动索和被动索的初始拉力,其余的刚性杆均设置成普通单元,点击开始找出结构的初始预应力分布。

Claims (7)

1.一种轮辐式张拉体系的结构确定方法,其特征在于,所述方法包括下列步骤:
1)根据轮辐式张拉体系的结构参数,确定轮辐式张拉体系的结构模型;
2)对步骤1)中建立的结构模型进行受力分析,得到相关的线性方程组,具体为:
21)接收外界输入的预应力分析初始参数;
22)根据步骤21)接收的预应力分析初始参数,计算步骤1)中建立的结构模型的预应力分布,具体为:通过矩阵平衡法计算,结合步骤21)接收的预应力分析参数,得到步骤1)中建立的结构模型的预应力分布;
23)根据步骤22)计算得到的结构模型的预应力分布,计算结构模型在初始态和荷载态下的受力性能,得到相关的线性方程组;
所述计算结构模型在初始态和荷载态下的受力性能具体为:通过牛顿-拉夫逊法将非线性的结构模型的预应力分布变换为线性方程组;
3)判断步骤2)得到的线性方程组是否得到最优解,若是则得到轮辐式张拉体系的结构,若否则修改轮辐式张拉体系的结构参数并返回步骤1)。
2.根据权利要求1所述的轮辐式张拉体系的结构确定方法,其特征在于,所述步骤1)具体为:
11)接收外界输入的轮辐式张拉体系的结构参数;
12)根据步骤11)接收的结构参数,通过参数化设计软件确定轮辐式张拉体系的结构模型。
3.根据权利要求2所述的轮辐式张拉体系的结构确定方法,其特征在于,所述轮辐式张拉体系的结构参数包括上压环参数、内拉环参数、下压环参数、索桁架榀数、承重索参数、吊索参数、稳定索参数、上压环环索参数和吊索数。
4.根据权利要求2所述的轮辐式张拉体系的结构确定方法,其特征在于,所述参数化设计软件通过Python语言实现。
5.根据权利要求1所述的轮辐式张拉体系的结构确定方法,其特征在于,所述预应力分析初始参数包括主动索初始拉力和被动索初始拉力。
6.根据权利要求1所述的轮辐式张拉体系的结构确定方法,其特征在于,所述判断步骤2)得到的线性方程组是否得到最优解具体为:根据结构模型的受力结果参数最小化原则判断步骤2)中得到相关的线性方程组是否得到最优解,若是则表明受力分析结果达到最优,若否则表明受力分析结果未达到最优。
7.根据权利要求6所述的轮辐式张拉体系的结构确定方法,其特征在于,所述受力结果参数包括应变能或预应力。
CN201710378618.5A 2017-05-25 2017-05-25 一种轮辐式张拉体系的结构确定方法 Active CN107194080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710378618.5A CN107194080B (zh) 2017-05-25 2017-05-25 一种轮辐式张拉体系的结构确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710378618.5A CN107194080B (zh) 2017-05-25 2017-05-25 一种轮辐式张拉体系的结构确定方法

Publications (2)

Publication Number Publication Date
CN107194080A CN107194080A (zh) 2017-09-22
CN107194080B true CN107194080B (zh) 2024-04-05

Family

ID=59874926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710378618.5A Active CN107194080B (zh) 2017-05-25 2017-05-25 一种轮辐式张拉体系的结构确定方法

Country Status (1)

Country Link
CN (1) CN107194080B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114781142B (zh) * 2022-04-11 2023-03-10 北京市建筑设计研究院有限公司 单环索轮辐式索桁架结构的形态分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622485A (zh) * 2012-03-16 2012-08-01 东南大学 确定索杆梁系空间结构放样态的形态分析逆迭代法
CN102720295A (zh) * 2012-04-04 2012-10-10 中国航空规划建设发展有限公司 一种基于索穹顶张拉和承载全过程分析的预应力确定方法
CN105243185A (zh) * 2015-09-10 2016-01-13 哈尔滨工业大学 一种基于蒙特卡洛法的张拉整体结构找形方法
CN106599509A (zh) * 2016-12-26 2017-04-26 河北建筑工程学院 一种模拟弦支穹顶结构施加预应力的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622485A (zh) * 2012-03-16 2012-08-01 东南大学 确定索杆梁系空间结构放样态的形态分析逆迭代法
CN102720295A (zh) * 2012-04-04 2012-10-10 中国航空规划建设发展有限公司 一种基于索穹顶张拉和承载全过程分析的预应力确定方法
CN105243185A (zh) * 2015-09-10 2016-01-13 哈尔滨工业大学 一种基于蒙特卡洛法的张拉整体结构找形方法
CN106599509A (zh) * 2016-12-26 2017-04-26 河北建筑工程学院 一种模拟弦支穹顶结构施加预应力的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
北京奥运会乒乓球馆大跨度预应力桁架壳体结构设计;何志军;丁洁民;巢斯;肖小凌;;建筑结构(第03期);全文 *

Also Published As

Publication number Publication date
CN107194080A (zh) 2017-09-22

Similar Documents

Publication Publication Date Title
CN103150460B (zh) 一种间接空冷塔的结构分析方法
CN101751376A (zh) 利用cpu和gpu协同工作对三角线性方程组求解的加速方法
CN106934179B (zh) 输电铁塔角钢主节点拉伸试验数据处理方法
CN106529719A (zh) 基于粒子群优化算法风速融合的风电功率预测的方法
CN107194080B (zh) 一种轮辐式张拉体系的结构确定方法
CN106096257A (zh) 一种非线性索单元分析方法及系统
CN109388858B (zh) 基于头脑风暴优化算法的非线性传感器校正方法
Zong-ling et al. The design of lightweight and multi parallel CNN accelerator based on FPGA
CN102004822B (zh) 基于模态的空间六自由度并联运动系统频率分析方法
CN111507027A (zh) 一种基于有限元软件的钢桁架塔式结构整体动力失效时刻的判别方法
Chen et al. A DNN optimization framework with unlabeled data for efficient and accurate reconfigurable hardware inference
CN115099522A (zh) 基于bp神经网络的专变用户有功无功碳排放预测方法
CN111597685B (zh) 一种面向多自由度波浪能装置获能体的形状优化方法
CN102279072A (zh) 步行器测力系统中的蚁群遗传融合标定方法
谢斌 et al. Structure optimization of frame for field vegetable power machinery based on measured load data
CN113761633A (zh) 一种剪力墙结构优化设计方法
CN105447263A (zh) 一种用于复杂连接结构非线性分析的螺栓预紧力施加方法
CN105334883B (zh) 一种振动控制系统中的智能前馈信号分析方法
Lei et al. Dynamic modeling of bicycle robot and nonlinear control based on feedback linearization of MIMO systems
CN111241731A (zh) 一种依据水平张力的输电导线找形方法
Bonafede et al. Compliant legs in biped robots for running optimization
Suzuki et al. ZMP reference trajectory generation for biped robot with inverted pendulum model by using virtual supporting point
CN204904325U (zh) 一种梯级水电站群多目标优化调度管理系统
CN107563612A (zh) 一种电力电站钛合金管板造价的分析方法
CN115033752A (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