基于区间参数摄动法的三维弹塑性区间有限元探讨.pdfVIP

第25卷增刊 计算力学学报 V01．25，Sup． 2008年12月 Mech痂cs December2008 CIline∞JournalofCo,npumdond 基于区间参数摄动法的 三维弹塑性区间有限元探讨 黄耀英¨， 吴中如2， 苏静波2 (1．三峡大学土木水电学院，宜昌443002；2．河海大学水利水电工程学院，南京210024) 摘要：首次探讨了弹塑性问题的区问分析，推导了三维弹塑性问题的区间参数摄动法有限元计算公式，研制 了相应的弹塑性区间有限元程序，通过算例分析得到结论：(1)采用增量法对弹塑性问题进行区间分析时，增 量位移离差随着增量位移均值的增大而增大．(2)单元高斯点应力没有达到屈服极限之前，高斯点的应力离差 很小：当单元高斯点应力达到屈服极限之后，高斯点应力离差逐渐增大．(3)是否考虑塑性矩阵的应力分量为 弹塑性参数的隐函数，对应力高差计算的影响较小． 关键词：不确定性：区问分析模型；区间参数；摄动法；弹塑性问题 中图分类号：034 文献标识码：A 考虑到不确定参数P，在小范围变化，令 1 引 言 r=置(，)，置j(P)=等l，-，，F=盂(，)， 近年来区间分析方法逐渐为人们所熟知并成 为一种新的工程结构不确定分析方法，它主要用 妒)=参．，(i=1∥2”，m) (1) 来研究具有区间特性的工程结构。目前，在区间 式中：置‘、F分别为均值整体劲度矩阵和均值荷 有限元方法用于工程结构分析的求解过程中，人 载列阵，一为区间参数向量，Pf、∥分别为P，的 们提出了很多的区间有限元求解方’法【卜71，如截断 均值和离差，m为区间参数个数。 方法、摄动方法、组合方法以及优化方法等用于 将含有区间参数的劲度矩阵置(∥)和荷载列 工程结构分析，其中采用参数子区间与摄动方法 相结合可以得到较好的精度18J。实际工程中的问题 量，可得 严格来说一般为非线性问题，将某些非线性问题 (2) 简化为线性问题来处理，计算得到的结构响应有 x(P’)=∥+∑纰7·《(P)=r+膳。 时与实际不符。而在实际工程中经常提供的力学 (3) R(P’)=∥+∑纰’·碍(P)=F+缱7 参数如弹性模量、泊松比、凝聚力、内摩擦角等 常在一定的区间范围内，此时需要进行弹塑性问 由区间有限元控制方程，有 题的区间分析。但目前弹塑性区间有限元的文献 (置‘+△置’)(口‘+血。)=(F+△∥) (4) 还未见报导。据此，本文对弹塑性问题的区间分 由线性方程的摄动法，有 析进行初步探讨，推导弹塑性问题的区间参数摄 lIf=(置‘)-l∥ (5) 动法的计算公式，研制三维弹塑性问题的区间有 △且’=(置‘)一1△露。一(置‘)一1△置’口‘ 限元程序，并进行算例分析。