平面弹性问题区间解及其在土木工程应用.pdfVIP

第十二届全国现代结构工程学术研讨会 平面弹性问题区间解及其在土木工程中的应用 王石瑛1聂慧1王日斌2 (1同济大学土来工程学院。上海200092-2．燕山大学里仁学院机械工程系，秦皇岛006004) 摘要；土木工程中存在大量不确定性的物理量，如岩石、土壤的相关物理参数，混凝土的弹性模量等。区间分析是可以处 理此类不确定性物理量的一种数学方法。使用区间分析…解决土木工程中不确定性问题时，很容易产生区间相关性 于重力坝的区间分析，并通过与较大样本数的MonteCarlo抽样法结果比较，初步表明了本文程序的正确性、有效性 和可行性。 关键词：不确定问题，区间分析，线弹性平面应变问题，EBE方法 一、引言 目前土木工程的结构分析和设计主要基于确定的结构参数及数学模型，然而实际工程中存在大量的不 确定因素，这些不确定因素会导致结构响应的不确定性。实际工程应用中的不确定性主要表现为：实际问题 被简化成理想模型时存在的模型不确定性：有限元进行离散时由单元数量、形函数选取造成的离散误差； 参数不确定性‘21即几何尺寸、材料的特性、边界条件的不确定性；舍入或截断误差。 二、基于区间方法的工程不确定性问题研究现状 不确定性问题的解决方法主要基于以下三种理论，即随机有限元理论、模糊集模型、区间模型。当结构 的不确定性缺乏统计特性，或缺乏足够的数据信息描述不确定参量的概率特性pJ，而只知道参数的界限时， 采用区间模型比较好。自90年代中期，很多学者便开始将区间理论与有限元方法相结合，建立起区间有限 元方法(IFEM)用以解决工程结构的不确定性问题。 引入区间量后，原来有限元中要求解的线性方程组在区间有限元中变成区间方程组，这样极大地增加 了求解的难度。现有区间有限元解法主要有如下几种：1)闭合解法，该方法适用于极简单问题。2)组合 解法，该方法对于简单问题结果较好，对于复杂问题误差较大且计算量庞大[41。3)直接优化法，当自由度 较多且有相关性存在时该法很难实施。4)摄动法，在计算时忽略方程中的高阶项，因此该法并不能保证真 实解在得到的解区间内。5)敏感性分析法，此法假定函数是单调的，利用敏感分析的概念来确定参数的上 足够大时，得到的最大及最小的区间会逼近真实的响应区间，但当样本数很庞大时会相当耗时。上述6种 方法用于求解区间有限元方程存在很多限制，在具体的实际问题中采用迭代法比较理想，但当存在相关性 问题时迭代法收敛比较困难。若采用EBE方法能很好的解决相关性的问题。本文基于论文【5】的EBE理论 研究，用基于MATLAB及Intlabt7】程序实现了该算法，通过算例验证该算法的可行性。 三、运用EBE技术的迭代法求解线性区间方程 工业建筑2012增刊 1273 第十_届个隔现代结构工程学术研讨会 (一)区间相关性问题 运用区间运算法则时，当同一区间变量出现的次数超过1次时会出现相关性问题，使得按照区间算法 函数值是0。这就是由于在区间运算时将函数的第一项和第二项X看成两个独立的变量，导致区间的扩张。 若要减少区间扩张，增加结果精确度就要保证减少区间量出现的次数，同时在实数和区间数的混合计算中 将区间量计算放到最后。 (二)Element Element技术 By 通用有限元在总刚进行组装时，因为总刚中的某些元素可能会来自予同一个单元刚度矩阵，且总刚中 的某个元素也可能由几个单元的贡献，因此总刚中的某些元素是相互关联的【6】。即刚度矩阵存在相关性问题。 计算时可以标记出这些相关的量，但是若对每个量都进行追踪则会相当麻烦，计算量也相当庞大。针对上 述问题张浩提出了ElementBy 相邻单元相同的节点都要分开编号，保证每个单元都是一个独立体。这样做会让节点数增加很多，但对于 已知的边界条件的节点就可以不集成该自由度的总刚元素，从而降低总刚的计算规模。将单元离散必定会 导致总刚矩阵的奇异性，相邻单元公共边编号不同的节点位移不协调，解决的办法是采用罚函数法【51在公共 边节点上施加一个约束使位移相同。 (三)基于EBE的固定点迭代法求解 采用EBE方法【5】组装的线弹性结构的区间方程的一般形式可表示为