非局部理论克服网格敏感性谱分析研究.pdfVIP

中国力学学会学术大会2009 （CCTAM2009 ） 非局部理论克服网格敏感性谱分析研究1) *，** 2) *，** 吕玺琳 黄茂松 * （同济大学地下建筑与工程系，上海 200092 ） ** （同济大学岩土与地下工程教育部重点实验室，上海 200092 ） 摘要：由于缺少一个特征长度，局部塑性软化模型使得控制偏微分方程呈现病态，导致有限元数值解具有网格敏感 性。通过对离散增量平衡方程刚度矩阵谱分析来研究导致网格敏感性的本质，结果表明，控制系统变形模式的主特 征向量极不规则且依赖于单元大小，局部化变形区域大小依赖于所引入的弱单元尺寸。通过引入一个特征长度，非 局部理论使变形模式变得光滑，局部化变形区域大小保持固定，从而使得塑性应变分布和荷载响应计算结果保持网 格客观性并随着网格细化趋近于理论解。将局部和非局部塑性软化模型用于双轴应变局部化试验数值模拟中，谱分 析表明，局部模型随着网格细化将有负特征值出现，解将出现分叉，边值问题出现病态。非局部化理论使得所有特 征值保持为正，得到了一致收敛解，这表明病态边值问题得到了正则化。 关键词：非局部理论，应变局部化，网格敏感性，谱分析 引 言 [1] 当用有限元模拟应变局部化现象时，常遇到计算结果存在网格敏感性的问题 。问题的本质在 于计算过程中部分区域内软化模型的采用使得控制偏微分方程变型，导致边值问题失去适定性而呈 现病态。要想克服网格敏感性，需引入某种正则化机制，当前普遍采用是通过在本构方程中引入特 征长度进行正则化，如粘塑性模型，微极(Cosserat)模型，梯度模型及积分非局部模型[2]等都属于这 一范畴。相对而言，非局部理论由于物理意义更明确，数值实现较容易，因而应用较广。非局部理 论虽在应变局部化问题中得到了一定应用，但其正则化机理目前尚不清楚，特别是如何克服数值解 的网格敏感性方面还需要深入探讨。本文将动力分析中广泛应用的谱分析方法[3]用于静力应变局部 化问题的研究中，通过分析变形模式来形象地揭示非局部理论克服网格敏感性的本质。 1 非局部塑性本构模型 在一维情形，局部塑性软化模型描述为： ⎧ p σ ε ε E ( − ) ⎪ ⎨ (1) σ ε σ σ σ σ F ε p H p =−( , ⎪⎩ ) − + ( y ) 0 p σ ε 式中： 为应力； 为应变；E 为弹性模量；ε 为塑性应变率；F 为屈服函数；σ 为当前屈 y 服应力；σ 为初始屈服应力；H 为塑性模量，当线性软化时，为一负的常数。 0 塑性应变演化规律为： ∂F p