配点型广义节点无网格法的基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。     

非等温广义牛顿流体的无网格模拟

基于原始变量法,将RPIM推广到非等温广义牛顿流动问题的求解.为了减少未知量的个数,压力和速度采用罚函数方法耦合;同时采用积分降阶技术以保证获得满意的数值解.数值实验结果表明,采用罚函数法处理速度和压力在无网格方法中同样适用,且RPIM用于非等温广义牛顿流动问题的求解时具有易施加本质边界条件、计算精度高和收敛性较快的优点.     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

局部彼得洛夫-伽辽金法分析各向异性板屈曲

基于Kirchhoff板理论和对挠度函数采用移动最小二乘近似函数进行插值，进一步研究无网格局部Petrov-Galerkin(MLPG)方法在各向异性板稳定问题中的应用．分析中，本质边界条件采用罚因子法施加，离散的特征值方程由板稳定控制方程的局部积分对称弱形式中得到．通过数值算例并与其他方法的结果进行比较，表明MLPG法求解各向异性薄板稳定问题具有收敛性好、精度高等一系列优点．     

固定网格下的特征线法求解溃坝问题

溃坝问题是典型的非线性双曲方程的Riemann问题,其数值求解的难点在于对间断面的捕捉以及避免间断面处在数值计算过程中产生数值色散,因而为求解此问题所产生的各种数值计算方法的优劣也体现在这两个方面。本文针对溃坝问题提出一种新的计算方法。该方法基于对偶变量推导的浅水波方程,根据方程的特点,从方程的特征值和黎曼不变量出发,采用高精度的激波捕捉方法计算黎曼不变量的位置随时间的变化,然后映射至不随时间变化的固定网格。根据黎曼不变量的位置,采用保形分段三次Hermite插值将物理量映射至网格节点。计算结果显示,该方法不仅操作简单,计算量小,而且结果准确。     

基于相场描述的拉压不对称强度结构拓扑优化

运用了基于相场描述的拓扑优化方法，来寻找在拉伸和压缩中表现出不对称强度行为的连续体结构的最优布局。依据Drucker-Prager屈服准则和幂率插值方案，优化问题可以描述为在局部应力约束下的最小化结构的体积。用qp放松法来解决应力约束的奇异性，并采用基于P-norm函数的聚合方法对应力约束进行凝聚，该方法实现了约束个数的降低，同时引入了稳定转化法来处理大量的局部应力约束和高度非线性的应力行为，以修正应力，提高优化收敛的稳定性。在优化问题求解时，使用拉格朗日乘子法对目标函数和应力约束进行处理。利用伴随变量法进行灵敏度分析，并通过求解Allen-Cahn方程更新相场函数设计变量。数值算例证明了该优化模型和相应数值技术的有效性，相关算例还揭示了考虑拉压不同强度和考虑同拉压强度约束时得到的结构优化拓扑构型具有显著的差异。     

固定网格下的特征线法求解溃坝问题

溃坝问题是典型的非线性双曲方程的Riemann问题，其数值求解的难点在于对间断面的捕捉以及避免间断面处在数值计算过程中产生数值色散，因而为求解此问题所产生的各种数值计算方法的优劣也体现在这两个方面。本文针对溃坝问题提出一种新的计算方法。该方法基于对偶变量推导的浅水波方程，根据方程的特点，从方程的特征值和黎曼不变量出发，采用高精度的激波捕捉方法计算黎曼不变量的位置随时间的变化，然后映射至不随时间变化的固定网格。根据黎曼不变量的位置，采用保形分段三次Hermite插值将物理量映射至网格节点。计算结果显示，该方法不仅操作简单，计算量小，而且结果准确。     

非稳态Hamilton-Jacobi方程的7阶加权紧致非线性格式

Hamilton-Jacobi (HJ) 方程是一类重要的非线性偏微分方程, 在物理学、流体力学、图像处理、微分几何、金融数学、最优化控制理论等方面有着广泛的应用. 由于HJ方程的弱解存在但不唯一, 且解的导数可能出现间断, 导致其数值求解具有一定的难度. 本文提出了非稳态HJ方程的7阶精度加权紧致非线性格式 (WCNS). 该格式结合了Hamilton函数的Lax-Friedrichs型通量分裂方法和一阶空间导数左、右极限值的高阶精度混合节点和半节点型中心差分格式. 基于7点全局模板和4个4点子模板推导了半节点函数值的高阶线性逼近和4个低阶线性逼近, 以及全局模板和子模板的光滑度量指标. 为避免间断附近数值解产生非物理振荡以及提高格式稳定性, 采用WENO型非线性插值方法计算半节点函数值. 时间离散采用3阶TVD型Runge-Kutta方法. 通过理论分析验证了WCNS格式对于光滑解具有最佳的7阶精度. 为方便比较, 经典的7阶WENO格式也被推广用于求解HJ方程. 数值结果表明, 本文提出的WCNS格式能够很好地模拟HJ方程的精确解, 且在光滑区域能够达到7阶精度; 与经典的同阶WENO格式相比, WCNS格式在精度、收敛性和分辨率方面更优, 计算效率略高.      

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

考虑耦合质量影响的饱和多孔介质动力响应分析的显式有限元法

根据Biot饱和多孔介质动力方程,采用解耦技术,提出了考虑耦合质量Pd影响的饱和多孔介质中动力响应分析的显式有限元法。文中建立并推导了显式有限元的公式,编制了相应的计算程序并进行了实例计算。计算结果与解析解进行了对比,两者符合很好,表明本文方法是处理饱和多孔介质动力问题的一种有效方法。文中还分析了耦合质量ρa对固相和液相动位移的影响。     

Velocity–pressure integrated versus penalty finite element methods for high-Reynolds-number flows

Velocity–pressure integrated and consistent penalty finite element computations of high-Reynolds-number laminar flows are presented. In both methods the pressure has been interpolated using linear shape functions for a triangular element which is contained inside the biquadratic flow element. It has been shown previously that the pressure interpolation method, when used in conjunction with the velocity-pressure integrated method, yields accurate computational results for high-Reynolds-number flows. It is shown in this paper that use of the same pressure interpolation method in the consistent penalty finite element method yields computational results which are comparable to those of the velocity–pressure integrated method for both the velocity and the pressure fields. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine grid finite difference computational results. Advantages and disadvantages of the two finite element methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow of Reynolds number 10 000, a laminar backward-facing step flow and a laminar flow through a nest of cylinders.     

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

LIQUID-SOLID COUPLED SYSTEM OF MICROPUMP

This paper employs the integral-averaged method of thickness to approximate the periodical flows in a piezoelectric micropump, with a shallow water equation including nonlinearity and viscous damp presented to characterize the flows in micropump. The finite element method is used to obtain a matrix equation of fluid pressure. The fluid pressure equation is combined with the vibration equation of a silicon diaphragm to construct a liquid-solid coupled equation for reflecting the interaction between solid diaphragm and fluid motion in a micropump. Numerical results of a mode analysis of the coupled system indicate that the natural frequencies of the coupled system are much lower than those of the non-coupled system. The influence of additional mass and viscous damp of fluid on the natural frequencies of the coupled system is more significant as the pump thickness is small. It is found that the vibration shape functions of silicon diaphragm of the coupled system are almost the same as those of the non-coupled system. This paper also gives the first-order amplitude-frequency relationship of the silicon diaphragm, which is necessary for the flow-rate-frequency analysis of a micropump.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.     

Superconvergence of nonconforming finite element penalty scheme for Stokes problem using L 2 projection method

A modified penalty scheme is discussed for solving the Stokes problem with the Crouzeix-Raviart type nonconforming linear triangular finite element. By the L ² projection method, the superconvergence results for the velocity and pressure are obtained with a penalty parameter larger than that of the classical penalty scheme. The numerical experiments are carried out to confirm the theoretical results.     

Variational asymptotic micromechanics modeling of heterogeneous porous materials

This paper presents a numerical technique to predict the effective elastic properties of heterogeneous fluid-filled porous media where the heterogeneity may result from dissimilar solid and fluid phase properties or due to mismatch in porous microstructure. The technique is based on the variational asymptotic method of homogenization where finite element method is employed for discretization. Biot’s theory of poroelasticity is used to describe porous media where both solid and fluid phase motions (u ? U formulation) are considered with associated strain measures. The method estimates the poroelastic constitutive law in single analysis which makes it very efficient compared to other finite element based homogenization techniques. The method is also general enough to compute all 28 elements of an anisotropic constitutive matrix. Other than estimating the effective properties the micro-stress/strain distribution is also obtained at no additional cost.The method is successfully applied for homogenization of porous media, fluid-filled cavity and finally for effective property estimation of bone lamella. In absence of any other direct method of porous media homogenization, the present technique is compared with classical homogenization methods with fluid approximated as solid of very high Poisson’s ratio. The suitability of this approximation and various other alternatives are also discussed. It is shown that the present homogenization method can be an efficient tool for bone property estimation where fluid-filled porous hierarchical micro-/nanostructure must be respected at all steps.     

流体饱和两相多孔介质拟静态问题的有限元解法

给出基于混合物理论的流体饱和两相多孔介质模型,该模型由一可变形固体 一流体相组成。采用Ｇａｌｅｒｋｉｎ加权残值法导出求解拟静态问题的有限元公式,并编制了二维有限元程序。用程序分析了一维和二维问题,得到合理的结果。     

Are fem solutions of incompressible flows really incompressible? (or how simple flows can cause headaches!)

It is generally accepted that mixed and penalty finite element methods can routinely solve the incompressible Navier–Stokes equations. This paper shows by means of simple examples that problems can arise even for the simpler Stokes equations. The causes of the problem fall in either of two categories: round-off and ill conditioning, or a poor choice of pressure discretization. Nonsensical solutions can be obtained. Computation of the discrete divergence of the flow field is a simple and powerful tool to diagnose such conditions. In the first part of the paper several simple techniques for minimizing the effect of round-off are reviewed. In the second part it is shown that, for coupled flow problems, care must be exercised in the choice of the pressure approximation. A unified treatment of various observations by different workers is presented. This should prove useful for general users of the finite element method.     

Finite element analysis of non-Newtonian fluid flow in 2-d branching channel

This paper presents finite element analysis of non-Newtonian fluid flow in 2-d branching channel. The Galerkin method and mixed finite element method are used. Here the fluid is considered as incompressible, non-Newtonian fluid with Oldyord differential-type constitutive equation. The non-linear algebraic equation system which is formulated with finite element method is solved by means of continuous differential method. The results show that finite element method is suitable for the analysis of non-Newtonian fluid flow with complex geometry.