有限到无限:基于离散静力平衡推导悬链线方程和梁位移方程

悬链线问题是一类非常经典的力学问题,推导出悬链线方程的方法从静力平衡到变分原理,非常之多。本文基于离散的静力平衡模型,采用从有限到无限的求解方法,建立并求解了悬链线方程。在此基础上,本文进一步建立了含有扭簧作用的离散静力平衡模型,推导并求解梁位移方程。最后在不同参数下求解梁位移方程得到其渐变解,并分析了多种平衡状态的存在情况。     

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

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

弹性悬链线方程参数变换法及其工程应用

建立了弹性锚泊线模型,该模型根据微段垂向力学平衡关系,通过引进一个变换参数u(sinhu=dy/dx),得到了弹性悬链线的参数方程.结合工程实际,将该模型用于海洋工程中的锚泊线,给出锚泊线拉力公式,与求解算法.该算法简单有效,易于编程实现.     

平面刚架弹塑性大位移分析的多刚体离散元法

本文基于多刚体－弹簧系统模型，给出了求解平面刚架结构弹塑性、大位移极限承载力分析的多刚体离散元法。文中首先推导了多刚体离散元法在总体坐标下的切线刚度阵，建立多刚体离散元法的增量平衡方程；而后推导了多刚体离散元的弹塑性弹簧系数矩阵，建立了多刚体离散元内力屈服面塑性铰法的增量求解格式，成功地进行了平面钢框架的弹塑性、大位移极限承载力分析。计算结果与其他数值方法或实验结果吻合良好，显示了多刚体离散元方法进行结构极限承载力分析这一复杂问题的优越性     

平面刚架弹塑性大位移分析的多刚体离散元法

本文基于多刚体－弹簧系统模型，给出了求解平面刚架结构弹塑性、大位移极限承载力分析的多刚体离散元法。文中首先推导了多刚体离散元法在总体坐标下的切线刚度阵，建立多刚体元法的增量平衡方程。     

水平悬链线中点集中力作用下的非线性分析及计算

首先,对均布力作用下悬链线理论进行了深入分析研究,推导出求解水平张力单参数和双参数形式的隐含超越方程,在算法上验证了两种方程的等价性。对于悬链线固有的广义倾角、曲率半径、约束反力和水平张力的嵌入关系等给出了具有新意的表达公式。然后,建立了水平悬链线在对称集中力作用下的平衡方程,对于不同垂向集中力作用下水平张力的变化情况和最大垂直距离差所在水平位置随集中力的变化等进行了非线性计算。本文在理论构建和分析方面具有创新性,在理论探索和工程应用上具有重要意义。     

悬索桥线形分析的悬链线单元法

从悬链线方程出发,以节点坐标和无应力索长作为未知量,推导了增量型索单元的刚度方程,提出了悬索桥在恒载作用下初始线形和无应力索长的计算公式及基于Newton-Raphson的求解方法.算例分析表明本文方法具有较高的精度.     

水平悬链线中点集中力作用下的非线性分析及计算

首先,对均布力作用下悬链线理论进行了深入分析研究,推导出求解水平张力单参数和双参数形式的隐含超越方程,在算法上验证了两种方程的等价性。对于悬链线固有的广义倾角、曲率半径、约束反力和水平张力的嵌入关系等给出了具有新意的表达公式。然后,建立了水平悬链线在对称集中力作用下的平衡方程,对于不同垂向集中力作用下水平张力的变化情况和最大垂直距离差所在水平位置随集中力的变化等进行了非线性计算。本文在理论构建和分析方面具有创新性,在理论探索和工程应用上具有重要意义。     

均布压力作用下悬链线波纹壳应力和位移的初参数积分方程分析

本文用初参数积分方程方法，对悬链线波纹壳在均布压力作用下的应力和位移进行了分析，并给出了一个算例。     

线性强化型弹塑性弯曲直梁挠曲线方程

针对材料在弹塑性阶段的应用不完全问题,本文用弹塑性分区最小势能原理,推导出线性强化模型下弯曲直梁的势能分区准则和欧拉方程.求解出集中载荷作用下悬臂梁和简支梁的挠曲线方程,将挠曲线方程代入MATLAB软件进行数值计算并将其结果与ANSYS对比分析.结果表明:数值解与有限元值均满足实际工程中允许的误差范围,给出的方法可为解...     

索牵引并联机器人静态平衡位置的正解研究

研究了一种大跨度但运动缓慢的索牵引并联机器人的位置正解问题,此时可认为系统处于静平衡状态,但必须考虑柔索垂度的影响。索牵引并联机器人的静态平衡位置的正解就是确定机器人末端执行器的位置和姿态,使得系统能在给定的柔索长度下处于静平衡状态。这涉及到系统的非线性力学分析而不仅仅是一个几何问题。本文采用悬链线方程描述柔索的力学特性,结合索端位移与机器人末端执行器位移之间的关系推导出机器人的末端执行器位移与索长变化之间的增量关系,据此提出一种迭代策略完成正解计算。数值算例验证了本文方法的有效性。     

STATIC STUDY OF CANTILEVER BEAM STICTION UNDER ELECTROSTATIC FORCE INFLUENCE

I. INTRODUCTION For capacitor-like microelectromechanical systems (MEMS) structure[1??6], the voltage betweenthe structure and substrate causes attractive force. The sources of the voltage can be an arti?ciallymounted device[2 , 7??9] or the temporar…     

具有动应力和动位移约束的离散变量结构形状优化设计方法

讨论了动应力、动位移约束下离散变量状优化设计问题。首先用拟静力算法,将结构惯性力极值作为静载荷施加到结构上,求得结构的动位移和动内力,然后将考虑动应力约束和动作移约束的离散变量结构优化设计问题化为静应力和静位移约束的优化问题。在求解过程中,将单元内力作了一阶近似,并将多约束问题转化为单约束问题,然后利用两类变量统一考虑的离散变量结构形状优化设计的综合算法进行求解。     

Modelling of dynamical crack propagation using time-domain boundary integral equations

We study dynamic antiplane cracks in the time domain by the boundary integral equation method (BIEM) based on the integral equation for displacement discontinuity (or crack opening displacement, COD) as a function of stress on the crack. This displacement discontinuity formulation presents the advantage, with respect to methods developed by Das and others in seismology, that it has to be solved only inside the crack. This BIEM is, however, difficult to implement numerically because of the hypersingularity of the kernel of the integral equation. Hence it is rewritten into a weakly singular form using a regularization technique proposed by Bonnet. The first step, following a method due to Sladek and Sladek, consists in converting the hypersingular integral equation for the displacement discontinuity into an integral equation for the displacement discontinuity and its tangential derivatives (dislocation density distribution); the latter involves a Cauchy type singular kernel. The second step is based on the observation that the hypersingularity is related to the static component of the kernel; the static singularity is then isolated and can be expressed in terms of weakly singular integrals using a result due to Bonnet. Although numerical applications discussed in this paper are all for the antiplane problem, the technique can be applied as well to in-plane crack dynamics.

The BIEM is implemented numerically using continuous linear space-time base functions to model the COD on the crack. In the present scheme the COD gradient interpolation is discontinuous at the element nodes while the integral equations are collocated at the element midpoints. This leads to an overdetermined discrete problem which is solved by standard least-squares methods. We use the dynamic BIEM to study a set of problems that appear in earthquake source dynamics, including the spontaneous dynamic crack propagation for a very simple rupture criterion. The numerical results compare favorably with the few exact solutions that are available. Then we demonstrate that difficulties experienced with finite difference simulations of spontaneous crack dynamics can be removed with the use of BIEM. The results are improved by the use of singular crack tip elements.     

Asymptotic solutions of coupled equations of supercritically axially moving beam

In supercritical regime, the coupled model equations for the axially moving beam with simple support boundary conditions are considered. The critical speed is determined by linear bifurcation analysis, which is in agreement with the results in the literature. For the corresponding static equilibrium state, the second-order asymptotic nontrivial solutions are obtained through the multiple scales method. Meantime, the numerical solutions are also obtained based on the finite difference method. Comparisons among the analytical solutions, numerical solutions and solutions of integro-partial-differential equation of transverse which is deduced from coupled model equations are made. We find that the second-order asymptotic analytical solutions can well capture the nontrivial equilibrium state regardless of the amplitude of transverse displacement. However, the integro-partial-differential equation is only valid for the weak small-amplitude vibration axially moving slender beams.     

DYNAMIC RESPONSE OF THE CATENARY RISER IN DEEP-SEA MINING SYSTEMS 1)

基于有限元模型,模拟、分析深海采矿系统悬链线立管在海流和水面船运动约束下的动力响应。结果表明：悬链线立管的最大等效应力和最大位移随时间呈周期性变化,且存在半个周期的相位差;当水面船运动到最高点时,悬链线立管位移达到最大值。悬链线立管最大 等效应力和最大位移随水面船运动位移幅值的增加而增加,随运动周期的增加而减小。进一步对比发现,水面船运动位移和周期对立管等效应力的影响大于对其位移的影响。     

关于Orr—Sommerfeld方程的Chebyshev谱方法的讨论

本文讨论了Ｏｒｒ－Ｓｏｍｍｅｒｅｌｄ方程的各种Ｃｈｅｂｙｓｈｅｖ谱离散方法，数值证明了Ｃｈｅｂｙｓｈｅｖ配置法离散Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程没有伪谱，并以此构造了适于任意平面平行速度剖面情形，对时间和时空稳定性模式一致有效的无伪谱的离离散方法，其中，对时空稳定性问题本文给出了一种新的迭代法可以快速有效地求出复频率的鞍点，对平面Ｐｏｉｓｅｕｉｌｌｅ流，Ｂｌａｓｉｕｓ边界层流和Ｇａｕｓｓ模型尾迹     

圆锥壳轴对称弯曲问题精确的四阶挠度微分方程和贝塞尔函数解

本文提出了轴对称圆锥壳精确的四阶挠度微分方程。和现行薄壳理论中常用的四阶剪力Q_1微分方程相比,挠度微分方程与其精度相同,阶数相同,而且满足边界条件简单,使圆锥壳的计算得到很大的简化。     

Response of an infinite beam resting on a tensionless elastic foundation subjected to arbitrarily complex transverse loads

This paper presents an investigation into the static response of an infinite beam supported on a unilateral (tensionless) elastic foundation and subjected to arbitrary complex loading, including self-weight. A new numerical method is developed to determine the initially unknown lengths that remain in contact. Based on the continuity conditions at the junctions of contact and non-contact segments, the response of the whole beam may be expressed through the displacement constants of the initial segment, reducing the contact problem to two nonlinear algebraic equations with two unknowns. The technique has been named the transfer displacement function method (TDFM). Comparison with the exact results of a particular limiting case shows the expected complete agreement. Finally, an example of a beam with several contact segments is presented and verified by the application of equilibrium conditions.