具有小周期孔隙复合材料弹性结构的双尺度有限元分析

对于具有小周期孔隙复合材料弹性结构,在双尺度渐近分析理论结果的基础上提出了双尺度有限元计算格式,并给出了严格的误差估计.     

整周期复合材料弹性结构的有限元计算

1．引言周期性复合材料与周期结构的弹性力学问题，由于其材料特征剧烈振荡，且周期很小，问题相当复杂，除了一些特殊的和简单的问题可以用解析法求解外，多数问题很难或不可能用解析法求解，需要采用数值方法计算，有限元法无疑是最有效的方法之一．用细观力学方法研究复合材料的力学问题时，需要涉及纤维的排列情况，纤维和基体的性能，界面的分布情况，以及细观的几何参数和物理参数等．由于复合材料细观构造的不均匀性和不规则性，损伤和缺陷的存在，以及许多难以精确测定的因素，使得复合材料的细观力学问题十分复杂，不作出一些简化…     

COUPLING OF FINITE ELEMENT AND BOUNDARY ELEMENT METHODS FOR THE SCATTERING BY PERIODIC CHIRAL STRUCTURES

Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In this paper, variational formulations coupling finite element methods in the chiral medium with a method of integral equations on the periodic interfaces are studied. The well-posedness of the continuous and discretized problems is established. Uniform convergence for the coupling variational approximations of the model problem is obtained.     

复合材料弹性结构的高精度多尺度算法与数值模拟

１．引言 由于复合材料结构物理参数的非均匀各向异性以及细部几何构形的复杂性,在计算它的位移场、应力、应变场时,传统的有限元法因网格生成困难和计算规模太大而难以实现．７０年代初,Ｉ．Ｂａｂｕｓｋａ,Ｊ．Ｌ．Ｌｉｏｎｓ等人针对复合材料弹性结构提出了均匀化方法,见文［１］,数值实验表明,均匀化方法对描述复合材料弹性结构的有效材料常数及刚度性质是有效的,但它不能刻画应力和应变场的局部变化,因而作为复合材料强度理论的判断依据,均匀化理论显然是不够的．为此,Ｊ．Ｌ．Ｌｉｏｎｓ,Ｏ．Ａ．Ｏｌｅｉｎｉｋ等分别就…     

复合材料拟周期结构的均匀化方法

1．预备知识和问题的提出文山讨论了复合材料周期结构的数学问题,即复合材料弹性体的物理参数,几何参数具有某种周期拓扑结构的问题．但在实际问题中,有时还会遇到这些情形,由于受湿热环境和结构设计等因素的影响,弹性体不具有整体周期性,而是具有局部周期性,这就是本文要研究的拟周期性．数学上,这些问题可以表述为下列边值问题：其中fiCR”是有界区域,具有LIPschitz连续边界,A‘巾,旬,＜二。‘x,八,k＝l,…,n．是nXn阶矩阵,满足下列的A一条件．J．L．Lions（见文门）等曾对拟周期的二阶椭圆型边值问题进行了讨论,…     

MULTIGRID ALGORITHM FOR THE COUPLING SYSTEM OF NATURAL BOUNDARY ELEMENT METHOD AND FINITE ELEMENT METHOD FOR UNBOUNDED DOMAIN PROBLEMS

In this paper, some V-cycle multigrid algorithms are presented for the coupling system arising from the discretization of the Dirichlet exterior problem by coupling the natural boundary element method and finite element method. The convergence of these multigrid algorithms is obtained even with only one smoothing on all levels. The rate of convergence is found uniformly bounded independent of the number of levels and the mesh sizes of all levels, which indicates that these multigrid algorithms are optimal. Some numerical results are also reported.     

基于双尺度渐近分析的有限元算法

1．引言正如文山所说，由于复合材料和周期结构的材料系数ail（x）在局部区域内间断且跳跃性很大，加上区域内含有周期性洞穴或裂缝，且周期长度很小．一般而言，直接采用有限元方法进行数值模拟，其计算量大得惊人，甚至难以实现．文山针对这种特征，提出了一种可计算的双尺度渐近分析模式，本文在此基础上给出了相应的有限元算法，它包括：1．周期解在一个基本构造上的有限元计算；2．边界层的有限元计算．同时，给出了相应的误差分析．2．周期解的有限元计算首先考虑下列形式的边值问题；其中把，代E尸（on叫，iii（0关于E—（EI，ZZ…     

Webster方程有限元法的超收敛性

<正>1前言由语音生成过程建立的数学模型是一个无界区域上的声压波动方程,引入人工边界可将这个无界区域分成一个内部有界区域和另一个无界区域,再将时间和空间进行变量分离,就可以得到一个如下一维简单的基于内部有界区间的微分方程边值问题     

求解美式回望期权的有限元方法

期权作为一种金融衍生产品,在欧美国家一直很受欢迎.由于其规避风险的特性,期权也吸引了中国投资者的兴趣.基于市场的需求,2015年初,上海证券交易所推出了中国首批期权产品,期权定价问题的研究热潮正席卷全球.本文研究的美式回望期权,是一种路径相关的期权,其支付函数不仅依赖于标的资产的现值,也依赖其历史最值.分析回望期权的特点,不难发现:1)这类期权空间变量的变化范围为二维无界不规则区域,难以应用数值方法直接求解;2)最佳实施边界未知,使得该问题变得高度非线性.本文的主要工作就是解决这两个困难,得到回望期权和最佳实施边界的数值逼近结果.现有的处理问题1)的有效方法是采用标准变量替换、计价单位变换以及Landau变换将定价模型化为一个[0,1]区间上的非线性抛物问题,本文也将沿用这些技巧处理问题1).进一步,采用有限元方法离散简化后的定价模型,并论证了数值解的非负性,提出了利用Newton法求解离散化的非线性系统.最后,通过数值模拟,验证了本文所提算法的高效性和准确性.     

GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS

1. IntroductionIn the numerical simulation of the Navier-Stokes equations one encounters three seriousdifficulties in the case of large Reynolds numbers f the treatment of the incomPressibility con-dition divu = 0, the treatment of the noIilinear terms and the large time integration. For thetreatment of the incoInPressibility condition, one use the penalty method in the case of finiteelemellts [1--2l and for the treatmen of the noulinar terms and the large tfor integration, oneuse the nonlin…     

一类波动方程有限元投影格式的误差分析

1.引 言 本文考虑如下不含阻尼项的波动方程的有限元逼近: 其中区域Ω Rd（d=2,3）是足够光滑的有界多边形区域,其边界为Γ=　Ω.初始条件为:当t=0时,u=u0,ut=u1.     

SUPERCONVERGENCE ANALYSIS FOR CUBIC TRIANGULAR ELEMENT OF THE FINITE ELEMENT SUPERCONVERGENCE ANALYSIS FOR CUBIC TRIANGULAR ELEMENT OF THE FINITE ELEMENT

1. IntroductionAlthough we bed proved the superconvergence of quadratic triangular elem6ntsbefore 1985, the superconvergence reseaxch of k(k 2 3)-degree triangulax elemeds Onlyhas a few advances, e.g.,.Lin, Yan and Zhou (see [15]) Prove that the three degreeHerlnite elements possess superconvergence and Walilbin (see [5-7]) obtains a roughresult by using a fine interinr estimation, that is, the placement fUnCtion or its gradestmay have weak superconvergence in the lOcal syUUnetric points,e…     

GLOBAL SUPERCONVERGENCE ESTIMATES OF FINITE ELEMENT METHOD FOR SCHRODINGER EQUATION

1.IntroductionConsiderthefollowinginitialboundaryvalueproblemofSchr6dingerequationwheren~[0,1]',at~On/Ot,T>0isaconstant.Theequivalentvariationalformof(1.1)is:foralltE[0,T],findu(t)6Hi(n)satisfiesthefollowingvariationalequation:where(w,v)~IwvdxdenotestheinnerproductofL'(fl)anda(u,v)~(Vu,Vv),ibetheimaginaryunit.Weassumethatthefunctionsarecomplex--valuedandHibertspacesarecomplexspaces.LetThbeaquasiuniformrectangulationoffiwithmeshsizeh>0andS'(O)CHi(fi)bethecorrespondingpiecewisebilinearpol…     

DISCONTINUOUS FINITE ELEMENT METHOD FOR CONVECTION-DIFFUSION EQUATIONS

1. IlltroductionThe finite element approximation of the convection--diffusin equations has been investigated using several different approaches (see e.g. [3] [4] and the references therein).Previous analysis in primal formulation of these problems was done for two types ofapproximation schemes: one which produces a continuous piecewise polynomial approximation and one which produces a piecewise polynomial approximation which arecontinuous for certain number of moments accross interelement edge…     

非定常的热传导-对流问题的混合有限元-Galerkin交替方向有限元法

1 引 言设Ω　R2为具有光滑边界的有界区域,考虑非定常的,无量纲化的,而且带有热传导的粘性不可压缩流体力学问题:问题（Ⅰ）:求u=（u1,u2）,p,T满足:     

CONVERGENCE OF AN EXPLICIT UPWIND FINITE ELEMENT METHOD TO MULTI-DIMENSIONAL CONSERVATION LAWS

1. IntroductionThe convergence problem for the numerical schemes to one dimensional cons~inn lawshas been edensively studied. By tensor product one dimensional schemes can be appliedto multi-dimensional equations. However the convergence of many of those schemes is stillchanknown even if it is true for one dimensional cases. Besides, for those Phys.ical domains withcomplicated geometry unstructured grids are more fiekible. In recent yeaes the convergellceproblem for unstructured grids has cal…     

A FINITE VOLUME ELEMENT METHOD FOR THERMAL CONVECTION PROBLEMS

Consider the finite volume element method for the thermal convection problem with the infinite Prandtl number. The author uses a conforming piecewise linear function on a fine triangulation for velocity and temperature, and a piecewise constant function on a coarse triangulation for pressure. For general triangulation the optimal order H^1 norm error estimates are given.     

有限元的一个局部超收敛结果

１．引 言 设 是一个有界开域，具充分光滑的边界 且设 是 上的一族拟一致的三角剖分，用 表示定义在Ｔｈ上的分片线性有限元空间，并置考虑模型问题 用 分别表示的有限元解及内插，那么有插值估计：（见［１］）一般地，如ｕ为问题（１．１）的解，我们有有限元逼近误差估计（见［３］） 命题１．设 并设 分别表示按定义的Ｇｒｅｅｎ函数及其有限元逼近，那么有其中 Ｃ与 ｚ，ｈ无关．（参见［３］） 注意．如 且 ，那么至少存在一个点 ，使即ｘ０是ｆ的奇点，例如其中 为常数， ，显然如果。，如果故我们假定 本文将证明，误差与ｆ的奇性…     

A FINITE ELEMENT METHOD WITH RECTANGULAR PERFECTLY MATCHED LAYERS FOR THE SCATTERING FROM CAVITIES

We develop a finite element method with rectangular perfectly matched layers （PMLs） for the wave scattering from two-dimensional cavities. The unbounded computational domain is truncated to a bounded one by using of a rectangular perfectly matched layer at the open aperture. The PML parameters such as the thickness of the layer and the fictitious medium property are determined through sharp a posteriori error estimates. Numerical experiments are carried out to illustrate the competitive behavior of the proposed method.