复杂表面问题的有限元计算与分析

摘要/Abstract

摘要： 碰撞、弹塑变形问题等交通运输器的安全试验通常采用实验室实物模型与计算机模型相结合进行模拟.对三维复杂区域的接触，国际标准方法采用接触变形算法模拟滑动控制.本课题采用流固耦合非牛顿流体方程初边值问题求解三维层结构特性，通过基于变分原理的摄动问题有限元方法，在高性能软件平台上实现数据的挖掘处理.由Sobolev空间嵌入原理，可将模型按接触区域进行分层单元剖分，将复杂区域剖分为若干相互连接、不重叠的六面体与空间平面四边形单元.同时，建立微观与宏观有限元双尺度计算模型进行模拟仿真对比，得到模型的能量与速度等一系列参数的变化曲线.此外，接触表面问题又可采用渐近摄动方法中的边界层理论进行研究，由此得到的微分方程特征函数空间，既可作为优化有限元基函数的解，又可用于建立一种新型的非线性特征值的渐近方法，也是估计材料特定参数的方法之一.最后，使用人工边界条件随机处理方法对求解结果的数据进行分析.

关键词: 复杂接触表面, 移动边界条件, 多层结构耦合

Abstract: The safety test of transportation devices, including impact test and plastoelastic deformation, is simulated both in laboratory and by computer. Contact deformation algorithm is an international standard method to simulate the skid control on complex 3-D areas. In this paper, coupled non-Newtonian fluid equations with initial boundary value are used to solve the 3-D layer structure. A finite element method (FEM) based on the variation principle is used to solve the perturbation problem, and data mining is processed by high performance software. According to the embedding principle, a stratified element division is processed, and the complex boundary is divided into several mutual connected yet not overlapped hexahedral and square elements. After building a FEM model of both macro scale and micro scale, varying curves of parameters including energy and velocity can be obtained. On the other hand, the boundary layer theory of asymptotic perturbation method is also a way to study the complex boundary problem. The characteristic function space gained can be used both to optimize the primary function of FEM, and to establish a new asymptotic method to solve the nonlinear eigen problem. It can be used to estimate the specific parameters of materials as well. The stochastic analysis of artificial boundary conditions is then used to study the resulting data.

Key words: complex contact boundary, moving boundary conditions, coupled multi-layer structure

中图分类号:

O242.1

侯磊，李涵灵，林德志，张敏， Daglish GEORGE. 复杂表面问题的有限元计算与分析[J]. 华东师范大学学报(自然科学版), 2012, 2012(4): 1-11,26.

HOU Lei, LI Han-ling, LIN De-zhi, ZHANG Min, Daglish GEORGE. Finite element computation and analysis on complex contact boundary[J]. Journal of East China Normal University(Natural Sc, 2012, 2012(4): 1-11,26.

参考文献

［1］HOU L. Failure modes analysis in the crash barrier simulation［C］// Proc 5th Euro LS-Dyna Conference. Birmingham: The ICC Birmingham, 2005. B-02.

［2］FRANK P. The work of Richard von Mises［J］. Science, 1954, 119(3102): 823-824.

［3］GIBSON L J, ASHBY M F. Cellular Solids, Structure and Properties［M］. Pergamon: Oxford, 1988.

［4］ZHOU Q, WIERZBICKI T. An incremental analysis of plane strain fully plastic crack growth in strain- hardening materials under extension［J］. International Journal of Fracture, 1996, 79: 27-48.

［5］CAO L Q, CUI J Z, ZHU D C. Multiscale asymptotic analysis and numerical simulation for the second order Helmholtz equations with rapidly oscillating coefficients over general convex domains［J］. SIAM J Numerical Analysis, 2002, 40(3): 543-577. 

［6］ZHANG Hui, ZHANG Pingwen. Local existence for the FENE-dumbbell model of polymeric fluids［J］. Archive for Rational Mechanics and Analysis, 2006, 181(2): 373-400.

［7］LIN Q, YAN N. Global superconvergence of mixed finite element methods for Maxwell’s equations［J］. Journal of Engineering Mathematics, 1996, 13: 1-10. 

［8］林群,谢和虎,罗福生,等. Stokes方程非协调混合元的特征值下界［J］. 数学的实践与认识, 2010, 40(19):157-168.

［9］PENG Shi Ge. Backward stochastic differential equations and application in optimal control［J］. Appl Math & Optima, 1993, 27: 125-144.

［10］ LIN Q, LU T. Asymptotic Expansions for Finite Element Eigenvalues and Finite Element Solution［M］. ［S.L.］: Bonn Math Schrift, 1984. 

［11］林群,周俊明,陈竤焘.椭圆形方程四面体线元的超逼近与外推［J］. 数学实践与认识, 2009, 39(15): 200-208.

［12］李开泰,黄爱香,黄庆怀.有限元方法及其应用［M］. 北京: 科学出版社, 2005.

［13］陈传淼.科学计算概论［M］. 北京: 科学出版社, 2007.

［14］ LIN Q, LIN J F. Finite Element Methods: Accuracy and Improvement［M］. Beijing: Science Press, 2006.

［15］ LIN Q, YAN N N. Construction and Analysis for Effective Finite Element Methods［M］. Hebei: Hebei University Press, 1996.

［16］ THOMAS J W. Numerical Partial Differential Equations［M］. New York: Springer-Verlag, 1995.

［17］ KELLER J B, GIVOLI D. Exact non-reflecting boundary conditions［J］. Comput Phys, 1989, 82: 172-192.

［18］ HAN H D. A new class of variational formulations for the coupling of finite and boundary element methods［J］. Comput Math, 1990, 8(3): 223-232.

［19］ GATICA G N, HSIAO G C. On the coupled BEM and FEM for a nonlinear exterior Dirichlet problem［J］. Numer Math, 1992, 61(2): 171-214.

［20］ XU J. Two-grid discretization techniques for linear and nonlinear PDEs［J］. SIAM J Numer Anal, 1996, 33(5): 1759-1777.〖HJ〗

［21］ MEDDAHI S, GONZALEZ M, PEREZ P. On a FEM-BEM formation for an exterior quasilinear problem in the plane［J］. SIAM J Numer Anal, 2000, 37(6): 1820-1837.

［22］ YU D H. Natural boundary integral method and its applications［M］. Beijing: Science Press/Kluwer Academic Publishers, 2002.

［23］ YING L A. Numerical Methods for Exterior Problems［M］. New Jersey: World Scientific, 2006.

［24］ HSIAO G C. Boundary element methods—an overview［J］. Appl Numer Math, 2006, 56(10-11): 1356-1369.

［25］ XU Z L, HAN H D, WU X N. Numerical method for the deterministic Kardar-Parisi-Zhang equation in unbounded domains［J］. Commum Comput Phys, 2006, 1(3): 481-495.

［26］ HOU L, QIU L. Computation and asymptotic analysis in the non-Newtonian impact problem［J］. ACTA Math Applic Sinica, 2009, 25(1): 117-126.

［27］ BURDEN R L, FAIRES J D. Numerical Analysis［M］. Beijing: Higher Education Press, 2001.

［28］ HOU L, PARIS R B, WOOD A D. Resistive interchange mode in the presence of equilibrium flow［J］. Physics of Plasmas, 1996, 3(2): 473-481.

［29］ HOU L, HARWOOD R. Nonlinear properties in the Newtonian and non-Newtonian equations［J］. Nonlinear analysis, 1997, 4(30): 2497-2505.

［30］ ADAMS R A.索伯列夫空间［M］.叶其孝译,北京: 人民教育出版社, 1981.

［31］吴其芬, 李桦. 磁耦合力学［M］. 北京: 国防科技大学出版社, 2007.

［32］陆金甫. 偏微分方程数值解法［M］2版.北京:清华大学出版社, 2003.

［33］王元明. 索伯列夫空间讲义［M］. 南京: 东南大学出版社, 2008.

［34］林群. 微分方程数值解法［M］. 北京: 科学出版社, 2003.

［35］张平文, 李铁军. 数值分析［M］. 北京: 北京大学出版社, 2007.

［36］黄世霖, 张金换, 王晓东, 等. 汽车碰撞与安全[M]. 北京: 清华大学出版社, 2000. ［37］MARCHAL J M, CROCHET M J. A new mixed finite element for calculating viscoelastic flow[J]. J NonNewt Fluid Mech, 1987, 26: 77-114. ［38］ HOU L, LI H L, ZHANG J, et al. Boundarylayer eigen solutions for multi-field coupled equations in the contact interface[J]. Appl Math Mech, 2010, 31(6):719-732. ［39］SUGENG F, PHANTHIEN N, TANNER R I. A study of non-isothermal non-newtonian extrudate swell by a mixed boundary element and finite element method[J]. J Rheology, 1987, 31(1): 37-58.

［40］FISHER R A. Theory of statistical estimation[C]// Proceedings of the Cambridge Phil Society, 1925, 22: 709-715. ［41］HOU L, LI H L, LIN D Z, et al. The stochastic boundary-layer in the non-Newtonian problem[C]// Proc WCE London, U K.: 2010, 3: 1401-1407. ［42］ HOU L, NASSEHI V. Evaluation of the stress effects flow in rubber mixing[J]. Nonlinear Analysis, 2001, 47(3): 1809-1820. ［43］HOU L, DING H Y. High performance computing and mathematical study for the viscoelastic impact problem[C]// Proc of HAPC. Shanghai: LNCSSpringer, 2010: 169-176.［44］ WEAIRE D L, HOU Lei. A statistic note on the mature 2D soap froth[J]. J Phil Mag Letters, 1990, 62: 427-430.