CN106407525A - 高拱坝三轴受力结构安全度简易算法 - Google Patents

高拱坝三轴受力结构安全度简易算法 Download PDF

Info

Publication number
CN106407525A
CN106407525A CN201610802280.7A CN201610802280A CN106407525A CN 106407525 A CN106407525 A CN 106407525A CN 201610802280 A CN201610802280 A CN 201610802280A CN 106407525 A CN106407525 A CN 106407525A
Authority
CN
China
Prior art keywords
concrete
point
state
stress
sigma
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
Application number
CN201610802280.7A
Other languages
English (en)
Other versions
CN106407525B (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.)
PowerChina Chengdu Engineering Co Ltd
Original Assignee
PowerChina Chengdu Engineering 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 PowerChina Chengdu Engineering Co Ltd filed Critical PowerChina Chengdu Engineering Co Ltd
Priority to CN201610802280.7A priority Critical patent/CN106407525B/zh
Publication of CN106407525A publication Critical patent/CN106407525A/zh
Application granted granted Critical
Publication of CN106407525B publication Critical patent/CN106407525B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及拱坝安全技术领域，尤其涉及是一种高拱坝三轴受力结构安全度简易算法。本发明公开的高拱坝三轴受力结构安全度简易算法能够较好地解决拱坝混凝土在考虑多轴受力状态下的如何确定坝体结构安全度的问题，具体包括步骤采集数据、拟合数据得到混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，以及将待分析拱坝的应力值σ1(i)、σ2(i)及σ3(i)带入关系曲线中进行插值从而确定拱坝三轴结构安全度K的步骤。本发明公开的方法混凝土σ3/fc与σ2/σ3、σ1/σ3的关系将拱坝三轴强度安全与单轴强度安全的表述联系起来，简化了拱坝三轴强度安全的描述。

Description

高拱坝三轴受力结构安全度简易算法
技术领域
本发明涉及拱坝安全技术领域，尤其涉及是一种高拱坝三轴受力结构安全度简易算法。
背景技术
在常规拱坝结构设计中，结构强度安全度的评价，均利用最大主应力与混凝土单轴强度之间的比值进行衡量，并不考虑第二、第三主应力对混凝土强度的影响，因此本质上属于单轴强度安全评价体系。而拱坝是一种受力结构相对复杂的坝型，是一种超静定结构，三个主应力之间相互影响，动态调整，且影响到混凝土的实际强度本身。因此，用混凝土单轴强度和最大主应力的比值来衡量拱坝的结构强度安全度仅仅是一个表观的设计安全度，并不真实反应结构的安全储备。
拱坝超静定结构在复杂受力状态下，较为充分地利用了混凝土天然的抗压能力，较大程度地发挥了材料的力学性能，然而由于其结构较为复杂，整个坝体受力状态复杂多变，到目前为止，尚没有一种很好的方法可以科学合理地解决拱坝混凝土在考虑多轴受力状态下的整体结构安全度。
发明内容
本发明所要解决的技术问题是：提供一种较好地解决拱坝混凝土在考虑多轴受力状态下的整体结构安全度难以计算的高拱坝三轴受力结构安全度简易算法。
本发明解决其技术问题所采用的高拱坝三轴受力结构安全度简易算法，包括如下步骤
A、采集数据，包括拱坝的混凝土单轴强度fc,及拱坝的点i对应的最大主应力值σ1(i)、中间主应力值σ2(i)和最小主应力值σ3(i)，其中i＝1，2，3，…；
B、通过采集的数据拟合出三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，及三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，其中σ1/σ3小于等于0.3；
C、用应力计算方法计算出待分析拱坝的应力值σ1(i)、σ2(i)及σ3(i)，其中i＝1，2，3，…，选取拱坝上的一点i，通过分析其σ1(i)、σ2(i)及σ3(i)的正负及大小判断该点所处状态，
若该点处于三轴压缩应力状态，则将当前应力点的σ1(i)/σ3(i)代入三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线进行插值计算，
若该点处于三轴拉压状态，则将当前应力点的σ1(i)/σ3(i)代入三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线进行插值计算，
从而获得相应的σ3(i)/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此获得该点在当前应力比状态下最大允许达到的[σ3(i)]；最后，通过公式K＝[σ3(i)]/σ3(i)计算出三轴结构安全度K。
进一步的是，步骤B中
三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线用插值公式表示为σ3/fc＝0.25+43×(σ1/σ3)1.3。
进一步的是，步骤B中
三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线用插值表格表示为
进一步的是，步骤B中
三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线用插值公式表示为
当σ2/σ3≤0.4时，
当σ2/σ3>0.4时，
进一步的是，步骤C中
若该点处于三轴拉压状态，利用当前应力点的σ2(i)/σ3(i)进行部位分区，如果σ2(i)/σ3(i)>0.4，则直接利用σ1(i)/σ3(i)进行插值，从而获取相应的σ3(i)/fc，如果σ2(i)/σ3(i)≤0.4，则利用σ1(i)/σ3(i)对斜率和截距同时进行插值，从而获取相应的σ3(i)/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此获得该应力点在当前应力比状态下最大允许达到的[σ3(i)]；最后，通过公式K＝[σ3(i)]/σ3(i)计算出拱坝三轴结构安全度K。
本发明的有益效果是：结合线弹性有限元或者拱梁分载法的计算成果，将拱坝结构任意一个点的三大主应力，用应力比的方式联系起来，从而反映出一个整体的受力形态；此外通过大量混凝土试验成果的总结，将混凝土多轴应力比与混凝土单轴强度的增减幅度联系起来，从而获得复杂受力状态下拱坝混凝土的三轴强度；最后利用本发明提供的算法，获得拱坝三轴受力结构安全度，因此，将拱坝复杂受力状态下结构的三轴结构强度安全问题，等效转换成单轴安全问题，从而用实际安全度一个系数获得拱坝在三轴受力状态下的安全度，一方面简化了真实三轴安全度评价的复杂程度，另一方面提供了比单轴强度更可靠的安全评价，从而较原规范安全度更接近真实的安全度；这种方法建立了大坝混凝土的强度与应力状态之间的联系，利用应力比的关系，将混凝土强度由单轴强度推广至多轴强度，且较为简单，快捷，难度较小，故容易为设计人员掌握。
附图说明
图1是三轴压缩情况下混凝土试验简化成果曲线；
图2是三轴拉压状态下混凝土试验简化成果曲线。
具体实施方式
下面结合附图对本发明作进一步说明。
本发明的高拱坝三轴受力结构安全度简易算法，包括如下步骤
A、采集数据，包括拱坝的混凝土单轴强度fc,及拱坝的点i对应的最大主应力值σ1(i)、中间主应力值σ2(i)和最小主应力值σ3(i)，其中i＝1，2，3，…；
B、通过采集的数据拟合出三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，及三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，其中σ1/σ3小于等于0.3；鉴于σ1/σ3小于等于0.3的试验数据较多，统计规律可靠，且大坝混凝土的应力比基本在这个范围内，故本发明要求σ1/σ3小于等于0.3，应力比超出这个范围的，可参考本发明的方法，数据可靠性无法保证，故不包括在本发明的权利要求之内。
C、用应力计算方法计算出待分析拱坝的应力值σ1(i)、σ2(i)及σ3(i)，其中i＝1，2，3，…，选取拱坝上的一点i，通过分析其σ1(i)、σ2(i)及σ3(i)的大小判断该点所处状态，
若该点处于三轴压缩应力状态，则将当前应力点的σ1(i)/σ3(i)代入三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线进行插值计算，从而获得相应的σ3(i)/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此获得该点在当前应力比状态下最大允许达到的[σ3(i)]；最后，通过公式K＝[σ3(i)]/σ3(i)计算出拱坝三轴结构安全度K，
若该点处于三轴拉压状态，则将当前应力点的σ1(i)/σ3(i)代入三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线进行插值计算，从而获得相应的σ3(i)/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此获得该点在当前应力比状态下最大允许达到的[σ3(i)]；最后，通过公式K＝[σ3(i)]/σ3(i)计算出拱坝三轴结构安全度K。
其中，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线的表现形式有多种，可以是如图1和图2所示的混凝土试验简化成果曲线，也可以是插值公式或者插值表格。
具体，经过统计分析，可以得出三轴压缩应力状下的混凝土试验简化成果曲线如图1所示，从图中可以看出，此种情况下，σ3/fc的值与σ1/σ3有关，而与σ2/σ3的大小相关性不大，可以忽略不计。也可以将σ3/fc与σ1/σ3之间的关系用插值公式表示：
σ3/fc＝0.25+43×(σ1/σ3)1.3 (1)
或者用插值表格表示：
在评判拱坝的安全强度前，首先基于多拱梁法或线弹性有限元的应力计算得到拱坝各点的应力值；在评判拱坝的安全强度时，先读取一点主应力状态，对主应力进行排序，使得σ1>σ2>σ3，拉应力为正，压应力为负，并根据主应力，判断该点处于三轴拉压应力状态还是三轴压缩应力状态，若σ1、σ2、σ3均为负，则此点处于三轴压缩状态，接下来首先读如图1所示的三轴压缩情况下混凝土试验简化成果曲线或者利用公式1，得到当前应力点所对应的σ1/σ3两个最接近的值，以及相应的σ3/fc的值，然后利用当前应力点的σ1/σ3进行插值，从而获取相应的σ3/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此可以获得该应力点在当前应力比状态下最大允许达到的[σ3]；最后，比较[σ3]和σ3，通过公式K＝[σ3]/σ3求得拱坝三轴结构安全度K。
具体，经过统计分析，可以得出三轴拉压状态下，混凝土试验简化成果曲线如图2所示，从图中可以看出，此种情况下，当σ2/σ3≤0.4时，σ3/fc的值与σ1/σ3和σ2/σ3均有关，当σ1/σ3的值一定时，σ3/fc的值与σ2/σ3几乎呈线性关系；当σ2/σ3>0.4时，σ3/fc的值仅与σ1/σ3有关。也可以将σ3/fc与σ1/σ3之间的关系用插值公式表示：
当σ2/σ3≤0.4时，
当σ2/σ3>0.4时，
同样，在评判拱坝的安全强度前，首先基于多拱梁法或线弹性有限元的应力计算得到拱坝各点的应力值；在评判拱坝的安全强度时，先读取一点主应力状态，对主应力进行排序，使得σ1>σ2>σ3，拉应力为正，压应力为负，并根据主应力，判断该点处于三轴拉压应力状态还是三轴压缩应力状态，若σ1为正，σ2和σ3中一个或者两个均为负，则此点处于三轴拉压状态，接下来首先读取如图2所示的三轴拉压状态下混凝土试验简化成果曲线，并利用当前应力点的σ2/σ3进行部位分区，如果σ2/σ3>0.4，则直接利用σ1/σ3进行插值，从而获取相应的σ3/fc，如果σ2/σ3≤0.4，则利用σ1/σ3对斜率和截距同时进行插值，从而获取相应的σ3/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此可以获得该应力点在当前应力比状态下最大允许达到的[σ3]；最后，比较[σ3]和σ3，通过公式K＝[σ3]/σ3求得拱坝三轴结构安全度K。
当然，在绘图时，也可以选择不同的σ1/σ3值及不同数量的σ1/σ3进行绘制，例如，在三轴压缩应力状态下，可以选择σ1/σ3＝[0.05 0.1 0.15 0.2 0.25 0.3]所对应的σ3/fc的值进行绘图；在三轴拉压应力状态下，可以选择σ1/σ3＝[-0.05-0.15-0.25-0.35-0.45-0.80-0.95]所对应的σ3/fc的值进行绘图。

Claims (5)

1.高拱坝三轴受力结构安全度简易算法，其特征在于：包括如下步骤
A、采集数据，包括拱坝的混凝土单轴强度fc,及拱坝的点i对应的最大主应力值σ1(i)、中间主应力值σ2(i)和最小主应力值σ3(i)，其中i＝1，2，3，…；
B、通过采集的数据拟合出三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，及三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线，其中σ1/σ3小于等于0.3；
C、用应力计算方法计算出待分析拱坝的应力值σ1(i)、σ2(i)及σ3(i)，其中i＝1，2，3，…，选取拱坝上的一点i，通过分析其σ1(i)、σ2(i)及σ3(i)的正负及大小判断该点所处状态，
若该点处于三轴压缩应力状态，则将当前应力点的σ1(i)/σ3(i)代入三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线进行插值计算，
若该点处于三轴拉压状态，则将当前应力点的σ1(i)/σ3(i)代入三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线进行插值计算，
从而获得相应的σ3(i)/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此获得该点在当前应力比状态下最大允许达到的[σ3(i)]；最后，通过公式K＝[σ3(i)]/σ3(i)计算出三轴结构安全度K。
2.根据权利要求1所述的高拱坝三轴受力结构安全度简易算法，其特征在于：步骤B中三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线用插值公式表示为
σ3/fc＝0.25+43×(σ1/σ3)1.3。
3.根据权利要求1所述的高拱坝三轴受力结构安全度简易算法，其特征在于：步骤B中三轴压缩状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线用插值表格表示为
4.根据权利要求1所述的高拱坝三轴受力结构安全度简易算法，其特征在于：步骤B中三轴拉压状态下，混凝土σ3/fc与σ2/σ3、σ1/σ3的关系曲线用插值公式表示为
当σ2/σ3≤0.4时，
当σ2/σ3>0.4时，
5.根据权利要求4所述的高拱坝三轴受力结构安全度简易算法，其特征在于：步骤C中若该点处于三轴拉压状态，利用当前应力点的σ2(i)/σ3(i)进行部位分区，如果σ2(i)/σ3(i)>0.4，则直接利用σ1(i)/σ3(i)进行插值，从而获取相应的σ3(i)/fc，如果σ2(i)/σ3(i)≤0.4，则利用σ1(i)/σ3(i)对斜率和截距同时进行插值，从而获取相应的σ3(i)/fc；其次根据该应力点所处的部位读取相应的混凝土单轴强度fc，由此获得该应力点在当前应力比状态下最大允许达到的[σ3(i)]；最后，通过公式K＝[σ3(i)]/σ3(i)计算出三轴结构安全度K。
CN201610802280.7A 2016-09-05 2016-09-05 高拱坝三轴受力结构安全度简易计算方法 Active CN106407525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610802280.7A CN106407525B (zh) 2016-09-05 2016-09-05 高拱坝三轴受力结构安全度简易计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610802280.7A CN106407525B (zh) 2016-09-05 2016-09-05 高拱坝三轴受力结构安全度简易计算方法

Publications (2)

Publication Number Publication Date
CN106407525A true CN106407525A (zh) 2017-02-15
CN106407525B CN106407525B (zh) 2019-07-23

Family

ID=57999445

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610802280.7A Active CN106407525B (zh) 2016-09-05 2016-09-05 高拱坝三轴受力结构安全度简易计算方法

Country Status (1)

Country Link
CN (1) CN106407525B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106906785A (zh) * 2017-03-13 2017-06-30 中国电建集团成都勘测设计研究院有限公司 基于结构强度、基础抗滑和整体稳定的拱坝安全评价方法
CN106996097A (zh) * 2017-05-24 2017-08-01 中国电建集团成都勘测设计研究院有限公司 一种拱坝结构的安全评价方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1196475A (zh) * 1997-04-15 1998-10-21 董伟 拱坝三维变形全自动追迹监测方法
CN104318030A (zh) * 2014-10-31 2015-01-28 中国电建集团成都勘测设计研究院有限公司 基于变形等效的拱坝坝基综合变形模量自动化计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1196475A (zh) * 1997-04-15 1998-10-21 董伟 拱坝三维变形全自动追迹监测方法
CN104318030A (zh) * 2014-10-31 2015-01-28 中国电建集团成都勘测设计研究院有限公司 基于变形等效的拱坝坝基综合变形模量自动化计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PENG LIN 等: "Failure and overall stability analysis on high arch dam based on DFPA code", 《ENGINEERING FAILURE ANALYSIS》 *
吴晓彬 等: "基于不同混凝土强度准则的高拱坝脆性破坏范围研究", 《水利水电技术》 *
钱向东 等: "混凝土的多轴强度安全系数及其应用", 《河海大学学报（自然科学版）》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106906785A (zh) * 2017-03-13 2017-06-30 中国电建集团成都勘测设计研究院有限公司 基于结构强度、基础抗滑和整体稳定的拱坝安全评价方法
CN106996097A (zh) * 2017-05-24 2017-08-01 中国电建集团成都勘测设计研究院有限公司 一种拱坝结构的安全评价方法
CN106996097B (zh) * 2017-05-24 2019-01-01 中国电建集团成都勘测设计研究院有限公司 一种拱坝结构的安全评价方法

Also Published As

Publication number Publication date
CN106407525B (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
Owen et al. Development of a criterion method to determine peak mechanical power output in a countermovement jump
CN105912772B (zh) 一种集频率、振型和频响函数为一体的模型修正方法
CN105004504B (zh) 复合绝缘子的起振风速评估方法及选型方法
CN105160067A (zh) 风机变桨轴承与轮毂和叶片连接螺栓强度性能的校核方法
CN104536941B (zh) 一种基于Tikhonov正则化的频域载荷识别方法
WO2008063527A3 (en) Detection of deception and truth-telling using fmri of the brain
Hull et al. Optimization of pedaling rate in cycling using a muscle stress-based objective function
CN109063310B (zh) 一种考虑振动模态偏角的高层建筑hffb风振分析方法
CN104899372B (zh) 结合仿真与可靠性分析高加速应力筛选试验剖面构造方法
CN106407525A (zh) 高拱坝三轴受力结构安全度简易算法
Holmes Along-and cross-wind response of a generic tall building: Comparison of wind-tunnel data with codes and standards
CN107290150A (zh) 发动机悬置总成垂向多级谱试验台架及其试验方法
CN110032776A (zh) 基于模态解耦的叶片动应变场与位移场重构方法及其系统
CN108984879A (zh) 一种串联多自由度系统的位移频率响应计算方法
CN116577051A (zh) 一种虑及损伤等效的随机振动疲劳寿命分析方法
CN115270296A (zh) 一种商用车驾驶室疲劳耐久性分析方法及系统
EP1928317A1 (en) Evaluation method of a sports performance
CN105445037B (zh) 一种汽车碰撞试验用假人颈部替代装置
CN107727350B (zh) 微纳卫星矢量振动试验方法
CN108489696B (zh) 一种消除频响函数中弹簧附加刚度的方法
Missirlis et al. Stress analysis of the aortic valve during diastole: important parameters
Pigatto et al. A new crank arm based load cell, with built-in conditioning circuit and strain gages, to measure the components of the force applied by a cyclist
CN106596131A (zh) 一种汽车白车身模态测试中激励力的位置、幅值以及方向的确定方法
CN115408909A (zh) 基于流固耦合和疲劳损伤的风机叶片形变分析方法及装置
CN110222405A (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