On an inverse problem for inhomogeneous thermoelastic rod

In recent years, different fields of engineering have been increasingly incorporating functionally graded materials with variable physical properties that significantly improve a quality of elements of designs. The efficiency of practical application of thermoelastic inhomogeneous materials depends on knowledge of exact laws of heterogeneity, and to define them it is necessary to solve coefficient inverse problems of thermoelasticity.In the present research a scheme of solving the inverse problem for an inhomogeneous thermoelastic rod is presented. Two statements of the inverse problem are considered: in the Laplace transform space and in the actual space. The direct problem solving is reduced to a system of the Fredholm integral equations of the 2nd kind in the Laplace transform space and an inversion of the solutions obtained on the basis of the theory of residues. The inverse problem solving is reduced to an iterative procedure, at its each step it is necessary to solve the Fredholm integral equation of the 1st kind; to solve it the Tikhonov method is used. Specific examples of a reconstruction of variable characteristics required are given.     

A fast volume integral equation method for the direct/inverse problem in elastic wave scattering phenomena

A fast method for solving the volume integral equation is introduced for the solution of forward and inverse multiple scattering problems in an elastic 3-D full space. For both forward and inverse scattering analysis, the volume integral equation in the wavenumber domain is used. By means of the discrete Fourier transform, the volume integral equation in the wavenumber domain can be dealt with as a Fredholm equation of the 2nd kind with respect to a non-Hermitian operator on a finite dimensional vector space. The Bi-CGSTAB method is employed to construct the Krylov subspace in the wavenumber domain. The current procedure establishes a fast and simplified method without requiring the derivation of a coefficient matrix. Several numerical results validate the accuracy and effectiveness of the current method for both forward and inverse scattering analysis. According to the numerical results, the reconstruction of inhomogeneities of the wave field is successful, even for multiple scattering of several cubes.     

Hydrodynamic methods in the inverse problem of gravitational prospecting

In [1, 2], a dynamical method is proposed for solving stationary inverse problems of potential theory, including the inverse problem of gravitational prospecting. It is based on analogy with the problem of establishing the interface of two immiscible fluids flowing in a porous medium. In the present paper, a system of two functional equations is derived from which one can obtain, as special cases, an equation corresponding to the method of [1, 2], and also a system of equations that enables one to propose a new and different method for solving the inverse problem of gravitational prospecting. Equations are derived in polar coordinates for plane Cauchy problems corresponding to both methods, and the results are also given of the solution of some model problems by these methods. Finally, ways of generating new methods of solution of the inverse problem of gravitational prospecting are considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 63–71, July–August, 1980.     

A generalized inverse procedure for bifurcation problems

We give a systematic procedure for handling the linear aspects of the investigation of a bifurcation problem for the solutions of a boundary-value problem with ordinary differential equations. The method is an alternative method. The required projections are uniquely specified. We also use the related generalized inverse of the linear mapping.

The technique is illustrated with two examples of periodic solutions of second-order nonlinear differential equations.     

A 2.5-D dynamic model for a saturated porous medium. Part II: Boundary element method

The 3-D boundary integral equation is derived in terms of the reciprocal work theorem and used along with the 2.5-D Green’s function developed in Part I [Lu, J.F., Jeng, D.S., Williams, S., submitted for publication. A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function. Int. J. Solids Struct.] to develop the 2.5-D boundary integral equation for a saturated porous medium. The 2.5-D boundary integral equations for the wave scattering problem and the moving load problem are established. The Cauchy type singularity of the 2.5-D boundary integral equation is eliminated through introduction of an auxiliary problem and the treatment of the weakly singular kernel is also addressed. Discretisation of the 2.5-D boundary integral equation is achieved using boundary iso-parametric elements. The discrete wavenumber domain solution is obtained via the 2.5-D boundary element method, and the space domain solution is recovered using the inverse Fourier transform. To validate the new methodology, numerical results of this paper are compared with those obtained using an analytical approach; also, some numerical results and corresponding analysis are presented.     

On inverse coefficient problems of poroelasticity

General approaches to the inverse coefficient problems of poroelasticity on the basis of a modified Biot model are considered. A generalized reciprocity relation is constructed, and an iteration process for determining the unknown coefficients is stated. By way of example, the problem of steady-state longitudinal vibrations of an inhomogeneous poroelastic layered system is considered, and integral equations for the direct and inverse problems are derived. The results of computational experiments where the elastic modulus and the Biot modulus were reconstructed for various laws of variation are given.     

An inverse problem formulation of the immersed-boundary method

We formulate the immersed-boundary method (IBM) as an inverse problem. A control variable is introduced on the boundary of a larger domain that encompasses the target domain. The optimal control is the one that minimizes the mismatch between the state and the desired boundary value along the immersed target-domain boundary. We begin by investigating a naïve problem formulation that we show is ill-posed: in the case of the Laplace equation, we prove that the solution is unique, but it fails to depend continuously on the data; for the linear advection equation, even solution uniqueness fails to hold. These issues are addressed by two complimentary strategies. The first strategy is to ensure that the enclosing domain tends to the true domain, as the mesh is refined. The second strategy is to include a specialized parameter-free regularization that is based on penalizing the difference between the control and the state on the boundary. The proposed inverse IBM is applied to the diffusion, advection, and advection-diffusion equations using a high-order discontinuous Galerkin discretization. The numerical experiments demonstrate that the regularized scheme achieves optimal rates of convergence and that the reduced Hessian of the optimization problem has a bounded condition number, as the mesh is refined.     

Solution of multiple crack problem in a finite plate using coupled integral equations

This paper studies a numerical solution of multiple crack problem in a finite plate using coupled integral equations. After using the principle of superposition, the multiple crack problem in a finite plate can be converted into two problems: (a) the multiple crack problem in an infinite plate and (b) a usual boundary value problem for the finite plate. For the former problem, the Fredholm integral equation is used. For the latter problem, a BIE based on complex variable is suggested in which a Cauchy singular kernel exists. For the proposed BIE, after using the inverse matrix technique, the dependence of the traction at a domain point from the boundary tractions is formulated indirectly. This is a particular advantage of the present study. Several numerical examples are provided and the computed results for stress intensity factor and T-stress at crack tips are given.     

Theory of non-linear wave propagation with application to the interaction and inverse problems

The propagation of non-linear deformation waves in a dissipativc medium is described by a unified asymptotic theory, making use of wave front kinematics and the concepts of progressive waves. The mathematical models are derived from the theories of thermoclasticity or viscoclasticity taking into account the geometric and physical non-linearities and dispersion. On the basis of eikonal equations for the associated linear problem the transport equations of the nth order are obtained. In the multidimensional case the method of matched separation of initial equations is proposed. The interaction problems which occur in head-on collisions and in reflection from boundaries or interfaces are analyzed. Conditions are also studied when the interaction of non-linear waves does not take place. The inverse problem of determining materials properties according to pulse shape changes is discussed.     

Large amplitude oscillations of thick hyperelastic cylindrical shells

We present numerical solutions to the problem of large amplitude oscillations of a thick-walled hyperelastic cylindrical shell employing the general theory of finite dynamic deformations of elastic bodies. The material of the shell is considered incompressible and of Mooney-Rivlin type rubbers.

We apply a fourth-order Runge-Kutta numerical technique to the governing equation which was originally derived by J.K. Knowles in 1960.

We consider the free as well as forced oscillations due to a Heaviside step load and display graphs for the variations of amplitude against time and frequencies for different thicknesses and material constants. Discussions are presented on the significances of the results obtained.     

Direct and inverse problems of oscillations of an elastoliquid layered medium

A plane problem of forced oscillations of an ideal compressible liquid bounded from above by an elastic layer with a rough lower surface and an inverse geometric problem of determining the shape of the rough lower surface of an elastic layer from the wave characteristics on the upper surface are considered. Three methods are used to solve the direct problem: the small parameter method, the boundary element method, and the Born approximation. Solving the inverse problem is reduced to solving the integral Fredholm equation of the first kind. Results of a numerical experiment are presented.     

General solutions to a class of time fractional partial differential equations

A class of time fractional partial differential equations is considered, which in- cludes a time fractional diffusion equation, a time fractional reaction-diffusion equation, a time fractional advection-diffusion equation, and their corresponding integer-order partial differential equations. The fundamental solutions to the Cauchy problem in a whole-space domain and the signaling problem in a half-space domain are obtained by using Fourier- Laplace transforms and their inverse transforms. The appropriate structures of the Green functions are provided. On the other hand, the solutions in the form of a series to the initial and boundary value problems in a bounded-space domain are derived by the sine- Laplace or cosine-Laplace transforms. Two examples are presented to show applications of the present technique.     

Biquaternion solution of the kinematic control problem for the motion of a rigid body and its application to the solution of inverse problems of robot-manipulator kinematics

The problem of reducing the body-attached coordinate system to the reference (programmed) coordinate system moving relative to the fixed coordinate system with a given instantaneous velocity screw along a given trajectory is considered in the kinematic statement. The biquaternion kinematic equations of motion of a rigid body in normalized and unnormalized finite displacement biquaternions are used as the mathematical model of motion, and the dual orthogonal projections of the instantaneous velocity screw of the body motion onto the body coordinate axes are used as the control. Various types of correction (stabilization), which are biquaternion analogs of position and integral corrections, are proposed. It is shown that the linear (obtained without linearization) and stationary biquaternion error equations that are invariant under any chosen programmed motion of the reference coordinate system can be obtained for the proposed types of correction and the use of unnormalized finite displacement biquaternions and four-dimensional dual controls allows one to construct globally regular control laws. The general solution of the error equation is constructed, and conditions for asymptotic stability of the programmed motion are obtained. The constructed theory of kinematic control of motion is used to solve inverse problems of robot-manipulator kinematics. The control problem under study is a generalization of the kinematic problem [1, 2] of reducing the body-attached coordinate system to the reference coordinate system rotating at a given (programmed) absolute angular velocity, and the presentedmethod for solving inverse problems of robotmanipulator kinematics is a development of the method proposed in [3–5].     

Performance evaluation of an inverse integral equation method applied to turbomachine cascades

An improved formulation of the inverse integral equation method proposed in Reference 1 is presented which allows, in particular, a well-posed problem to be ensured. The corresponding computation code is tested in an exhaustive manner for axial and radial compressor and turbine cascades. The agreement between the velocity field obtained with the inverse method and that resulting from a direct calculation is examined for subsonic, transonic and supersonic flows. Accuracy and reliability of the solution to the boundary condition problem are excellent for the subsonic and transonic flows. However, for the supersonic flow, the application of the method seems to be limited by the use of elementary solutions of the Laplace operator.     

Reflection and transmission of elastic waves by the spatially periodic interface between two solids (Theory of the integral-equation method)

The linear theory of two-dimensional reflection and transmission of time-harmonic, elastic waves by the spatially periodic interface between two perfectly elastic media is developed. A given phase progression of the incident wave in the direction of periodicity induces a modal structure in the elastodynamic field and leads to the introduction of the so-called spectral orders. The main tools in the analysis are the elastodynamic Green-type integral relations. They follow from the two-dimensional form of the elastodynamic field reciprocity theorem, where in the latter a Green state adjusted to the periodicity of the structure at hand is used. One of these relations is a vectorial integral equation from which the elastodynamic field quantities can be determined.

The consequences of field reciprocity in the structure and of the conservation of energy are developed in view of their serving as a check on numercal results to be obtained from the relevant integral equations.

The formalism thus developed applies to profiles, if periodic, of arbitrary shape and size and can without too serious difficulties be implemented on a computer. The major difficulty in this respect is the relevant Green function, the series representation of it being slowly convergent. Its evaluation becomes tractable after an appropriate technique for accelerating the convergence. The only practical limitations are then put by the speed of the computer and its storage capacity.     

An inverse for the Jaumann derivative and some applications to the rheology of viscoelastic fluids

Summary By using a generalization of the matrizant of matrix calculus, it is shown how one can construct formally an inverse, or integral, for the well-knownJaumann derivative of continuum mechanics. Some applications to fluid rheology are then considered. First, it is shown that this integral provides, via theBoltzmann super-position principle, a generalization of Oldroyd's quasi-linear fluid model, which is related to the molecular model ofBueche. Explicit expressions for the stresses arising in a general laminar shear flow are then derived for this model. Secondly, it is indicated how the operation can be used with rheological equations which are nonlinear in the deformation-rate, but quasi-linear in stress, to solve explicitly for the stress in terms of kinematic quantities. As an example, a rheological equation for suspensions of viscoelastic spheres in aNewtonian fluid is treated.     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

Inverse problem for the viscoelastic medium with discontinuous wave impedance

Ｉｎｖｅｒｓｉｏｎｏｆｔｈｅｍａｔｅｒｉａｌｆｕｎｃｔｉｏｎｓｗｉｔｈｔｈｅｍｅａｓｕｒｅｍｅｎｔｄａｔａｏｂｔａｉｎｅｄｆｒｏｍｔｈｅｓｃａｔｔｅｒｉｎｇｅｘｐｅｒｉｍｅｎｔｉｓｏｆｃｅｎｔｒａｌｉｍｐｏｒｔａｎｃｅｉｎｓｏｍｅｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｆｏｒｅｘａｍｐｌｅ,ｓｃａｔｔｅｒｅｄｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｗａｖｅｓｃａｎｂｅｕｓｅｄｔｏｒｅｃｏｎｓｔｒｕｃｔｔｈｅｐｅｒｍｉｔｔｉｖｉｔｙ,ｃｏｎｄｕｃｔｉｖｉｔｙａｎｄｓｕｓｃｅｐｔｉｂｉｌｉｔｙｋｅｒｎｅｌｏｆ…     

硬脑膜受压膨出挠度的力学分析

硬脑膜是一种粘弹性材料,为控制硬脑膜在脑压作用下的膨出度,对粘弹性薄膜受压膨出挠度作力学分析。以位移为未知量,从粘弹性材料的分型本构关系出发将Foepple薄膜大挠度理论从弹性推广到粘弹性膜,得到一组非线性积分偏微分方程。先在空间上运用Galerkin方法将积分偏微分方程组化为积分常微分方程组。然后,在时间域上运用数值积分和有限差分将方程离散为非线性代数方程组。本文对四周固定夹紧的圆形、椭圆形和矩形薄膜进行了求解,并将求解结果用于颅底缺损重建膜的膨出量计算,计算值与实验值吻合,为颅底外科提供一个理论分析方法。     

GLOBAL SOLUTION OF THE INVERSE PROBLEM FOR ACLASS OF NONLINEAR EVOLUTION EQUATIONSOF DISPERSIVE TYPE

ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｗｅｌｌ＿ｋｎｏｗｎｔｈａｔｐｓｅｕｄｏ＿ｐａｒａｂｏｌｉｃｅｑｕａｔｉｏｎｗｉｔｈｐｒｉｎｃｉｐａｌｐａｒｔｕｔ -ｕｘｘｔｈａｓｂｅｅｎｓｔｕｄｙｉｎｇｒｅｃｅｎｔｌｙ ,ｂｅｃａｕｓｅｔｈｅｒｅｅｘｉｓｔｓｗｉｄｅｐｈｙｓｉｃａｌｂａｃｋｇｒｏｕｎｄｆｏｒｔｈｉｓｃｌａｓｓｏｆｅｑｕａｔｉｏｎｓ.[1 ]ｓｔｕｄｉｅｄｔｈｅｍｕｌｔｉ＿ｄｉｍｅｎｓｉｏｎｉｎｖｅｒｓｅｐｒｏｂｌｅｍｆｏｒｔｈｅｆｏｌｌｏｗｉｎｇｃｌａｓｓｏｆｎｏｎｌｉｎｅａｒｅｖｏｌｕｔｉｏｎｅｑｕ…