CN105608291A - 预应力混凝土桥梁调束快速方法 - Google Patents
预应力混凝土桥梁调束快速方法 Download PDFInfo
- Publication number
- CN105608291A CN105608291A CN201610044218.6A CN201610044218A CN105608291A CN 105608291 A CN105608291 A CN 105608291A CN 201610044218 A CN201610044218 A CN 201610044218A CN 105608291 A CN105608291 A CN 105608291A
- Authority
- CN
- China
- Prior art keywords
- effect
- tendon
- load
- result
- calculating
- 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
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
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)
- Bridges Or Land Bridges (AREA)
Abstract
预应力混凝土桥梁调束快速方法,以有效提高预应力混凝土桥梁的设计效率。按如下步骤进行:1)确定结构的几何尺寸、单元节点组成、边界条件及施工阶段划分,组建有限元模型;2)计算成桥阶段影响线加载的活载效应,只改变荷载向量,重复利用总刚度三角分解后结果计算结构活载效应,计算机内存中保存结构的活载效应;3)调整预应力钢束;4)依次精确计算各施工阶段的结构效应,保存各施工阶段的刚度矩阵三角分解结果供下一次调束计算时使用;5)将最后一个施工阶段的结构效应与步骤2的结构活载效应进行荷载组合;6)根据荷载组合结果计算结构的内力、应力和位移;7)重复循环步骤3、4、5、6,直到步骤6得到的计算结果满足规范要求,计算结束。
Description
技术领域
本发明是预应力混凝土桥梁设计方法,具体为预应力混凝土桥梁快速调束方法。
背景技术
近年来我国在土木工程投资和建设规模方面均居世界前列。在混凝土工程和预应力技术应用方面取得了巨大进步。近二三十年来,我国预应力混凝土桥梁发展很快,无论桥型、跨度以及施工方法与技术方面都有了突破性发展。
现有的预应力混凝土桥梁结构分析一般要经过:组建有限元模型、调整预应力钢束、施工阶段计算、成桥活载效应计算和荷载组合分析等过程。其中预应力钢束布置往往需要经过大量的调束尝试来确定,现有的调束计算方法每调整一次预应力钢束,就要重新进行各施工阶段和成桥阶段活载效应计算。其中计算成桥阶段活载效应,由于需要进行影响线计算和动态规划法加载,所以最耗费计算机时间,这就必然导致计算效率低下。为了提高预应力混凝土桥梁的计算效率,首先要能够快速的获得调束结果,这就需要从计算方法上做出发明和创新。
发明内容
本发明所要解决的技术问题是提供一种预应力混凝土桥梁调束快速方法,以有效提高预应力混凝土桥梁的设计效率。
本发明解决其技术问题所采用的技术方案如下:
本发明的预应力混凝土桥梁快速调束的方法,按如下步骤进行:
1)确定结构的几何尺寸、单元节点组成、边界条件及施工阶段划分,组建有限元模型;
2)计算成桥阶段影响线加载的活载效应,对结构总刚度矩阵进行一次三角分解,计算活载时只改变荷载向量,重复利用总刚度三角分解后结构计算结构活载效应,计算机内存中保存结构的活载效应,供后续的荷载组合使用;
3)调整预应力钢束;
4)依次精确计算各施工阶段的结构效应,粗略计算时忽略钢束变化对截面几何特性的影响,保存各施工阶段的刚度矩阵三角分解结果供下一次调束计算时使用;
5)将最后一个施工阶段的结构效应与步骤2的结构活载效应进行荷载组合;
6)根据荷载组合结果计算结构的内力、应力和位移;
7)重复循环步骤3、4、5、6,直到步骤6得到的计算结果满足规范要求,计算结束。
本发明的有益效果是,将最耗费计算机时间的结构活载效应计算前置并保存计算结果,结构活载效应计算结果可以为每次调束使用;计算活载效应时,保存结构总刚度矩阵的三角分解结果,当单位荷载移位时只需要改变荷载向量,从而提高了矩阵求解速度,有利于提高桥梁设计效率。
附图说明
本说明书包括如下两幅附图:
图1是现有通用预应力混凝土桥梁调束快速方法的流程图;
图2是本发明预应力混凝土桥梁调束快速方法的流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
参照图2,本发明的预应力混凝土桥梁快速调束的方法。
本发明目的的实现方式为,将结构的活载效应计算前置,按如下步骤进行:
1)确定结构的几何尺寸、单元节点组成、边界条件及施工阶段划分,组建有限元模型;
2)计算成桥阶段影响线加载的活载效应,为了提高计算效率只对结构总刚度矩阵进行一次三角分解,计算活载时只改变荷载向量,重复利用总刚度三角分解后结果计算结构活载效应,计算机内存中保存结构的活载效应,供后续的荷载组合使用;
3)调整预应力钢束;
4)依次精确计算各施工阶段的结构效应。粗略计算时可忽略钢束变化对截面几何特性的影响,保存各施工阶段的刚度矩阵三角分解结果供下一次调束计算时使用;
5)将最后一个施工阶段的结构效应与步骤2的结构活载效应进行荷载组合;
6)根据荷载组合结果计算结构的内力、应力和位移;
7)重复循环步骤3、4、5、6,直到步骤6得到的计算结果满足规范要求,计算结束。
所述步骤(2)中活载效应计算如下:
成桥阶段的结构总刚度矩阵为[K],用改进平方根法将结构总刚度矩阵[K]分解为三个矩阵的乘积:
[K]=[L][D][L]T(1)
式(1)中:[L]为主对角元为1的下三角矩阵,[D]为仅有主对角元的矩阵。
可详细表示如下:
将矩阵[L]和[D][L]T保存在内存中,为了后续的计算使用。
矩阵[K]分解之后,结构总体节点平衡方程为:
[L][D][L]T{x}={p}(3)
式中:{x}为节点位移向量,{p}为节点荷载向量。
求解式(3)可得到节点的位移,进一步计算可得到单元节点内力和截面应力。
以上公式及公式中各参数的含义为本领域技术人员所常用,详见《桥梁结构分析及程序系统》(肖汝诚人民交通出版社),其他公式推导则和常规推导公式相类似,在此不再赘述。
如果桥面影响线的点数为n,则需要求解n次刚度矩阵。由于成桥阶段的结构刚度矩阵秩最大,矩阵三角分解和求解方程与其他阶段相比都耗时较多。
通过对预应力混凝土桥梁计算原理的分析,可以发现最耗费计算机时间的步骤在结构刚度矩阵的三角分解和求解矩阵,活载效应要计算结构每个单元、节点的内力影响线,势必进行很多次的刚度矩阵求解,从而耗费大量的计算机时间。预应力钢束对结构而言是内力,调束对相关单元的截面几何特性影响不大,对结构总刚度的影响更是非常小。桥梁结构荷载组合是将不同状态时取用的荷载及其代表值进行组合,所以可以将活载效应计算前置,并保存活载效应结果,供每次调束后的荷载组合使用,避免了每次调束都进行大量的刚度矩阵三角分解和求解计算,从而提高了调束效率。
以铁路某预应力混凝土连续梁(72.8+136+72.8m)调束为例:
建立该桥的全桥有限元模型,详细信息见下表:
表1
类别 | 数量 |
单元个数 | 90 |
节点个数 | 91 |
截面类型个数 | 46 |
约束个数 | 18 |
施工阶段个数 | 48 |
荷载组合个数 | 3 |
采用的计算机配置和操作系统:
操作系统:MicrosoftWindowsXPProfessionalServicePack3
系统硬件:i7-2600CPU3.40GHz4GB内存
通用调束方法耗时:9078.5898毫秒
本发明调束方法耗时:2453.2506毫秒
减少耗时比例:73%
由此可见,本发明方法在调整预应力钢束时将结构活载效应计算前置,计算预应力混凝土连续梁时与通用调束方法相比节省时间达到73%,从而极大的提高了调束效率。
Claims (1)
1.预应力混凝土桥梁快速调束的方法,按如下步骤进行:
1)确定结构的几何尺寸、单元节点组成、边界条件及施工阶段划分,组建有限元模型;
2)计算成桥阶段影响线加载的活载效应,对结构总刚度矩阵进行一次三角分解,计算活载时只改变荷载向量,重复利用总刚度三角分解后结果计算结构活载效应,计算机内存中保存结构的活载效应,供后续的荷载组合使用;
3)调整预应力钢束;
4)依次精确计算各施工阶段的结构效应,粗略计算时忽略钢束变化对截面几何特性的影响,保存各施工阶段的刚度矩阵三角分解结果供下一次调束计算时使用;
5)将最后一个施工阶段的结构效应与步骤2的结构活载效应进行荷载组合;
6)根据荷载组合结果计算结构的内力、应力和位移;
7)重复循环步骤3、4、5、6,直到步骤6得到的计算结果满足规范要求,计算结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044218.6A CN105608291B (zh) | 2016-01-22 | 2016-01-22 | 预应力混凝土桥梁调束快速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044218.6A CN105608291B (zh) | 2016-01-22 | 2016-01-22 | 预应力混凝土桥梁调束快速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105608291A true CN105608291A (zh) | 2016-05-25 |
CN105608291B CN105608291B (zh) | 2019-01-04 |
Family
ID=55988223
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610044218.6A Active CN105608291B (zh) | 2016-01-22 | 2016-01-22 | 预应力混凝土桥梁调束快速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105608291B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111291510A (zh) * | 2020-01-20 | 2020-06-16 | 中车齐齐哈尔车辆有限公司 | 有限元工况加载方法及分析方法 |
CN111985084A (zh) * | 2020-07-24 | 2020-11-24 | 华设设计集团股份有限公司 | 基于有限变量的桥梁力学实时分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101587046A (zh) * | 2009-06-19 | 2009-11-25 | 同济大学 | 基于动测法的偏心直线预应力混凝土梁桥承载力评估方法 |
CN103065035A (zh) * | 2012-09-29 | 2013-04-24 | 长沙理工大学 | 一种变截面预应力混凝土连续梁桥纵向预应力束设计方法 |
CN103268371A (zh) * | 2013-04-26 | 2013-08-28 | 重庆交通大学 | 一种基于影响矩阵的桥梁实时荷载识别方法 |
CN104933285A (zh) * | 2015-03-05 | 2015-09-23 | 西南交通大学 | 一种桥梁现场静载试验评定方法 |
-
2016
- 2016-01-22 CN CN201610044218.6A patent/CN105608291B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101587046A (zh) * | 2009-06-19 | 2009-11-25 | 同济大学 | 基于动测法的偏心直线预应力混凝土梁桥承载力评估方法 |
CN103065035A (zh) * | 2012-09-29 | 2013-04-24 | 长沙理工大学 | 一种变截面预应力混凝土连续梁桥纵向预应力束设计方法 |
CN103268371A (zh) * | 2013-04-26 | 2013-08-28 | 重庆交通大学 | 一种基于影响矩阵的桥梁实时荷载识别方法 |
CN104933285A (zh) * | 2015-03-05 | 2015-09-23 | 西南交通大学 | 一种桥梁现场静载试验评定方法 |
Non-Patent Citations (1)
Title |
---|
邢兵 等: "初始应力状态对加固后混凝土梁力学性能的影响", 《桥梁建设》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111291510A (zh) * | 2020-01-20 | 2020-06-16 | 中车齐齐哈尔车辆有限公司 | 有限元工况加载方法及分析方法 |
CN111291510B (zh) * | 2020-01-20 | 2022-04-29 | 中车齐齐哈尔车辆有限公司 | 有限元工况加载方法及分析方法 |
CN111985084A (zh) * | 2020-07-24 | 2020-11-24 | 华设设计集团股份有限公司 | 基于有限变量的桥梁力学实时分析方法 |
CN111985084B (zh) * | 2020-07-24 | 2023-10-27 | 华设设计集团股份有限公司 | 基于有限变量的桥梁力学实时分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105608291B (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101706741B (zh) | 一种基于负载平衡的cpu和gpu两级动态任务划分方法 | |
CN108595728A (zh) | 一种蜂窝材料的铺层等效有限元模型构建方法 | |
CN109359434B (zh) | 模板支架体系衍生式设计方法 | |
CN112784359B (zh) | 薄壁梁约束扭转极限承载力迭代计算方法 | |
JP2019537079A (ja) | 大規模再生可能エネルギーのデータについて確率モデルを構築する方法 | |
CN106096257A (zh) | 一种非线性索单元分析方法及系统 | |
CN105608291A (zh) | 预应力混凝土桥梁调束快速方法 | |
CN102087678A (zh) | 通过时域划分实现并行计算的电路仿真方法 | |
CN101788425A (zh) | 一种结构件复合型裂纹前缘应力强度因子分离和分布的确定方法 | |
CN103915840A (zh) | 一种基于Givens正交增量行变换的大电网状态估计方法 | |
CN106096183B (zh) | 一种基于特征线法的多重并行方法 | |
CN106650155A (zh) | 一种大规模振动工程控制装置优化布置方法 | |
CN113378355A (zh) | 一种考虑弯矩作用的薄壁梁剪切极限强度计算方法 | |
CN111274624B (zh) | 一种基于rbf代理模型的多工况异形节点拓扑优化设计方法 | |
CN104732043A (zh) | 变电架构的设计方法 | |
CN115034015A (zh) | 一种办公塔楼结构系统多级优化方法 | |
CN102254245A (zh) | 电力系统调度日计划分时段并行安全稳定校核方法 | |
CN114626127A (zh) | 一种输电塔主材弯矩调整系数的计算方法、系统及设备 | |
CN112784333B (zh) | 一种现浇混凝土箱梁横梁的设计方法 | |
CN104166593A (zh) | 一种计算多应用功能异步并发调度方法 | |
CN111985024A (zh) | 一种调整建筑结构的构件的方法、装置以及电子设备 | |
CN110135038B (zh) | 一种应用于客车快速轻量化的分析方法 | |
CN104281494A (zh) | 一种基于插值算法面向计算通信联合优化的负载均衡方法 | |
CN109033572B (zh) | 一种箱型截面框架结构极限承载力分析的线弹性方法 | |
CN114936393B (zh) | 一种用于确定斜拉桥最优梁高的方法、装置及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |