CN111967176B - 一种外闪木框架的纠偏荷载预测方法 - Google Patents
一种外闪木框架的纠偏荷载预测方法 Download PDFInfo
- Publication number
- CN111967176B CN111967176B CN202010829983.5A CN202010829983A CN111967176B CN 111967176 B CN111967176 B CN 111967176B CN 202010829983 A CN202010829983 A CN 202010829983A CN 111967176 B CN111967176 B CN 111967176B
- Authority
- CN
- China
- Prior art keywords
- frame
- displacement
- rebound
- loading
- deviation
- 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
- 239000002023 wood Substances 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000006073 displacement reaction Methods 0.000 claims abstract description 60
- 230000003014 reinforcing effect Effects 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 239000000463 material Substances 0.000 claims abstract description 8
- 238000011835 investigation Methods 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 16
- 230000002787 reinforcement Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 2
- 230000035882 stress Effects 0.000 abstract description 7
- 230000032683 aging Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000010276 construction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 240000004307 Citrus medica Species 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种外闪木框架的纠偏荷载预测方法,属于建筑结构分析技术领域,步骤包括:S1.现场调研,获得待纠偏框架的几何与材料参数、约束边界条件以及框架的初始偏移;S2.建立待纠偏框架的有限元模型;S3.对外闪框架加载试探位移,求解结构变形和单元应力;S4.对加载后结构进行加固,建立有限元模型,释放加载位移,使结构在内力作用下回弹;S5.测算回弹后的框架倾斜度,若达不到纠偏要求,则修正得到新的试探位移,重复步骤S3‑S5,直至框架在回弹后倾斜度为零;S6.当框架倾斜度达到纠偏要求,加载的试探位移所对应施加的纠偏力即为预测得到的木框架纠偏荷载。本发明易操作、计算效率高,解决了由地震、强风、地基沉陷或老化作用下发生倾斜变形的木框架的纠偏问题。
Description
技术领域
本发明属于建筑结构分析技术领域,涉及一种纠偏荷载预测方法,特别是一种外闪木框架的纠偏荷载预测方法。
背景技术
中国传统木框架建筑大多以横向梁柱结构为承重结构,在构件之间多采用榫卯进行连接,而构架之间由檩、橼连成整体。由于榫卯节点为非整体性连接,木框架的抗变形刚度较弱,在地震、强风或地基沉陷的作用下易发生节点松脱和框架外闪,即向外发生了倾斜。我国现存的古建筑木框架大部分存在不同程度的外闪变形,对外闪木框架进行纠偏加固,是古建筑木结构修缮加固的重点工作。
传统的木框架纠偏方法是利用张拉装置或者千斤顶,根据框架外闪的情况以经验为主进行现场加载,纠偏效果难以保证。目前对于纠偏方案进行合理规划,并拥有较好纠偏效果的方法较少,因此在具有理论依据和可操作性的前提下,解决古建筑木结构纠偏加固的问题显得较为重要。
发明内容
本发明的目的是针对上述现有技术中存在的不足,提出一种外闪木框架的纠偏荷载预测方法,可在施工前对木框架的纠偏方案及加载后果做出预测,能以较高的精度确保框架最终回弹到理想的位置,以解决由地震、强风、地基沉陷或老化作用下发生倾斜变形的木框架的纠偏问题。
本发明的技术方案:一种外闪木框架的纠偏荷载预测方法,其特征在于:所述方法包括如下步骤:
S1.由现场调研,获得待纠偏框架的几何与材料参数、约束边界条件以及框架的初始偏移;
S2.建立待纠偏框架的有限元模型,并预先设定结构的纠偏目标;
S3.对倾斜框架加载一个试探位移,然后求解结构的变形和单元应力;
S4.对加载后的结构进行加固,建立加固后的有限元模型,然后释放加载位移,使结构在内力作用下回弹;
S5.测算回弹后框架的倾斜度,若回弹后的结构仍然倾斜,则修正获得新的试探位移,然后重复步骤S3-S5,直至框架回弹后倾斜度为零;
S6.当框架倾斜度达到纠偏要求后,所加载试探位移对应施加的纠偏力即为预测得到的纠偏荷载。
步骤中S3中所述对倾斜结构加载一个试探位移,然后求解结构的变形和单元应力,具体步骤如下:
S31.对于代表真实倾斜框架的有限元模型,在预先设定的加载节点,该模型与垂直位置的初始偏移等于Δ0;如果将试探位移Δ*应用于加载节点,那么获得的节点位移向量U将分为两部分:第一部分ua为纠偏时非加载节点的位移,第二部分ub为纠偏时加载点的位移;对于该过程,有限元方程为:
其中,α是一个大数,取1010,用于乘以刚度矩阵中加载自由度的对角线元素;
S32.根据求解结果,向量ub等于Δ*,ua是的非加载节点位移向量;各变形单元的应力向量σe在单元局部坐标系内计算:
其中,De为材料的本构矩阵,Be为构件的应变位移矩阵,ue为单元的位移向量,ue,a和ue,b分别是ua和ub中的单元位移分量。
步骤中S4中所述对加载后的结构进行加固,建立加固后的有限元模型,然后释放加载位移,使结构在内力作用下回弹,具体步骤如下:
S41.在加载后,结构的几何形态发生变化,为改善结构的性能需采用加固措施,因此需建立一个加固后的新模型:
其中,上标“′”表示加固后的结构矩阵,以表示与加固前矩阵的不同;
S42.结构的内力向量F通过组合各个变形单元的等效节点力得到:
其中,B′e为加固后结构的应变位移矩阵;
S43.求解获得加固后的结构在内力下的回弹位移u′a和u′b。
步骤中S5中所述测算回弹后框架的倾斜度,若结构仍然倾斜,达不到纠偏目标,则修正获得新的试探位移,然后重复步骤S3-S5,直至框架回弹后倾斜度为零,具体步骤如下:
S51.纠偏时消除加载节点上预先存在的偏移Δ0,效果进行如下检测:
ub+u′b=-Δ0 (5)
若检测不满足要求,则修正获得一个新的试探位移:
Δ*=(1-ω)ub+ω(-Δ0-u′b) (6)
其中,ω是比例因子,取1.3~1.6;
S52.由于一次计算无法满足纠偏要求,需重复步骤S3-S5进行迭代计算,以比例因子ω连续调整试探位移Δ*,以保证倾斜框架在加固并释放荷载后能够回弹到无倾斜的理想形态。
本发明的有益效果为:本发明提出的一种外闪木框架的纠偏荷载预测方法,与现有的外闪木框架纠偏方法相比,本发明能够在施工前对木框架的纠偏方案及加载后果做出合理的预测和规划,能以较高的精度确保框架最终回弹到理想位置,方法新颖且科学合理,利用修正算法,将外闪木框架的纠偏问题转换为结构加载最优解的估计问题,开发了适用于计算机分析的迭代流程,大大降低了木框架纠偏问题的分析难度,具有易操作和计算效率高的特点,解决了由地震、强风、地基沉陷或老化作用下发生倾斜变形的木框架的纠偏问题。
附图说明
图1为本发明逻辑流程框图。
图2为本发明中外闪木框架的初始状态示意图。
图3为本发明中外闪木框架的纠偏加载示意图。
图4为本发明中框架弹簧节点的加固示意图。
图5为本发明中框架卸载、回弹后的变形示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
实施例:
对图2所示外闪木框架利用本发明的方法进行纠偏。结构高度为h,宽度为l,顶部存在初始水平偏差Δ0。梁与柱、柱与基础之间采用非线性的转动弹簧进行连接。
利用本发明的技术对外闪木框架进行纠偏荷载计算,具体包括以下步骤:
S1.由现场调研,获得待纠偏框架的几何与材料参数、约束边界条件以及框架的初始偏移;在本例中,结构高度h和宽度l均为4m,梁柱的横截面高度为0.3m,宽度为0.2m,材料的弹性模量为1.13×1010N/m2,泊松比为0.45,结构的初始水平偏差Δ0=0.4m。连接弹簧的转角(rad)—弯矩(M·m)关系如下:
M=1.691×107θ4-6.619×106θ3+5.991×105θ2+3.072×104θ+84.59 (1)
S2.建立待纠偏框架的有限元模型,并预先设定结构的纠偏目标;本例中,拟在框架的横梁处进行纠偏加载,纠偏目标为使初始偏差Δ0消失为零。
S3.对倾斜框架加载一个试探位移,然后求解结构的变形和单元应力,包括以下步骤:
S31.对结构加载水平向右的试探位移Δ*=0.4m,如图3所示,然后采用有限元方法求解结构的位移向量ua和ub:
根据求解结果,ub约等于Δ*。
S32.根据结构的整体变形,求解各个单元内部的应力向量:
其中,De为材料的本构矩阵,Be为单元的应变位移矩阵,ue为单元的位移向量,ue,a和ue,b分别是ua和ub中的单元位移分量。
S4.对加载后的结构进行加固,建立加固后的有限元模型,然后释放加载位移,使结构在内力作用下回弹,包括以下步骤:
S41.对加载后的结构进行加固,在本例中是对弹簧节点进行加固,如图4所示。加固后的弹簧的转动承载力为原来的2倍,然后根据结构几何形态的变化以及加固情况,建立新的有限元模型:
其中,上标“′”表示加固后的结构向量和矩阵,以区分加固前的模型。
S42.结构内力向量F由集成各个变形单元的等效节点力得到:
其中,B′e为加固后结构单元的应变位移矩阵。
S43.释放荷载,求解加固后的结构在内力作用下的回弹位移u′a和u′b。
S5.测算回弹后框架的倾斜度,如果回弹后的结构仍然倾斜,则修正获得新的试探位移,然后如图1所示重复步骤S3-S5,直至框架回弹后倾斜度为零,包括以下步骤:
S51.结构纠偏应消除加载节点上预先存在的位移0.4m,纠偏效果应进行如下检测:
ub+u′b=-0.4 (7)
如果不满足要求,则通过修正得到新的试探位移:
Δ*=(1-ω)ub+ω(-0.4-u′b) (8)
其中,ω是比例因子,本例中取1.3。
S52.由于一次纠偏计算无法满足要求,需重复步骤S3-S5进行迭代计算,以比例因子1.3连续调整试探位移Δ*,以保证倾斜框架在加固并释放荷载后能够回弹到无倾斜的理想形态。当达到以上目标时,此时的试探位移就是最终的加载位移。
在本例中,经过6次迭代计算,得到最终的纠偏加载位移ub为向左0.644m,弹簧节点经过加固并释放荷载后,结构回弹位移ub′为向右0.244m,恰好能够消除初始偏移0.4m。
S6.当框架倾斜度达到纠偏要求后,所加载试探位移对应施加的纠偏力即为预测得到的纠偏荷载。
在本例中,纠偏加载位移为向左0.644m,对应施加的纠偏荷载为向左4.36kN;在纠偏加载后,弹簧节点的变形值为0.052rad,梁柱中的最大等效应力为1.35MPa,结构较为安全。
由上可知,本发明的外闪木框架纠偏荷载预测方法,可在施工前对木框架的纠偏方案及加载后果做出预测,能以较高的精度确保框架最终回弹到理想的位置,方法科学,具有易操作和计算效率高的特点,可较好解决由地震、强风、地基沉陷或老化作用下发生倾斜变形的木框架的纠偏问题。
Claims (2)
1.一种外闪木框架的纠偏荷载预测方法,其特征在于:所述方法包括如下步骤:
S1.现场调研,获得待纠偏框架的几何与材料参数、约束边界条件以及框架的初始偏移;
S2.建立待纠偏框架的有限元模型,并预先设定结构的纠偏目标;
S3.对倾斜框架加载一个试探位移,然后求解结构的变形和单元应力;
S4.对加载后的结构进行加固,建立加固后的有限元模型,然后释放加载位移,使结构在内力作用下回弹;
S5.测算回弹后框架的倾斜度,若回弹后的结构仍然倾斜,则修正获得新的试探位移,然后重复步骤S3-S5,直至框架回弹后倾斜度为零;
S6.当框架倾斜度达到纠偏要求后,所加载试探位移对应施加的纠偏力即为预测得到的纠偏荷载;
步骤中S3中所述对倾斜框架加载一个试探位移,然后求解结构的变形和单元应力,具体步骤如下:
S31.对于代表真实倾斜框架的有限元模型,在预先设定的加载节点,该模型与垂直位置的初始偏移等于Δ0;如果将试探位移Δ*应用于加载节点,那么获得的节点位移向量U将分为两部分:第一部分ua为纠偏时非加载节点的位移,第二部分ub为纠偏时加载点的位移;对于该过程,有限元方程为:
其中,α是一个大数,取1010,用于乘以刚度矩阵中加载自由度的对角线元素;
S32.根据求解结果,向量ub等于Δ*,ua是的非加载节点位移向量;各变形单元的应力向量σe在单元局部坐标系内计算:
其中,De为材料的本构矩阵,Be为构件的应变位移矩阵,ue为单元的位移向量,ue,a和ue,b分别是ua和ub中的单元位移分量;
步骤中S4中所述对加载后的结构进行加固,建立加固后的有限元模型,然后释放加载位移,使结构在内力作用下回弹,具体步骤如下:
S41.在加载后,结构的几何形态发生变化,为改善结构的性能需采用加固措施,因此需建立一个加固后的新模型:
其中,上标“'”表示加固后的结构矩阵,以表示与加固前矩阵的不同;
S42.结构的内力向量F通过组合各个变形单元的等效节点力得到:
其中,Be'为加固后结构的应变位移矩阵;
S43.求解获得加固后的结构在内力下的回弹位移u'a和u'b。
2.根据权利要求1所述的一种外闪木框架的纠偏荷载预测方法,其特征在于,步骤中S5中所述测算回弹后框架的倾斜度,若结构仍然倾斜,达不到纠偏目标,则修正获得新的试探位移,然后重复步骤S3-S5,直至框架回弹后倾斜度为零,具体步骤如下:
S51.纠偏时消除加载节点上预先存在的偏移Δ0,效果进行如下检测:
ub+u'b=-Δ0 (5)
若检测不满足要求,则修正获得一个新的试探位移:
Δ*=(1-ω)ub+ω(-Δ0-u'b) (6)
其中,ω是比例因子,取1.3~1.6;
S52.由于一次计算无法满足纠偏要求,需重复步骤S3-S5进行迭代计算,以比例因子ω连续调整试探位移Δ*,以保证倾斜框架在加固并释放荷载后能够回弹到无倾斜的理想形态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010829983.5A CN111967176B (zh) | 2020-08-18 | 2020-08-18 | 一种外闪木框架的纠偏荷载预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010829983.5A CN111967176B (zh) | 2020-08-18 | 2020-08-18 | 一种外闪木框架的纠偏荷载预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111967176A CN111967176A (zh) | 2020-11-20 |
CN111967176B true CN111967176B (zh) | 2024-03-01 |
Family
ID=73389200
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010829983.5A Active CN111967176B (zh) | 2020-08-18 | 2020-08-18 | 一种外闪木框架的纠偏荷载预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111967176B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113806942B (zh) * | 2021-09-22 | 2024-07-23 | 福州大学 | 一种获取倾斜框架有损变形的方法 |
CN115048703B (zh) * | 2022-06-27 | 2023-08-29 | 中建三局第一建设工程有限责任公司 | 迭代纠偏结构变形的高层建筑施工方法及高层建筑 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102995670A (zh) * | 2012-09-28 | 2013-03-27 | 上海宝冶工程技术有限公司 | 一种既有建筑物掏土纠偏施工监控方法 |
JP2016121475A (ja) * | 2014-12-25 | 2016-07-07 | 株式会社巴コーポレーション | 鉄骨部材の局部座屈抑制構造及び同構造を適用した構造骨組 |
CN109299540A (zh) * | 2018-09-25 | 2019-02-01 | 重庆大学 | 基于刚体准则的平面框架结构静力非线性分析方法 |
CN109753744A (zh) * | 2019-01-14 | 2019-05-14 | 哈尔滨工业大学 | 一种基于简化模型的框架结构模型修正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007052784A1 (ja) * | 2005-11-07 | 2007-05-10 | Keio University | 荷重伝達経路法に基づく数値構造解析装置 |
-
2020
- 2020-08-18 CN CN202010829983.5A patent/CN111967176B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102995670A (zh) * | 2012-09-28 | 2013-03-27 | 上海宝冶工程技术有限公司 | 一种既有建筑物掏土纠偏施工监控方法 |
JP2016121475A (ja) * | 2014-12-25 | 2016-07-07 | 株式会社巴コーポレーション | 鉄骨部材の局部座屈抑制構造及び同構造を適用した構造骨組 |
CN109299540A (zh) * | 2018-09-25 | 2019-02-01 | 重庆大学 | 基于刚体准则的平面框架结构静力非线性分析方法 |
CN109753744A (zh) * | 2019-01-14 | 2019-05-14 | 哈尔滨工业大学 | 一种基于简化模型的框架结构模型修正方法 |
Non-Patent Citations (1)
Title |
---|
邵满秋等.列车交汇过桥时预应力混凝土轨道箱形桥梁动力可靠度研究.《工程建设》.2019,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111967176A (zh) | 2020-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111967176B (zh) | 一种外闪木框架的纠偏荷载预测方法 | |
CN101615212B (zh) | 一种排桩-内支撑支护体系简化设计计算方法 | |
CN112081285B (zh) | 一种索结构预应力态拉索长度确定方法 | |
CN107274490A (zh) | 一种考虑材料损伤和钢筋作用的自由曲面结构形态创建方法 | |
CN102890739A (zh) | 一种用于建筑物钢桁架的应力卸除加固方法 | |
CN106503305B (zh) | 一种考虑损伤的自由曲面形态创建方法 | |
CN114330019B (zh) | 体内无粘结预应力波形钢腹板组合梁抗弯承载力计算方法 | |
Yuan et al. | Experimental and theoretical studies on the seismic performance of CFST battened built-up column piers | |
CN110837690A (zh) | 蜂窝结构非线性本构关系的建立方法、介质和设备 | |
CN109241629A (zh) | 一种确定除尘器箱体立柱轴压稳定承载力的方法 | |
CN106599489A (zh) | 空间圆管结构极限承载力分析的一次线弹性估算方法 | |
Wong et al. | Geometrically nonlinear analysis of elastic framed structures | |
JP7348509B2 (ja) | 連続梁の評価方法及び連続梁の評価プログラム | |
CN105781126B (zh) | 一种张弦结构被动建立预应力施工方法 | |
CN112502527A (zh) | 一种杆塔原位升高方法 | |
CN116226986A (zh) | 一种拉索无应力长度计算方法 | |
CN114945923A (zh) | 连续梁的评价方法、连续梁的评价程序及组合梁的评价方法 | |
CN115758812A (zh) | 一种劲性骨架混凝土拱桥拱圈外包混凝土工作段浇筑顺序的计算方法 | |
CN110472376B (zh) | 一种异形塔固结系统支撑件刚度识别的方法 | |
CN111241727B (zh) | 一种利用单跨有限元模型进行豪华邮轮极限强度计算的方法 | |
Koo et al. | A study on the stability of the single-layer latticed dome during erection using the step-up method | |
Zhang et al. | Structural Stability Analysis of Steel Tubular Scaffold with Couplers Based on Direct Analysis Method | |
CN102094459B (zh) | 整体倾斜相互链接结构体系及其施工方法 | |
CN114135001B (zh) | 一种超高层拉索结构建造方法 | |
CN104850532B (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 |