正交各向异性弹性问题的规则化边界元法

本文致力于平面正交各向异性弹性问题的规则化边界元法研究,提出了新的规则化边界元法的理论和方法。对问题的基本解的特性进行了研究,确立基本解的积分恒等式,提出一种基本解的分解技术,在此基础上,结合转化域积分方程为边界积分方程的极限定理,建立了新颖的规则化边界积分方程。和现有方法比,本文不必将问题变换为各向同性的去处理,从而不含反演运算,也有别于Galerkin方法,无需计算重积分,因此所提方法不仅效率高,而且程序设计简单。特别是,所建方程可计算任何边界位移梯度,进而可计算任意边界应力,而不仅限于面力。数值实施时,采用二次单元和椭圆弧精确单元来描述边界几何,使用不连续插值逼近边界函数。数值算例表明,本文算法稳定、效率高,所取得的边界量数值结果与精确解相当接近。     

New numerical method for volterra integral equation of the second kind in piezoelastic dynamic problems

The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly simplifies the solving procedure for such elastodynamic problems. Meanwhile, it becomes very important to find a way to solve the second kind Volterra integral equation effectively and quickly. By using an interpolation function to approximate the unknown function, two new recursive formulae were derived, based on which numerical solution can be obtained step by step. The present method can provide accurate numerical results efficiently. It is also very stable for long time calculating.     

Solution of the laminar duct flow problem by a boundary integral equation method

Summary A boundary integral equation method is proposed for approximate numerical and exact analytical solutions to fully developed incompressible laminar flow in straight ducts of multiply or simply connected cross-section. It is based on a direct reduction of the problem to the solution of a singular integral equation for the vorticity field in the cross section of the duct. For the numerical solution of the singular integral equation, a simple discretization of it along the cross-section boundary is used. It leads to satisfactory rapid convergency and to accurate results. The concept of hydrodynamic moment of inertia is introduced in order to easily calculate the flow rate, the main velocity, and the fRe-factor. As an example, the exact analytical and, comparatively, the approximate numerical solution of the problem of a circular pipe with two circular rods are presented. In the literature, this is the first non-trivial exact analytical solution of the problem for triply connected cross section domains. The solution to the Saint-Venant torsion problem, as a special case of the laminar duct-flow problem, is herein entirely incorporated.     

Solution of the problem of a flow of rarefied gas of constant density around a semiinfinite plate by the integral diffusion method

Several papers [1–4] have proposed approximate diffusion models which can be used to examine the transport process in a rarefied gas where the mean free path is large and transport is not determined by the local gradient of the particular quantity.In this paper the integral diffusion model [2] is used to solve the problem of determination of the friction stress and velocity of a flow of an incompressible gas around a plane semi-infinite plate in the whole range of Knudsen numbers. The obtained solution is compared with published solutions and experimental data [9].     

A domain decomposition method for modelling Stokes flow in porous materials

An algorithm is presented for solving the Stokes equation in large disordered two‐dimensional porous domains. In this work, it is applied to random packings of discs, but the geometry can be essentially arbitrary. The approach includes the subdivision of the domain and a subsequent application of boundary integral equations to the subdomains. This gives a block diagonal matrix with sparse off‐block components that arise from shared variables on internal subdomain boundaries. The global problem is solved using a biconjugate gradient routine with preconditioning. Results show that the effectiveness of the preconditioner is strongly affected by the subdomain structure, from which a methodology is proposed for the domain decomposition step. A minimum is observed in the solution time versus subdomain size, which is governed by the time required for preconditioning, the time for vector multiplications in the biconjugate gradient routine, the iterative convergence rate and issues related to memory allocation. The method is demonstrated on various domains including a random 1000‐particle domain. The solution can be used for efficient recovery of point velocities, which is discussed in the context of stochastic modelling of solute transport. Copyright © 2002 John Wiley & Sons, Ltd.     

压电桁架作动器/传感器优化配置算法研究

针对自适应压电桁架结构振动控制,建立了作动器/传感器优化配置数学模型,并提出一种优化配置的新方法。为了减少结构分析次数,该方法将近似概念、对偶法和遗传算法相结合,首先采用多点近似技术建立原问题的序列近似问题,再对近似问题中的作动器/传感器位置离散变量和控制增益连续变量采用遗传算法和对偶方法分别寻优的分层优化策略。为了提高近似问题对原问题的逼近程度,本文提出一种适于离散变量结构优化的分段多点近似函数。算例表明本文方法能够以很少的结构分析次数得到最优解。     

The symmetrical adhesive contact problem for circular elastic cylinders

The rolling contact problem involving circular cylinders is at the heart of numerous industrial processes, and critical to any elastohydrodynamic lubrication analysis is an accurate knowledge of the associated contact pressure for the static dry problem. In a recent article [1] the authors have obtained new horizontal pressure distributions, both exact and approximate for various problems involving the symmetrical contact of circular elastic cylinders. In [1] it is assumed that only the circumferential horizontal displacement is prescribed in the contact region while the vertical circumferential displacement is left arbitrary and is assumed to take on whatever value is predicted by the deformation. The advantage of this assumption is that the problem reduces to a single singular integral equation which by transformations can be simplified to an integral equation involving the standard finite Hilbert transform. Here we consider the more general displacement boundary value problem within the contact region, and to be specific we examine the problem with zero vertical circumferential displacement and prescribed horizontal circumferential displacement. The solution of this problem involves two coupled singular integral equations for the horizontal and vertical pressure distributions. Basic equations and some approximate analytical solutions are obtained for symmetrical contact of circular elastic cylinders by both parallel plates and circular cylinders which are either rigid or elastic. Numerical results for the approximate analytical solutions are given for contact by rigid parallel plates and rigid circular cylinders.     

Solutions of thermoelastic crack problems in bonded dissimilar media or half-plane medium

The two-dimensional thermoelastic crack problem in bonded dissimilar media or in a half-plane medium is considered. The proposed method for solving this problem consists of two parts. In the first part, complex potential functions are derived which are enforced to satisfy the continuity conditions across the interface, while the second part consists of the derivation of singular integral equations by introducing the dislocation functions along the crack border which are solved numerically. For both half-plane and two bonded half-plane problems associated with an insulated crack, the thermal stress intensity factors are computed numerically by using the appropriate interpolation formulae. The results compared with those of the homogeneous case given in the literature show that the method proposed here is effective, simple and general.     

弹性体表面分层压缩失稳问题的数学弹性力学解

采用数学弹性力学的稳定平衡方程并结合富氏积分变换的方法研究了含表面平行裂纹的弹性体在压缩载荷下的表面分层失稳问题。导出了一级显式的精确齐次奇异积分方程组,然后.通过Gauss-Chebyshev积分公式,得到一组齐次代数方程组,从而求出临界压缩载荷。并将结果与经典的材料力学梁板稳定的研究方法所得结果进行了比较,指出经典方法误差太大而不适于求解此问题。最后,利用数学弹性力学解求出的等效弹性支承常数给出一个简单精确的临界压缩载荷计算公式。     

Finding minimum energy configurations for constrained beam buckling problems using the Viterbi algorithm

In this work, we present a novel technique to find approximate minimum energy configurations for thin elastic bodies using an instance of dynamic programming called the Viterbi algorithm. This method can be used to find approximate solutions for large deformation constrained buckling problems as well as problems where the strain energy function is non-convex. The approach does not require any gradient computations and could be considered a direct search method. The key idea is to consider a discretized version of the set of all possible configurations and use a computationally efficient search technique to find the minimum energy configuration. We illustrate the application of this method to a laterally constrained beam buckling problem where the presence of unilateral constraints together with the non-convexity of the energy function poses challenges for conventional schemes. The method can also be used as a means for generating “very good” starting points for other conventional gradient search algorithms. These uses, along with comparisons with a direct application of a gradient search and simulated annealing, are demonstrated in this work.     

Diffraction of antiplane shear waves by two cracks in the presence of a magnetic field

The problem of scattering of horizontal polarized shear waves by the two cracks in a uniform magnetostatic field is considered. The magnetic field is assumed to be parallel to the crack surfaces as well as perpendicular to the crack surfaces. The elastic medium under consideration is a homogeneous, isotropic and infinitely conducting one. The solution of the problem is reduced into a pair of triple integral equations having trigonometrical kernels. Using the finite Hilbert transform technique, solution of the pair of triple integral equations is obtained for the low frequencies. Finally, approximate formulae are derived for the stress intensity factors.     

桁架拓扑优化的多点逼近遗传算法

提出一种基于多点逼近函数和遗传算法的桁架拓扑优化方法。该方法建立了包含连续尺寸和离散拓扑两类变量的优化模型，并通过构造多点逼近函数建立了结构优化问题的第一级序列显式近似，然后采用分层优化方法：在外层对拓扑变量采用遗传算法进行优化；在内层对尺寸变量通过可由对偶法求解的第二级序列近似问题进行优化。几个经典的桁架拓扑优化算例表明该方法能以较少的结构分析次数获得比较理想的概率意义上的最优解。     

Explicit solutions for shearing and radial stresses in curved beams

In this paper the formulae for the shearing and radial stresses in curved beams are derived analytically based on the solution for a Volterra integral equation of the second kind. These formulae satisfy both the equilibrium equations and the static boundary conditions on the surfaces of the beams. As some applications, the resulting solutions are used to calculate the shearing and radial stresses in a cantilevered curved beam subjected to a concentrated force at its free end. The numerical results are compared with other existing approximate solutions as well as the corresponding solutions based on the theory of elasticity. The calculations show a better agreement between the present solution and the one based on the theory of elasticity. The resulting formulae can be applied to more general cases of curved beams with arbitrary shapes of cross-sections.     

A formalism for the solution of problems involving chemical reaction and concentration-dependent diffusion coefficients

By using a previously proposed novel integral transform approach, approximate analytical solutions for a certain class of concentration-dependent diffusion and reaction problems can be obtained. Three example problems of increasing complexity are presented to demonstrate the versatility of this technique. Perturbation techniques are also employed to determine asymptotic behaviour of the solutions in the limit of large or small values of some parameters.     

Calculation of contact temperatures generated by the rotation of a shaft in a bearing

The problem of the distribution of contact stresses resulting from the interaction between a journal and its bearing was considered in [1]. This paper deals with the problem of temperature distribution in the area of contact of a rotating cylindrical shaft and a bearing. The process is assumed to be stabilized.The problem reduces to an integral equation with respect to the contact temperature at the shaft surface.An approximate method is proposed for solving the integral equation which had permitted the derivation of a simple approximate formula for the contact temperature within any range of variation of the parameters of this problem.     

物性值随温度变化热弹性问题的摄动—边界元分析

本文将摄动法和边界元法相结合求解物性值随温度变化的热弹性问题,简述了基本方程和积分方程的建立,导出了有关计算公式。算例表明本文方法简便、有效。     

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题.已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算.本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法.该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解.值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差.辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.     

SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR DYNAMIC ANALYSIS OF FOUNDATION

A semi-analytical and semi-numerical method is proposed for the dynamic analysis of foundations. The Lamb's solution and the approximate formulae were used to establish the relation of the contact force and deflection between the foundation and soil. Therefore, the foundation can be separated from soil and analyzed by FEM as for the static cases. The plate can be treated as that the known forces are acting on the upper surface, and the contact pressure from soil can be represented as the deflection. So that only the plate needs to be divided into elements in the analysis. By this method, a series of vibration problems, including various shapes and rigidities of foundations, different excitation frequencies, were analyzed. Furthermore, it can be used for the embedded foundation. The numerical examples show that this method has simplicity, highly accurate and versatile. It is an effective method for the dynamic analysis of foundations.     

Interaction of a die with a layered elastic foundation

Plane and axisymmetric contact problems for a three-layer elastic half-space are considered. The plane problem is reduced to a singular integral equation of the first kind whose approximate solution is obtained by a modified Multhopp-Kalandiya method of collocation. The axisymmetric problem is reduced to an integral Fredholm equation of the second kind whose approximate solution is obtained by a specially developed method of collocation over the nodes of the Legendre polynomial. An axisymmetric contact problem for an transversely isotropic layer completely adherent to an elastic isotropic half-space is also considered. Examples of calculating the characteristic integral quantities are given. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 165–175, May–June, 2006.     

A least-squares procedure for the solution of transport problems

In this paper a least-squares formulation associated with a conjugate gradient algorithm is proposed for the solution of transport problems. In this procedure the advection–diffusion equation is first discretized in time using an implicit scheme. At each time step the resulting partial differential equation is replaced by an optimal control problem. This minimization problem involves the minimization of a functional defined via a state equation. This functional is chosen in order to force the numerical solution of the advection–diffusion equation to be equal to the hyperbolic advective part of this equation. The effectiveness of the method is shown through a one-dimensional example involving advective and diffusive transport. No oscillation and high accuracy have been obtained for the entire range of Peclet numbers with a Courant number well in excess of unity.