基于对偶混合变分原理的Signorini问题的数值模拟

计算物理 ›› 2002, Vol. 19 ›› Issue (2): 149-154.

基于对偶混合变分原理的Signorini问题的数值模拟

王光辉¹, 王烈衡²

1. 清华大学计算机科学与技术系, 北京 100084;
2. 中科院计算数学与科学工程计算研究所, 科学与工程计算国家重点试验室, 北京 100080

收稿日期:2000-07-24 修回日期:2001-03-05 出版日期:2002-03-25 发布日期:2002-03-25
作者简介:王光辉(1962-),男,湖楠益阳,博士,从事应用数值分析研究.
基金资助:
国家自然科学基金(19672064)资助项目

NUMERICAL MODELING OF SIGNORINI PROBLEM BASED ON DUAL MIXED VARIATIONAL PRINCIPLE

WANG Guang-hui¹, WANG Lie-heng²

1. Department of Computer Science and Technology, Tsinghua University, Beijing 100084, P R China;
2. Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, State Key Labora

Received:2000-07-24 Revised:2001-03-05 Online:2002-03-25 Published:2002-03-25

摘要/Abstract

摘要： 基于Signorini问题的对偶混合变分形式,提出了一种非协调有限元逼近格式,证明了离散的B-B条件,获得了Raviart-Thomas(k=0)有限元逼近的误差界O(h^3/4),并且Uzawa型算法对协调与非协调有限元逼近格式进行了数值求解.根据数值结果的分析和比较,表明应用非协调有限元逼近格式求解更有效.

关键词: Signorini问题, 对偶混合变分形式, Raviart-Thomas元, 非协调有限元, Uzawa算法

Abstract: Based on the dual mixed variational formulation for the Signorini problem,a nonconforming finite element method is proposed.The discrete B-B condition is confirmed and the error estimation O(h^3/4)for Raviart-Thomas(k=0)finite element is achieved.A Uzawa type algorithm is used for solving the Signorini problem which is discretized by conforming finite element method as well as nonconforming finite element method.The accuracy and efficiency of both are demonstrated by numerical results.By contrast to the conforming one,nonconforming method is more cost-effective.

Key words: Signorini problem, dual mixed variational formulation, Raviart-Thomas element, nonconforming finite element, Uzawa algorithm

中图分类号:

O753

王光辉, 王烈衡. 基于对偶混合变分原理的Signorini问题的数值模拟[J]. 计算物理, 2002, 19(2): 149-154.

WANG Guang-hui, WANG Lie-heng. NUMERICAL MODELING OF SIGNORINI PROBLEM BASED ON DUAL MIXED VARIATIONAL PRINCIPLE[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2002, 19(2): 149-154.

[1]	邢红玉, 叶文江, 周璇, 孙玉宝, 张志东, 宣丽. 挠曲电效应对向列相液晶盒响应时间的影响[J]. 计算物理, 2013, 30(5): 727-732.
[2]	张艳君, 王超, 叶文江, 张志东. Gruhn-Hess两体势修正模型向列相液晶微滴的Monte Carlo模拟[J]. 计算物理, 2013, 30(2): 244-250.
[3]	周璇, 张志东. 基于修正的Gruhn-Hess两体势模型研究弯曲形变向列相液晶盒[J]. 计算物理, 2010, 27(4): 598-602.
[4]	杨培培, 刘红(小), 刘红. 双轴向列相液晶的Monte Carlo模拟[J]. 计算物理, 2009, 26(4): 597-602.
[5]	叶文江, 邢红玉, 张志东, 陈国鹰. 挠曲电效应对向列相液晶盒电光效应的影响[J]. 计算物理, 2009, 26(1): 114-120.
[6]	叶文江, 邢红玉, 杨国琛. 栅状表面液晶盒的光学属性[J]. 计算物理, 2008, 25(5): 591-596.
[7]	张艳君, 张志东, 朱礼智. Monte Carlo模拟空间各向异性势向列相液晶微滴[J]. 计算物理, 2008, 25(3): 315-322.
[8]	叶文江, 邢红玉, 杨国琛. 向列相液晶盒中的挠曲电体积效应[J]. 计算物理, 2007, 24(3): 337-341.
[9]	张志东, 姜丽, 魏怀鹏. 反扭曲向列相液晶盒中的引流[J]. 计算物理, 2006, 23(2): 199-203.
[10]	张志东, 叶文江, 邢红玉. 液晶盒中挠曲电效应的计算[J]. 计算物理, 2004, 21(2): 156-160.
[11]	向兰, 张志东. 向列相液晶沿面到垂面排列转变的分子场理论[J]. 计算物理, 2003, 20(3): 210-214.
[12]	张志东, 卢遵铭. 液晶胆甾相中的分子短程关联[J]. 计算物理, 2003, 20(3): 215-218.
[13]	张志东, 孔浩, 薛山岭. 摩擦基板间液晶薄层的分子场理论[J]. 计算物理, 2001, 18(5): 443-446.
[14]	孙玉宝, 张志东, 马红梅. 手性沿面排列IPS模式液晶盒中的动态响应[J]. 计算物理, 2001, 18(5): 447-450.
[15]	张志东, 刘金伟, 马红梅. 弱锚泊条件下超扭曲向列相液晶显示的光学阈值电压[J]. 计算物理, 2000, 17(S1): 39-44.