无拉力Winkler地基上自由边矩形薄板的弯曲

本文应用Fourier级数加补充项的方法求解了无拉力Winkler地基上自由边矩形薄板的弯曲问题.通过适当设定满足可导条件的Fourier级数加补充项形式的挠度函数,把给定边界条件下的微分方程化成一个无穷代数方程组.因接触区的边界预先不能定出,故这组方程为弱非线性方程.使用迭代法获得解答.     

Winkler地基上固支薄板自由振动问题的准Green函数方法

将准Green函数方法应用于求解Winkler地基上固支薄板的自由振动问题.即利用问题的基本解和边界方程构造一个准Green函数,这个函数满足了问题的齐次边界条件.采用Green公式,将Winkler地基上固支薄板自由振动问题的振型控制微分方程化为第二类Fredholm积分方程.通过边界方程的适当选择,积分方程核的奇异性被克服了.数值算例表明,该方法具有较高的精度,是一种有效的数学方法.     

锥壳的位移解及应用*

本文提出求锥壳方程通解的另一种方法——位移法.文中根据文献[1]给出的壳体基本关系,导出锥壳一般弯曲问题的位移方程组,然后通过引入一个位移函数U(s,θ)(在极限情况下,就变为对于圆柱壳所引入的位移函数),从而将锥壳基本方程组化成关于位移函数U(s,θ)的8阶可解偏微分方程(控制微分方程).对于一般弯曲问题,该方程的一般解以广义超几何函数给出;对于轴对称弯曲问题,用Bessel函数给出其一般解.作为锥壳位移解法的应用,讨论了Winkler地基模式上的锥壳的轴对称弯曲问题,给出数值结果.     

基于凝聚函数的互补问题的光滑化算法

对于不可微的"极大值"形式的函数,可以利用凝聚函数对其进行光滑逼近.借助这个技术,给出了求解线性互补问题的光滑方程组算法.首先是将互补问题转化为等价的非光滑方程组,再利用凝聚函数进行光滑逼近,从而转化为光滑方程组的求解问题.通过一些考题对这个算法进行了数值试验,结果显示了该算法的有效性和稳定性.     

电力系统中地网腐蚀诊断的一种新的数学模型及其仿真计算

对电力系统中具有重大应用价值的地网腐蚀诊断问题抽象出仿真求解的一种新的数学模型:即求解带约束的非线性隐式方程组模型.但由于问题本身的物理特性决定了所建立的数学模型具有以下特点:一是非线性方程组为欠定方程组,而且非线性程度非常高;二是方程组的所有函数均为隐函数;三是方程组附加若干箱约束条件.这种特性给模型分析与算法设计带来巨大困难.对于欠定方程组的求解,文中根据工程实际背景,尽可能地扩充方程的个数,使之成为超定方程组,然后对欠定方程组和超定方程组分别求解并进行比较.将带约束的非线性隐函数方程组求解问题,转化为无约束非线性最小二乘问题,并采用矩阵求导等技术和各种算法设计技巧克服隐函数的计算困难,最后使用拟牛顿信赖域方法进行计算.大量的计算实例表明,文中所提出的数学模型及求解方法是可行的.与目前广泛采用的工程简化模型相比较,在模型和算法上具有很大优势.     

不可压磁流体力学方程组的高精度数值解法

为数值求解低雷诺数下不可压流体在电磁场作用下的流动,提出一种四阶紧致差分方法.由二维原始变量的MHD方程组出发,推导出具有较少未知量的电流密度-涡量-流函数形式MHD方程组.建立了求解二维非定常不可压MHD方程组的电流密度-涡量-流函数形式的四阶精度紧致差分格式.为验证本文提出的高精度紧致差分方法的精确性和可靠性,对有...     

薄板弯曲自由振动问题的高精度近似解析解及改进研究

对于薄板弯曲自由振动问题,已有如下方法:在Hamilton(哈密顿)体系下基于分离变量法得到挠度的解析形式,并建立自振频率联立方程组,给出求解振动频率和振型函数的方法.笔者指出该方法中所用挠度函数的解析式实际上是一种满足位移边界条件的高精度近似解,基于Rayleigh-Ritz(瑞利-里茨)法再次求近似频率后发现,原方法的近似解的精度很高.另外,对于含有固支、简支等不同的边界形式,恰当地选取不同位置作为坐标系的原点,得到含有频率的方程组的统一形式,且较为简洁.这些形式可基于四边固支、四边简支等边界条件的矩形板研究,依照板变形的对称性可验证频率方程组形式的正确性,并得到不同边界条件下频率方程形式之间的联系与转化.     

一个求解广义圆锥互补问题的无导数光滑算法

本文研究了一个求解广义圆锥互补问题的无导数光滑算法.利用光滑函数将广义圆锥互补问题等价转化成一个光滑方程组,然后再利用牛顿法求解此方程组.该算法采用了一种新的非单调无导数线搜索技术,并且在适当条件下具有全局和局部二次收敛性质.数值实验结果表明算法是非常有效的.     

一个求解互补问题的光滑Newton方法

１．引言 考虑非线性互补问题ＮＣＰ（Ｆ）：其中 Ｆ： 是连续可微函数．目前比较流行的求解ＮＣＰ（Ｆ）的方法之一是首先把它转化为一个方程组,然后通过求解方程组的方法［１］间接求解,这样的方法通常是通过Ｆｉｓｃｈｅｒ函数来完成的［２］容易验证所以求解ＮＣＰ（Ｆ）可以等价求解一个ｎ维方程组 然而函数φ有一个缺点,即它在零点不可微．这就导致Φ在某些点不可微．因此传统的求解方程组的方法并不能直接应用到Φ上．为克服这个缺点,可使用它的光滑形式［４］： 我们注意到,只要μ＞０,φμ就是可微的,而且对任意μ有所以可…     

变系数李方程组的精确行波解

利用修正的简单方程法对变系数李方程组进行求解,给出了变系数李方程组的双曲函数形式的行波解,当参数取特殊值时,便可以得到该方程组的精确孤波解.     

Differential quadrature method based on the highest derivative and its applications

The numerical differentiation is often used when dealing with the differential equations. Using the numerical differentiation, the differential equations can be transformed into algebraic equations. Then we can get the numerical solution from the algebraic equations. But the numerical differentiation process is very sensitive to even a small level of errors. In contrast, it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, we provide a new method using the DQ method based on the interpolation of the highest derivative (DQIHD) for the differential equations. The original function is then obtained by integration. In this paper, the DQIHD method was applied to the buckling analysis of thin isotropic plates and Winkler plates, the numerical results agree well with the analytic solutions, and the results show that our method is of high accuracy, of good convergence with little computational efforts. And it is easy to deal with the boundary conditions.     

Development of a higher order finite element on a Winkler foundation

Analyzing thick plates as a construction component has been of interest to structural engineering research for several decades. In particular, thick plates resting on elastic foundations are more specific. Mindlin's plate theory for thick plate analysis and the Winkler theory for elastic foundation analyses have wide applications. The current research considers analysis of isotropic plates on a Winkler foundation according to Mindlin's plate theory. The analysis uses a higher order plate element to avoid shear locking phenomena in the plate. The main features of this element are representation of real displacement functions of the plate perfect and shear locking do not occur at the plates modeled with this element. Derivation of the equations for finite element formulation for thick plate theory uses fourth-order displacement shape functions. A computer program using the finite element method, coded in C++, analyzes the plates resting on an elastic foundation. The analysis involves a 17-noded finite element. The study's graphs and tables assist engineers' designs of thick plates resting on elastic foundations. The study concludes with the computer-coded program, which allows effective use for the shear locking-free analysis of thick Mindlin plates resting on elastic foundations.     

关于一类有限格蕴涵代数方程的解法

研究一类含有有限个未知量且含有蕴涵算子的格蕴涵代数方程,给出方程有解的充分必要条件。在方程的解集非空时,讨论解集的一些结构性质。并刻画出方程在具有某种限制条件下的整个解集。最后,相关结论被应用到一个数值例子。     

A nonlinear model for footings on granular bed–stone column reinforced earth beds

In the present study, an attempt has been made to model and analyze a combined footing supporting column, which is to be constructed on very soft soil. In view of small bearing capacity and very large deflections, foundations are constructed after improving the original ground. Here, the ground has been improved by providing the stone columns in the soft soil and on top of this improved ground; a granular fill layer has been placed just below the footing. The footing has been modeled as a beam having finite flexural rigidity. Granular fill layer, soft soil and stone columns have been represented by Pasternak shear layer, Kelvin–Voigt body and the Winkler springs, respectively. Nonlinear behavior of these has been considered by means of hyperbolic constitutive relationships. Governing differential equations for response of the system have been derived and presented in non-dimensional form. These equations have been solved using appropriate boundary conditions by means of an iterative Gauss Elimination technique.     

On One-Parameter Families of Painlevé III

Albrecht, Mansfield, and Milne developed a direct method with which one can calculate special integrals of polynomial type (also known as one-parameter family conditions, Darboux polynomials, eigenpolynomials, or algebraic invariant curves) for nonlinear ordinary differential equations of polynomial type. We apply this method to the third Painlevé equation and prove that for the generic case, the set of known one-parameter family conditions is complete.     

Free vibration of functionally graded rectangular plates using first-order shear deformation plate theory

The main objective of this research work is to present analytical solutions for free vibration analysis of moderately thick rectangular plates, which are composed of functionally graded materials (FGMs) and supported by either Winkler or Pasternak elastic foundations. The proposed rectangular plates have two opposite edges simply-supported, while all possible combinations of free, simply-supported and clamped boundary conditions are applied to the other two edges. In order to capture fundamental frequencies of the functionally graded (FG) rectangular plates resting on elastic foundation, the analysis procedure is based on the first-order shear deformation plate theory (FSDT) to derive and solve exactly the equations of motion. The mechanical properties of the FG plates are assumed to vary continuously through the thickness of the plate and obey a power law distribution of the volume fraction of the constituents, whereas Poisson’s ratio is set to be constant. First, a new formula for the shear correction factors, used in the Mindlin plate theory, is obtained for FG plates. Then the excellent accuracy of the present analytical solutions is confirmed by making some comparisons of the results with those available in literature. The effect of foundation stiffness parameters on the free vibration of the FG plates, constrained by different combinations of classical boundary conditions, is also presented for various values of aspect ratios, gradient indices, and thickness to length ratios.     

正规升列在参数代数方程组求解中的应用

本文提出一种利用多项式系统的正规零点分解的算法来求解代数方程组以及带有参数的代数方程组的方法．对于给定的的代数方城组,通过正规分解,可以得到一组具有三角形式的分解．根据这种三角形式,我们可以给出代数方程组的所有解．而对于带有参数方程组,将给出方程组有解时参数需满足的条件．进一步,对于给定的参数值,正规分解中得到三角形式仍然保持,通过求解三角形式的方程组从而得出原参数方程组的解．     

On The Navier-Stokes Equations with Constant Total Temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.     

Differential algebraic equations with after-effect

In this paper, we are concerned with the solution of delay differential algebraic equations. These are differential algebraic equations with after-effect, or constrained delay differential equations. The general semi-explicit form of the problem consists of a set of delay differential equations combined with a set of constraints that may involve retarded arguments. Even simply stated problems of this type can give rise to difficult analytical and numerical problems. The more tractable examples can be shown to be equivalent to systems of delay or neutral delay differential equations. Our purpose is to highlight some of the complexities and obstacles that can arise when solving these problems, and to indicate problems that require further research.