一种非对称大型有限元方程组的有效分块波前解法

提出分块波前法求解非对称大型有限元方程组，该方法综合分块解法和波前法的优点，根据计算机容量和方程的阶数优化计算过程，以在经济效率和存储效率之间达到平衡。本文利用该方法求解了轮胎稳态滚动有限元方程组，表明该方法是有效、可靠的。     

边界元法中线性方程组的拟波阵分块求解

以边界单元法形成的线性方程组为例提出了一种分块求解具有非线称满系数矩阵方程组的拟波阵方法。这种在于高斯消去法基本原理的分块求解法,采用了拟波阵求解技术,使常驻计算机内存的仅为二行分块的系数子矩阵。     

充气结构与流场的耦合求解方法

针对一种充气前缘(inflatable leading edge, ILE)增升技术,建立了其充气结构与流场耦合作用的运动方程. 将方程写成状态空间形式,采用时域推进方法求解. 对使用了变前缘增升技术的NACA63-212翼型进行了充气结构静变形的数值计算,结果表明充气结构的刚度对翼型的气动特性有明显影响. 与原翼型相比,在不考虑充气结构变形时, 该增升技术大约能使翼型的失速迎角增加30{\%},最大升力系数增加22{\%};考虑结构变形后增升效果有所降低. 刚度较低的薄膜在前缘吸力峰的作用下会隆起形成鼓包,容易引起流动分离.      

非对称荷载下内撑式基坑支护结构参数分析及设计优化探讨

非对称荷载作用下,内撑式基坑会向荷载小侧整体偏移,这将不利于基坑的稳定,因此有必要对支护结构进行设计优化。本研究以某偏压综合管廊基坑为例,通过有限元软件PLAXIS 2D对该基坑进行了仿真模拟,模拟结果与实测数据吻合良好,验证了模型的有效性。进一步对该基坑进行参数分析,研究了支护结构长度、等效厚度和弹性模量对支护结构变形和基坑安全性的影响。研究结果表明,在两侧支护结构均按偏压侧进行设计的基础上,增大某侧支护结构长度,会使该侧支护结构水平位移最大值增大,分别增大两侧支护结构等效厚度、弹性模量,均可以减小该侧支护结构水平位移最大值和非偏压侧支护结构顶部逆向位移。基于参数分析,提出一种针对非对称荷载基坑的简单有效、经济可行的优化思路,并通过数值模拟对优化效果进行了分析,证明了该优化思路的可行性。     

大型边界元方程组的并行直接分块求解算法

针对大型边界元方程组和网络微机机群环境提出了一种并行直接分块求解算法，算法基于分块高斯-若当消去法的原理，采用内外存交互技术，并行分块消去方法，节点超行的卷帘存储方案和并行环状循环逐次修正策略，增大了解题规模，提高了计算速度。算例计算结果表明该算法具有较高的并行加速比和并行效率，适用于大型问题的边界元法求解。     

弹性地基上矩形薄板问题的Hamilton正则方程及解析解

利用辛算法求出弹性地基上矩形薄板问题的解析解,将弹性地基视为双参数弹性地基,直接从弹性矩形薄板的控制方程推导出了问题的Hamilton正则方程,为求出任意边界条件下问题的理论解奠定了基础,并且通过算例验证了文中所采用方法的正确性.     

一种求解大型结构动力响应的快速计算方法

针对大型工程结构动力响应求解效率较低的问题,提出了一种基于减基法的快速求解方法。该方法对动力学方程进行时间域积分构建减基空间,利用Galerkin映射向减基空间进行投影得到减缩方程,利用减缩方程快速求得原系统的逼近解,从而大大提高了动力学方程的求解效率。该方法还考虑了Galerkin映射下的奇异性计算,文中算例验证了该...     

多刚体系统动力学的正则方程及数值解法

本文研究了多刚体系统动力学的正则方程,并由此论证了多刚体系统动力学方程的运动稳定性态.在比较不同数值积分方法的基础上,得出结论:基于正则方程的数值积分方法能较好地保证计算精度.     

有限元—传递矩阵法求解链式结构非线性系统的动力响应

本文提出一分析多自由度非线性振动系统时域响应的数值计算法,并就一三自由度非线性系统,分别计算了此系统具有一个非线性弹簧、两个非线性弹簧和三个非线性弹簧三种情况下的动态响应;通过与对应的线性系统响应结果的比较,验证了此法的有效性、可靠性和多种优点.     

约束多体系统动力学正则方程的约束变尺度方法

对约束多体系统动力学的正则方程建立了最优化问题模型,采用辛差分格式,提出相应的最优化问题的约束变尺度法,该方法有效地提高了计算稳定性,算例说明了方法的有效性。     

On the canonical equations of Kirchhoff-Love theory of shells

The paper outlines a procedure to derive the canonical system of equations of the classical theory of thin shells using Reissner’s variational principle and partial variational principles. The Hamiltonian form of the Reissner functional is obtained using Lagrange multipliers to include the kinematical conditions that follow from the Kirchhoff-Love hypotheses. It is shown that the canonical system of equations can be represented in three different forms: one conventional form (five equilibrium equations) and two forms that are equivalent to it. This can be proved by reducing them to the same system of three equations. For problems with separable active and passive variables, partial variational principles are formulated __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 10, pp. 99–107, October 2007.     

Mixed variational principle in elasticity theory and canonical systems of equations

The paper proposes a modification of the mixed variational principle from which stationarity conditions are derived in the form of a mixed system of equations resolved for the first derivatives of the displacement and stress components acting in a plane perpendicular to one of the coordinate axes. The variational principle allows decreasing the dimension of the problem of elasticity thus reducing the system of equations to a canonical form. The modified mixed principle helps immediately obtain a canonical system of equations for various applied theories. This possibility is demonstrated with the example of the Timoshenko theory of plates __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 5, pp. 55–62, May 2007.     

The new solution for the axially symmetrical laminar boundary layer equations between two parallel sherical surfaces

The transformations, which are similar to Mangler’s transformations, are given in this paper. They change the entrance region flow of axially symmetrical laminar boundary layer between two parallel spherical surfaces into the flow of two-dimensional boundary layer, and simplify the problems. The simplified equations can be solved by the two-dimensional boundary layer theory and numerical methods. Therefore, a new way is opened up to solve the diffusive laminar flow in the entrance region between two parallel spherical surfaces.     

Application of canonical coordinates for solving single-freedom constraint mechanical systems

This paper introduces the canonical coordinates method to obtain the first integral of a single-degree freedom constraint mechanical system that contains conserva-tive and non-conservative constraint homonomic systems. The definition and properties of canonical coordinates are introduced. The relation between Lie point symmetries and the canonical coordinates of the constraint mechanical system are expressed. By this re-lation, the canonical coordinates can be obtained. Properties of the canonical coordinates and the Lie symmetry theory are used to seek the first integrals of constraint mechanical system. Three examples are used to show applications of the results.     

结构可靠度分析的响应面法及其Matlab实现

对功能函数不能明确表达的问题进行可靠度分析,常采用响应面法。其中以求得验算点为目的迭代的二次多项式序列响应面法应用较为广泛,本文给出了该方法的Mattab源程序。提出了基于Mattab的插值响应面法和BP神经网络响应面法,介绍了其在Mattab环境下的实现方法,并进行了三种方法的对比分析。Mattab语言基本元素是矩阵,提供了各种矩阵的运算和操作,其中包含结构可靠度计算中常用的各种数值计算方法工具箱。采用Mattab语言构造响应面函数,进行结构可靠度计算,可充分发挥其矩阵运算功能及各种工具箱的作用,使编程效率大大提高,且语法简便,易于掌握。Mattab语言在可靠度计算中的应用,会对结构可靠性理论的推广使用起到积极推进作用。     

Eigenoscillations of a fluid in a canonical domain and functional difference equations

In this work we construct and discuss special solutions of a homogeneous problem for the Laplace equation in a domain with cone-shaped boundaries. The problem at hand is interpreted as that describing oscillatory linear wave movement of a fluid under gravity in such a domain. These solutions are found in terms of the Mellin transform and by means of the reduction to some new functional-difference equations solved in an explicit form (by quadrature). The behavior of the solutions at large distances is studied by use of the saddle point technique. The corresponding eigenoscillations of a fluid are then interpreted as generalized eigenfunctions of the continuous spectrum.     

The asymptotic solving equations of thick ring shell and its solution under moment Mo

In this paper, from the fundamental equations of three dimensional elastic mechanics, we have found a sequence of asymptotic solving equations of thick ring shell (or body) applied arbitrary loads by the perturbation method based upon a geometric small parameter a=r_o/R_o, which may be divided into two independent equation groups which are similar to the equation groups for plane strain and torsional problems. Using these equations, we have also found first order and second order approximate solutions of thick ring shell applied moment M_o.     

Differential characteristic set algorithm for the complete symmetry classification of partial differential equations

In this paper, we present a differential polynomial characteristic set algorithm for the complete symmetry classification of partial differential equations （PDEs） with some parameters. It can make the solution to the complete symmetry classification problem for PDEs become direct and systematic. As an illustrative example, the complete potential symmetry classifications of nonlinear and linear wave equations with an arbitrary function parameter are presented. This is a new application of the differential form characteristic set algorithm, i.e., Wu＇s method, in differential equations.     

Linearization criteria for a system of second-order quadratically semi-linear ordinary differential equations

Conditions are derived for the linearizability via invertible maps of a system of n second-order quadratically semi-linear differential equations that have no lower degree lower order terms in them, i.e., for the symmetry Lie algebra of the system to be sl(n + 2, ℝ). These conditions are stated in terms of the coefficients of the equations and hence provide simple invariant criteria for such systems to admit the maximal symmetry algebra. We provide the explicit procedure for the construction of the linearizing transformation. In the simplest case of a system of two second-order quadratically semi-linear equations without the linear terms in the derivatives, we also provide the construction of the linearizing point transformation using complex variables. Examples are given to illustrate our approach for two- and three-dimensional systems.