New hybrid Laplace transform/finite element method for three-dimensional transient heat conduction problem

The paper presents results obtained by the implementation of a new hybrid Laplace transform/finite element method developed by the authors. The present method removes the time derivatives from the governing differential equation using the Laplace transform and then solves the associated equation with the finite element method. Previously reported hybrid Laplace transform/finite element methods¹ have been confined to one nodal solution at a time. When applied to many nodes it takes an excessive amount of computer time. By using a similarity transform method on the matrix of the complex number coefficients this restriction is removed and the reported new method provides a more useful tool for the solution of linear transient problems. Test examples are used to show that the basic accuracy is comparable to that obtainable by analytical, finite difference and finite element methods.     

Numerical operational methods for time-dependent linear problems

A general and systematic discussion on the use of the operational method of Laplace transform for numerically solving complex time-dependent linear problems is presented. Application of Laplace transform with respect to time on the governing differential equations as well as the boundary and initial conditions of the problem reduces it to one independent of time, which is solved in the transform domain by any convenient numerical technique, such as the finite element method, the finite difference method or the boundary integral equation method. Finally, the time domain solution is obtained by a numerical inversion of the transformed solution. Eight existing methods of numerical inversion of the Laplace transform are systematically discussed with respect to their use, range of applicability, accuracy and computational efficiency on the basis of some framework vibration problems. Other applications of the Laplace transform method in conjunction with the finite element method or the boundary integral equation method in the areas of earthquake dynamic response of frameworks, thermaliy induced beam vibrations, forced vibrations of cylindrical shells, dynamic stress concentrations around holes in plates and viscoelastic stress analysis are also briefly described to demonstrate the generality and advantages of the method against other known methods.     

A numerical treatment to the solution of quasiparabolic partial differential equations

This paper presents results obtained by the implementation of a hybrid Laplace transform finite element method to the solution of quasiparabolic problem. The present method removes the time derivatives from the quasiparabolic partial differential equation using the Laplace transform and then solves the associated equation with the finite element method. The numerical inverse of the Laplace transform is realized by solving linear overdetermined systems and a polynomial equation of the kth order. Test examples are used to show that the numerical solution is comparable to the exact solution of the initial-boundary value problem at the given grid points.     

Hybrid method for transient response of circular pins

The combined application of the Laplace transform and the finite element method is used to analyse the transient response of circular pins. The present method removes the time-dependence terms from the governing differential equations and boundary conditions using the Laplace transform and then the eigenfunction expansion method is applied to reduce the two-dimensional boundary value problem to that of one dimension. Accordingly, the final transformed equation can easily be solved by the finite element method. The transformed temperature is inverted to the physical quantity numerically. The present results agree well with analytical solutions. In addition, it is seen that the results of axisymmetric transient heat conduction problems with the central node at r = 0 can accurately be obtained using the present method.     

On the Laplace-transformed-based local meshless method for fractional-order diffusion equation

In this work, we formulate a local meshless method based on Laplace transform to estimate the solution of a time-fractional diffusion equation. The collocation is constructed over small subdomains and combined with Laplace transform for a temporal variable. In this approach, the differentiation matrices are constructed by solving small systems over small local domains instead of a large global collocation matrix. The application of Laplace transform avoids the classical time-stepping procedure. This method is capable of solving fractional differential equations in multidimensions with higher accuracy.     

Numerical modeling of micro- and macro-behavior of viscoelastic porous materials

This paper presents a numerical method for solving the two-dimensional problem of a polygonal linear viscoelastic domain containing an arbitrary number of non-overlapping circular holes of arbitrary sizes. The solution of the problem is based on the use of the correspondence principle. The governing equation for the problem in the Laplace domain is a complex hypersingular boundary integral equation written in terms of the unknown transformed displacements on the boundaries of the holes and the exterior boundaries of the finite body. No specific physical model is involved in the governing equation, which means that the method is capable of handling a variety of viscoelastic models. A truncated complex Fourier series with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the boundaries of the holes. A truncated complex series of Chebyshev polynomials with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the straight boundaries of the finite body. A system of linear algebraic equations is formed using the overspecification method. The viscoelastic stresses and displacements are calculated through the viscoelastic analogs of the Kolosov–Muskhelishvili potentials, and an analytical inverse Laplace transform is used to provide the time domain solution. Using the concept of representative volume, the effective viscoelastic properties of an equivalent homogeneous material are then found directly from the corresponding constitutive equations for the average field values. Several examples are given to demonstrate the accuracy of the method. The results for the stresses and displacements are compared with the numerical solutions obtained by commercial finite element software (ANSYS). The results for the effective properties are compared with those obtained with the self-consistent and Mori–Tanaka schemes.     

New implementation of MLBIE method for heat conduction analysis in functionally graded materials

The meshless local boundary integral equation (MLBIE) method with an efficient technique to deal with the time variable are presented in this article to analyze the transient heat conduction in continuously nonhomogeneous functionally graded materials (FGMs). In space, the method is based on the local boundary integral equations and the moving least squares (MLS) approximation of the temperature and heat flux. In time, again the MLS approximates the equivalent Volterra integral equation derived from the heat conduction problem. It means that, the MLS is used for approximation in both time and space domains, and we avoid using the finite difference discretization or Laplace transform methods to overcome the time variable. Finally the method leads to a single generalized Sylvester equation rather than some (many) linear systems of equations. The method is computationally attractive, which is shown in couple of numerical examples for a finite strip and a hollow cylinder with an exponential spatial variation of material parameters.     

纤维增强塑料齿轮轮齿弯曲强度分析

本文由纤维增强塑料齿轮轮齿的细观结构建立了各向异性粘弹性力学的计算模型,采用Laplace变换和有限单元法相结合求解轮齿的粘弹性应力场。通过引入时间折算因子,改进了Schapery的数值反演过程,给出了一个可行的复合材料齿轮弯曲强度实用分析方法。      

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.     

The boundary element method applied to the analysis of two-dimensional nonlinear sloshing problems

This paper deals with an application of the boundary element method to the analysis of nonlinear sloshing problems, namely nonlinear oscillations of a liquid in a container subjected to forced oscillations. First, the problem is formulated mathematically as a nonlinear initial-boundary value problem by the use of a governing differential equation and boundary conditions, assuming the fluid to be inviscid and incompressible and the flow to be irrotational. Next, the governing equation (Laplace equation) and boundary conditions, except the dynamic boundary condition on the free surface, are transformed into an integral equation by employing the Galerkin method. Two dynamic boundary condition is reduced to a weighted residual equation by employing the Galerkin method. Two equations thus obtained are discretized by the use of the finite element method spacewise and the finite difference method timewise. Collocation method is employed for the discretization of the integral equation. Due to the nonlinearity of the problem, the incremental method is used for the numerical analysis. Numerical results obtained by the present boundary element method are compared with those obtained by the conventional finite element method and also with existing analytical solutions of the nonlinear theory. Good agreements are obtained, and this indicates the availability of the boundary element method as a numerical technique for nonlinear free surface fluid problems.     

Boundary integral equation method for linear porous-elasticity with applications to soil consolidation

For physical phenomena governed by the Biot model of porous-elasticity, a reciprocal relation, similar to the Betti's recoprocal theorem in elasticity, is constructed in Laplace transformed space. Integrating the reciprocal relation enables one to formulate boundary integral equations. The fundamental kernels for the integral equations are solved in closed forms for the case of isotropic material. Numerical implementation of two-dimensional problems includes finite element ideas of discretization and polynomial interpolation, and numerical inversion of a Laplace transform. Practical applications of the method are found in consolidation problems in soils which contain compressible as well as incompressible pore fluids. Also, as a numerical experiment, consolidation of partially saturated soil is simulated and interesting phenomena are observed. The currently developed boundary integral equation method (BIEM) for porous-elasticity may be viewed as an efficient and accurate alternative of existing finite element and finite difference methods. For linear consolidation problems, application of BIEM is always preferred to the other numerical methods whenever possible.     

An electro-magneto-thermo-visco-elastic problem in an infinite medium with a cylindrical hole

A magneto-thermo-visco-elastic problem of an infinite isotropic medium is considered in this problem. The analysis encompasses Lord–Shulman and Green–Lindsay theories of generalized thermo-elasticity to account for the finite velocity of heat equation. The problem is solved by taking Laplace transform techniques and the inversion of the Laplace transform is carried out using a numerical approach.     

Green's functions for the stress intensity factor evolution in finite cracks in orthotropic materials

The elastodynamic response of an infinite orthotropic material with finite crack under concentrated loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some example materials are obtained. This solution can be used as a Green's function to solve dynamic problems involving fini te cracks.     

任意激励作用下变厚度Reissner板的动力响应

本文提出了一个分析中厚度矩形板动力响应的一般方法。文中借助Laplace 变换和传递矩阵法,并引进节线的相关矩阵,较简便地求得了响应在象空间中的一般解,然后再利用拉氏逆变换的数值解求出板的瞬时位移场和内力场。列举的实例和以往的结果作了比较。     

Time-dependent fundamental solutions for homogeneous diffusion problems

This paper describes the applications of the method of fundamental solutions (MFS) for 1-, 2- and 3-D diffusion equations. The time-dependent fundamental solutions for diffusion equations are used directly to obtain the solution as a linear combination of the fundamental solution of the diffusion operator. The proposed scheme is free from the conventionally used Laplace transform or the finite difference scheme to deal with the time derivative of the governing equation. By properly placing the field points and the source points at a given time level, the solution is advanced in time until steady state solutions are reached. Test results obtained for 1-, 2- and 3-D diffusion problems show good comparisons with the analytical solutions and some with the MFS based on the modified Helmholtz fundamental solutions, thus the demonstration present numerical scheme of MFS with the space–time unification has been demonstrated as a promising mesh-free numerical tool to solve homogeneous diffusion problem.     

On the comparison of the methods used for the solutions of the governing equation for unsteady unidirectional flows of second grade fluids

In order to solve the governing equation for unsteady unidirectional flows of second grade fluids, the use of the Laplace and the Fourier transform methods are discussed. Three characteristic examples which are unsteady flow between two parallel plates, unsteady pipe flow and unsteady flow over a plane wall are considered. It is found that the solution for unsteady flow in bounded regions obtained by the Laplace and the Fourier transform methods are exactly the same as the case of the unsteady flow of a Newtonian fluid. It is shown that the Laplace transform method for small values of time is useful for flows of Newtonian fluids but it is not convenient for flows in unbounded regions of second grade fluids. Furthermore, it is explained that for some unsteady flows of second grade fluids, the solution obtained by using the Laplace transform does not satisfy the initial condition and therefore the Fourier transform method is used.     

Finite integration method for nonlocal elastic bar under static and dynamic loads

The finite integration method is proposed in this paper to approximate solutions of partial differential equations. The coefficient matrix of this finite integration method is derived and its superior accuracy and efficiency is demonstrated by making comparison with the classical finite difference method. For illustration, the finite integration method is applied to solve a nonlocal elastic straight bar under different loading conditions both for static and dynamic cases in which Laplace transform technique is adopted for the dynamic problems. Several illustrative examples indicate that high accurate numerical solutions are obtained with no extra computational efforts. The method is readily extendable to solve more complicated problems of nonlocal elasticity.     

Magnetothermoelastic transient response of a multilayered hollow cylinder subjected to magnetic and vapor fields

Summary This paper deals with one-dimensional axisymmetric quasi-static coupled magnetothermoelastic problems subjected to magnetic and vapor fields. Laplace transform and finite difference methods are used to analyze the problems. Using the Laplace transform with respect to time, the general solutions of the governing equations are obtained in the transform domain. The solution is obtained by using the matrix similarity transformation and inverse Laplace transform. We obtain solutions for the temperature and thermal deformation distributions for a transient and steady state. It is demonstrated that the computational procedures established in this paper are capable of solving the generalized magnetothermoelasticity problem of a hollow cylinder with nonhomogeneous layers.     

Dynamic Stress Intensity Factor of Interfacial Finite Cracks in Orthotropic Materials

The elastodynamic response of an infinite non-homogeneous orthotropic material with an interfacial finite crack under distributed normal and shear impact loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some materials are obtained. Interfacial cracks between two different materials and between two pieces of the same material but different fiber orientation are considered. Bimaterial formulation of a crack problem is shown to converge to the mono-material formulation, derived independently, in the limiting case when both materials are the same.     

Analysis for Microscopic Hyperbolic Two-Step Heat Transfer Problems

This work analyzes the heat transfer problems in thin metal films using the microscopic hyperbolic two-step model. It is necessary in dealing with such problems to solve a set of the coupled energy equations or an equation containing higher-order mixed derivatives in both time and space. This present numerical scheme eliminates the coupling between energy equations with the Laplace transform technique and leads to a second-order governing differential equation in the transform domain. Afterward, the transformed second-order governing differential equation is discretized by the control volume scheme. To demonstrate the efficiency and accuracy of the present numerical scheme, a comparison between the present numerical results and the analytical solution is made. Theoretical insight into the hyperbolic two-step heat conduction is provided. Results show that the thermal propagation velocity is finite and is independent of the coupling factor and the volumetric heat capacity ratio between electrons and the lattice.