Unconditionally Superclose Analysis of a New Mixed Finite Element Method for Nonlinear Parabolic Equation

This paper develops a framework to deal with the unconditional superclose analysis of nonlinear parabolic equation. Taking the finite element pair $Q_{11}/Q_{01} × Q_{10}$ as an example, a new mixed finite element method (FEM) is established and the $τ$ -independent superclose results of the original variable $u$ in $H^1$-norm and the flux variable $\mathop{q} \limits ^{\rightarrow}= −a(u)∇u$ in $L^2$-norm are deduced ($τ$ is the temporal partition parameter). A key to our analysis is an error splitting technique, with which the time-discrete and the spatial-discrete systems are constructed, respectively. For the first system, the boundedness of the temporal errors is obtained. For the second system, the spatial superclose results are presented unconditionally, while the previous literature always only obtain the convergent estimates or require certain time step conditions. Finally, some numerical results are provided to confirm the theoretical analysis, and show the efficiency of the proposed method.     

NUMERICAL SIMULATION OF FIBROBLAST DUROTAXIS MIGRATION INDUCED BY STIFFNESS GRADIENT OF SUBSTRATE1)

建立了一种细胞趋硬性迁移的理论模型和有限元分析框架,为连续变刚度人工基质的试验设计提供理论依据。考虑了细胞体的黏弹性属性,以及细胞与基质间的配受体动态反应过程,并以配受体合成时间为时间步长,将细胞运动方程化为静力学形式进行求解。对有限元过程提出一种动约束,便于消除其结构矩阵的奇异性。结果表明,模型能够模拟黏着斑内部力的快速波动现象,细胞的运动速度与观测数据一致,可有效模拟20,h以上的长时程问题。     

A Low Order Nonconforming Mixed Finite Element Method for Non-Stationary Incompressible Magnetohydrodynamics System

A low order nonconforming mixed finite element method (FEM) is established for the fully coupled non-stationary incompressible magnetohydrodynamics (MHD) problem in a bounded domain in 3D. The lowest order finite elements on tetrahedra or hexahedra are chosen to approximate the pressure, the velocity field and the magnetic field, in which the hydrodynamic unknowns are approximated by inf-sup stable finite element pairs and the magnetic field by $H^1(\Omega)$-conforming finite elements, respectively. The existence and uniqueness of the approximate solutions are shown. Optimal order error estimates of $L^2(H^1)$-norm for the velocity field, $L^2(L^2)$-norm for the pressure and the broken $L^2(H^1)$-norm for the magnetic field are derived.     

有限元/边界元耦合方法在海洋生物目标强度特性研究的应用*

目标强度特性是海洋生物声学识别与资源量评估的重要依据,其中,基于近似几何体和声阻抗特性的理论模型法是研究海洋生物目标强度的重要手段。由于对几何形态近似处理以及数值求解方法的限制,传统理论模型对声波频率、入射方位以及目标声阻抗、形态尺寸等均有各自不同的适用范围,单一模型难以满足不同种类或同一种类但不同尺寸海洋生物的目标强度求解。本文尝试将逐渐见诸应用的有限元/边界元耦合方法用于海洋生物目标强度特性研究,分别以球形生物、纺锤形鱼类尾明角灯鱼(Ceratoscopelus warmingii)和细长形浮游动物南极大磷虾(Euphausia superba)为例进行仿真计算,并与相适应的经典理论模型进行对比分析。结果表明,对于球形生物,有限元/边界元耦合方法与解析模型的目标强度频响曲线完全吻合;对于纺锤形鱼类,有限元/边界元耦合方法可有效弥补基于模态级数解的形变圆柱体模型在中低频和两端入射时的准确性问题;对于细长形浮游动物,有限元/边界元耦合方法与畸变波玻恩近似模型高度吻合。综上,有限元/边界元耦合方法对多种海洋生物目标强度求解均具有较好的适用性,未来有待进一步结合实验测定进行验证。     

考虑膨胀力的非饱和介质热-水-应力耦合二维有限元分析

从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,使用Galerk in方法,将各控制方程分别在空间域和时间域进行离散,开发出了一个可考虑膨胀力的用于分析非饱和介质中热-水-应力耦合弹塑性问题的二维有限元程序.通过对一个假定的核废料地下处置库的热-水-应力耦合问题的数值计算,比较了无、有膨胀力时的情况,在定性上验证了该程序的正确性.     

刚柔耦合建模理论的实验验证

传统的混合坐标建模理论是零次近似方法，对刚柔耦合问题的描述存在缺陷，研究以一个由中心刚体，柔性梁及末端质量组成的刚柔耦合系统对象，建立了精确的一次近似的刚柔耦合动力学方程，在该模型中计及了结构阻尼及风阻的影响，利用单轴气浮台动力学实验平台，通过与实验数据的比较，说明了传统零次近似方法在某些条件下已不能下确描述柔性的刚性正确性和可靠性。     

Dynamics of Multi-Body Rotors: Numerical and Experimental FEM Analysis of the Scientific Earth Experiment Galileo Galilei Ground

Some small scientific satellites can be regarded as free, multi-body, supercritical and statorless rotors. Classical rotordynamics does not cover the dynamic behaviour and the stability analysis of such rotating systems, being focused on fixed rotors. Attitude dynamics usually deals with single rigid spacecrafts, sometimes equipped with flexible appendices like solar arrays or antennas. The case of the Galileo Galilei Ground (GGG) test facility is herewith analyzed through a numerical and experimental investigation on the dynamic behaviour of a fixed multi-body fast-spinning rotor in order to validate the design approach proposed for Galileo Galilei (GG) spacecraft.     

自适应有限元局部P加密的面向对象程序实现

本文采用了面向对象的程序设计方法，以C 作为编程语言，构造了一个局部自适应有限元P加密程序。研究内容涉及从Serendipity函数的构造、单元的形成、自升阶算法以及误差估计和判断。最后采用C 的一系列类及成员函数实现了这些功能，并构造了相应的输入输出模块，使其成为一个完整的自适应有限元分析程序。文中给出的计算结果表明，在P加密有限元分析程序中，误差随着阶数的升高而下降。而且，局部P加密的区域也集中在应力奇异的部位。这验证了文中的理论推导及程序实现的正确性。     

ANALYSIS ON ACOUSTICAL SCATTERING BY A CRACKED ELASTIC STRUCTURE

The acoustical scattering by a cracked elastic structure is studied. The mixed method of boundary element and fractal finite element is adopted to solve the cracked structure-acoustic coupling problem. The fractal two-level finite element method is employed for the cracked structure, which can reduce the degree of freedoms (DOFs) greatly, and the boundary element method is used for the exterior acoustic field which can automatically satisfy Sommerfeld‘s radiation condition. Numerical examples show that the resonance frequency is lower with the crack‘s depth increase, and that the effect on the acoustical field by the crack is particularly pronounced in the vicinity of the crack tip. This mixed method of boundary element and finite element is effective in solving the scattering problem by a cracked structure.     

空间几何构造分析的有限单元法

提出空间杆系几何构造分析的有限单元法，构造了两种单元(链杆单元和准梁单元)的几何约束矩阵，集成为整体矩阵并引入承条件后，通过对其阶数与秩的比较分析确定体系的几何可变性及静定性．本法原理简单，便于计算机实施，结果完备：对于几何不变体系，可指出多余约束的数目；对于几何可变体系，可给出体系的自由度数及相应的运动模态，并确定自由度的常变瞬变性质．