Fracture analysis of a weak-discontinuous interface in a symmetrical functionally gradient composite strip loaded by anti-plane impact

The anti-plane impact fracture analysis was performed for a weak-discontinuous interface in a symmetrical functionally gradient composite strip. A new bi-parameter exponential function was introduced to simulate the continuous variation of material properties. Using Laplace and Fourier integral transforms, we reduced the problem to a dual integral equation and obtained asymptotic analytical solution of crack-tip stress field. Based on the numerical solution of the second kind of Fredholm integral equation transformed from the dual integral equation, the effects of the two non-homogeneity parameters on DSIF were discussed. It was indicated that the relative stiffness of the interface and the general stiffness of the whole structure are two important factors affecting the impact fracture behavior of the weak-discontinuous interface. The greater the relative stiffness of the interface is, the higher the value of the dynamic stress intensity factor will be. The greater the general stiffness of the whole structure is, the shorter the time for DSIF to arrive at the peak value and then to stabilize to the steady one. If the general stiffness of the whole structure is great enough, there will be an oscillation between the peak and steady values of DSIF.     

变刚度梁挠度曲线的Green函数法求解

在变刚度梁的线弹性问题中,求解梁受静力荷载的挠度曲线常用解法有积分法与单位荷载法.本文从变刚度梁挠度曲线的微分方程出发,给出了变刚度梁挠度曲线的Green函数法解答,并分析了该解法的优点.从推导结果可以看到,本文提出的公式具有统一、精确、简洁、适合电算的特点,在编制杆系结构计算软件中将具有重要应用价值.     

Nonlinear free vibration of nanotube with small scale effects embedded in viscous matrix

In this paper, the nonlinear free vibration of the nanotube with damping effects is studied. Based on the nonlocal elastic theory and Hamilton principle, the governing equation of the nonlinear free vibration for the nanotube is obtained. The Galerkin method is employed to reduce the nonlinear equation with the integral and partial differential characteristics into a nonlinear ordinary differential equation. Then the relation is solved by the multiple scale method and the approximate analytical solution is derived. The nonlinear vibration behaviors are discussed with the effects of damping, elastic matrix stiffness, small scales and initial displacements. From the results, it can be observed that the nonlinear vibration can be reduced by the matrix damping. The elastic matrix stiffness has significant influences on the nonlinear vibration properties. The nonlinear behaviors can be changed by the small scale effects, especially for the structure with large initial displacement.     

带支撑广义Maxwell粘弹性阻尼耗能结构风振响应分析

工程中粘弹性阻尼器的安装通过支撑与结构进行连接,但在安装粘弹性阻尼器的耗能结构随机响应分析中,为了简化计算过程,将支撑的刚度看作无穷大或者忽略支撑刚度的影响。实际上,对支撑刚度的影响加以考虑更能符合工程实际。针对考虑支撑刚度影响的粘弹性阻尼耗能结构风振响应分析过程复杂的问题,提出了一种求解考虑支撑影响的广义Maxwell粘弹性阻尼耗能结构基于Davenport谱风振响应的简明解析解法。在广义Maxwell粘弹性阻尼器微分型本构模型基础上,给出了考虑支撑刚度的粘弹性阻尼器等效本构关系。将粘弹性阻尼器等效本构关系与结构运动方程联立,采用复模态法将其解耦,获得结构风振响应的统一表达式。将耗能结构在Davenport风速谱下的系列响应功率谱密度函数分解为频域响应函数与Davenport功率谱密度函数的乘积形式,基于随机振动理论中谱矩的定义,对响应功率谱密度函数积分后获得无积分项的系列响应谱矩表达式。在算例中通过与虚拟激励法计算结果对比,验证了所提解法的准确性,并分析了支撑刚度在耗能系统中的影响。     

超静定梁?柱的解析解研究

本文采用渐进积分法研究了超静定梁?柱的弯曲问题. 首先建立超静定梁?柱的四阶挠度微分方程, 考虑到边界条件和连续光滑条件, 采用连续分段独立一体化积分法求解得到了挠度的精确解析解. 为了满足工程设计需要, 构造了超静定梁?柱的四阶挠度微分迭代方程, 选取无轴向力作用时超静定梁的挠曲线作为梁的初函数, 将初函数代入梁的四阶挠度微分迭代方程进行积分, 利用边界条件和连续光滑条件确定积分常数, 得到下一次迭代挠度函数, 依次进行迭代积分运算. 计算出了最大挠度、最大转角和最大弯矩等用轴向力放大系数表示的多项式解析函数解. 本文选取了两种边界条件下受分布力作用的超静定梁?柱进行分析, 计算结果表明, 当超静定梁?柱所受的轴向力小于欧拉临界力的1/2时, 迭代六次误差就可以控制在1%以内; 不仅梁?柱最大位移和最大内力的大小随轴向力的增大而增大, 而且其位置也随轴向力的增大而发生迁移. 本文的研究对揭示轴向力对超静定梁?柱变形和内力的影响有重要意义, 为超静定梁?柱的实际设计提供了一定的理论基础.      

空间网格结构动力分析的非线性模态方法

空间网格结构因自由度数多且无简化的力学模型,非线性动力分析通常要耗费大量时间．传统的非线性模态方法用于求解多高层结构的局部非线性问题已获得良好的效果,但对系统非线性问题的应用尚缺少研究．对比分析多高层结构和空间网格结构动力性能差异,指出网格结构动力非线性分析存在的问题．以主振型理论和切线刚度分离法为基础,将非线性模态方法用于几何非线性效应显著的空间网格结构动力分析．通过对运动方程的非线性恢复力进行拆分,形成线性表达形式,然后解耦到主振型所在的广义坐标系,以达到缩减自由度数量的目的．并通过实例验证非线性模态方法的高效性与适用性．     

A generalized model of elastic foundation based on long-range interactions: Integral and fractional model

The common models of elastic foundations are provided by supposing that they are composed by elastic columns with some interactions between them, such as contact forces that yield a differential equation involving gradients of the displacement field. In this paper, a new model of elastic foundation is proposed introducing into the constitutive equation of the foundation body forces depending on the relative vertical displacements and on a distance-decaying function ruling the amount of interactions. Different choices of the distance-decaying function correspond to different kind of interactions and foundation behavior. The use of an exponential distance-decaying function yields an integro-differential model while a fractional power-law decay of the distance-decaying function yields a fractional model of elastic foundation ruled by a fractional differential equation. It is shown that in the case of exponential-decaying function the integral equation represents a model in which all the gradients of the displacement function appear, while the fractional model is an intermediate model between integral and gradient approaches. A fully equivalent discrete point-spring model of long-range interactions that may be used for the numerical solution of both integral and fractional differential equation is also introduced. Some Green’s functions of the proposed model have been included in the paper and several numerical results have been also reported to highlight the effects of long-range forces and the governing parameters of the linear elastic foundation proposed.     

The integral bifurcation method combined with factoring technique for investigating exact solutions and their dynamical properties of a generalized Gardner equation

It is well known that it is difficult to obtain exact solutions of some partial differential equations with highly nonlinear terms or high order terms because these kinds of equations are not integrable in usual conditions. In this paper, by using the integral bifurcation method and factoring technique, we studied a generalized Gardner equation which contains both highly nonlinear terms and high order terms, some exact traveling wave solutions such as non-smooth peakon solutions, smooth periodic solutions and hyperbolic function solutions to the considered equation are obtained. Moreover, we demonstrate the profiles of these exact traveling wave solutions and discuss their dynamic properties through numerical simulations.     

非圆截面杆中非线性扭转波方程的解析求解

研究了非圆截面杆中非线性扭转波动方程的精确求解问题. 利用直接积分与微分变换相结合的方法,得到了该方程的隐式通解. 通过对积分常数和方程系数的不同情形的讨论, 给出了该方程的三角函数、双曲函数、椭圆函数、指数函数以及它们的组合形式的解,分别对应于的非线性扭转波的孤立波、周期波以及冲击波等多种传播形式.     

功能梯度材料涂层半空间的轴对称光滑接触问题

求解了功能梯度材料涂层半空间的轴对称光滑接触问题,其中梯度层剪切模量按照线性变化,利用Hankel积分变换方法求解微分方程,将问题化为具有Cauchy型奇异核的积分方程.数值方法求解表明:功能梯度材料涂层半空间在刚性柱形压头和球形压头作用下,接触表面分布应力,接触半径以及最大压痕受材料梯度效应的影响较大.     

Stress Analysis for an Elastic Semispace with Surface and Graded Layer Coatings under Induced Torsion

In this paper, the axisymmetric torsional problem of a coating structure consisting of a surface coating, a functionally graded layer and a substrate under a rigid cylindrical punch is investigated. The coating and substrate are homogeneous materials with distinct physical properties while the intermediate layer is inhomogeneous with its shear modulus changing exponentially along the thickness direction. The Hankel integral transform technique is employed to reduce the torsional problem to a singular integral equation with a Cauchy kernel. The circumferential shear stress and displacement fields in the coating structure are calculated by solving the integral equation numerically. The results show that the stiffness ratio has significant effect on the distribution of the circumferential stress and displacement at the interface.     

Quasi-Green＇s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green＇s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi- Green＇s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob- lem. The mode shape differential equations of the free vibration problem of a simply- supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa- tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green＇s function method.     

基于胡聿贤谱的带支撑广义Maxwell阻尼隔震结构随机响应分析

工程中通过设置支撑将阻尼器和建筑结构连接, 但常为了简化分析, 将支撑的水平刚度看成无穷大, 即不考虑支撑对耗能结构随机地震响应的影响. 实际上, 考虑有限水平刚度的支撑对耗能结构响应的影响更加符合工程实际, 为考虑支撑影响的广义Maxwell耗能隔震结构在胡聿贤谱地震激励下的响应分析, 提出一种求解随机地震响应的简明解析解法. 将带支撑广义Maxwell阻尼器等效本构关系、隔震结构运动方程以及胡聿贤谱滤波方程联合组成非经典阻尼系统, 运用复模态法对该非经典阻尼系统解耦, 通过不同响应模态获得耗能隔震系统系列响应基于白噪声激励的Duhamel积分表达式; 利用Dirac函数的性质, 将系统系列响应协方差简化为无积分运算的表达式, 根据功率谱密度函数与其协方差函数的Wiener-Khinchin关系, 得到耗能隔震系统系列响应功率谱和地面加速度功率谱, 基于随机振动理论中谱矩的定义, 得到耗能隔震系统系列响应0 ~ 2阶谱矩. 算例通过与虚拟激励法对比分析, 验证所提方法在该耗能隔震系统分析的正确性和高效性, 并讨论了不同支撑刚度对阻尼器减震效果的影响.      

The uniqueness and existence of solution of the characteristic problem on the generalized KdV equation

The generalized KdV equationu ₁+auu_a+μu_a3+eu_a5=0^[1] is a typical integrable equation. It is derived studying the dissemination of magnet sound wave in cold plasma^[2], the isolated wave in transmission line^[3], and the isolated wave in the boundary surface of the divided layer fluid^[4]. For the characteristic problem of the generalized KdV equation, this paper, based on the Riemann function, designs a suitable structure, then changes the characteristic problem to an equivalent integral and differential equation whose corresponding fixed point, the above integral differential equation has a unique regular solution, so the characteristic problem of the generalized KdV equation has a unique solution. The iteration solution derived from the integral differential equation sequence is uniformly convegent in .     

GREEN QUASIFUNCTION METHOD FOR FREE VIBRATION OF SIMPLY-SUPPORTED TRAPEZOIDAL SHALLOW SPHERICAL SHELL ON WINKLER FOUNDATION

The idea of Green quasifunction method is clarified in detail by considering a free vibration problem of simply-supported trapezoidal shallow spherical shell on Winkler foundation.A Green quasifunction is established by using the fundamental solution and boundary equation of the problem.This function satisfies the homogeneous boundary condition of the problem.The mode shape differential equation of the free vibration problem of simply-supported trapezoidal shallow spherical shell on Winkler foundation is reduced to two simultaneous Fredholm integral equations of the second kind by Green formula.There are multiple choices for the normalized boundary equation.Based on a chosen normalized boundary equation, the irregularity of the kernel of integral equations is avoided.Finally, natural frequency is obtained by the condition that there exists a nontrivial solution in the numerically discrete algebraic equations derived from the integral equations.Numerical results show high accuracy of the Green quasifunction method.     

硬夹心矩形夹层板的整体稳定性分析

摘要：本文在Reissner型理论给出的位移模式基础上,修正其软夹心假设,考虑夹心层面内刚度,给出了硬夹心夹层板的几何方程、物理方程,建立了硬夹心夹层板结构在面内纵向载荷作用下的平衡微分方程,并对方程进行了简化,通过理论计算得到了四边简支条件下硬夹心矩形夹层板整体失稳临界载荷的解析解,并分别计算了夹心层材料的弹性模量 、厚度 、泊松比 对硬夹心夹层板临界载荷的影响,结果证明,对于硬夹心夹层结构,夹心层面内刚度对硬夹心夹层板整体失稳临界载荷的影响较大,考虑其面内刚度是必要的。     

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation

The quasi-Green’s function method is used to solve the free vibration problem of clamped thin plates on the Winkler foundation. Quasi-Green’s function is established by the fundamental solution and the boundary equation of the problem. The function satisfies the homogeneous boundary condition of the problem. The mode-shape differential equation of the free vibration problem of clamped thin plates on the Winkler foundation is reduced to the Fredholm integral equation of the second kind by Green’s formula. The irregularity of the kernel of the integral equation is overcome by choosing a suitable form of the normalized boundary equation. The numerical results show the high accuracy of the proposed method.     

求解多层弹性半空间轴对称问题的精确刚度矩阵法

本文首先从弹性力学的基本方程出发，利用Hankel积分变换等数学手段，推导出了单层弹性半空问轴对称问题的刚度矩阵，然后按传统的有限元方法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程和Hankel积分逆变换就可解出静荷载作用下多层弹性半空间轴对称问题的精确解。由于刚度矩阵的元素中不含有正指数项，计算时不会出现溢出的现象，从而克服了传递矩阵法的缺点。由于在推导过程中摒弃了应力函数的选择，使得问题的求解更加理论化和合理化，同时也为进一步研究这类问题如温度场，动力学等方向奠定了理论基础。最后，文中还给出了计算实例来证明推导结果的准确性。     

Quasi-Green’s function method for free vibration of simply-supported trapezoidal shallow spherical shell

The idea of quasi-Green’s function method is clarified by considering a free vibration problem of the simply-supported trapezoidal shallow spherical shell. A quasi-Green’s function is established by using the fundamental solution and boundary equation of the problem. This function satisfies the homogeneous boundary condition of the prob-lem. The mode shape differential equations of the free vibration problem of a simply-supported trapezoidal shallow spherical shell are reduced to two simultaneous Fredholm integral equations of the second kind by the Green formula. There are multiple choices for the normalized boundary equation. Based on a chosen normalized boundary equa-tion, a new normalized boundary equation can be established such that the irregularity of the kernel of integral equations is overcome. Finally, natural frequency is obtained by the condition that there exists a nontrivial solution to the numerically discrete algebraic equations derived from the integral equations. Numerical results show high accuracy of the quasi-Green’s function method.     

Recursive geometric integrators for wave propagation in a functionally graded multilayered elastic medium

The differential equations governing transfer and stiffness matrices and acoustic impedance for a functionally graded generally anisotropic magneto-electro-elastic medium have been obtained. It is shown that the transfer matrix satisfies a linear 1st order matrix differential equation, while the stiffness matrix satisfies a nonlinear Riccati equation. For a thin nonhomogeneous layer, approximate solutions with different levels of accuracy have been formulated in the form of a transfer matrix using a geometrical integration in the form of a Magnus expansion. This integration method preserves qualitative features of the exact solution of the differential equation, in particular energy conservation. The wave propagation solution for a thick layer or a multilayered structure of inhomogeneous layers is obtained recursively from the thin layer solutions. Since the transfer matrix solution becomes computationally unstable with increase of frequency or layer thickness, we reformulate the solution in the form of a stable stiffness-matrix solution which is obtained from the relation of the stiffness matrices to the transfer matrices. Using an efficient recursive algorithm, the stiffness matrices of the thin nonhomogeneous layer are combined to obtain the total stiffness matrix for an arbitrary functionally graded multilayered system. It is shown that the round-off error for the stiffness-matrix recursive algorithm is higher than that for the transfer matrices. To optimize the recursive procedure, a computationally stable hybrid method is proposed which first starts the recursive computation with the transfer matrices and then, as the thickness increases, transits to the stiffness matrix recursive algorithm. Numerical results show this solution to be stable and efficient. As an application example, we calculate the surface wave velocity dispersion for a functionally graded coating on a semispace.