CN110866300B - 大型建筑的裂缝预测方法 - Google Patents
大型建筑的裂缝预测方法 Download PDFInfo
- Publication number
- CN110866300B CN110866300B CN201911117493.6A CN201911117493A CN110866300B CN 110866300 B CN110866300 B CN 110866300B CN 201911117493 A CN201911117493 A CN 201911117493A CN 110866300 B CN110866300 B CN 110866300B
- Authority
- CN
- China
- Prior art keywords
- large building
- modulus
- displacement
- model
- velocity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 230000008859 change Effects 0.000 claims abstract description 19
- 239000013598 vector Substances 0.000 claims description 25
- 238000013016 damping Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 5
- 238000009966 trimming Methods 0.000 claims description 4
- 238000010008 shearing Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 5
- 230000007423 decrease Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种大型建筑的裂缝预测方法,包括以下步骤:步骤S1:建立大型建筑的模型;步骤S2:设定条件参数;步骤S3:设置边界条件和初始条件;步骤S4:计算预定时间内的应力变化值及位置变化值;步骤S5:计算整合位移。如此能够预测大型建筑的裂缝的变化趋势。
Description
技术领域
本发明涉及建筑技术领域,特别是一种大型建筑的裂缝预测方法。
背景技术
大型建筑,如水坝、城墙等,由于建筑面积及体积大、跨度较大,在设定使用年限时除了评估建筑本身质量外,还需要评估地壳运动给其带来的影响,如由于非均匀的地壳运动产生的剪切应力产生的裂缝。经研究发现,大型建筑的裂缝往往在靠近底部的位置最大,往上逐渐减小,这种裂缝就是非均匀的地壳运动产生的,远比风晒雨淋或者普通外力破坏产生的裂缝更为严重,也对建筑影响更大。本方案即是通过研究非均匀的地壳运动随着时间变化在大型建筑上产生的裂缝的趋势及变化,得出大型建筑裂缝的产生规律及变化趋势,预测裂缝的产生及变化程度,为大型建筑的使用年限的设定及保护提供参考或依据。
发明内容
有鉴于此,本发明提供了一种能够预测大型建筑的裂缝的变化趋势的大型建筑的裂缝预测方法,以解决上述问题。
一种大型建筑的裂缝预测方法,包括以下步骤:
步骤S1:建立大型建筑的模型,根据大型建筑的实际尺寸建立一个等比例的三维均质各向同性的模型;
步骤S2:设定条件参数;
步骤S3:设置边界条件和初始条件;
应变张量εij满足:
其中,E是弹性模量,v是泊松比,σkk为主应力之和;δij是克罗内克变量;
步骤S5:计算整合位移。
进一步地,所述条件参数包括弹性体积模量、内聚力、弹性切边模量、抗拉强度、内摩擦角、剪胀角、局部阻尼及密度。
进一步地,所述弹性体积模量为2.0E+7Pa,内聚力为1.0E+5Pa,弹性切边模量为1.0E+7Pa,抗拉强度为1.0E+5Pa,内摩擦角为25°,剪胀角为0°,局部阻尼为0.05,密度为1950Kg/m3。
与现有技术相比,本发明的大型建筑的裂缝预测方法包括以下步骤:步骤S1:建立大型建筑的模型,根据大型建筑的实际尺寸建立一个等比例的三维均质各向同性的模型;步骤S2:设定条件参数;步骤S3:设置边界条件和初始条件;步骤S4:计算预定时间内的应力变化值及位置变化值;应力张量σij满足:其中,xj是坐标向量,ρ是体积密度,gi是重力加速度,是速度矢量,ui是位移矢量,t为时间,i是坐标值;应变速度张量满足:其中,xi是不同方向的坐标向量,是不同方向的速度矢量;应变张量εij满足:其中,E是弹性模量,v是泊松比,okk为主应力之和;δij是克罗内克变量;步骤S5:计算整合位移。如此能够预测大型建筑的裂缝的变化趋势,为大型建筑的使用年限的设定及保护提供参考或依据。
附图说明
以下结合附图描述本发明的实施例,其中:
图1为本实施例中的城墙的概念模型示意图。
图2为图1中的城墙的不同位置的剖面随着时间变化的Smax分布图。
图3为图1中的城墙的不同位置的剖面随着时间变化的应力状态分布图。
图4为图1中的城墙的不同位置的剖面随时间变化的y向位置等值线图。
图5为图1中的城墙的不同位置的剖面随时间变化的x向位置等值线图。
图6为本发明的大型建筑的裂缝预测方法的流程示意图。
具体实施方式
以下基于附图对本发明的具体实施例进行进一步详细说明。应当理解的是,此处对本发明实施例的说明并不用于限定本发明的保护范围。
请参考图6,本发明提供的大型建筑的裂缝预测方法包括以下步骤:
步骤S1:建立大型建筑的模型。
建立一个三维均质各向同性的模型,可根据大型建筑的实际尺寸设定一个等比例的模型,如x轴,y轴和z轴方向的长度分别为Lx=1000m,Ly=10m,Lz=15m,如图1所示。
步骤S2:设定后续计算过程中所需的条件参数,如弹性体积模量、内聚力、弹性切边模量、抗拉强度、内摩擦角、剪胀角、局部阻尼及密度。本实施方式中,以城墙为例,各种条件参数如表1所示:
表1设定条件参数
步骤S3:设置边界条件和初始条件。如:
1.z方向的位置矢量的边界条件:
uz=0在z=0m
来自模型的底部在垂向上的位移受约束,因为与水平方向的变形相比,底部垂向上的变形是可以忽略的。
2.x方向的位置矢量的边界条件:
ux=0在x=0m和x=Lx
模型在x方向的边界,即x=0和x=Lx,水平位移预设是零。模型的顶部边界,即城墙的顶部,默认为自由边界。
3.y方向的速度矢量的局部边界条件:
模型在y轴方向上承受的非均匀地壳运动则通过GPS测定结果。在模型y方向上y=0m的边界,即模型的前边界,于边界的右半部分即x=500m至x=1000m的边界,且z=0处,施加顺y轴方向的速度场,即即在y=0和x在Lx和Lx/2的范围内的速度场。如图1所示。
4.y方向的速度矢量的局部边界条件:
5.应力张量的初始条件:
σij=σo当t=0
σo是初始状态下的初始应力场。
上述边界条件和初始条件也可根据实际采用不同的值。
步骤S4:计算预定时间内的应力变化值及位置变化值。
大型建筑的模型的力学性质受控于弹性力学基本方程,包括几何方程,动力方程和本构方程。根据弹性理论和运动方程,模型的力学性质主要受控于动力方程:
其中σij是应力张量(单位面积上的应力)[N/L2],xj是坐标向量[L],ρ是体积密度,gi是重力加速度,是速度矢量[L/a],ui是位移矢量[L],t为时间[a],i(x,y,z)是坐标值。L是长度单位,N是力,a是时间单位:年。
速度矢量如上述边界条件所示,体积密度如上述条件参数所示。重力加速度已知。
速度引起的应变通过几何方程描述:
采用线性孔隙弹塑性理论描述模型变形的应力-应变关系:
其中,εij是应变张量,E是弹性模量[N/L2],v是泊松比,σkk为主应力之和(σkk=σx+σy+σz);δij是克罗内克变量(Kronecker),i=j时,δij=0,i≠j时,δtf=1。
具体计算时,可运用FLAC3D软件,结合摩尔-库伦本构模型,摩尔-库伦本构模型常被用来模拟材料的变形,至达到变形极限,最后破裂的过程。本实施例中模型网格为60*20*30,共50375个节点;动力计算时间步设置为0.001年。局部阻尼根据FLAC3D手册是地质材料中最小临界阻尼比的一般值,将其值定为0.05。FLAC3D是基于三维有限差分显示算法的软件,在静力计算、动力计算、固液耦合方面广泛应用。
图2展示了y方向上y=0,y=2.5m,y=5m处的三个xz剖面的最大主应力(Smax)随时间的变化,yr表示年份,图中展示了时间为0年、1年和2年的变化趋势,最大主应力以颜色深浅表示。初始应力场主要受控于重力,因此Smax在x方向上是常值,但是在z方向上是变化的(图2中第一列);顶部是正值,即拉张状态,底部逐渐为负值,即压缩状态,但是由于模型的底部垂向位移和模型的两个侧向边界的x方向的水平位移是固定的,因此应力在y方向是变化的。随着时间的增加,由于连续的速度场加载在模型底部的不同边界,Smax越来越在剖面的中下部集中(图2的中间行及底部行)。
由于模型的边界(y=0和y=Ly)受速度场影响最大,且这种影响随着远离边界而逐渐减弱,因此Smax随着y的增加逐渐减低,在y=5m处最小。可以看到,最大拉张应力在剖面底部的中间区域形成(图2的左下角的底部)。根据摩尔-库伦本构模型破坏准则判定,该区域可能发生变形,并随着z的增加,这种变形效应逐渐减弱,正如图3模拟的应力状态结果。这种变形主要是来自于模型边界速度场的施加导致的剪切变形。如此使得Smax分布及塑性变形区分布,都分布于建筑的底部,且随着高度的增加而减少。产生的剪切变形以较浅的颜色表示。
图4主要展示y方向的y=0,y=2.5m,y=5m处的三个不同xz剖面的y向位移随时间的变化。y向的位置初始位移值设为0,随着速度场施加时间的增加,位移由0增至近1cm。这种剪切变形出现在剖面中部,并随高度逐渐减小(图4的中间行和底部行)。在0≤x≤Lx/2部分的位移是负值,主要是因为模型左半部分的速度场方向与y轴相反,在Lx/2≤x≤Lx部分是正值,因为模型的右半部分的速度场与y轴是一致的。产生的位移量以颜色深浅表示。
图5展示的是相同剖面的x向位移随时间的变化。同样地,x方向的初始位移设置为0,随着速度场加载时间的增加,x向位移由0增至0.7cm。拉张变形区的出现并随高度的减弱的原因与上述原因一致。与上述情况一致,位移的正负主要看位移的方向与坐标轴是否一致。0≤x≤Lx/2部分是负值,是由于模型左半部分的x向位移与x方向相反;Lx/2≤x≤Lx是正值,是由于模型右半部分的x向位移与x方向一致,并且模型侧面的x向位移是固定的,但是模型内部的x向位移是自由的。产生的位移量以颜色深浅表示。
步骤S5:计算整合位移。不同方向的位移,如y向位移和x向位移,应该整合为总位移才是非均匀地壳运动产生的剪切应力导致的实际变形量,通过角度函数整合为整合位移即实际位移。模型底部的中间部分的位移最大且随高度减小,且由于速度场的施加,有一个很明显的剪切变形面。
如此能够模拟并预测大型建筑的裂缝的产生及变化趋势,能够预测未来一段时间裂缝的尺寸及延伸方向。用于可根据预测的结果为大型建筑的使用年限的设定及保护提供参考或依据。如裂缝的最大宽度超过预定值时,该大型建筑应报废拆除。
上述实施例仅作为示例,实际应用时各种数值可根据实际情况设定,通过本方法预测出未来一段时间的裂缝的变化趋势,从而为大型建筑的使用年限的设定及保护提供参考或依据。
与现有技术相比,本发明的大型建筑的裂缝预测方法包括以下步骤:步骤S1:建立大型建筑的模型;步骤S2:设定条件参数;步骤S3:设置边界条件和初始条件;步骤S4:计算预定时间内的应力变化值及位置变化值;步骤S5:计算整合位移。如此能够预测大型建筑的裂缝的变化趋势。
以上仅为本发明的较佳实施例,并不用于局限本发明的保护范围,任何在本发明精神内的修改、等同替换或改进等,都涵盖在本发明的权利要求范围内。
Claims (2)
1.一种大型建筑的裂缝预测方法,其特征在于:包括以下步骤:
步骤S1:建立大型建筑的模型,根据大型建筑的实际尺寸建立一个等比例的三维均质各向同性的模型;
步骤S2:设定条件参数,所述条件参数包括弹性体积模量、内聚力、弹性切边模量、抗拉强度、内摩擦角、剪胀角、局部阻尼及密度;
步骤S3:设置边界条件和初始条件;
步骤S4:计算预定时间内的应力变化值及位置变化值;应力张量σij满足:
应变张量εij满足:
其中,E是弹性模量,v是泊松比,okk为主应力之和;δij是克罗内克变量;步骤S5:计算整合位移。
2.如权利要求1所述的大型建筑的裂缝预测方法,其特征在于:所述弹性体积模量为2.0E+7Pa,内聚力为1.0E+5Pa,弹性切边模量为1.0E+7Pa,抗拉强度为1.0E+5Pa,内摩擦角为25°,剪胀角为0°,局部阻尼为0.05,密度为1950Kg/m3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911117493.6A CN110866300B (zh) | 2019-11-15 | 2019-11-15 | 大型建筑的裂缝预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911117493.6A CN110866300B (zh) | 2019-11-15 | 2019-11-15 | 大型建筑的裂缝预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110866300A CN110866300A (zh) | 2020-03-06 |
CN110866300B true CN110866300B (zh) | 2022-11-25 |
Family
ID=69654577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911117493.6A Active CN110866300B (zh) | 2019-11-15 | 2019-11-15 | 大型建筑的裂缝预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110866300B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111651816A (zh) * | 2020-05-25 | 2020-09-11 | 上海建工一建集团有限公司 | 裂缝识别方法及建筑结构安全性监测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE0602417L (sv) * | 2006-11-14 | 2008-05-15 | Ragnar Slunga | Sätt att bestämma den spänningstensor som har utlöst ett jordskalv |
CN105093309A (zh) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | 基于构造应力分析进行裂缝检测及有效性验证的方法 |
CN105319603A (zh) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | 致密砂岩储层复杂网状裂缝的预测方法 |
CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730596B (zh) * | 2015-01-25 | 2017-11-07 | 中国石油大学(华东) | 一种基于多尺度因素约束的离散裂缝建模方法 |
-
2019
- 2019-11-15 CN CN201911117493.6A patent/CN110866300B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE0602417L (sv) * | 2006-11-14 | 2008-05-15 | Ragnar Slunga | Sätt att bestämma den spänningstensor som har utlöst ett jordskalv |
CN105093309A (zh) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | 基于构造应力分析进行裂缝检测及有效性验证的方法 |
CN105319603A (zh) * | 2015-11-06 | 2016-02-10 | 中国石油大学(华东) | 致密砂岩储层复杂网状裂缝的预测方法 |
CN108868748A (zh) * | 2018-04-28 | 2018-11-23 | 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 |
Non-Patent Citations (2)
Title |
---|
地应力与构造裂缝参数间的定量关系;季宗镇等;《石油学报》;20100115(第01期);全文 * |
裂缝宽度预测的有限元数值模拟;练章华等;《天然气工业》;20010528(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110866300A (zh) | 2020-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A continuous contact force model for impact analysis in multibody dynamics | |
Ma et al. | Enhancement factors for grounded ice and ice shelves inferred from an anisotropic ice-flow model | |
CN102720295B (zh) | 一种基于索穹顶张拉和承载全过程分析的预应力确定方法 | |
Huétink et al. | Progress in mixed Eulerian‐Lagrangian finite element simulation of forming processes | |
Tvergaard | Effect of stress-state and spacing on voids in a shear-field | |
Wu et al. | On forming limit stress diagram analysis | |
Dong et al. | Solid-mechanics finite element simulations of the draping of fabrics: a sensitivity analysis | |
CN110866300B (zh) | 大型建筑的裂缝预测方法 | |
Kukielka et al. | Incremental modelling and numerical solution of the contact problem between movable elastic and elastic/visco-plastic bodies and application in the technological processes | |
Chakraborty et al. | Experimental investigation of surface pressure on ‘+’plan shape tall building | |
Burchitz | Improvement of springback prediction in sheet metal forming | |
Robert et al. | Comparison between incremental deformation theory and flow rule to simulate sheet-metal forming processes | |
Marvi-Mashhadi et al. | Surrogate models of the influence of the microstructure on the mechanical properties of closed-and open-cell foams | |
Xiao et al. | Estimation of peening effects of random and regular peening patterns | |
CN115169181A (zh) | 一种考虑局部劣化的土体边坡稳定性分析方法及系统 | |
Liu et al. | Springback behaviors of extruded 6063 aluminum profile in subsequent multi-stage manufacturing processes | |
CN106777484A (zh) | 一种表面裂纹扩展形状的计算方法 | |
CN113239427A (zh) | 既有桩基抵抗侧向变形力学性能增强方法 | |
Gao | New expanding cavity model for indentation hardness including strain-hardening and indentation size effects | |
Wang et al. | Numerical simulation of flexible micro-bending processes with consideration of grain structure | |
Nakamachi et al. | Process metallurgy design of aluminum alloy sheet rolling by using two-scale finite element analysis and optimization algorithm | |
Smith et al. | An approach towards the reduction of sink marks in sheet molding compound | |
CN112883477B (zh) | 基于双向受力模型的墙体加固性能评价方法及装置 | |
Utsunomiya et al. | Evolution of matt surface topography in aluminium pack rolling. Part I: model development | |
Han et al. | Design of wood-like metallic material using metal sheet architecture |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20221102 Address after: 5 / F, 277 Huqingping Road, Minhang District, Shanghai, 201105 Applicant after: SHANGHAI HUANLIAN ECOLOGICAL TECHNOLOGY Co.,Ltd. Address before: 201417 Room 1601, Building 5, No. 215, Lianhe North Road, Fengxian District, Shanghai Applicant before: Shanghai Pinshi Environmental Technology Co.,Ltd. Applicant before: SHANGHAI HUANLIAN ECOLOGICAL TECHNOLOGY Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |