复合材料反平面切口端部奇性研究

基于复合材料切口尖端位移场的渐近展开,将切口的反平面平衡控制方程转化为关于切口奇性指数的特征微分方程,采用一种变换将其化为线性特征微分方程组,引入插值矩阵法计算相应边界条件下方程组的特征值以获取切口尖端的应力奇性指数.研究单相材料切口、双相材料切口及止于异质界面切口的奇异性,算例表明该方法可以一次性计算出多阶奇性指数.对所取得的非奇异指数尽管切口不表现出奇性状态,但却是描述切口尖端完整应力场必不可少的参量.     

正交各向异性材料切口尖端热弹奇性特征分析

研究正交各向异性平面V形切口,计算其热弹奇性特征.通过引入切口尖端物理场的渐近级数展开式,将应力和热流平衡方程转化为关于奇性指数的特征常微分方程组,采用插值矩阵法求解,获取切口尖端的热流、应力奇性指数和对应的特征角函数.算例表明,该法精度高适应性强.     

基于平面问题的位移压力混合配点法

引入压力变量,将弹性力学控制方程表达为位移和压力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式.采用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法求解过约束方程组,得到平面问题位移数值解.数值算例验证了所提方法的有效性和计算精度.     

标量波传播问题的双互易精细积分法

为避免使用计算多种特征频率下的声场响应,采用双互易方法将边界积分方程中时间二次导数项的域积分转化为边界积分.首先,将计算场点配置在边界上并考虑边界条件,可以获得由内部节点上声压量线性表示的边界节点上的物理量;其次,将计算场点配置于域内离散节点上,将所得边界积分方程组中关于边界物理量用内部节点的声压量线性表示,获得关于声压量的二阶常微分方程组;第三,引入声压变化速度作为未知量,将二阶常微分方程组转化为一阶常微分方程组;最后,采用精细积分法精确求解常微分方程组.数值算例验证了双互易精细积分法的正确性和稳定性.     

粒子输运离散纵标方程基于界面修正的并行计算方法

为了改造粒子输运方程求解的隐式格式,研究设计适应大型并行计算机的并行计算方法,介绍一类求解粒子输运方程离散纵标方程组的基于界面修正的源迭代并行计算方法.应用空间区域分解,在子区域内界面处首先采用迎风显式差分格式进行预估,构造子区域的入射边界条件,然后,在各个子区域内部进行源迭代求解隐式离散纵标方程组.在源迭代过程中,在内界面入射边界处采用隐式格式进行界面修正.数值算例表明该并行计算方法在精度、并行度、简单性诸方面均具有良好的性质.     

含双时间步法的化学非平衡流解耦算法

发展基于隐式双时间步法的化学非平衡流解耦型计算方法.采用算子分裂法对流动和反应进行解耦处理,流动方程组通过双时间步方法求解;源项方程组采用二阶梯形公式迭代求解;提出"源项消去"法,以消除化学反应源项对流动求解引入的误差,从而保证流动方程组求解的时间精度.理论分析和计算结果表明,方法既可以保证双时间步法的求解效率,又可以获得比较精确的非定常计算结果.     

正交各向异性板液/固界面的声反射与声透射

采用Legendre正交多项式法,对液/固界面声波反射和透射系数进行求解。利用Legendre正交多项式对正交各向异性板中的位移解进行展开,推导出板中的应力和波动控制方程。联立液/固界面的边界条件和波动控制方程,建立线性无关方程组,用以计算液/固界面的反射和透射系数及Legendre多项式的展开系数,计算所得铝板液/固界面的反射和透射系数与传递矩阵法的计算结果吻合良好。以单向纤维增强复合材料板为例,在不同的方位角下,分析了反射和透射系数随斜入射角度、入射波频率的变化关系,以及板中声场的位移分布。所取Legendre正交多项式截止阶数越大,可用来计算的频厚积范围越大。研究拓展了Legendre正交多项式法的适用范围,为材料力学性能的声学测量提供了理论基础。     

正交各向异性板液/固界面的声反射与声透射

采用Legendre正交多项式法,对液/固界面声波反射和透射系数进行求解。利用Legendre正交多项式对正交各向异性板中的位移解进行展开,推导出板中的应力和波动控制方程。联立液/固界面的边界条件和波动控制方程,建立线性无关方程组,用以计算液/固界面的反射和透射系数及Legendre多项式的展开系数,计算所得铝板液/固界面的反射和透射系数与传递矩阵法的计算结果吻合良好。以单向纤维增强复合材料板为例,在不同的方位角下,分析了反射和透射系数随斜入射角度、入射波频率的变化关系,以及板中声场的位移分布。所取Legendre正交多项式截止阶数越大,可用来计算的频厚积范围越大。研究拓展了Legendre正交多项式法的适用范围,为材料力学性能的声学测量提供了理论基础。      

横向各向同性柱状复合结构对超声波的散射

首先描述了横向各向同性复合圆柱结构的入射声场、散射声场及内部的驻波声场，然后利用转移矩阵方法导出了求解散射声场的方程组，计算了铝／各向异性界面层／纤维复合结构对斜入射声波的背向散射谱和散射截面积。将柱状复合结构中各向异性界面薄层相应的转移矩阵作渐近展开，建立了模拟这种界面薄层的弹簧模型及界面处广义边界条件。结果表明，模型中劲度常数仅依赖于界面薄层厚度及界面层媒质的弹性常数Ｃｌｌ，Ｃ１２和Ｃ４４，而振子质量与Ｃｌｌ，Ｃ４４，Ｃ１３和Ｃ３３有关。     

流体力学方程的间断有限元方法

在二维区域三角形网格上应用一阶、二阶和三阶精度间断有限元方法,对流体力学方程和方程组进行了数值模拟.计算结果与差分方法计算结果比较,认为间断有限元方法在求解复杂边界条件和区域问题上有一定的优势.     

Symplectic model for piezoelectric wedges and its application in analysis of electroelastic singularities

In this paper, a symplectic model, based on the Hamiltonian system, is developed for analyzing singularities near the apex of a multi-dissimilar piezoelectric wedge under antiplane deformation. The derivation is based on a modified Hellinger–Reissner generalized variational principle or a differential equation approach. The study indicates that the order of singularity depends directly on the non-zero eigenvalue of the proposed Hamiltonian operator. Using the coordinate transformation technique and continuity conditions on the interface between two dissimilar materials, the orders of singularity for multi-dissimilar piezoelectric and piezoelectric–elastic composite wedges are determined. Numerical examples are considered to show potential applications and validity of the proposed method. It is found that the order of singularity also depends on the piezoelectric constant, in addition to the geometry and shear modulus.     

Melting phenomenon in magneto hydro-dynamics steady flow and heat transfer over a moving surfacein the presence of thermal radiation

The Lie group method is applied to present an analysis of the magneto hydro-dynamics(MHD) steady laminar flow and the heat transfer from a warm laminar liquid flow to a melting moving surface in the presence of thermal radiation.By using the Lie group method,we have presented the transformation groups for the problem apart from the scaling group.The application of this method reduces the partial differential equations(PDEs) with their boundary conditions governing the flow and heat transfer to a system of nonlinear ordinary differential equations(ODEs) with appropriate boundary conditions.The resulting nonlinear system of ODEs is solved numerically using the implicit finite difference method(FDM).The local skin-friction coefficients and the local Nusselt numbers for different physical parameters are presented in a table.     

An Alternative Algorithm for the Symmetry Classification of Ordinary Differential Equations

This is the first paper on symmetry classification for ordinary differential equations (ODEs) based on Wu’s method. We carry out symmetry classification of two ODEs, named the generalizations of the Kummer-Schwarz equations which involving arbitrary function. First, Lie algorithm is used to give the determining equations of symmetry for the given equations, which involving arbitrary functions. Next, differential form Wu’s method is used to decompose determining equations into a union of a series of zero sets of differential characteristic sets, which are easy to be solved relatively. Each branch of the decomposition yields a class of symmetries and associated parameters. The algorithm makes the classification become direct and systematic. Yuri Dimitrov Bozhkov, and Pammela Ramos da Conceição have used the Lie algorithm to give the symmetry classifications of the equations talked in this paper in 2020. From this paper, we can find that the differential form Wu’s method for symmetry classification of ODEs with arbitrary function (parameter) is effective, and is an alternative method.     

Eigenvalues of structures with uncertain elastic boundary restraints

The eigenvalue problems of structures with random elastic boundary supports are studied in this paper. Using the perturbation method, the differential equations including stochastic distributed parameters and random boundary conditions that govern the eigenproblems are transformed to a series of deterministic differential equations and boundary conditions. Then the differential equations and boundary conditions are discretized utilizing the finite element method (FEM). This process is easy to be implemented since the resulting perturbation equations with different orders own the same FEM meshes. The first-order and second-order sensitivities of eigenvalues are derived through solving the deterministic algebraic equations. Upon determining these sensitivities of eigenvalues, the approximate statistic expressions of random eigenvalues considering uncertain elastic boundary supports are established. At the end, several numerical examples are given to illustrate the application and effectiveness of the present method. Comparison of the present numerical results with those from the Monte-Carlo simulation method verifies the accuracy of the developed method.     

A search for the simplest chaotic partial differential equation

A numerical search for the simplest chaotic partial differential equation (PDE) suggests that the Kuramoto-Sivashinsky equation is the simplest chaotic PDE with a quadratic or cubic nonlinearity and periodic boundary conditions. We define the simplicity of an equation, enumerate all autonomous equations with a single quadratic or cubic nonlinearity that are simpler than the Kuramoto-Sivashinsky equation, and then test those equations for chaos, but none appear to be chaotic. However, the search finds several chaotic, ill-posed PDEs; the simplest of these, in the discrete approximation of finitely many, coupled ordinary differential equations (ODEs), is a strikingly simple, chaotic, circulant ODE system.     

The analysis of the flow of a micropolar fluid between two orthogonally moving porous disks with counter rotating directions

The present work investigates the unsteady, imcompressible flow of a micropolar fluid between two orthogonally moving porous coaxial disks. The lower and upper disks are rotating with the same angular speed in counter directions. The flows are driven by the contraction and the rotation of the disks. An extension of the Von Kármán type similarity transformation is proposed and is applied to reduce the governing partial differential equations (PDEs) to a set of non-linear coupled ordinary differential equations (ODEs) in dimensionless form. These differential equations with appropriate boundary conditions are responsible for the flow behavior between large but finite coaxial rotating disks. The analytical solutions are obtained by employing the homotopy analysis method. The effects of some various physical parameters like the expansion ratio, the rotational Reynolds number, the permeability Reynolds number, and micropolar parameters on the velocity fields are observed in graphs and discussed in detail.     

Interaction of the near field with a gyrotropic sectorial structure

A matrix method is proposed to obtain the characteristic equations for quantities governing the degree of the field singularity on the edge of sectorial gyrotropic structures. Exact analytic solutions are analyzed for giving mixed boundary conditions on the outer faces of a wedgelike domain.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 10, pp. 64–69, October, 1987.