压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

???????????????????????

采用谱单元方法推导带阻尼梁的传递函数矩阵. 采用一维连续梁的解析解作为动力形状函 数. 与有限元法相比,谱单元方法的自由度和单元数目小且计算精度高. 传递函数表示为梁 的几何和物理参数的超越隐函数,不同于用模态参数表示的传统传递函数. 提出了采用遗传 算法的结构物理参数识别方法. 以变截面悬臂梁为数值算例,显示该方法的有效性.     

框架结构屈曲的精确有限元求解

基于屈曲微分控制方程的一般解,构造了Euler梁在轴力作用下的精确形函数,建立了用于框架结构屈曲分析的精确有限单元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,并提出了基于常规特征值计算的迭代算法以确定屈曲载荷及相应失稳模态的精确解. 研究表明, 对于线性稳定性分析而言,常规框架有限单元可视为精确有限单元的一种近似. 若采用精确单元,无需进行网格细分就可以获得精确的屈曲载荷和失稳模态. 数值算例证明了新单元和算法的效率和精度.      

The scattering of SH waves by numerous inhomogeneities in an anisotropic body

The fundamental solution of antiplane strain motion in an anisotropic infinite medium is derived. On the basis of this, an integral equation for the scattered displacement field associated with many inhomogeneities of arbitrary shapes is proposed with the help of an equivalent body force and two-dimensional Heaviside functions. By using the Born approximation theory the scattered far-field is discussed when only two flaws (or voids) are considered. At the same time the differential cross line is defined and deduced. In the end, numerical examples are given.     

均布荷载作用下功能梯度悬臂梁弯曲问题的解析解

采用弹性力学半逆解法,假设所有材料常数沿梁厚度方向按同一函数规律变化,求得了功能梯度悬臂梁在均布载荷作用下的解析解.该解退化到各向同性均匀弹性情况时与已有的理论解相一致.对弹性模量按指数函数梯度变化的算例进行了分析.所得到的解对任意梯度函数均成立,可作为数值解以及简化理论的检验依据.     

FREE VIBRATION ANALYSIS OF SIMPLY SUPPORTED BEAM WITH BREATHING CRACK USING PERTURBATION METHOD

In this paper, a new analytical method for vibration analysis of a cracked simply supported beam is investigated. By considering a nonlinear model for the fatigue crack, the governing equation of motion of the cracked beam is solved using perturbation method. The solution of the governing equation reveals the superhaxmonics of the fundamental frequency due to the nonlinear effects in the dynamic response of the cracked beam. Furthermore, considering such a solution, an explicit expression is also derived for the system damping changes due to the changes in the crack parameters, geometric dimensions and mechanical properties of the cracked beam. The results show that an increase in the crack severity and approaching the crack location to the middle of the beam increase the system damping. In order to validate the results, changes in the fundamental frequency ratios against the fatigue crack severities are compared with those of experimental results available in the literature. Also, a comparison is made between the free response of the cracked beam with a given crack depth and location obtained by the proposed analytical solution and that of the numerical method. The results of the proposed method agree with the experimental and numerical results.     

均布荷载作用下悬臂磁电弹性梁的解析解

对磁电弹性平面问题进行了研究，给出了用拟调和位移函数表达的通解，进而以试凑法按平面应力问题推导出了均布荷载作用下悬臂磁电弹性粱的解析解，所得解有易于理解、便于校对、形式统一简洁的特点。本文还将计算结果与压电材料和弹性材料相应结果进行了分析、比较，为验证各种数值计算方法提供了参考依据。     

Buckling reliability evaluation of a clamped beam with random imperfections subjected to axial impact using probability density evolution method

Based on the generalized probability density evolution theory, the joint probability density function of impact responses and random factors, that is, random initial imperfection and load intensity for a clamped beam, is derived. And a Total Variation Diminishing (TVD)-based numerical scheme is employed to obtain the numerical solutions. With given threshold values, the dynamic reliability is assessed in terms of a bilateral displacement criterion. Two numerical examples are then presented to demonstrate the proposed method. One is a clamped beam with random initial perfections subjected to a deterministic impulse load, and another is a clamped beam with a deterministic initial imperfection subjected to impact loading with a random amplitude. The predicted dynamic reliability is in good agreement with the result by Monte Carlo method.     

基于无网格法的Timoshenko曲梁动力分析

曲梁具有外形美观、受力性能良好的优点,故在工程中得到广泛应用。本文基于移动最小二乘近似和一阶剪切变形理论,提出一种对Timoshenko曲梁自由振动和受迫振动进行分析的无网格方法。通过一系列离散点离散曲梁,建立曲梁无网格模型,然后推导曲梁势能和动能方程,通过哈密顿原理给出曲梁自由振动和受迫振动的控制方程,因为本文方法不能直接施加边界条件,所以使用完全转换法处理本质边界条件,最后求解方程得到频率和振动模态。文末通过算例验证了本文方法的有效性,且通过收敛性分析表明本文方法具有较好的收敛性,并进一步分析了不同边界条件、跨高比和变截面变曲率对曲梁自由振动和受迫振动的影响,将计算结果与文献解或ABAQUS解进行对比分析,表明本文方法具有较高的精度,且适用于实际工程情况。     

Plastic modal approximations in analyzing beam-on-beam collisions

Dynamic behavior of a moving free–free beam striking the tip of a cantilever beam, as a typical example of collision between two deformable structures, is analysed by employing modal approximation techniques. The applicability of both rigid-plastic and elastic–plastic mode approximations is examined in predicting the energy partitioning between the two colliding beams.Three rigid-plastic modes (RP-Modes) are considered and the Lee’s functional is applied to select the appropriate mode. It is found that one of the beams would absorb all the initial kinetic energy, unless a higher-order RP-mode is adopted.To incorporate the effect of elastic deformation into the modal solution, an elastic, perfectly plastic mode (EP-Mode) approximation for the same problem is proposed. By replacing each of the plastic hinges in the RP-Mode with a nature hinge and an elastic–plastic rotational spring, the fundamental features of the dynamic elastic–plastic behavior of the two colliding beams are revealed. Both beams participate in energy dissipation, while the structural and geometrical parameters greatly influence the energy partitioning. It is shown from numerical examples that the EP-Mode solution provides a fairly good approximation compared with the RP complete solution and finite element simulation.     

Static analysis of anisotropic functionally graded magneto-electro-elastic beams subjected to arbitrary loading

This paper considers the analytical and semi-analytical solutions for anisotropic functionally graded magneto-electro-elastic beams subjected to an arbitrary load, which can be expanded in terms of sinusoidal series. For the generalized plane stress problem, the stress function, electric displacement function and magnetic induction function are assumed to consist of two parts, respectively. One is a product of a trigonometric function of the longitudinal coordinate (x) and an undetermined function of the thickness coordinate (z), and the other a linear polynomial of x with unknown coefficients depending on z. The governing equations satisfied by these z-dependent functions are derived. The analytical expressions of stresses, electric displacements, magnetic induction, axial force, bending moment, shear force, average electric displacement, average magnetic induction, displacements, electric potential and magnetic potential are then deduced, with integral constants determinable from the boundary conditions. The analytical solution is derived for beam with material coefficients varying exponentially along the thickness, while the semi-analytical solution is sought by making use of the sub-layer approximation for beam with an arbitrary variation of material parameters along the thickness. The present analysis is applicable to beams with various boundary conditions at the two ends. Two numerical examples are presented for validation of the theory and illustration of the effects of certain parameters.     

基于线性互补模型的梯度塑性连续体无网格方法

对梯度塑性连续体提出了一个归结为线性互补问题的数值分析方法. 塑性乘子与位移均为主要未知变量,并采用基于移动最小二乘的无网格方法分别在积分点与节点上插值. 联立弱形式下的平衡方程与积分点上逐点满足的非局部本构方程和屈服准则可以导出一个线性互补问题,并通过Lexico-Lemke算法求解. 构造了一个基于N-R方法的迭代方案,使得不需要形成一致性切线刚度矩阵而仍保持二阶收敛性. 一维和二维的数值算例证明了所提出的方法处理由应变软化引起的应变局部化问题的有效性.      

基于渐近均匀化方法的准周期梁结构等效刚度数值实现方法研究

相比周期梁结构,准周期梁结构沿轴向梯度变化,具有更大的设计自由度,能够获得更好的结构性能。由于其非均质性,一般将其均匀化为具有等效性质的均质梁结构,但现有工作很少涉及准周期梁结构等效性质的计算。本文针对由周期梁结构映射而成的准周期梁结构,通过引入雅可比矩阵,基于渐近均匀化方法推导的单胞方程及其等效性质计算列式,并建立了其单胞方程及等效刚度的有限元求解列式。该方法可以处理沿轴向变形的任意微单胞构型,数值算例验证了其正确性和有效性。     

Galerkin technique based on beam functions in application to the parametric instability of thermal convection in a vertical slot

A Fourier–Galerkin spectral technique for solving coupled higher‐order initial‐boundary value problems is developed. Conjugated systems arising in thermoconvection that involve both equations of fourth and second spatial orders are considered. The set of so‐called beam functions is used as basis together with the harmonic functions. The necessary formulas for expressing each basis system into series with respect to the other are derived. The convergence rate of the spectral solution series is thoroughly investigated and shown to be fifth‐order algebraic for both linear and nonlinear problems. Though algebraic, the fifth‐order rate of convergence is fully adequate for the generic problems under consideration, which makes the new technique a useful tool in numerical approaches to convective problems. An algorithm is created for the implementation of the method and the results are thoroughly tested and verified on different model examples. The spatial and temporal approximation of the scheme is tested. To further validate the scheme, a singular asymptotic expansion is derived for small values of the modulation frequency and amplitude and the numerical and analytic results are found to be in good agreement. The new technique is applied to the G‐jitter flow, and the Floquet stability diagrams are produced. We obtain the expected alternating isochronous and subharmonic branches and find that stable motions are always isochronous while unstable motions can be either isochronous or subharmonic. The numerical investigation also leads to novel conclusions regarding the dependence of the amplitude of the solutions on some of the governing parameters. Copyright © 2008 John Wiley & Sons, Ltd.     

Geometrically nonlinear,steady state vibration of viscoelastic beams

The problem of geometrically non-linear steady state vibrations of beams excited by harmonic forces is considered in this paper. The beams are made of a viscoelastic material defined by the classic Zener rheological model - the simplest model that takes into account all the basic properties of real viscoelastic materials. The constitutive stress-strain relationship for this type of material is given as a differential equation containing derivatives of both stress and strain. This significantly complicates the solution to the problem. The von Karman theory is applied to describe the effects of geometric nonlinearities of beam deformations. The equations of motions are derived using the finite element methodology. A polynomial approximation of bending moments is used. The order of basis functions is set so as to obtain a coherent approximation of moments and displacements. In the steady-state solution of equations of motion, only one harmonic is taken into account. The matrix equations of amplitudes are derived using the harmonic balance method and the continuation method is applied for solving them. The tangent matrix of equations of amplitudes is determined in an explicit form. The stability of steady-state solution is also examined. The resonance curves for beams supported in a different way are shown and the results of calculation are briefly discussed.     

采用有限单元法计算梁任意形状截面特性

采用将梁截面离散化的方式,用数值积分计算截面的几何特性,并根据梁剪切变形和扭转理论,利用变分原理建立截面的有限元法方程,求解任意形状截面的扭转常数、剪切中心以及剪切面积修正系数等特性.本方法适用于各种形式的截面,具有计算精度高及适应性强的特点.根据上述理论编制了相应程序,按照不同的单元划分方式,分别计算出矩形截面截面特性,与理论解进行比较;又对舟山市定海长峙至岙山预应力混凝土连续箱梁截面进行了计算,并与Ansys结果进行比较,均证明采用本文的计算方法能得到满意的结果,且该方法适用于各种形状的截面形式.     

考虑液面波动和液体压缩时一侧受液作用悬臂梁的横振分析

本文分析一侧部分受有液体作用悬臂梁的横向自由振动，同时考虑了液面波动和液体可压缩性对悬臂梁自振特性的影响。利用一组广义三角级数的正交完备性，求得了悬臂梁与液体耦联振动的振型函数和频率方程的精确解析解，最后给出了几个数值算例。     

无约束修正Timoshenko梁的冲击问题

介绍了修正后的Timoshenko梁运动方程，并比较了修正Timoshenko梁与经 典Timoshenko梁的运动方程. 推导了考虑剪切变形引起的转动惯量的修正Timoshenko 梁的正交条件，推导了集中质量对无约束修正Timoshenko梁的正碰撞对梁所引起的瞬态冲 击响应公式，并用算例进行了分析，且与集中质量对经典的无约束Timoshenko梁的正碰撞 对梁所引起的冲击响应进行了比较，另外还用算例分析了梁的刚度的变化和冲击质量比对其 冲击响应产生的影响.     

内共振条件下轴向运动梁磁弹性超谐共振研究

以载流导线激发的磁场中轴向运动梁为研究对象,同时考虑外激励力作用,推导出梁的磁弹性非线性振动方程．通过位移函数的设定和伽辽金积分法,将非线性振动方程离散为常微分方程组．采用多尺度法得到系统的近似解析解．应用Matlab 和Mathematica 软件求解幅频响应方程,并对稳态解进行稳定性判定．通过具体算例得到前两阶假设模态的响应幅值随不同参数的变化规律．结果发现：系统在内共振条件下发生超谐波共振时,二阶假设模态幅值明显小于一阶;随着外激励的增大,多值解区域范围明显缩小;随着电流强度增加,振动幅值减小,表明载流导线能够起到控制共振的作用．     

AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS

The hybrid boundary node method （HBNM） is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares（MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.