石英晶体板非线性高频振动的拓展伽辽金法分析

针对考虑几何和材料非线性的石英晶体板厚度剪切振动和弯曲振动的方程组,利用扩展伽辽金法对该方程组进行转化和求解,分别获得了强烈耦合的厚度剪切振动模态和弯曲振动模态的频率响应关系,绘制了不同振幅比和不同驱动电压影响下的频率响应曲线图。数值计算结果表明可以选取石英晶片的最佳长厚比尺寸来避免两种模态的强烈耦合。驱动电压的变化将引起石英晶体谐振器厚度剪切振动频率的明显改变,必须将振动频率的漂移值控制在常用压电声波器件的允许值之内。扩展伽辽金法对石英晶体板非线性振动方程组的求解为非线性有限元分析和偏场效应分析奠定了基础。     

Nonlinear elastic constitutive relations of residually stressed composites with stiff curved fibres

Applied Mathematics and Mechanics - Mechanical models of residually stressed fibre-reinforced solids, which do not resist bending, have been developed in the literature. However, in some residually...     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.     

A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Euler Characteristics on virtually free products

We define Euler characteristics on classes of residually finite and virtually torsion free groups and we show that they satisfy certain formulas in the case of amalgamated free products and HNN extensions over finite subgroups. These formulas are obtained from a general result which applies to the rank gradient and the first L²?Betti number of a finitely generated group.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

Handling contacts in an Eulerian frame: a finite element approach for fluid structures with contacts

ABSTRACT

Eulerian variational formulations for deformable solids, with or without fluids around them, end up, after implicit time discretisation, as large non-linear systems for the velocities in the moving domains. Handling moving domains and moving boundaries requires careful meshing procedures; on the other hand, the detection of contact is particularly simple with a distance function. Then at every time step, a variational inequality can be used to update the velocities. This article gives new implementation details and two new complex simulations: a very soft bouncing ball in an axisymmetric flow and a disk hit by a club.     

围压与径向荷载共同作用下巴西盘裂纹应力强度因子的解析解

利用权函数法推导了围压和径向荷载共同作用下，考虑裂纹面摩擦的预制裂纹巴西盘应力强度因子计算公式，从理论上分析了围压、径向荷载和裂纹面摩擦对巴西盘应力强度因子的影响。结果表明，围压对I型应力强度因子有很大影响，I型应力强度因子随围压的增大而减小。当裂纹面闭合后围压和摩擦系数对II型应力强度因子同样具有显著影响，考虑裂纹面有效剪应力的权函数法理论解与有限元数值解相吻合，表明理论分析的正确性。     

Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

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

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