含铰接杆系结构几何非线性分析子结构方法

将细长杆系结构按长度方向划分为多个子结构,由于在子结构坐标系下的节点位移均是小位移,可以将子结构内部自由度凝聚到边界. 考虑到子结构端面在变形过程中保持为刚性截面,将端面节点自由度进一步凝聚到端面形心点,这样每一个子结构就减缩成形式上只有两个节点的广义梁单元,大大减缩了自由度. 大位移大转动是细长杆系结构产生几何非线性效应的一个重要原因,基于共旋坐标法,建立了随单元一起运动的随动坐标系,推导了子结构单元的节点力平衡方程及其切线刚度阵. 同时,考虑到工程机械中细长杆系结构含有相互铰接的刚体加强块,给出了非独立自由度节点力转换到独立参数下的广义节点力及其导数. 最后,通过履带式起重机的副臂工况算例,给出了其在不同载荷下的臂架结构位移,验证了方法的正确性.     

一种考虑剪切变形的平行四边形厚/薄板通用单元

根据Timoshenko二广义位移梁理论，构造了深梁位移场的插值函数。利用斜坐标系与直角坐标系的变换关系、有限条带思想和深梁位移插值函数，构造了一种考虑剪切变形的平行四边形厚／薄板弯曲通用单元的位移(曲率、剪应变、转角、横向位移)插值函数，导出了刚度矩阵和非结点荷载等效力。并对简支阍支方板、Razzaque斜板、四边简支斜交板弯曲进行了数值计算。算例表明此单元有较好的精度，对于薄板不出现剪切闭锁，可适应于目前桥梁建设中大量采用的斜交板桥结构分析。     

考虑切向摩擦力时弹性地基梁的振动

假定切向摩擦力与梁底面的纵向位移成正比,通过引入广义剪力,得到了梁的位移型平衡方程。将位移及荷载展开为带附加项的Fourier级数,利用平衡方程和边界条件研究了弹性地基梁的自由振动和简谐振动。通过算例结果分析表明:纵向摩擦力对梁的固有频率、位移和内力均有影响。梁的最大挠度、转角、弯矩及剪力随着地基纵向反力系数的增大而减小;梁的固有频率、轴向位移和轴力则随着地基纵向反力系数的增大而增大;同时轴力引起的轴向位移和转角引起的梁底面纵向位移具有同一数量级。     

回传波矩阵法在杆系结构静力分析中的应用

回传波矩阵法最初是由Pao等人分析二维框架结构动力响应时提出的。对于三维杆系结构的静力分析，为了确定结构的位移和内力，先要建立传递分配矩阵和载荷源向量，这可通过列出所有节点的静力平衡方程和位移协调方程来实现。同时，通过分析每根杆近端位移和远端位移的关系，建立结构的回传波矩阵（重分配矩阵）。在此基础上求解线性方程组，就可以得到结构的位移和内力。本文推导了空间杆系结构的有关矩阵方程式，并给出了一固定梁的两端弯矩求解算例。     

NODAL DISPLACEMENT ACCURACY OF ONE-DIMENSIONAL FINITE ELEMENTS

研究高次杆单元和梁单元的节点位移精度问题.首先求出一端固支均匀杆和悬臂梁在任意次多项式形式分布载荷作用下的位移精确解,然后用二次杆单元、五次欧拉梁单元和三次铁木辛柯梁单元求得了节点位移.通过比较有限元解与精确解以及利用静力凝聚方法,发现一次以上杆单元、三次以上欧拉梁单元以及三次以上铁木辛柯梁单元都可以给出精确的端点位移.     

一维有限单元的节点位移精度

研究高次杆单元和梁单元的节点位移精度问题.首先求出一端固支均匀杆和悬臂梁在任意次多项式形式分布载荷作用下的位移精确解,然后用二次杆单元、五次欧拉梁单元和三次铁木辛柯梁单元求得了节点位移.通过比较有限元解与精确解以及利用静力凝聚方法,发现一次以上杆单元、三次以上欧拉梁单元以及三次以上铁木辛柯梁单元都可以给出精确的端点位移.     

建立基于力的变截面梁单元质量矩阵的新方法

基于位移的有限梁单元中三次Hermite插值函数不能有效地描述变截面梁单元内部位移变化,只能通过加密网格增加单元数解决,会造成计算量增大。基于力的有限梁单元由于使用的力插值函数不受截面形状变化的影响,在处理变截面梁时有很大优势,可以得到精确的位移插值函数,利用较少的单元可以达到很高的精度,解决了基于位移的有限梁单元在处理变截面梁时的不足。本文得到了考虑剪切变形的位移插值函数和考虑转动惯量的一致质量矩阵。利用算例验证了本文理论的正确性和高效性。     

平面刚架(梁)单元杆端力向量公式推导

本文通过位移法的转角位移方程，导出平面刚架（梁）单元杆端力向量公式和单元刚度矩阵．     

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

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

管道流固耦合振动的行波方法研究

应用波动方法分析载流管道流固耦合振动问题。采用直梁模型，推导了管道流固耦合系统的轴向、横向波导方程及平面管系的坐标转换矩阵，通过节点位移连续和力平衡条件建立了耦合系统的散射模型。最终分析了两种耦合对载流管道振动的影响机理及流速变化对管道振动稳定性的影响。数值算例表明本文方法的正确性及有效性。     

考虑分布轴向力的细长杆横向振动与失稳分析

为量化梁、杆、柱的自重(下称分布轴向力)对静力失稳和动力横向振动的影响,在《材料力学》和《机械系统动力学》教材的基础上,建立了分布轴向力下的杆柱失稳和横向振动的力学、数学模型。采用有限差分法、伽辽金法和数值积分法获取计算结果。结果表明:考虑分布轴向力的杆柱横向振动固有频率随杆长增加而减小,杆柱失稳时一阶固有频率为0;分布轴向力对较短杆柱的临界载荷影响较小,对较长杆柱影响较大;伽辽金法不适用计算超长杆柱的失稳临界载荷。     

Buckling of short beams with warping effect included

A beam theory for the stability analysis of short beam that includes shear deformation and warping of the cross-section is developed. The warping of the cross-section is taken to be an independent kinematics quantity and corresponding force resultants are defined. For the beam subjected to the external loading only at the ends of the beam, equilibrium equations have been obtained by the principle of virtual work. The variations of lateral displacement, rotational angle of the cross-section and the multiplier of the warping shape along the beam axis are solved in closed form and expressed in terms of deformation quantities at the ends of the beam. Based on this beam theory, the lateral stiffness of the beam sustained an axial compression force and a lateral shear force at one end is explicitly derived, from which the equation of the buckling load is established and the buckling load can be solved. When the effect of cross-section warping is neglected, the derived lateral stiffness and buckling load converge to the solutions of the Haringx theory.     

Surface and non-local effects for non-linear analysis of Timoshenko beams

In this paper, we present a non-local non-linear finite element formulation for the Timoshenko beam theory. The proposed formulation also takes into consideration the surface stress effects. Eringen׳s non-local differential model has been used to rewrite the non-local stress resultants in terms of non-local displacements. Geometric non-linearities are taken into account by using the Green–Lagrange strain tensor. A C⁰ beam element with three degrees of freedom has been developed. Numerical solutions are obtained by performing a non-linear analysis for bending and free vibration cases. Simply supported and clamped boundary conditions have been considered in the numerical examples. A parametric study has been performed to understand the effect of non-local parameter and surface stresses on deflection and vibration characteristics of the beam. The solutions are compared with the analytical solutions available in the literature. It has been shown that non-local effect does not exist in the nano-cantilever beam (Euler–Bernoulli beam) subjected to concentrated load at the end. However, there is a significant effect of non-local parameter on deflections for other load cases such as uniformly distributed load and sinusoidally distributed load (Cheng et al. (2015) [10]). In this work it has been shown that for a cantilever beam with concentrated load at free end, there is definitely a dependency on non-local parameter when Timoshenko beam theory is used. Also the effect of local and non-local boundary conditions has been demonstrated in this example. The example has also been worked out for other loading cases such as uniformly distributed force and sinusoidally varying force. The effect of the local or non-local boundary conditions on the end deflection in all these cases has also been brought out.     

再论从一个怪例看细长梁理论的基本假定

针对均布力偶作用下两端固支欧拉梁挠度为零的状态, 通过欧拉梁挠度的四 阶微分方程分析欧拉梁挠度为零的本质. 再考虑6种边界下均布力偶作用下的欧拉梁, 发现 边界条件中剪力是否为零起决定性作用.     

Symmetry breaking in an initially curved pre-stressed micro beam loaded by a distributed electrostatic force

The symmetric and asymmetric buckling of an initially curved micro beam subjected to an axial pre-stressing load and transversal distributed electrostatic force is studied. The analysis is based on a reduced order (RO) model resulting from the Galerkin decomposition with buckling modes of a straight beam used as the base functions. The criteria of symmetric limit point buckling and of non-symmetric bifurcation are derived in terms of the geometric parameters of the beam and the axial load. Two symmetry breaking conditions, defining the relations between the axial load and the geometric parameters of beams for which an asymmetric response bifurcates from the symmetric one, are obtained. The necessary criterion establishes the conditions for the appearance of bifurcation points on the unstable branch of the symmetric response curve; the sufficient criterion assures a realistic asymmetric buckling bifurcating from the stable branches of the symmetric response curve. A comparison between the RO model results and those obtained by direct numerical analysis shows good agreement between the two and indicates that the obtained criteria can be used to predict symmetric and non-symmetric buckling in electrostatically actuated curved pre-stressed micro beams. It is shown that while the symmetry breaking conditions are affected by the nonlinearity of the electrostatic force, its influence is less pronounced than in the case of the symmetric snap-through criterion. The nature of the latter and the relations between it and the symmetry breaking criteria are found to go through a prominent qualitative change as the initial distance between the beam and the electrode, characterizing the electrostatic force, changes.     

Symmetry breaking in an initially curved micro beam loaded by a distributed electrostatic force

The asymmetric buckling of a shallow initially curved stress-free micro beam subjected to distributed nonlinear deflection-dependent electrostatic force is studied. In order to highlight the symmetry breaking phenomenon and the approach to its analysis, the subsidiary simplified problem of a curved beam attached to a linearly elastic foundation, and subjected to uniformly distributed “mechanical” load, which is independent of deflections, is addressed first. The analysis is based on a two degrees of freedom reduced order (RO) model resulting from the Galerkin decomposition with linear undamped eigenmodes of a straight beam used as the base functions. Simple approximate expressions are derived defining the geometric parameters of beams for which an asymmetric response bifurcates from the symmetric one. The necessary criterion establishes the conditions for the appearance of bifurcation points on the unstable branch of the symmetric limit point buckling curve; the sufficient criterion assures a realistic asymmetric buckling bifurcating from the stable branches of the curve. It is shown that while the symmetry breaking conditions are affected by the nonlinearity of the electrostatic force, its influence is less pronounced than in the case of the symmetric snap-through criterion. A comparison between the RO model results and those obtained by direct numerical analysis shows good agreement between the two and indicates that the obtained criteria can be used to predict non-symmetric buckling in electrostatically actuated bistable micro beams.     

Stability and optimal shape of Pflüger micro/nano beam

This paper deals with optimal shapes against buckling of an elastic nonlocal small-scale Pflüger beams with Eringen’s model for constitutive bending curvature relationship. By use of the Pontryagin’s maximum principle the optimality condition in form of a depressed quartic equation is obtained. The shape of the lightest (having the smallest volume) simply supported beam that will support given uniformly distributed follower type of load and axial compressive force of constant intensity without buckling, is determined numerically. A special attention is paid to the influence of the characteristic small length scale parameter of the nonlocal constitutive law to both critical load and optimal shape of the analyzed beams. For the case when distributed follower type of load is zero, our results reduce to those obtained recently for compressed nonlocal beam. Also the post buckling shape of the optimally shaped rod is studied numerically.     

基于均布荷载挠度曲率的梁结构损伤识别方法

为运用荷载挠度曲线进行损伤识别，通过力法推导了梁结构在均布荷载作用下的挠度曲率理论公式，提出通过损伤前后的挠度曲率差进行损伤定位，针对多跨连续梁均布荷载下挠度曲率指标存在的损伤识别漏判问题，提出逐跨均布荷载挠度曲率指标以避免其影响，并建立了挠度曲率与损伤程度的理论关系式，可对损伤程度进行较精确的定量描述．通过一简支梁和一三跨连续梁算例，考虑多种损伤工况，分析了指标在不同程度噪声水平下的抗干扰能力，验证了挠度曲率损伤指标应用于实际的可行性．     

Plastic collapse of general frames

A structural theory is presented for the large static plastic deformation of space frames composed of thin walled members. Displacements comparable to the overall structural dimensions are contemplated. The frame is considered to consist of an arbitrary number of beam elements connected at node points. The analysis assumes that plastic deformation is confined to idealized hinges located at the node points. As a basis for a general frame computer program, the equations for a beam element are derived as a relationship between appropriate generalized force and deformation rates. The structural constitutive theory employed for the plastic hinge includes biaxial bending, torsion, and axial extension. It accounts for reduction in the load carrying capacity of the hinge due to local deformation. Predicted force-deformation curves for a space frame are in good agreement with experimental results.     

Shear and Bending Response of a Rigid-Plastic Beam to Partly Distributed Blast-Type Loading

ABSTRACT

A study is undertaken on dynamic response of a simply supported rigid perfectly plastic beam that is subjected to partly distributed blast-type pressure loading. The beam material has finite shear strength and obeys a square yield criterion relating bending moment and transverse shear force. The transverse dynamic load is uniformly and symmetrically distributed over a middle portion of the span. Various patterns of deformation, which combine plastic bending and shear sliding, are obtained for a wide range of parameters, and the effects of transverse shear forces and time dependence of the dynamic pressure are examined.