分数阶理论下实心球体受热冲击作用的渐进分析

基于分数阶广义热弹性理论,针对实心球体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征,借助于Laplace 正、反变换技术及柱函数的渐进性质,推导了热冲击作用周期内位移场、温度场和应力场的渐进表达式. 通过计算,得到了不同传热能力下受热冲击作用时热波、热弹性的传播规律以及位移场、温度场及应力场的分布规律. 结果表明:分数阶参数取值的不同,热波、热弹性波的传播以及各物理场的分布均有所不同,分数阶参数可视为延迟时间的影响因子,通过改变延迟效应对热弹性行为的影响来改变热冲击的作用效果.      

Static Model Describing the Inhomogeneous Deformation of Unidirectional Fibrous Composites under Plane Strain

A mathematical model is proposed to describe the static deformation of a unidirectional fiber-reinforced composite under plane-strain conditions in a plane perpendicular to the fiber direction. The model is based on the stochastic static equations for an elastic homogeneous two-component medium with nonzero body forces. Applying the method of conditional moments and double Fourier transform yields a nonlocal model in the form of integro-differential equations. Expanding the integral kernels into series in expansion coefficients yields differential equations whose order depends on the number of terms in the series. The zero-order approximation leads to the theory of effective moduli, the first-order approximation to the refined theory of effective moduli, and the second-order approximation to fourth-order equilibrium equations for mean displacements (mathematical expectations) and formulas for displacements, strains, and stresses averaged over the composite and its components. All the coefficients in theses formulas are expressed in terms of the elastic constants of the components and the geometric parameters of the structure. The model is valid for heavy stress gradients. In a particular case, the static theory of two-component elastic mixtures follows from the model__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 2, pp. 41–51, February 2005.     

Boundary-layer analysis of propagating mode II cracks in porous elastic media

The steady advance of semi-infinite cracks in porous elastic media is examined in the limiting cases of very slow and very rapid propagation. The low-speed limit is the elastic solution representing fully-drained behavior, with no boundary layers present. In the high-speed limit, two boundary layers are identified: one along the crack face and another centered at the crack tip. The zero-order outer solution is the elastic solution representing fully-undrained behavior, with both boundary layers acting as zones of transition to the effectively drained conditions at the crack itself. Closedform expressions are derived for the stresses and pore pressure within the crack-tip boundary layer. The true stress intensity factor is found to be reduced by an amount related to the disparity between the drained and undrained Poisson's ratios of the medium.     

基于不同板理论的嵌于柔性层中薄膜的皱曲分析

论文对于柔性层-薄膜-柔性层三层结构系统,基于经典板理论、一阶剪切变形理论和高阶剪切变形理论,分别推导给出薄膜皱曲的控制方程.对于两个柔性层,则把它们处理成具有有限厚度的平面应变弹性体.针对上下柔性层固支边界或自由边界条件,利用线性扰动方法得到柔性层对薄膜的横向压力差,最终获得确定薄膜具有周期性正弦型皱曲的临界载荷方程...     

SOLUTION AND ANALYSIS OF CIRCULAR TUNNEL FOR THE STRAIN-SOFTENING ROCK MASSES CONSIDERING THE AXIAL IN SITU STRESS AND SEEPAGE FORCE

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明：在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.     

弹性力学问题的一阶有限元解法

求解弹性力学问题的应力时,如果采用常规的位移有限元法,需要先求得单元的节点位移,再经过求导运算得到。为了解决这种求解方式引起的应力精度下降的问题,提出了弹性力学问题的一阶多变量形式,使得应力与位移精度同阶,并推导了弱形式。采用有限元方法,对弹性力学问题给出了一阶解法的二维、三维数值算例,并且将一阶解法的结果与常规位移有限元法的解进行了比较。数值计算结果表明,一阶解法有效提高了应力的精度,并且应力的误差和节点位移的误差具有相同的收敛阶,验证了本文方法的有效性,为提高有限元法的应力精度提供了新的思路。     

Cohesive zone modeling of interfacial stresses in plated beams

The elastic analysis of interfacial stresses in plated beams has been the subject of several investigations. These studies provided both first-order and higher-order solutions for the distributions of interfacial shear and normal stresses close to the plate end in the elastic range. The notable attention devoted to this topic was driven by the need to develop predictive models for plate end debonding mechanisms, as the early models of this type adopted debonding criteria based on interfacial stresses. Currently, approaches based on fracture mechanics are becoming increasingly established. Cohesive zone modeling bridges the gap between the stress- and energy-based approaches. While several cohesive zone analyses of bonded joints subjected to mode-II loading are available, limited studies have been conducted on cohesive zone modeling of interfacial stresses in plated beams. Moreover, the few available studies present complex formulations for which no closed-form solutions can be found. This paper presents an analytical cohesive zone model for the determination of interfacial stresses in plated beams. A first-order analysis is conducted, leading to closed-form solutions for the interfacial shear stresses. The mode-II cohesive law is taken as bilinear, as this simple shape is able to capture the essential properties of the interface. A closed-form expression for the debonding load is proposed, and the comparison between cohesive zone modeling and linear-elastic fracture mechanics predictions is discussed. Analytical predictions are also compared with results of a numerical finite element model where the interface is described with zero-thickness contact elements, using the node-to-segment strategy and incorporating decohesion and contact within a unified framework.     

Non-axisymmetric thermal stress of a functionally graded coated circular inclusion in an infinite matrix

This paper is to study the non-axisymmetric two-dimensional problem of thermal stresses in an infinite matrix with a functionally graded coated circular inclusion based on complex variable method. With using the method of piece-wise homogeneous layers, the general solution for the functionally graded coating having radial arbitrary elastic properties is derived when the matrix is subjected to uniform heat flux at infinity, and then numerical results are presented for several special examples. It is found that the existence of the functionally graded coating can change interfacial thermal stresses, and choosing proper change ways of the radial elastic properties in the coating can obviously reduce the thermal stresses.     

Modified equations of finite-size layered plates made of orthotropic material. Comparison of the results of numerical calculations with analytical solutions

This paper describes the modified bending equations of layered orthotropic plates in the first approximation. The approximation of the solution of the equation of the three-dimensional theory of elasticity by the Legendre polynomial segments is used to obtain differential equations of the elastic layer. For the approximation of equilibrium equations and boundary conditions of three-dimensional theory of elasticity, several approximations of each desired function (stresses and displacements) are used. The stresses at the internal points of the plate are determined from the defining equations for the orthotropic material, averaged with respect to the plate thickness. The construction of the bending equations of layered plates for each layer is carried out with the help of the elastic layer equations and the conjugation conditions on the boundaries between layers, which are conditions for the continuity of normal stresses and displacements. The numerical solution of the problem of bending of the rectangular layered plate obtained with the help of modified equations is compared with an analytical solution. It is determined that the maximum error in determining the stresses does not exceed 3 %.     

EXPERIMENTS ON EVOLUTION CHARACTERISTICS FOR THE MUSHROOM-LIKE VORTEX STRUCTURE GENERATED BY A SUBMERGED LAMINAR ROUND JET

利用溢流恒压装置产生具有稳定出流速度的层流圆管潜射流,结合染色液流态显示方法,在多种射流时间和雷诺数组合下,实验研究了该射流动量在密度均匀黏性流体中的演化机理及其表现特征,定量分析了蘑菇型涡结构的无量纲射流长度L^*、螺旋型涡环半径R^*及其包络外形长度d^*等几何特征参数随无量纲时间t^*的变化规律。系列实验结果表明,蘑菇型涡结构的形成与演化过程可分为3个不同的阶段,分别为启动阶段、发展阶段和衰退阶段。在启动阶段,L^*和d^*随t^*线性变化,而R^*则近似为一个常数。在发展阶段,蘑菇型涡结构的演化具有自相似性,在各种射流时间和雷诺数组合下,L^*,R^*和d^*与t^*1/2均为同一正比关系,而且实验结果与基于斯托克斯（Stokes）近似的理论解结果一致。在衰退阶段,蘑菇型涡结构会发生两类形式的演化,第1类衰退出现在射流结束之后,其间L^*和R^*与t^*1/5相关,而d^*近似为一个常数;第2类衰退出现在射流结束之前,当射流动量达到某个临界值后,蘑菇型涡结构就会发生破碎等现象。     

A HIGH-ORDER RIGID-FLEXIBLE COUPLING MODEL OF A ROTATING FLEXIBLE BEAM UNDER LARGE DEFORMATION

对在平面内做大范围转动的中心刚体-柔性梁系统的刚柔耦合建模理论进行了深入研究,建立了系统的高次耦合动力学模型. 该动力学模型考虑了柔性梁横向弯曲变形和纵向伸长变形,且在纵向位移中计及由于横向变形而引起的纵向缩短项,即非线性耦合变形项,并保留了与非线性耦合项相关的一些高阶项,最终得到了系统的高次刚柔耦合动力学方程. 由此得到的动力学方程不仅能适用于柔性梁的小变形问题,也同样适用于大变形问题,弥补了一次近似耦合模型在处理柔性梁大变形问题上的不足. 通过与绝对节点坐标法以及一次近似耦合模型的对比验证了高次耦合模型的正确性.     

Analysis of nonlinear acoustoelastic effect of surface acoustic waves in laminated structures by transfer matrix method

The propagation of surface acoustic waves in a layered half-space is investigated in this paper, where a thin cubic Ge film is perfectly bonded to an isotropic elastic Si half-space. Application of the transfer matrix and by solving the coupled field equations, solutions to the mechanical displacements are obtained for the film and elastic substrate, respectively. The phase velocity equations for surface acoustic waves are obtained. Effects of the homogeneous initial stresses induced by the mismatch of the film and substrate are discussed in detail. The results are useful for the design of acoustic surface wave devices.     

Analytical Solution for Radial Deformations of Functionally Graded Isotropic and Incompressible Second-Order Elastic Hollow Spheres

We analytically analyze radial expansion/contraction of a hollow sphere composed of a second-order elastic, isotropic, incompressible and inhomogeneous material to delineate differences and similarities between solutions of the first- and the second-order problems. The two elastic moduli are assumed to be either affine or power-law functions of the radial coordinate R in the undeformed reference configuration. For the affine variation of the shear modulus μ, the hoop stress for the linear elastic (or the first-order) problem at the point R=(R _ou R _in (R _ou +R _in )/2)^1/3 is independent of the slope of the μ vs. R line. Here R _in and R _ou equal, respectively, the inner and the outer radius of the sphere in the reference configuration. For μ(R)∝R ⁿ , for the linear problem, the hoop stress is constant in the sphere for n=1. However, no such results are found for the second-order (i.e., materially nonlinear) problem. Whereas for the first-order problem the shear modulus influences only the radial displacement and not the stresses, for the second-order problem the two elastic constants affect both the radial displacement and the stresses. In a very thick homogeneous hollow sphere subjected only to pressure on the outer surface, the hoop stress at a point on the inner surface depends upon values of the two elastic moduli. Thus conclusions drawn from the analysis of the first-order problem do not hold for the second-order problem. Closed form solutions for the displacement and stresses for the first-order and the second-order problems provided herein can be used to verify solutions of the problem obtained by using numerical methods.     

Thermo-mechanical solution of film/substrate systems under local thermal load and application to laser lift-off of GaN/sapphire structures

Film/substrate structures may undergo a localized thermal load, which can induce stresses, deformation and defects. In this paper, we present the solutions of temperature and stresses in a film/substrate structure under a local thermal load on the film surface. Then, the generalized Stoney formula, which connects the curvature of deformation and the stress field is obtained. The present solution takes into account the non-uniformity of the temperature field both in the width and thickness directions of the film. The thermo-mechanical solution is applied to the analysis of the temperature distribution, stresses, and damage of a GaN/sapphire system during the laser lift-off (LLO) process. It is shown that the laser with the Gaussian distribution of energy density causes much smaller tensile stresses at the edge of the heated area in the film than the laser with the uniform distribution of energy density, and thus can avoid damage to the GaN films separated from the substrate.     

Elastic analyses of heterogeneous hollow cylinders

Two different kinds of heterogeneous elastic hollow cylinders are studied in the present paper. One is a multi-layered cylinder with different values in different layers for both elastic modulus and Poisson’s ratio. Another is an elastic hollow cylinder with continuously graded material properties. By introducing two recursive algorithms, the extrusion stresses between two neighbor layers in the multi-layered cylinder submitted to uniform pressures on the inner and outer surfaces can be simply determined. Then the exact solutions of the multi-layered structure can be found based on Lamé’s solution. For the hollow cylinder with continuously graded properties, the displacement method is used. Both Whittaker equation and hyper-geometric equation are derived and successfully solved, and then the exact solutions are found. The results obtained in the present paper are compared with the numerical solutions and good agreements are found. At the end of the present paper, some inherent properties of these two different kinds of heterogeneous elastic hollow cylinders are presented and discussed. The results obtained in the present paper are useful in the design and analysis for composites reinforced by unidirectional fiber layers.     

A NEW STATE SPACE SOLUTION FOR RECTANGULAR THICK LAMINATES WITH CLAMPED EDGES Hu Wenfeng Liu Yihua

通过把固支边上的边界位移函数作为状态变量引入状态方程,得到了含固支边矩形单层与叠层厚板的精确解. 在求解过程中,将非齐次状态方程的求解变为齐次状态方程的求解,省略了求解待定常数的中间过程,使求解过程变得简单. 所得到的解能够严格满足固支边界条件,在同一材料层内不需作分层处理,因而更加精确. 此外,对固支边的应力提出了新的计算方法,能够得到更精确的边界应力. 算例表明,本文解比现有精确解收敛快,与有限元解吻合的更好,尤其是在固支边处体现得更加明显.     

Dynamic thermoelastic effects for half-planes and half-spaces with nearly-planar surfaces

The effects of non-planarity on the dynamic surface temperature changes induced for plane-strain and 3D problems on the nearly-planar surfaces of, respectively, coupled thermoelastic half-planes and half-spaces by surface heat fluxes are treated. The nearly-planar nature of the surfaces allows the problem solutions to be written, following a standard perturbation scheme, as series expansions in a dimensionless surface contour amplitude parameter. The first, or zero-order, terms represent the ideal (planar) surface solutions, while the second, or first-order, terms represent corrections for non-planarity.Because the characteristic thermoelastic time is of O(10^–7)s, large-time asymptotic forms of the exact integral transform solutions can be used. These can be inverted exactly and used in Green's function operations to yield analytic, or integrals of analytic, expressions. Two types of thermal loading for the half-plane and yet a third type of thermal loading for the half-space are considered. Comparison of the zero- and first-order surface temperature changes for each case indicate that non-planarity gives rise for large times to changes in surface regions beyond those predicted by an ideal surface analysis. Moreover, the magnitudes of these changes can be more significant than the ideal surface results.     

Free vibration analysis of elastically supported functionally graded annular plates subjected to thermal environment

Free vibration analysis of functionally graded (FG) thin-to-moderately thick annular plates subjected to thermal environment and supported on two-parameter elastic foundation is investigated. The material properties are assumed to be temperature-dependent and graded in the thickness direction. The equations of motion and the related boundary conditions, which include the effects of initial thermal stresses, are derived using the Hamilton’s principle based on the first order shear deformation theory (FSDT). The initial thermal stresses are obtained by solving the thermoelastic equilibrium equations. Differential quadrature method (DQM) as an efficient and accurate numerical tool is adopted to solve the thermoelastic equilibrium equations and the equations of motion. The formulations are validated by comparing the results in the limit cases with the available solutions in the literature for isotropic and FG circular and annular plates. The effects of the temperature rise, elastic foundation coefficients, the material graded index and different geometrical parameters on the frequency parameters of the FG annular plates are investigated. The new results can be used as benchmark solutions for future researches.     

Thermomagnetic viscoelastic responses in a functionally graded hollow structure

This paper presents an analytical solution for the interaction of electric potentials,electric displacements,elastic deformations,and thermoelasticity,and describes electromagnetoelastic responses and perturbation of the magnetic field vector in hollow structures(cylinder or sphere),subjected to mechanical load and electric potential.The material properties,thermal expansion coefficient and magnetic permeability of the structure are assumed to be graded in the radial direction by a power law distribution.In the present model we consider the solution for the case of a hollow structure made of viscoelastic isotropic material,reinforced by elastic isotropic fibers,this material is considered as structurally anisotropic material.The exact solutions for stresses and perturbations of the magnetic field vector in FGM hollow structures are determined using the infinitesimal theory of magnetothermoelasticity,and then the hollow structure model with viscoelastic material is solved using the correspondence principle and Illyushin’s approximation method.Finally,numerical results are carried out and discussed.     

Analytical elasto-creep model of interfacial thermal stresses and strains in trilayer assemblies

A two-dimensional model has been developed for thermal stresses, elastic strains, creep strains, and creep energy density at the interfaces of short and long trilayer assemblies under both plane stress and plane strain conditions. Both linear (viscous) and non-linear creep constitutive behavior under static and cyclic thermal loading can be modeled for all layers. Interfacial stresses and strains are approximated using a combination of exact elasticity solutions and elementary strength of materials theories. Partial differential equations are linearized through a simple finite difference discretization procedure. The approach is mathematically straightforward and can be extended to include plastic behavior and problems involving external loads and a variety of geometries. The model can provide input data for thermal fatigue life prediction in solder or adhesive joints. For a typical solder joint, it is demonstrated that the predicted cyclic stress–strain hysteresis shows shakedown and a rapid stabilization of the creep energy dissipation per cycle in agreement with the predictions of finite element analysis.