动态弹塑性扭转问题的区域分解方法.pdfVIP

第24卷第4期 苏 州 大 学 学 报 (自然科学版) Vo1．24 No．4 2008年 1O月 JOURNALOFSUZHOUUNIVERSITY(NATURALSCIENCEEDITION) 0ct．2o08 动态弹塑性扭转问题的区域分解方法 吴 娟 ，丁 睿h ，李荣军2 (1．苏州大学数学科学学院，江苏苏州215006；2．兰飞中学，甘肃兰州730000) 摘 要：研究了动态弹塑性扭转问题描述的发展型变分不等式的区域分解方法，给出了方法的实现步骤及收 敛性结果，实现了一维、二维数值算例，说明了方法的有效性． 关键词：动态弹塑性扭转问题；区域分解法；发展型变分不等式 中图分类号：0175．2 文献标识码：A 文章编号：1000—2073(2008)04—0024—40 0 介 绍 动态弹塑性扭转问题 ¨在物理、力学及工程应用中是一个很重要的问题．目前，求解此类问题的抛物变 分不等式的常用数值方法有，采用时间半离散化为椭圆变分不等式，再利用有限差分或有限元方法 ]、对偶 方法 、罚方法 ，等． 区域分解方法 (DDM) 是一种新兴的求解偏微分方程的数值方法，具有高效、并行的特点．目前，用于求 解变分不等式的区域分解方法主要局限于椭圆变分不等式 J，发展型变分不等式的区域分解方法还不多 见．在文[6]中，Zhou等研究了双侧障碍问题的抛物型变分不等式的区域分解方法．本文将区域分解方法推广 到由动态弹塑性扭转问题描述的发展型变分不等式中，给出了该方法的算法及收敛性结果．全文结构如下：第 一 部分给出了动态弹塑性扭转问题的Euler—Galerkin形式；第二部分陈述了动态弹塑性扭转问题的区域分解 算法及收敛性结果；第三部分通过一维、二维数值算例说明了方法的有效性． 讨论如下动态弹塑性扭转问题 的发展型变分不等式： 求u∈L(0，T；K)， ∈L(0，T；L(力))，使得 f(“，一M)+n(u，一 )≥ ，一M)，V ∈K，t∈(0，T]， ， ．、 u【(，0)=u()， 在 _j二． 、 其中 是R (d：1，2，3)中的有界区域,fEL(0，；L2())， K={∈ ()：l l≤1，在 上}， n。∈K，(．，．)是 上的内积，n(“，)：J u·Vvdx． JJ』 注：由发展型变分不等式解的存在唯一性定理 可知：问题(0．1)解的存在唯一性． 1 Euler—Galerkin方法 对 (O．1)中时间采用向后差分格式，记 t0=0t…tⅣ= ，r=t一t ，则问题为： 求u ∈K，使得 ff【。“_n+。un‘，一n“≥gn，一un，V∈， (1．1) o=0(戈)， =1，2，…， 收稿 日期：2007—10—24 基金项目：国家自然科学基金10672111)资助项 目 作者简介：吴 娟(1982一)，女，江苏南通人，硕士研究生，主要从事变分不等式数值解的研究 通讯联系人 第4期 吴 娟 ，丁 睿，李荣军：动态弹塑性扭转问题的区域分解方法 其中g= +f = ，t)． 设 为 上网格宽度为h的拟一致三角剖分， 为相应的有限元空间，： 为线性插值算子， = {∈ ：l l≤1，在 节点上)，则(1．1)的Euler—Galerkin方法如下 ： 求 ∈ ={∈Vh：l I≤1，在 节点上)，使得 (‘u ，一 ^)+0(Unh，一 )≥(g，一un^)，V ∈ ， ，1 ，)、 【11,o^=厶M0()