随机平面弹性问题的协调有限元分析.pdfVIP

2014年 9月 高 等 学 校 计 算 数 学 学 报 第 36卷第 3期 随机平面弹性问题的协调有限元分析 许小静 范文文 谢小平 (四川大学数学学院，成都610064) CoNF0RM ING FIN ITE ELEM ENT ANALYSIS FOR STOCHASTIC PLANE ELASTICITY PR0BLEM S XuXiaojing FanWenwen XieXiaoping (SchoolofMathematics，SichuanUniversity，Chengdu610064) Abstract Thispaperconsidersthenumericalsolutionofplaneelasticity equa- tionswithstochasticYoung’Smodulusandstochasticloads．Thestochasticfields areapproximatedbyKarhunen--LofiveexpansionandW ienerpolynomialchaosex-． pansion alongsamplepaths，and in spacecontinuousLagrangian finiteelements areused．Errorestimatesarederived．Numericalexperimentsaredonetoverify thetheoreticalresults Keywords stochasticelasticityequation ，Karhunen—Lo~veexpansion，W iener— Chaosexpansion ，Lagrangianfiniteelement AMS(2000)subjectclassifications 65N30 中图法分类号 0242．21 国家自然科学基金 ，四川大学优秀青年学者基金 (2011SCU04B28) 收稿 日期：2012—04—06． 208 · 许小静等：随机平面弹性问题的协调有限元分析 第 3期 1 引 言 微分方程的随机性通常源于材料性能参数、物体所受力以及区域，这些随机项一般不 能显式表出，其处理方法可分为统计方法和非统计方法．统计方法，即MonteCarlo模拟 法 [9】，通过产生随机项的样本函数，将随机方程确定化．非统计方法需要对随机项的信息 有足够的了解，如其期望值、二阶矩等，运用 Taylor级数展开、KarhunenLo~ve展开 5『1． Stochastic等方法将随机项显式表出，从而将随机方程化成确定方程． 随机有限元结合随机分析方法和有限元方法，对随机方程的空间部分或整个确定化 的随机方程运用有限元法求解，如摄动方法、Neumaim展开方法 5『]以及 P×h型、k×h 型有限元法 [2]．对于二阶随机椭圆型微分方程，文献 4【]采用P×h型有限元法针对只有 参数或所受力是随机的情况进行了求解和误差分析；文献 2『]证明了参数和力都是随机时 解的存在唯一性和分别用蒙特卡洛方法、P×h型和 ×h型有限元方法求解时解的收敛 性；文献 [8]考虑了非线性随机椭圆方程的求解；文献 1【]针对随机区域的Neumann边值 问题进行求解；针对随机平面弹性问题，文献 6『1中用P×h型有限元法求解，但无方程解 的存在唯一性讨论及误差估计． 设D∈R。为有界区域，rd和 r 为OD 的两个开子集，满足OD=一FdUF— ，rdnr =0， 且meas(Fd)0．(2【，，P)是概率空间，其中Q为样本空间， 为事件域，P为概率．考虑 如下随机平面弹性问题 ：对几乎所有的0∈Q div~(u(-，))=f(·，)， inD (u(．，))：2pc(u(．，))+Adivu(