Analytical solution for bending problem of moderately thick composite annular sector plates with general boundary conditions and loadings using multi-term extended Kantorovich method

An analytical solution for bending of composite sector plates is presented using multi-term extended Kantorovich method (MTEKM). The governing equations are derived using the displacement field of the first-order shear deformation theory and converted into two sets of coupled ordinary differential equations (ODEs) utilizing MTEKM. Next, an analytical iterative procedure is presented for solving the derived sets of ODEs based on state-space method. Numerous examples are studied by the present method, and as special cases, solid sector and rectangular plates are also investigated. Next, the results obtained by the present method are compared to those of finite element method and other results available in the literature. It is found that the present method has a high convergence rate as well as good accuracy in all cases.     

Dynamic solution of a multilayered orthotropic piezoelectric hollow cylinder for axisymmetric plane strain problems

The dynamic solution of a multilayered orthotropic piezoelectric infinite hollow cylinder in the state of axisymmetric plane strain is obtained. By the method of superposition, the solution is divided into two parts: one is quasi-static and the other is dynamic. The quasi-static part is derived by the state space method, and the dynamic part is obtained by the separation of variables method coupled with the initial parameter method as well as the orthogonal expansion technique. By using the obtained quasi-static and dynamic parts and the electric boundary conditions as well as the electric continuity conditions, a Volterra integral equation of the second kind with respect to a function of time is derived, which can be solved successfully by means of the interpolation method. The displacements, stresses and electric potentials are finally obtained. The present method is suitable for a multilayered orthotropic piezoelectric infinite hollow cylinder consisting of arbitrary layers and subjected to arbitrary axisymmetric dynamic loads. Numerical results are finally presented and discussed.     

Non-Equilibrium Mixing of Turbulence Scales Using a Continuous Hybrid RANS/LES Approach: Application to the Shearless Mixing Layer

A new subgrid-scale (SGS) model based on partially integrated transport method (PITM) is applied to the case of a turbulent spectral non-equilibrium flow created by the mixing of two turbulence fields of differing scales: the shearless mixing layer. The method can be viewed as a continuous hybrid RANS/LES approach. In this model the SGS length scale is no longer given by the size of the discretization step, but is dynamically estimated using an additional transport equation for the dissipation rate. The results are compared to those corresponding to the classical model of Smagorinsky and to the experimental data of Veeravalli and Warhaft. A method for creating an anisotropic analytical pseudo-random field for inflow conditions is also proposed. This approach based on subgrid-scale transport modelling combined with anisotropic inlet conditions gives better results for the prediction of the shearless mixing layer.     

Dynamic analysis of functionally graded plates using the hybrid Fourier-Laplace transform under thermomechanical loading

In this study, the analytical solution is presented for dynamic response of a simply supported functionally graded rectangular plate subjected to a lateral thermomechanical loading. The first-order and third-order shear deformation theories and the hybrid Fourier-Laplace transform method are used. The material properties of the plate, except Poisson’s ratio, are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is subjected to a heat flux on the bottom surface and convection on the upper surface. A third-order polynomial temperature profile is considered across the plate thickness with unknown constants. The constants are obtained by substituting the profile into the energy equation and applying the Galerkin method. The obtained temperature profile is considered along with the equations of motion. The governing partial differential equations are solved using the finite Fourier transformation method. Using the Laplace transform, the unknown variables are obtained in the Laplace domain. Applying the analytical Laplace inverse method, the solution in the time domain is derived. The computed results for static, free vibration, and dynamic problems are presented for different power law indices for a plate with simply supported boundary conditions. The results are validated with the known data reported in the literature. Furthermore, the results calculated by the analytical Laplace inversion method are compared with those obtained by the numerical Newmark method.     

Transient responses of a multilayered spherically isotropic piezoelectric hollow sphere

The dynamic solution of a multilayered spherically isotropic piezoelectric hollow sphere subjected to radial dynamic loads is obtained. By the method of superposition, the solution is divided into two parts: one is quasi-static and the other is dynamic. The quasi-static part is derived by the state-space method, and the dynamic part is obtained by the method of separation of variables coupled with the initial parameter method as well as the orthogonal expansion technique. By using the quasi-static and dynamic parts, the electric boundary conditions as well as the electric continuity conditions, a Volterra integral equation of the second kind with respect to a function of time is derived, which can be solved successfully by means of the interpolation method. The displacements, stresses and electric potentials are finally obtained. The present method is suitable for a multilayered spherically isotropic piezoelectric hollow sphere consisting of arbitrary layers and subjected to arbitrary spherically symmetric dynamic loads. Finally, numerical results are presented and discussed.     

Wave motion in a thick cylindrical rod undergoing longitudinal impact

The paper presents a new formulation and the comprehensive analytical solution to longitudinal impact of thick elastic rods. In contrast to previously published works, the solution is derived based on the exact three-dimensional theory without using the classic Skalak’s decomposition. This new direct approach makes the analytical solution more transparent and much easier to obtain. The resulting formulas for basic mechanical quantities are derived using the residue theorem and their evaluation is made in such a way that the accuracy of presented results is significantly higher than those previously published. Based on these results, the transient wave phenomena occurring in the rods are discussed in detail. Additionally, the solution in time domain is obtained also by semi-analytical approach making use of numerical inverse Laplace transform. It is shown that the selected FFT based algorithm is accurate and robust enough, such that the analysis of wave motion in spatial and time domain can be done effectively preserving the results precision. Presented solution can be used as a benchmark for verification of numerical and experimental methods applied to elastodynamics problems.     

层合球面各向同性热释电空心球的瞬态响应

运用叠加原理，将层合球面各向同性热释电空心球的球对称动力学问题的解分成准静态和动 态两部分，准静态部分首先运用状态空间法给出了显式表达式，然后运用分离变量法、初参 数法和特征函数展开技术，给出了动态部分的表示式，再结合内外表面上的电学边界条件和 界面上的电学连续条件，导出一个关于时间函数的第二类Volterra积分方程，运用插值法 可成功地给出此积分方程的高精度数值解，最终可求得原问题的位移、应力、电位移以及电 势的响应. 此方法适用任意层数且各层是任意厚度的层合热释电空心球作用随时间以任意形 式变化的球对称温度场. 文中还给出了数值结果.     

Modelling of Advection-Dominated Transport in a Porous Column with a Time-Decaying Flow Field

This paper examines the problem of the advective-dispersive movement of a non-decaying, inert chemical dye solution through the pore space of a fluid saturated porous column. The objective of the paper is to present a complete study of the one-dimensional advective-dispersive transport problem by considering certain analytical solutions, experimental results and their comparisons with specific computational simulations. Dye concentrations obtained by means of an image processing method are used in conjunction with an analytical solution to identify the hydrodynamic dispersion coefficient that governs the advective-dispersive transport problem. The experimental results and identified parameters are also used to assess the computational estimates derived from several stabilized computational schemes available in the literature, for examining advection-dominated transport processes in porous media.     

水平圆管中受压扭作用管桩屈曲后的解析解

首先建立了考虑自重、受水平圆管约束的管柱在压担组合作用下的屈曲方程：此方程为四阶强非线性常微分方程,难以直接精确求解;通过管柱的正弦屈由分析,找到了方程中的小参数;利用小参数摄动法求解了强非线性屈曲方程,得到了管柱螺旋屈曲后的实际屈曲构形解析解,由此可进一步得到管柱发生螺旋屈曲的临界载荷;所得解析结果与数值计算结果及本文结果退化后与文献中的相关结果比较都有良好的一致性．     

An efficient method to solve the strongly coupled nonlinear differential equations of impact dampers

In this paper, ongoing studies to solve nonlinear differential equations are extended by combining the Newmark-beta integration method and the piecewise linearization approach. The discussed method is illustrated with a practical example. In doing so, the coupled nonlinear differential equations of an impact oscillator, which incorporates the Hertzian contact, are derived. To investigate this problem, an object-oriented computer code, based on the presented method, is written in MATLAB. Furthermore, the discussed problem is solved numerically using the Runge–Kutta commercial code. To verify the calculated results, the contact durations, which are obtained using the discussed methods, are compared with the previous analytical results. In this study, accuracy of solution and the process time (cost) are selected as two main parameters of the solution method. The so-called adequacy factor is presented to combine the two main parameters of solution. Finally, it is shown that in the case of Hertzian contact, the presented method can be more adequate than the Runge–Kutta method.     

Analytical solution to the unsteady one-dimensional conduction problem with two time-varying boundary conditions: Duhamel’s theorem and separation of variables

An analytical resolution of the time-dependent one-dimensional heat conduction problem with time-dependent boundary conditions using the method of separation of variables and Duhamel’s theorem is presented. The two boundary conditions used are a time-dependent heat flux at one end and a varying temperature at the other end of the one-dimensional domain. It is put forth because the author found that the prescribed resolution method using separation of variables and Duhamel’s theorem presented in heat conduction textbooks is not directly applicable to problems with more than one time-dependent boundary condition. The analytical method presented in this paper makes use of one of the property of the heat conduction equation: the apparent linearity of the solutions. For that reason, in order to solve a problem with two time-dependent boundary conditions, the author first separates the initial problem into two independent but complementary problems, each with only one time-dependent boundary condition. Doing that, both simpler problems can be solved independently using a prescribed method that is known to work and the final solution can be obtained by joining the two independent solutions from the simpler separated problems. Every step of the resolution method is presented in this paper, along with a numerical validation of the final solution of three test case problems.     

Chaotic Oscillation of Saddle Form Cable-Suspended Roofs under Vertical Excitation Action

The purpose of this paper is to investigate the dynamic behaviour of saddle form cable-suspended roofs under vertical excitation action. The governing equations of this problem are system of nonlinear partial differential and integral equations. We first establish a spectral equation, and then consider a model with one coefficient, i.e., a perturbed Duffing equation. The analytical solution is derived for the Duffing equation. Successive approximation solutions can be obtained in likely way for each time to only one new unknown function of time. Numerical results are given for our analytical solution. By using the Melnikov method, it is shown that the spectral system has chaotic solutions and subharmonic solutions under determined parametric conditions.     

An analytic solution to the electric analog simulation of the regenerative heat exchanger with time-varying fluid inlet temperatures

Direct simulation on an electric analog of heat transfer problems, has been extensively used, for obtaining approximate solutions. The electric simulation is based on discrete space continuous time differential equations, using electric capacitances and resistances with special devices for simulating convective terms and boundary conditions. (See ref. [16] in the text). In this paper, an analytic solution of the discretized space continuous time differential equations, describing the performance of the periodic-flow regenerative heat exchanger, is presented. With this method the cyclic periodic operation and the starting-up periods of the regenerator, having arbitrary fluid inlet and initial matrix temperatures, can be determined. Several cases with constant inlet fluid temperatures are given and the results are in a good agreement with those previously reported in the literature; the computations for the case of a starting-up period having linearly varying fluid inlet temperature are also presented, which compares well with the results given in a previous paper. Using the present method Hausen's simplified theory is tested and it is found to be justifiable.     

Thermal stresses in infinite circular cylinder subjected to rotation

The present investigation is concerned with the effect of rotation on an infinite circular cylinder subjected to certain boundary conditions.An analytical procedure for evaluation of thermal stresses,displacements,and temperature in rotating cylinder subjected to thermal load along the radius is presented.The dynamic thermal stresses in an infinite elastic cylinder of radius a due to a constant temperature applied to a variable portion of the curved surface while the rest of surface is maintained at zero temperature are discussed.Such situation can arise due to melting of insulating material deposited on the surface cylinder.A solution and numerical results are obtained for the stress components,displacement components,and temperature.The results obtained from the present semi-analytical method are in good agreement with those obtained by using the previously developed methods.     

Transient heat conduction in functionally graded thick hollow cylinders by analytical method

In this article, transient heat conduction in a cylindrical shell of functionally graded material is studied by using analytical method. The shell is assumed to be in axisymmetry conditions. The material properties are considered to be nonlinear with a power law distribution through the thickness. The temperature distribution is derived analytically by using the Bessel functions. To verify the proposed method the obtained numerical results are compared with the published results. The comparisons of temperature distribution between various time and material properties are presented.     

Elastodynamic solution for plane-strain response of functionally graded thick hollow cylinders by analytical method

An elastodynamic solution for plane-strain response of functionally graded thick hollow cylinders subjected to uniformly-distributed dynamic pressures at boundary surfaces is presented. The material properties, except Poisson’s ratio, are assumed to vary through the thickness according to a power law function. To achieve an exact solution, the dynamic radial displacement is divided into two quasi-static and dynamic parts, and for each part, an analytical solution is derived. The quasi-static solution is obtained by means of Euler’s equation, and the dynamic solution is derived using the method of the separation of variables and the orthogonal expansion technique. The radial displacement and stress distributions are plotted for various functionally graded material (FGM) hollow cylinders under different dynamic loads, and the advantages of the presented method are discussed. The proposed analytical solution is suitable for analyzing various arrangements of hollow FGM cylinders with arbitrary thickness and arbitrary initial conditions, which are subjected to arbitrary forms of dynamic pressures distributed uniformly on their boundary surfaces.     

A study on flow resistance in the entrance region of isosceles triangular ducts

In the present paper the variational solution of velocity profile for anincompressible laminar and fully developed flow in isosceles triangular ducts is derivedby applying the Kantovorich method.The theoretical and experimental results ofpressure loss are also given.The velocity distribution model,additional pressure losscoefficient and calculating method of inlet length in the entrance region of isoscelestriangular ducts are also derived,which are suitable for various kinds of vertex angles.The calculations and experiments are also performed for two models:the isoscelestriangular channels with vertex angles2α=45.1°and2α=60°.Comparisons aremade between the theoretical analysis in this paper and those of the other authors.Itcan be seen that the present analytical result is of high.accuracy and widepracticability,and agrees well with the authors’experiment.     

Elastodynamic solution for plane-strain response of functionally graded thick hollow cylinders by analytical method

An elastodynamic solution for plane-strain response of functionally graded thick hollow cylinders subjected to uniformly-distributed dynamic pressures at boundary surfaces is presented. The material properties, except Poisson’s ratio, are assumed to vary through the thickness according to a power law function. To achieve an exact solution, the dynamic radial displacement is divided into two quasi-static and dynamic parts, and for each part, an analytical solution is derived. The quasi-static solution is obtained by means of Euler’s equation, and the dynamic solution is derived using the method of the separation of variables and the orthogonal expansion technique. The radial displacement and stress distributions are plotted for various functionally graded material (FGM) hollow cylinders under different dynamic loads, and the advantages of the presented method are discussed. The proposed analytical solution is suitable for analyzing various arrangements of hollow FGM cylinders with arbitrary thickness and arbitrary initial conditions, which are subjected to arbitrary forms of dynamic pressures distributed uniformly on their boundary surfaces.     

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.      

Analytical solution to a mixed boundary value elastic problem of a roller-guided panel of laminated composite

This study presents an analytical solution to elastic field in a roller-guided panel of symmetric cross-ply laminated composite material. The mixed boundary value two-dimensional plane stress elasticity problem is formulated in terms of a single displacement potential function. This reduces the problem to the solution of a single fourth order partial differential equation of equilibrium as the other equilibrium equation is satisfied automatically. The solution is obtained in terms of an infinite Fourier series. To present some numerical results, a panel of glass/epoxy laminated composite is considered and different components of stress and displacement at different sections of the panel are presented graphically. To justify the present analytical solution, it is compared with the finite element solution obtained by using the commercial software ANSYS. It is found that the two solutions agree well with each other. This ensures that the formulation developed in this study based on the displacement potential approach can be used to obtain analytical solution of an elastic field in structural elements of laminated composite under any mode of boundary conditions prescribed in terms of either stress, displacement or any combination of these.