Exact analytic method for solving variable coefficient differential equation

Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations under arbitrary boundary condition. By this method, the general computation format is obtained. Its convergence is proved. We can get analytic expressions which converge to exact solution and its higher order derivatives unifornuy. Four numerical examples are given, which indicate that satisfactory results can be obtained by this method.     

APPROXIMATE SOLUTION FOR BENDING OF RECTANGULAR PLATES Kantorovich-Galerkin’s Method

This paper derives the cubic spline beam function from the generalized beam differential equation and obtains the solution of the discontinuous polynomial under concentrated loads, concentrated moment and uniform distributed by using delta function. By means of Kantorovich method of the partial differential equation of elastic plates which is transformed by the generalized function (δ function and σ function), whether concentrated load, concentrated moment, uniform distributed load or small-square load can be shown as the discontinuous polynomial deformed curve in the x-direction and the y-direction. We change the partial differential equation into the ordinary equation by using Kantorovich method and then obtain a good approximate solution by using Glerkin’s method. In this paper there ’are more calculation examples involving elastic plates with various boundary-conditions, various loads and various section plates, and the classical differential problems such as cantilever plates are shown.     

微分求积区域分裂法在裂缝问题上的应用

微分求积法DQM在处理裂缝问题时,会产生很大的误差。因此,本文用微分求积法结合不带重叠的区域分裂法DQDDM来求解。通过本文的讨论,可以看到DQDDM在处理裂缝问题时,在节点数目不多的条件下获得比较精确的解,同时计算量又不大。     

On one class of differential equations of fractional order

We consider the Darboux problem for a differential equation of fractional order that contains a regularized mixed derivative. Sufficient conditions for the existence and uniqueness of a solution of this problem are obtained in the class of continuous functions. We also propose a method for finding an approximate solution of this problem and prove the convergence of this method.     

变厚度板弯曲的边界元分析法

由于变厚度板弯曲问题的控制分方程复杂，直接求解其基本解推导边界积分方程建立边界元分析法较为困难，本文通过引入等效荷载，等效刚度，将此问题的控制微分方程化成与普通薄板弯曲问题基本方程相同的形式，利用求解通板弯曲问题的边界元迭代求解，建立了分析变厚度板弯曲问题的蛤法，算例表明本方法理正确，精度良好。     

Isospectral Euler-Bernoulli beams via factorization and the Lie method

We obtain isospectral Euler-Bernoulli beams by using factorization and Lie symmetry techniques. The canonical Euler-Bernoulli beam operator is factorized as the product of a second-order linear differential operator and its adjoint. The factors are then reversed to obtain isospectral beams. The factorization is possible provided the coefficients of the factors satisfy a system of non-linear ordinary differential equations. The uncoupling of this system yields a single non-linear third-order ordinary differential equation. This ordinary differential equation, called the principal equation, is analyzed, reduced or solved using Lie group methods. We show that the principal equation may admit a one-dimensional or three-dimensional symmetry Lie algebra. When the principal system admits a unique symmetry, the best we can do is to depress its order by one. We obtain a one-parameter family of invariant solutions in this case. The maximally symmetric case is shown to be isomorphic to a Chazy equation which is solved in closed form to derive the general solution of the principal equation. The transformations connecting isospectral pairs are obtained by numerically solving systems of ordinary differential equations using the fourth-order Runge-Kutta method.     

Exponential dichotomies of nonlinear discrete systems and its application to numerical analysis and computation

In this pepar we consider the upwind difference scheme of a kind of boundary value problems for nonlinear, second order, ordinary differential equations. Singular perturbation method is applied to construct the asymptotic approximation of the solution to the upwind difference equation. Using the theory of exponential dichotomies we show that the solution of an order-reduced equation is a good approximation of the solution to the upwind difference equation except near boundaries. We construct correctors which yield asymptotic approximations by adding them to the solution of the order-reduced equation. Finally, some numerical examples are illustrated.     

Motion of a strong shock front in a two-dimensional gas layer with moving internal boundary

We examine the problem of planar one-dimensional motion of a strong shock wave with moving internal boundary in which the initial position of the front, its intensity, the mass of the gas involved in the motion, and the energy contained in this gas are known. The problem is not self-similar and its exact solution, which involves working with partial differential equations, presents serious difficulties. In the following we determine the law of shock-front motion in this problem via the method of [1], which makes it possible to find a system of ordinary differential equations for the problem. The method is based on an initial specification of the power-law coupling between the dimensionless Lagrangian and Eulerian variables and replacement of the energy equation by this coupling and the energy integral. The solution is sought in the first approximation.     

Differentiator series solution of linear differential ordinary equation

１Ｄｉｆｅｒｅｎｔｉａｔｏｒ,ＩｎｖｅｒｓＯｐｅｒａｔｏｒａｎｄＴｈｅｉｒＰｒｏｐｅｒｔｉｅｓ１．１ＤｉｆｅｒｅｎｔｉａｔｏｒａｎｄｉｎｖｅｒｓｏｐｅｒａｔｏｒＳｕｐｐｏｓｅｔｈｅｌｉｎｅａｒｏｒｄｉｎａｒｙｄｉｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｏｆ...     

Uniform convergence for the difference scheme in the conservation form of ordinary differential equation with a small parameter

In this paper, we consider a singular perturbation boundary problem for a self-adjoint ordinary differential equaiton. We construct a class of difference schemes with fitted factors, and give the sufficient conditions under which the solution of difference scheme converges uniformly to the solution of differential equation. From this we propose several specific schemes under weaker conditions, and give much higher order of uniform convergence.     

A new semi-analytical method for the non-linear static analysis of an infinite beam on a non-linear elastic foundation: A general approach to a variable beam cross-section

We propose a new non-linear method for the static analysis of an infinite non-uniform beam resting on a non-linear elastic foundation under localized external loads. To this end, an integral operator equation is newly formulated, which is equivalent to the original differential equation of non-uniform beam. By using the integral operator equation, we propose a new functional iterative method for static beam analysis as a general approach to a variable beam cross-section. The method proposed is fairly simple as well as straightforward to apply. An illustrative example is presented to examine the validity of the proposed method. It shows that just a few iterations are required for an accurate solution.     

The uniformly convergent difference schemes for a singular perturbation problem of a self-adjoint ordinary differential equation

In this paper, we construct a class of difference schemes with fitted factors for a singular perturbation problem of a self-adjoint ordinary differential equation. Using a different method from [1], by analyzing the truncation errors of schemes, we give the sufficient conditions under which the solution of the difference scheme converges uniformly to the solution of the differential equation. From this we propose several specific schemes under weaker conditions, and give much higher order of uniform convergence, and applying them to example, obtain the numerical results.     

Exact finite element method

In this paper,a new method,exact element method for constructing finite element,ispresented.It can be applied to solve nonpositive definite or positive definite partialdifferential equation with arbitrary variable coefficient under arbitrary boundarycondition.Its convergence is proved and its united formula for solving partial differentialequation is given.By the present method,a noncompatible element can be obtained and thecompatibility conditions between elements can be treated very easily.Comparing the exactelement method with the general finite element method with the same degrees of freedom,the high convergence rate of the high order derivatives of solution can be obtained.Threenumerical examples are given at the end of this paper,which indicate all results canconverge to exact solution and have higher numerical precision.     

Application of equation regulation method to the solution of inverse problems in multiphase materials

In this paper, we propose a new method of determining local material properties of multiphase composites given the experimentally measured displacements and known traction boundary conditions. In the proposed method, an “observation” term is added to the original differential equation, and the modified equation is solved in terms of a regulation parameter. We call this approach the equation regulation (ER) method. By appropriately adjusting the value of the regulation parameter based on the noise level in the input data, we get faster convergence and improved stability than prevailing methods of solving the inverse problem in elliptic ordinary differential equations. Several numerical examples to the solution of this non-linear problem with continuous and discontinuous coefficient functions are given to show the accuracy and reliability of the proposed method.     

边界节点法的计算稳定性研究

边界节点法利用满足控制方程的非奇异通解作为基函数,半解析边界数值离散偏微分方程,具有精度高、收敛快、易编程等优点,是一种纯无网格配点方法.但是在求解具体问题时,随着节点数的增加,边界节点法经常得到严重病态的插值矩阵.本文利用有效条件数评价边界节点法求解Helmholtz问题线性方程组的计算稳定性;然后利用三种正则化方法处理其病态的线性方程组,并与高斯消元法比较计算精度和收敛性.通过数值实验,本文研究了有效条件数、误差和正则化方法之间的关系.     

THE EXPLICIT SOLUTION OF HOMOGENEOUS LINEAR ORDINARY DIFFERENTIAL EQUATIONS WITH CONSTANT COEFFICIENTS

In this paper by using tensor analysis we give the explicit expressions of the solution of the initial-value problem of homogeneous linear differential equations with constant coefficients and the n th-order homogeneous linear differential equation with constant coefficients. In fact, we give the general formula for calculating the elements of the matrix exp[At]. We also give the results when the characteristic equation has the repeated roots. The present method is simpler and better than the other methods.     

人行荷载下梁式结构振动分析的DQ-IQ混合法

为了评估人行荷载作用下梁式结构的振动舒适度，利用微分求积-积分求积，即DQ-IQ混合法求解移动荷载作用下梁的振动响应。人行荷载作用下梁式结构的振动控制方程是含Dirac函数的偏微分方程，首先利用IQ法离散与时间相关的Dirac函数，再利用DQ法把控制方程转化为二阶常系数微分方程，最后利用Newmark算法求解微分方程。以某钢结构连廊为例，利用DQ法计算结构自振频率并与解析解进行对比，结果验证了节点选取和边界条件施加的合理性，再利用DQ-IQ混合法和振型叠加法分别计算了不同行走步频下连廊的响应，计算结果表明，DQ-IQ混合法具有较高的可靠性和精确性。DQ-IQ混合法也可以推广到诸如车辆荷载作用下路面或桥梁的动力响应等其他移动荷载下结构的振动分析。     

Parabolic method of solving problems of the flow of a sonic gas stream around a thin symmetric profile

A parabolic method consisting of replacement of the stream acceleration ?_xx in the non-linear member of (1.1) by a specially chosen constant has been proposed [1] for the solution of the mixed-type transonic equation with boundary conditions on the profile, and the solution of the linear parabolic-type equation obtained can be considered as a certain approximation to the solution of the initial problem. An improvement of the parabolic method is the method of local linearization [2] (see [3] also), in which the acceleration ?_xx fixed from the beginning is replaced by a function of the coordinate x which satisfies some condition. An ordinary first-order differential equation is obtained for the velocity distribution along the profile in [2]. Another method of “defrosting” the acceleration ?_xx “frozen” from the beginning is proposed in this paper; a second-order ordinary differential equation is obtained for the velocity on the profile, which permits getting rid of some disadvantages of the local linearization method. Several solutions of (2.3) are presented in comparison to known exact solutions.     

The explicit solution of homogeneous linear oridinary differential equations with constant coefficients

In this paper by using tensor analysis we give the explicit expressions of the solution of the initial-value problem of homogeneous linear differential equations with constant coefficients and the n th-order homogeneous linear differential equation with constant coefficients. In fact, we give the general formula for calculating the elements of the matrix exp [At]. We also give the results when the characteristic equation has the repeated roots. The present method is simpler and better than, the other methods.     

Solving Huxley equation using an improved PINN method

Although many effective methods for solving partial differential equations (PDEs) have been proposed, there is no universal method that can solve all PDEs. Therefore, solving partial differential equations has always been a difficult problem in mathematics, such as deep neural network (DNN). In recent years, a method of embedding some basic physical laws into traditional neural networks has been proposed to reveal the dynamic behavior of equations directly from space-time data [i.e., physics-informed neural network (PINN)]. Based on the above, an improved deep learning method to recover the new soliton solution of Huxley equation has been proposed in this paper. As far as we know, this is the first time that we have used an improved method to study the numerical solution of the Huxley equation. In order to illustrate the advantages of the improved method, we use the same network depth, the same hidden layer and neurons contained in the hidden layer, and the same training sample points. We analyze the dynamic behavior and error of Huxley’s exact solution and the new soliton solution and give vivid graphs and detailed analysis. Numerical results show that the improved algorithm can use fewer sample points to reconstruct the exact solution of the Huxley equation with faster convergence speed and better simulation effect.