基于单位分解积分的伽辽金无网格方法研究

数值积分是伽辽金无网格方法实施的一个重要环节,提出了一种适合于伽辽金无网格方法的单位分解积分技术．该积分技术建立在有限覆盖和单位分解基础之上,不需要对积分区域进行分解,具有较高的积分精度．并以无单元伽辽金方法为例,详细说明了基于单位分解积分的伽辽金无网格方法的实现过程．这样,在近似函数建立和数值积分过程中都不需要进行网格划分,从而形成一种“真正的”无网格方法．     

Galerkin meshless methods based on partition of unity quadrature

Introduction Meshlessmethodsarenewmethodsofnumericalcomputationwhichhavebeendeveloped rapidlyinrecentyears.Inthesemethods,onlynodesareneeded,meshinformationistotally unnecessary.Thiscanavoidorpartlyavoidthedifficultyofmeshgeneration.Duetohigh accuracyandstability,Galerkinmeshlessmethodsareappliedbroadly,butitisunavoidable tocomputetheintegrationoverthewholephysicaldomaininGalerkinweakform,whichisa greatchallengeforGalerkinmeshlessmethodsbecauseoftheabsenceofmesh.TocarryouttheintegrationinGal…     

A(SLq(2)) at Roots of Unity is a Free Module over A(SL(2))

It is shown that when q is a primitive root of unity of order not equal to 2 mod 4, A(SL_q(2)) is a free module of finite rank over the coordinate ring of the classical group SL(2). An explicit set of generators is provided.     

Extended finite element modeling of deformable porous media with arbitrary interfaces

In this paper, an enriched finite element method is presented for numerical simulation of saturated porous media. The arbitrary discontinuities, such as material interfaces, are encountered via the extended finite element method (X-FEM) by enhancing the standard FEM displacements. The X-FEM technique is applied to the governing equations of porous media for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. In X-FEM, the material interfaces are represented independently of element boundaries and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of elements. Finally, several numerical examples are analyzed, including the dynamic analysis of the failure of lower San Fernando dam, to demonstrate the efficiency of the X-FEM technique in saturated porous soils.     

Partition of unity refinement for local approximation

In this article, we propose a Partition of Unity Refinement (PUR) method to improve the local approximations of elliptic boundary value problems in regions of interest. The PUR method only needs to refine the local meshes and hanging nodes generate no difficulty. The mesh qualities such as uniformity or quasi‐uniformity are kept. The advantages of the PUR include its effectiveness and relatively easy implementation. In this article, we present the basic ideas and implementation of the PUR method on triangular meshes. Numerical results for elliptic Dirichlet boundary value problem on an L‐shaped domain are shown to demonstrate the effectiveness of the proposed method. The extensions of the PUR method to multilevel and higher dimension are straightforward. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 803–817, 2011     

Modeling quasi-static crack growth with the extended finite element method Part I: Computer implementation

The extended finite element method (X-FEM) is a numerical method for modeling strong (displacement) as well as weak (strain) discontinuities within a standard finite element framework. In the X-FEM, special functions are added to the finite element approximation using the framework of partition of unity. For crack modeling in isotropic linear elasticity, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are used to account for the crack. This enables the domain to be modeled by finite elements without explicitly meshing the crack surfaces, and hence quasi-static crack propagation simulations can be carried out without remeshing. In this paper, we discuss some of the key issues in the X-FEM and describe its implementation within a general-purpose finite element code. The finite element program Dynaflow™ is considered in this study and the implementation for modeling 2-d cracks in isotropic and bimaterial media is described. In particular, the array-allocation for enriched degrees of freedom, use of geometric-based queries for carrying out nodal enrichment and mesh partitioning, and the assembly procedure for the discrete equations are presented. We place particular emphasis on the design of a computer code to enable the modeling of discontinuous phenomena within a finite element framework.     

一般n阶特征量泛函的Euler-Lagrange方程及与定量因果原理、相对性原理和广义牛顿三定律的统一

本文获得了有各种相互作用的一般n阶特征量泛函,其耦合系数反映了不同特征量泛函之间的耦合强度.依据定量因果原理,导出了一般n阶特征量泛函的变分原理,获得了一般n阶特征量泛函的Euler-Lagrange方程,它的不同系数可拟合不同的物理现实,如从线性到任意n阶非线性物理系统,使复杂难解的任意n阶非线性物理系统变得具体可解.并获得了该对称变换下不变的m个的守恒量,以及它们之间的关系和统一描述.依据定量因果原理导出了相对性原理,证明了绝对加速参考系、牵连参考系和相对参考系的力都有来自加速度和质量变化的贡献.利用定量因果原理自然导出了广义牛顿第一定律和广义牛顿第二定律,而且还导出了一个新定律,即广义牛顿第三定律,亦即平移不变性系统合力为零定理.进而将研究结论应用于对银河系的修正引力势、分子势、夸克禁闭势等,且其结果与物理实验一致.     

网络化靶场综合模拟试验网时间统一系统研究

为了在靶场网络化测试模拟试验网中统一多台测试设备的本地时间,以便组织网络化多台设备组网综合数据融合模拟试验,需要设计时统专用标定设备。在分析了几种常用授时方式的基础上,提出了采用单脉冲硬件驱动定标机制来统一标定各测试设备的本地时间。此硬件设备以微控制器为控制核心,通过统一启动时间基准标定信号,来对各测试设备的主机本地时间进行统一标定。在试验中,数据融合系统能够在新的统一时间基准下工作,融合后的数据具有实时驱动其他测试设备的能力。试验结果表明此时间基准设备能够达到系统测试功能要求,可用作在靶场网络化综合测试模拟试验中对多台测试设备组网的时间的统一标定。     

A greedy algorithm for partition of unity collocation method in pricing American options

A greedy algorithm in combination with radial basis functions partition of unity collocation (GRBF‐PUC) scheme is used as a locally meshless method for American option pricing. The radial basis function partition of unity method (RBF‐PUM) is a localization technique. Because of having interpolation matrices with large condition numbers, global approximants and some local ones suffer from instability. To overcome this, a greedy algorithm is added to RBF‐PUM. The greedy algorithm furnishes a subset of best nodes among the points X. Such nodes are then used as points of trial in a locally supported RBF approximant for each partition. Using of greedy selected points leads to decreasing the condition number of interpolation matrices and reducing the burdensome in pricing American options.     

Enriched meshless manifold method for two-dimensional crack modeling

The meshless manifold method is based on the partition of unity method and the finite cover approximation theory which provides a unified framework for solving problems dealing with both continuum with and without discontinuities. The meshless manifold method employs two cover systems. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the partition of unity functions. And the physical cover system describes geometry of the domain and the discontinuous surfaces in the domain. The shape functions are derived by the partition of unity and the finite covers approximation theory. In meshless manifold method, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the meshless manifold method is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the meshless manifold method, the enriched methods are introduced in this work for crack problems.