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.     

弹性薄板分析的条形传递函数方法

提出一种用于矩形弹性薄板变形分析的条形传递函数方法．一个矩形区域首先沿某一个方向被剖分成若干个条形子域,分割这些子域的直线称为结线,在结线上定义位移函数,它是结线坐标的一维函数,结线的两个端点称为结点．为适应复杂边界条件,在边界结线上定义若干结点,该结线的位移函数用结点位移参数插值表示．每个条形子域的变形用结线位移函数和适当的插值函数（形函数）表示．结线位移函数和结点位移参数满足的平衡微分方程及代数方程由变分原理给出     

Method of constructing estimates of the solution of equations of filtration of an aerated liquid

The exact solutions of the nonlinear equations of filtration of an aerated liquid have been obtained in [1–3]. In [4] the system of equations of an aerated liquid have been reduced to the heat-conduction equation under certain assumptions. An approximate method of computing the nonsteady flow of an aerated liquid is given in [5], where the real flow pattern is replaced by a computational scheme of successive change of stationary states. In [6] the same problem is solved by the method of averaging. In the present article estimates of the solution of the equations for nonstationary filtration of an aerated liquid in one-dimensional layer are constructed under certain conditions imposed on the desired functions. These estimates can be used as approximate solutions with known error or for the verification of the accuracy of different approximate methods. We note that the use of comparison theorem for the estimate of solutions of equations of nonlinear filtration is discussed in [7–9]. The methods of constructing estimates of solutions of various problems of heat conduction are given in [10, 11]     

On solutions of the problem in stresses with the use of Maxwell stress functions

The spatial problems of elasticity are mainly solved in displacements [1, 2], i.e., the Lamé equations are taken as the initial equations. This is related to the lack of general solutions for the system of basic equations of elasticity expressed in stresses. In this connection, a new variational statement of the problem in stresses was developed in [3, 4]; this statement consists in solving six generalized equations of compatibility for six independent components of the stress tensor, while the three equilibrium equations are transferred to the set of boundary conditions. This method is more convenient for the numerical solution of problems in stresses and has been tested when solving various boundary value problems. In the present paper, analyzing the completeness of the Saint-Venant identities and using the Maxwell stress functions, we obtain a new resolving system of three differential equations of strain compatibility for the three desired stress functions φ, ξ, and ψ. This system is an alternative to the three Lamé equilibrium equations for three desired displacement components u, v, w and is simpler in structure. Moreover, both of these systems of resolving equations can be solved by the new recursive-operator method [5, 6]. In contrast to well-known methods for constructing general solutions of linear differential equations and their systems, the solutions obtained by the recursive-operator method are constructed as operator-power series acting on arbitrary analytic functions of real variables (not necessarily harmonic), and the series coefficients are determined from recursive relations (matrix in the case of systems of equations). The arbitrary functions contained in the general solution can be determined directly either from the boundary conditions (the obtained system of inhomogeneous equations with a right-hand side can also be solved by the recursive-operator method [6]) or by choosing them from various classes of analytic functions (elementary, special); a complete set of particular solutions can be obtained in the same function classes, and the coefficients of linear combinations of particular solutions can be determined by the Trefftz method, the least-squares method, and the collocation method.     

NONLINEAR DYNAMICAL BIFURCATION AND CHAOTIC MOTION OF SHALLOW CONICAL LATTICE SHELL

The nonlinear dynamical equations of axle symmetry are established by the method of quasi-shells for three-dimensional shallow conical single-layer lattice shells. The compatible equations are given in geometrical nonlinear range. A nonlinear differential equation containing the second and the third order nonlinear items is derived under the boundary conditions of fixed and clamped edges by the method of Galerkin. The problem of bifurcation is discussed by solving the Floquet exponent. In order to study chaotic motion, the equations of free oscillation of a kind of nonlinear dynamics system are solved. Then an exact solution to nonlinear free oscillation of the shallow conical single-layer lattice shell is found as well. The critical conditions of chaotic motion are obtained by solving Melnikov functions, some phase planes are drawn by using digital simulation proving the existence of chaotic motion.     

An exact analytic method applied to nonhomogeneous ring- and stringer-stiffened cylindrical shells

In this paper, the nonlinear axial symmetric deformation problem of nonhomogeneous ring- and stringer-stiffened shells is first solved by the exact analytic method. An analytic expression of displacements and stress resultants is obtained and its convergence is proved. Displacements and stress resultants converge to exact solution uniformly. Finally, it is only necessary to solve a system of linear algebraic equations with two unknowns. Four numerical examples are given at the end of the paper which indicate that satisfactory results can be obtained by the exact analytic method.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

THE GENERAL SOLUTION FOR DYNAMIC RESPONSE OF NONHOMOGENEOUS BEAM WITH VARIABLECROSSSECTION

ＴＨＥＧＥＮＥＲＡＬＳＯＬＵＴＩＯＮＦＯＲＤＹＮＡＭＩＣＲＥＳＰＯＮＳＥＯＦＮＯＮＨＯＭＯＧＥＮＥＯＵＳＢＥＡＭＷＩＴＨＶＡＲＩＡＢＬＥＣＲＯＳＳＳＥＣＴＩＯＮＪｉＺｈｅｎ－ｙｉ（纪振义）（ＡｎｈｕｉＡｒｃｈｉｔｅｃｔｕｒａｌＩｎｄｕｓｔｒｙＣｏｌｌ...     

On quadruple integral equations related to a certain crack problem

An exact solution of a four part mixed boundary value problem representing a three colinear crack system connected with specified crack opening displacements between the cracks is obtained. The three cracks thus become one with pressure and/or opening displacement prescribed on the crack face. From considerations of dual symmetry and a formulation based on Papkovich-Neuber harmonic functions, the boundary value problem is reduced to solving a quadruple set of integral equations. An exact solution of these equations is derived using a modified finite Hilbert transform technique. The closed form results for the stress distributions and the crack-tip stress intensity factors are presented. Limiting cases of the solution yield results which agree with well known solutions.     

用离散速度方法计算浅水长波方程

用离散速度法计算浅水波方程，将空气动力学方程和浅水波方程作了比较，用Nadiga提出的近平衡流动方法模拟浅水波方程的连续和间断解。计算了一维的溃坝波问题和Thacker提出的连续解问题，结果与精确解作了比较，并且计算了水流跃过障碍物的问题。     

扁球面网壳的混沌运动研究

在圆形三向网架非线性动力学基本方程的基础上,用拟壳法给出了圆底扁球面三向网壳的非线性动力学基本方程．在固定边界条件下,引入了异于等厚度壳的无量纲量,对基本方程和边界条件进行无量纲化,通过Galerkin作用得到了一个含二次、三次的非线性动力学方程．为求Melnikov函数,对一类非线性动力系统的自由振动方程进行了求解,得到了此类问题的准确解．在无激励情况下,讨论了稳定性问题．在外激励情况下,通过求Melnikov函数,给出了可能发生混沌运动的条件．通过数字仿真绘出了平面相图,证实了混沌运动的存在．     

Passage of a supersonic flow over intersecting surfaces

Using the linear formulation, the problem of passage of a supersonic flow over slightly curved intersecting surfaces whose tangent planes form small dihedral angles with the incident flow velocity at every point is considered. Conditions on the surfaces are referred to planes parallel to the incident flow forming angles 0≤γ≤2π at their intersection [1]. The problem reduces to finding the solution of the wave equation for the velocity potential with boundary conditions set on the surfaces flowed over and the leading characteristic surface. The Volterra method is used to find the solution [2]. This method has been applied to the problem of flow over a nonplanar wing [3] and flow around intersecting nonplanar wings forming an angle γ=π/n (n=1, 2, 3, ...) with consideration of the end effect on the wings forming the angle [4]. In [5] the end effect was considered for nonplanar wings with dihedral angle γ=m/nπ. In the general case of an arbitrary angle 0≤γ≤2π the problem of finding the velocity potential reduces to solution of Volterra type integrodifferential equations whose integrands contain singularities [1]. It was shown in [6] that the integrodifferential equations may be solved by the method of successive approximation, and approximate solutions were found differing slightly from the exact solution over the entire range of interaction with the surface and coinciding with the exact solution on the characteristic lines (the boundary of the interaction region, the edge of the dihedral angle). The solution of the problem of flow over intersecting plane wings (the conic case) for an arbitrary angle γ was obtained in terms of elementary functions in [7], which also considered the effect of boundary conditions set on a portion of the leading wave diffraction. In [8, 9] the nonstationary problem of wave diffraction at a plane angle π≤γ≤2π was considered. On the basis of the wave equation solution found in [8], this present study will derive a solution which permits solving the problem of supersonic flow over nonplanar wings forming an arbitrary angle π≤γ≤2π in quadratures. The solutions for flow over nonplanar intersecting surfaces for the cases 0≤γ≤π [6] and π≤γ≤2π, found in the present study, permit calculation of gasdynamic parameters near a wing with a prismatic appendage (fuselage or air intake). The study presents a method for construction of solutions in various zones of wing-air intake interaction.     

Plate cascade in subsonic unsteady gas flow

Accounting for fluid compressibility creates serious difficulties in solving the problem of oscillations of a grid of thin, slightly curved profiles in a subsonic stream. The problem has been solved in [1–3] for a widely-spaced cascade without stagger whose profiles oscillate in phase opposition. The phenomenon of aerodynamic (acoustic) resonance, which may arise in a grid in the direction transverse to the stream for definite values of the stream velocity and profile oscillation frequency, was discovered in [2]. An approximate solution of the problem in which account is not taken of the effect of the vortex trails on the gas flow has been obtained in [4]. In [5, 6] Meister studied in the exact linear formulation the problem of unsteady gas motion through an unstaggered cascade of semi-infinite plates. In [7] Meister considered a grid of profiles with finite chords, but the problem solution was not completed. The problem of subsonic gas flow through a staggered lattice whose profiles oscillate following a single law with constant phase shift was solved most completely in the studies of Kurzin [8, 9] using the method of integral equations. A method of solving the problem for the case of arbitrary harmonic oscillation laws for the lattice profiles was indicated in [10]. The results of the calculation of the unsteady aerodynamic forces for the particular case of a plate cascade without stagger are presented in [9,11], and the possibility of the occurrence of aerodynamic resonance in the cascade in the directions transverse to and along the stream is indicated.Another method of solving the problem is given in [12], in which the more general problem of unsteady subsonic gas flow through a three-dimensional cascade of plates is solved. In the present study this method is applied to the solution of the problem of oscillations of staggered plate cascades in a two-dimensional subsonic gas flow. The results are presented of an electronic computer calculation of the unsteady aerodynamic characteristics of the cascade profiles, which show the essential influence of fluid compressibility on these characteristics. In particular, a sharp decrease of the aerodynamic damping in the acoustic resonance regimes is obtained.     

A Mode- II Griffith Crack in Decagonal Quasicrystals

By using the method of stress functions, the problem of mode-II Griffith crack in decagonal quasicrystals was solved. First, the crack problem of two-dimensional quasicrystals was decomposed into a plane strain state problem superposed on anti-plane state problem and secondly, by introducing stress functions, the18 basic elasticity equations on coupling phonon-phason field of decagonal quasicrystals were reduced to a single higher-order partial differential equations. The solution of this equation under mixed boundary conditions of mode-II Griffith crack was obtained in terms of Fourier transform and dual integral equations methods. All components of stresses and displacements can be expressed by elemental functions and the stress intensity factor and the strain energy release rate were determined. Biography: GUO Yu-cui (1962-), Associate professor, Doctor     

Incompatible curved quadrilateral plate bending element with 12-degrees of freedom

This paper presents a new curved quadrilateral plate element with12-degree freedom by the exact element method.The method can be used to arbitrary non-positive and positive definite partial differential equations without variation principle.Using this method,the compatibility conditions between element can be treated very easily,if displacements and stress resultants are continuous at nodes between elements.The displacements and stress resultants obtained by the present method can converge to exact solution and have the second order convergence speed.Numerical examples are given at the end of this paper,which show the excellent precision and efficiency of the new element.     

Determining linear vibration frequencies of a ferromagnetic shell

The problems of determining the roots of dispersion equations for free bending vibrations of thin magnetoelastic plates and shells are of both theoretical and practical interest, in particular, in studying vibrations of metallic structures used in controlled thermonuclear reactors. These problems were solved on the basis of the Kirchhoff hypothesis in [1–5]. In [6], an exact spatial approach to determining the vibration frequencies of thin plates was suggested, and it was shown that it completely agrees with the solution obtained according to the Kirchhoff hypothesis. In [7–9], this exact approach was used to solve the problem on vibrations of thin magnetoelastic plates, and it was shown by cumbersome calculations that the solutions obtained according to the exact theory and the Kirchhoff hypothesis differ substantially except in a single case. In [10], the equations of the dynamic theory of elasticity in the axisymmetric problem are given. In [11], the equations for the vibration frequencies of thin ferromagnetic plates with arbitrary conductivity were obtained in the exact statement. In [12], the Kirchhoff hypothesis was used to obtain dispersion relations for a magnetoelastic thin shell. In [5, 13–16], the relations for the Maxwell tensor and the ponderomotive force for magnetics were presented. In [17], the dispersion relations for thin ferromagnetic plates in the transverse field in the spatial statement were studied analytically and numerically.In the present paper, on the basis of the exact approach, we study free bending vibrations of a thin ferromagnetic cylindrical shell. We obtain the exact dispersion equation in the form of a sixth-order determinant, which can be solved numerically in the case of a magnetoelastic thin shell. The numerical results are presented in tables and compared with the results obtained by the Kirchhoff hypothesis. We show a large number of differences in the results, even for the least frequency.     

Determination of an Exact Solution of the Integral Gelfand-Levitan-Marchenko Equation for the Sturm-Liouville Operators with Step-Type Potential

We propose a method for solving the integral Gelfand-Levitan-Marchenko (GLM) equation for the Sturm-Liouville operator with step-type potential. The scattering function is found in explicit form. An associated system of infinite recurrence equations is solved. The integral operator kernel is presented in exact form using Bessel functions. A series of new integral representations for Bessel functions is obtained for the first time.     

On singular perturbation boundary-value problem of coupling type system of convection-diffusion equations

In this paper we consider the singular perturbation boundary-value problem of thefollowing coupling type system of convection-diffusion equationsWe advance two methods:the first one is the initial value solving method,by which theoriginal boundary-value problem is changed into a series of unperturbed initial-valueproblems of the first order ordinary differential equation or system so that an asymptoticexpansion is obtained;the second one is the boundary-value solving method,by which theoriginal problem is changed into a few boundary-value problems having no phenomenon ofboundary-layer so that the exact solution can be obtained and any classical numericalmethods can be used to obtain the numerical solution of consismethods can be used to obtainthe numerical solution of consistant high accuracy with respect to the perturbationparameterε     

Investigation of a griffith crack subject to uniform tension using the non-local theory by a new method

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｓｅｖｅｒａｌｐｒｅｖｉｏｕｓｐａｐｅｒｓ［１,２,３］,Ｅｒｉｎｇｅｎｄｉｓｃｕｓｓｅｄｔｈｅｓｔａｔｅｏｆｓｔｒｅｓｓｎｅａｒｔｈｅｔｉｐｏｆａｓｈａｒｐｌｉｎｅｃｒａｃｋｉｎａｎｅｌａｓｔｉｃｐｌａｔｅｓｕｂｊｅｃｔｔｏｕｎｉｆｏｒｍｔｅｎｓｉｏｎ,ｓｈｅａｒａｎｄａｎｔｉ＿ｐｌａｎｅｓｈｅａｒ．Ｔｈｅｆｉｅｌｄｅｑｕａｔｉｏｎｓｅｍｐｌｏｙｅｄｉｎｔｈｅｓｏｌｕｔｉｏｎｏｆｔｈｅｓｅｐｒｏｂｌｅｍｓａｒｅｔｈｏｓｅｏｆｔｈｅｔｈｅｏｒｙｏｆｎｏｎ＿ｌｏｃａｌｅ…