层状压电半空间非轴对称问题的状态变量解

在柱坐标下 ,通过引入状态变量 ,建立了横观各向同性压电介质空间非轴对称问题的状态变量方程 .利用Fourier级数和Hankel变换 ,将文中提出的状态变量方程转换为一阶常微分方程组 .采用求解常微分方程组的矩阵方法 ,得到以状态变量和传递矩阵的乘积的形式表示的单层压电介质的解析解 .给出了状态变量解的应用 ,即利用状态变量解求解半无限压电体在坐标原点作用着垂直集中力Pz,x方向的水平集中力Px 和集中点电荷的解 .由层间完全接触的条件 ,给出了N层压电体非轴对称问题的一般解析表达式 .     

层状压电半空间非轴对称问题的状态变量解

在柱坐标下,通过引入状态变量,建立了横观各向同性压电介质空间非轴对称问题的状态变量方程.利用Fourier级数和Hankel变换,将文中提出的状态变量方程转换为一阶常微分方程组.采用求解常微分方程组的矩阵方法,得到以状态变量和传递矩阵的乘积的形式表示的单层压电介质的解析解.给出了状态变量解的应用,即利用状态变量解求解半无限压电体在坐标原点作用着垂直集中力Pz,x方向的水平集中力Px和集中点电荷的解.由层间完全接触的条件,给出了N层压电体非轴对称问题的一般解析表达式.     

非线性波方程的精确孤立波解

立了一种求解非线性波方程精确孤立波解的双曲函数方法,并在计算机代数系统上加以实现,推导出了一大批非线性波方程的精确孤立波解.方法的基本原理是利用非线性波方程孤立波解的局部性特点,将方程的孤立波解表示为双曲函数的多项式,从而将非线性波方程的求解问题转化为非线性代数方程组的求解问题.利用吴消元法或Gröbner基方法在计算机代数系统上求解非线性代数方程组, 最终获得非线性波方程的精确孤立波解,其中有很多新的精确孤立波解.     

Hamilton体系下含弱粘接复合材料层合板的灵敏度分析研究

基于径向基点插值函数(RPIM),在Hamilton体系下研究了含弱粘接复合材料层合板的灵敏度分析问题．利用弹簧层模型和修正H-R(Hellinger-Reissner)变分原理,推导了可用于含弱粘接复合材料层合板响应和灵敏度分析的混合控制方程,给出了基于该混合控制方程进行灵敏度分析的解析法(AM)、半解析法(SA)和有限差分法(FD)．该混合控制方程的主要优点是可以在进行灵敏度分析过程中避免卷积运算．另外,利用该混合控制方程进行灵敏度分析不仅能够同时得到响应结果和灵敏度系数,而且还考虑了层合板的层间弱粘接问题．     

复合材料力学的Hamilton体系和辛几何方法(Ⅰ)——一般原理*

首次把用于动态体系的Hamilton系统引入到静力学中,建立了与原控制方程相对应的Hami-lton方程,可以对全状态向量分离变量,求出解析解和半解析解,特别适合于求解矩形域平面问题和柱形域空间问题.本文建立了一种求解偏微分方程的新方法,并对复合材料力学中的层合板的弯曲和平面应力问题的求解做了详细说明.     

一类求行波解的线性方法

基于齐次平衡法和李志斌的 tanh函数法 ,本文得到一类简单有效的求解非线性发展方程的线性方法 .这类方法利用非线性发展方程孤立波的局部性特点 ,适当地选取函数 f 和 g,将孤波表示为 f,g的多项式 ,从而将非线性发展方程求解问题转化为非线性代数方程组的求解问题 ,再利用吴消元法求解方程组从而得到非线性发展方程的行波解     

边界节点法计算二维瞬态热传导问题

采用边界节点法(BKM)结合双重互易法(DRM)求解二维瞬态热传导问题.采用差分格式处理时间变量,可将原瞬态热传导方程转化为一系列非齐次修正的Helmholtz方程.随后,方程的解可分为特解和齐次解两部分计算,引入双重互易法在区域内部配点求解方程的特解,采用边界节点法仅需边界配点求解方程的齐次解.给出的数值算例显示该方法计算精度高,适用性好,具有很好的稳定性和收敛性,适合求解瞬态热传导问题.     

非线性波方程准确孤立波解的符号计算

该文将机械化数学方法应用于偏微分方程领域，建立了构造一类非线性发展方程孤立波解的一种统一算法，并在计算机数学系统上加以实现，推导出了一批非线性发展方程的精确孤立波解．算法的基本原理是利用非线性发展方程孤立波解的局部性特点，将孤立波表示为双曲正切函数的多项式．从而将非线性发展方程（组）的求解问题转化为非线性代数方程组的求解问题．利用吴文俊消元法在计算机代数系统上求解非线性代数方程组，最终获得非线性发展方程（组）的准确孤立波解.     

复合材料层合板的分析

在复合材料层合板静力和动力分析方面,本文提出了一个层合板理论.此理论从板的总挠度中分开了由弯曲所产生的挠度W_b与由剪切产生的挠度W_s,因而使得求控制方程的解变得容易了.而且便于讨论和分析横向剪切变形对层合板弯曲,振动和稳定性的影响.     

应用(1/G)-展开法求解五阶KdV方程

本篇论文首次提出(1/G) -展开法,用于求解非线性演化方程的行波解.将该法应用于五阶KdV方程的求解,当参数满足一定条件时,该方程可化为Sawada-Kotera (SK)方程、Caudrey-Dodd-Gibbon(CDG)方程、Kaup-Kupershmidt (KK)方程、Lax方程和Ito方程.其解可被表示为...     

合成展开法求解圆薄板大挠度问题

本文用合成展开摄动法,把外场解和内层解结合起来,求解圆薄板大挠度问题.本文把Hencky的薄膜解当作外场解的一级近似解,并求出了外场解的二级近似解.利用边界内层坐标,求得了相应的各级内层解,即边界层解.本文采用最大位移和板厚之比的倒数作为小参数,所得结果大大改进了1948年作者所得的结果.     

Outer solution for elastic torsion by the method of boundary layer residual states

By the method of boundary layer residual state (BLRS), it is possible to specify the unknown parameters in the general form of the outer asymptotic solution of the governing differential equations for linear boundary value problems (BVP) without any reference to the inner asymptotic solutions of the same problem and the matching procedure. The method accomplishes this task by rationally assigning a portion of the prescribed boundary data to the outer solution. Specifically, the method requires certain weighted averages of the outer solution to be equal to the same averages of the data over the (localized) boundary where the data is prescribed. These weighted averages are consequences of a reciprocity relation inherent in the BVP and the stipulation that the difference between the outer solution and the exact solution (called the residual solution) of the BVP be a boundary layer phenomenon.¶The weighted average requirements are only necessary conditions for the residual state to be a boundary layer. Unfortunately, there are generally countably infinite number of (2) states, many more than the available degrees of freedom in the outer solution to satisfy them. We must show that there is no over-determination or non-uniqueness of the outer asymptotic solution, the abundance of necessary conditions notwithstanding. The present note describes an approach to assuring a well-specified outer solution (up to the expected accuracy) by way of the problem of Saint-Venant torsion. The same approach also also applies to other linear BVP, deducing the appropriate outer solution whenever the determination of the relevant inner solutions is not practical.     

A PRACTICAL PARALLEL DIFFERENCE SCHEME FOR PARABOLIC EQUATIONS

A practical parallel difference scheme for parabolic equations is constructed as follows: to decompose the domain Ω into some overlapping subdomains, take flux of the last time layer as Neumann boundary conditions for the time layer on inner boundary points of subdomains, solve it with the fully implicit scheme on each subdomain, then take correspondent values of its neighbor subdomains as its values for inner boundary points of each subdomain and mean of its neighbor subdomain and itself at overlapping points. The scheme is unconditionally convergent. Though its truncation error is O(τ h), the convergent order for the solution can be improved to O(τ h2).     

Viscoelastic flow in a curved channel: A similarity solution for the Oldroyd-B fluid

A similarity solution is used to analyse the flow of the Oldroyd fluid B, which includes the Newtonian and Maxwell fluids, in a curved channel modelled by the narrow annular region between two circular concentric cylinders of large radius. The solution is exact, including inertial forces. It is found that the non-Netonian kinematics are very similar to the Newtonian ones, although some stress components can become very large. At high Reynolds number a boundary layer is developed at the inner cylinder. The structure of this boundary layer is asymptotically analysed for the Newtonian fluid. Non-Newtonian stress boundary layers are also developed at the inner cylinder at large Reynolds numbers.     

On the shear boundary layer of tension plates

The present paper deals with a generic class of problems for plates subjected to loadings combining a high in-plane tension and a small lateral pressure. It develops the governing differential equations in the singular pertubation form, through the postulation of retaining only one of the Kirchhoff's assumptions, that the plate thickness in the boundary layer region is invariant. The solution by using the standard perturbation method is discussed. The postulation is justified when it is demonstrated that in the shear boundary layer the plate thickness is of higher-order smallness. The general method of solution by the standard perturbation technique is applied to an annular plate problem. Problems of different combinations of supports at the inner and the outer boundaries are solved. The case in which both edges are simply supported is presented as an illustration of the solution technique. In other cases results only are presented. The effect of support on the boundaries is also discussed. The shear effect is found to be most significant at a clamped edge. In the special geometry, it is possible to demonstrate that, when the condition on membrane force is not met as required in the general theory, thagnitude of the boundary layer changes. Specifically, the paper presents a case in which the membrhich the membrane force is zero at the inner edge.     

On the use of a Pohlhausen method in three dimensional boundary layers

Summary The paper deals with the problem of a three dimensional boundary layer flow on the inner surface of a cone when the flow external to the boundary layer is due to a potential line vortex on the axis of the cone together with a radial flow parallel to the cone generators.The momentum integral equations are solved using a Pohlhausen method. An asymptotic solution of the boundary layer equations is used to overcome the arbitrariness which commonly arises in this type of approach and, together with a discussion of certain other boundary conditions, lays a foundation for the use of the method.

---

Zusammenfassung Die Arbeit behandelt das dreidimensionale Grenzschichtproblem der Strömung entlang des Inneren an einem Kegel, wobei die Strömung ausserhalb der Grenzschicht durch einen Potentialwirbel entlang der Kegelachse und durch Radialbewegung entlang der Kegel-Generatoren gebildet wird.Die Momenten-Integral-Gleichung wird mit einer Pohlhausen-Methode gelöst. Eine asymptotische Lösung der Grenzschichtgleichung wird dazu benützt, eine gewisse Willkür zu beseitigen, die gewöhnlich bei solchen Methoden auftritt, und gibt zusammen mit der Diskussion von gewissen Grenzbedingungen die Grundlage der Methode.

     

奇摄动问题中阶梯状解的缝接法

In view of singularly perturbed problems with complex inner layer phenomenon,including contrast structures(step-step solution and spike-type solution),corner layer behavior and right-hand side discontinuity,we carry out the process with sewing connection.The presented method of sewing connection for singularly perturbed equations is based on the two points singularly perturbed simple boundary problems.By means of sewing orbit smoothness,we get the uniformly valid solution in the whole interval.It is easy to prove the existence of solutions and deal with the high dimensional singularly perturbed problems.     

Analytical approximate solution of the cooling problem by Adomian decomposition method

The Adomian decomposition method (ADM) can provide analytical approximation or approximated solution to a rather wide class of nonlinear (and stochastic) equations without linearization, perturbation, closure approximation, or discretization methods. In the present work, ADM is employed to solve the momentum and energy equations for laminar boundary layer flow over flat plate at zero incidences with neglecting the frictional heating. A trial and error strategy has been used to obtain the constant coefficient in the approximated solution. ADM provides an analytical solution in the form of an infinite power series. The effect of Adomian polynomial terms is considered and shows that the accuracy of results is increased with the increasing of Adomian polynomial terms. The velocity and thermal profiles on the boundary layer are calculated. Also the effect of the Prandtl number on the thermal boundary layer is obtained. Results show ADM can solve the nonlinear differential equations with negligible error compared to the exact solution.     

两互相垂直的平面间的层流边界层

本文得到了两互相垂直的平面间的层流边界层的三级近似解.在边界层中,边界层方程中的粘性项和惯性项具有相同的数量级[3].本文则首先假定惯性项大于粘性项去求解边界层方程;然后,令粘性项大于贯性项.最后,取二者的平均值作为边界层方程的真实解.本文所得一级及二级近似解和文献[1]的结果相同.本文的三级近似解则较[1]的结果更精确.     

Approximate method for the numerical solution of singular perturbation problems

We present an approximate method for the numerical solution of linear singularly perturbed two point boundary value problems in ordinary differential equations with a boundary layer on the left end of the underlying interval. It is motivated by the asymptotic behavior of singular perturbation problems. The original problem is divided into inner and outer region problems. The reduced problem is solved to obtain the terminal boundary condition. Then, a new inner region problem is created and solved as a two point boundary value problem. In turn, the outer region problem is also modified and the resulting problem is efficiently treated by employing the trapezoidal formula coupled with discrete invariant imbedding algorithm. The proposed method is iterative on the terminal point. Some numerical experiments have been included to demonstrate its applicability.