反分析计算法.pdfVIP

为 工黼 杆 反 分 析 计 算 法 6) 马选荣 高效伟 (宁夏大学) 摘 要 ： 反分析计算法是 目前工程计算中一种实用且经济的计算法。本文 对工程 中常用的几种反分析计算作了舟绍 ，并 比较 了它们在实际使用中的优缺 点 ． 一 、 概 述 在岩土工程问题 中，一般都是已知边界条件、初始条件和材料的力学参数来求解区域内 的应力分布、位移分布等，这种分析过程通常称之为正分析。但在某些实际问题中，边界条件 和力学参数的确定却是很棘手的，例如，地下工程问题中，地应力的测量费时、费力、费财、技 术要求高 ，而测量的结果往往是离散性很大，要得到比较可靠的地应 力值 ，就必须进行大量 的测量，这是一般工程投资所担负不起的；因而近年来人们一直在寻找既经济又实用的确定 地应 力和力学参数的方法。反分析计算法就是这种方法中较好的一种 。 所谓反分析计算，就是由在施工中实测到的一组位移或应力，反过来计算原始地应力或 力学参数。由于位移的测量既经济又方便，而且能反映具体问题的实际情况，所以近年来位 移反分析计算法得到了广泛的应用 “ 。 反分析计算法主要分为两大类，一类是直接逼近法 ，另一类称为逆算法。直接逼近法的 计算过程与正分析中的计算过程相似，即先假设一组地应力和力学参数的值，然后按照正分 析法的步骤求解位移值 ，再利用最小二乘法原理 ，由计算值与实测值的偏差对所设的地应力 和力学参数进行修正 ，重复上述步骤直到修正前后的地应力值和力学参数相差很小为止。直 接逼近法具有应用方便，编制程序简单等优点，但计算时间长、甚至会出现不收敛的现象 。 逆算法是把基本方程倒过来求解，即将地应力和力学参数作为未知量写出基奉方程。这 种方法的特点是计算速度 ’，收敛性好，缺点是编制程序时工作量大一些。 除了直接逼近法和逆算方法外，属于反分析的还有位移联图法，多元回归法等。以下就 几种常见的反分析计算法作一介绍。 二、 有限元反分析法 1．一般方法： 有限单元的基本方程可表示为 OK3{n)一 {F) (1) 其 中[K]为总刚矩阵(对于非线性问题，是应力或应变的函数)，{u)为结点位移矩阵，{F)为 29l 等效结点力矩 阵。 将 (1)式按测量的已知位移 u1和未知位移u_分开表示时，可写成； ft ]··CK·] ·’l= -’1 ) ． CK．] [ ] 【 )J 【{，)J (3) 由(2)、(3)可得 ： {F．卜CK．]CK ] {F，)一 (CK ]一 CK,]CK ] [ ]) ) (●) (4)式中．右端为已知量 ，左端为含有未知量地应力和力学参量的项。求解方程(4)，即可求出 未知量来 (对一非线性 问题则要进行迭代)。应当说明，这里方程的数 目应大于未知量的个 数 ，所 以还要进行优化计算。 2．弹性矩阵分解法 对(4)所示的方程，反算出地应力和弹性模量E，比较容易，但要同时反算出E和泊松 比 则比较困难，原因是 u不能直接从剐度矩阵中提出来，但用弹性矩阵分解法可 以解决此问 题 。 单元刚度矩阵可 以表示为 O O ．