弹性参数反演(基本)演示文件修改版.pptVIP

反演该弹性模量的逆问题可归结为下式约束条件下，以上式的综合残差Re为极小化目标函数的n维数学优化问题。 所选优化方法：复合形直接搜索法 复合形法是解决有约束条件的非线性优化问题的有效方法之一，是一种直接法。 由于该法可以在全部可行域内寻找，性能较稳定、可靠、编程简单方便，对解决工程问题有较强的适应性，所以，在工程中应用较多。 复合形法的基本思想来自单纯形法，这种方法是在n维受非线性约束的设计空间内，由Kn＋1个顶点（常取K=2n）构成多面体，称之为复合形。 然而对复合形的各顶点函数值逐一进行比较，不断丢掉使目标函数值最大的顶点，代之既能使目标函数有所改进，又满足约束条件的断点，逐步调向最优点。 根据复合形法原理编制相应反演分析程序，其中的非线性有限元计算子程序考虑随水位的不断上升，荷载的逐级加大来模拟蓄水期水荷载对大坝的真实的加载过程，期望能更好的反演出坝体的真实弹模 反演分析流程图见下图。 参照坝体各材料的室内实验提供的结果，坝体浆砌块石弹模的大小参考设计的弹模E0取搜索区间为（0.8e10，1.2e10）Pa，ε0 的搜索范围取在（0.002，0.004）。 5、反演成果及分析 ?（1）.利用位移资料反演参数结果 基于上述参照混凝土本构关系而假设的本构关系，对浆砌石坝体弹模E0、及ε0反演。 反演计算结果为该混凝土浆砌块石重力坝的初始弹模估计值为E0=1.11E10Pa，本构关系为 大坝初始弹模及ε0反演估计值与迭代次数T之间的收敛过程线分别如下图所示 初始弹模、ε0反演值与迭代次数关系图 计算收敛速度较快，大约经过24步的迭代计算，问题得到了解决 此时的E0=1.11E10Pa，ε0＝0.0034，此时目标函数F的值为0.462305E-6m2。 利用该本构关系对坝体砌石材料采用非线性弹性，对大坝其他材料采用弹性反演出蓄水期各个部分的初始弹性模量值。 反演结果见下表。 从表中可以看出，在－20％～＋20％的区间内搜索得到的结果和设计采用值相比较变化不大，由此可见复合形反分析结果是令人满意的。 计算的收敛速度是较快的，大约经过24步的迭代计算，问题得到解决，得到比较稳定的值，此时ΔF=0.4626126×10-6mm2，E1=2.63×1010Pa，E2=0.85×1010Pa，E3=1.53×1010Pa，E4=2.13×1010Pa。 （2）.利用应变资料反演参数成果 本坝内部埋设众多应变计，用以监测大坝的变形。 故可对埋设于4#坝段混凝土坝体内部的S4和S5两组五向应变计的监测资料进行分析处理，并分离出水压分量所产生的数值，用以反演计算。 待反演混凝土坝体弹模E取上、下20％的搜索区间，即弹模E在(2.0e10,3.0e10)Pa取值。 除坝体混凝土参数待反演外，其它各部分的参数值采用前面位移参数反演的结果，反演结果为2.64e10，反演结果与位移资料反演结果很接近。 6、非线性弹性与线弹性模型反演精度的比较 由于非线性理论本身复杂性，在以往很多工程中，人们很难都用非线性的理论来考虑问题，而是把问题作了简化。 用弹性问题来代替非线性的问题，对很多问题采用了理性模型考虑，与实际情况存在一定的差异，结果自然也不太一样。 按前面非线性弹性反演计算的思路，只是把非线性弹性的本构关系部分换成线弹性的本构关系对大坝各物理参数进行优化反演，将二者反演结果对比，对比结果如下。 （1）.在相同情况下，非线性弹性与线弹性反演收敛时，目标函数值分别为ΔF1=0.501506E-6和ΔF2＝0.5429E-6。 （2）.按非线性弹性和线弹性反演出的参数值，通过有限元计算，在不同水位情况下产生的位移值与监测资料值对比如下： 从对比结果可以看出，利用非线性弹性本构关系反演出来的结果与按线弹性反演的结果，无论是从目标函数还是从水压与位移关系都可以看出非线性弹性反演结果比较理想。 7、小结 基础：位移观测资料和应变计观测资料 利用该资料对研究对象的本构关系、分区物理力学参数进行反演分析； 按非线性弹性反演比线弹性模型反演精度高、效果好 岩土数值分析岩土工程参数反分析方法 三、物理力学参数反演分析 如前所述，岩土工程优化反分析最常用，也是岩土工程反分析研究中一个最重要、最实用的研究方向。 岩土工程优化反分析本质上看是一个典型的复杂非线性函数优化问题，该目标优化函数是一个高度复杂的非线性多峰函数。 而普通的优化反分析中多次调用正分析的特点使得整个算法的计算效率很低。 要解决这类复杂优化函数，采用传统优化方法往往难以奏效。采用全局优化算法是解决这个问题的理想途径。 结合有限元 研究对象不同，反分析的参数不同，具体方法有所差异 优化方法 以某坝及岩基的弹性参数反演为例 前期要求：充足的监测