具有分数导数本构关系的非线性粘弹性Timoshenko梁动力学行为分析

本文利用分数导数型本构关系建立了在有限变形情况下Timoshenko梁的控制方程并利用Galerkin方法进行简化。然后利用一种存储部分历史数据的分数积分的计算方法对梁的控制方程进行求解。考察了载荷参数和分数导数参数对梁振动的影响,并采用非线性动力学中的各种数值方法,如时程曲线、功率谱、相图、Poincare截面等,揭示了非线性粘弹性Timoshenko梁丰富的动力学行为。     

在有限变形条件下损伤粘弹性梁的动力学行为

本文在有限变形条件下，根据损伤粘弹性材料的一种卷积型本构关系和温克列假设，建立了粘弹性基础上损伤粘弹性Timoshenko梁的控制方程。这是一组非线性积分——偏微分方程。为了便于分析，首先利用Galerkin方法对该方程组进行简化，得到一组非线性积分一常微分方程。然后应用非线性动力学中的数值方法，分析了粘弹性地基上损伤粘弹性Timoshenko梁的非线性动力学行为，得到了简化系统的相平面图、Poincare截面和分叉图等。考察了材料参数和载荷参数等对梁的动力学行为的影响。特别，考察了基础和损伤对粘弹性梁的动力学行为的影响。     

非线性粘弹性Timoshenko梁动力学行为的分布

本文讨论了有限变形粘弹性Timoshenko梁的动力学行为。首先由Timoshenko梁的理论和分数导数型本构关系给出了梁的控制方程。其次为了便于求解,采用Galerkin方法对系统进行了简化,并比较了1阶和2阶截断系统的动力学性质,它们具有相同的定性性质,说明Galerkin方法的合理性。给出了求解包含分数积分的积分－微分方程的一种新方法,以便求解系统的长时间的解。综合利用非线性动力系统中的经典方法,揭示了梁在有限变形情况下丰富的动力学行为,并分别考察了载荷参数的材料参数对结构的动力学行为的影响。     

复合载荷作用下具有光滑中心波纹膜片的非线性分析

采用轴对称旋转壳体的简化Reissner方程，研究了在复合载荷作用下具有光滑中心波纹膜片的非线性弯曲问题。应用积分方程方法，可以获得膜片的特征关系（载荷－中心挠度关系）。文末给出了实例计算的数值结果。     

二自由度结构受迫非线性振动响应的数值计算

本文研究了由飞机尾部简化而成的二自由度结构受迫非线性振动,由于对复杂的运动方程不作进一步简化,方程中含有强非线性耦合项,首先采用数字滤波技术对系统受确定性激励的响应进行了计算;然后在数字模拟计算的基础上用嫁接的ITO积分方法对系统受随机白噪音激励的响应进行了计算。结果表明,ITO积分的综合嫁接方法有更高的计算效率,结果并表明系统的非线性耦合部分对响应有重大影响,是不能加以忽略的,本文的研究为多自由度非线性结构的响应计算提供了有用的思想和方法。     

精细积分的非线性动力学积分方程及其解法

给出了非线性动力学积分方程的表达式,针对该方程提出了一个显式预测－校正的单步四阶精度的积分算法,适用于多自由度、强非线性,非保守系统。算例表明该方法精度高、计算量较少。     

剪切可变形梁非线性静态响应的精确解

本文给出了纵横向载荷作用下,梁非线性静态问题的精确解。基于非线性一阶剪切变形梁理论,导出了梁非线性静态问题的基本方程。将三个非线性方程化简为一个关于横向挠度的非齐次四阶非线性积分-微分方程,当只有轴向载荷作用时,该方程和相应的边界条件构成微分特征值问题。直接求解该方程,得到了梁非线性静态变形闭合形式的解,这个解显式地给出了梁的变形与外载荷之间的非线性关系,描述了梁变形后的非线性平衡路径。利用这个解,得到了梁临界屈曲载荷的一阶结果与经典结果。为考察载荷、长高比以及边界条件的影响,根据得到的解析解给出了一些数值算例,并讨论了梁不同阶屈曲模态下非线性静态响应的一些性质。结果表明:对应于方程特征参数λ的不同取值区间,梁的轴向载荷-挠度曲线有不同的解支;而对应于参数λ的同一取值区间,梁分别对应两个不同的屈曲模态。     

NONLINEAR VIBRATIONS OF AXIALLY ACCELERATING VISCOELASTIC TIMOSHENKO BEAMS

研究了轴向加速黏弹性Timoshenko梁的非线性参数振动。参数激励是由径向变化张力和轴向速度波动引起的。引入了取决于轴向加速度的径向变化张力,同时还考虑了有限支撑刚度对张力的影响。应用广义哈密尔顿原理建立了Timoshenko梁耦合平面运动的控制方程和相关的边界条件。黏弹性本构关系采用Kelvin模型并引入物质时间导数。耦合方程简化为具有随时间和空间变化系数的积分-偏微分型非线性方程。采用直接多尺度法分析了Timoshenko梁的组合参数共振。根据可解性条件得到了Timoshenko梁的稳态响应,并应用Routh-Hurvitz判据确定了稳态响应的稳定性。最后通过一系列数值例子描述了黏弹性系数、平均轴向速度、剪切变形系数、转动惯量系数、速度脉动幅值、有限支撑刚度参数以及非线性系数对稳态响应的影响。     

分数导数型本构关系描述粘弹性梁的振动分析

本文研究粘弹性梁在周期激励作用下的受迫振动问题。梁的材料满足Kelvin-Volgt分数导数型本构关系。基于动力学方程、本构关系和应变－位移关系建立了小变形粘弹性梁的振动方程。采用分离变量法分析粘弹性梁的自由振动,导出模态坐标满足的常微分－积分方程和模态函数满足的常微分方程,对于两端简支的截面梁给出了固有频率和模态函数。对于简谐激励作用下粘弹性梁的受迫振动,利用模态叠加得到了稳态响应。最后给出数值算例说明本文方法的应用。     

非线性MEMS微梁的重心有理插值迭代配点法分析

文章利用重心有理插值迭代配点法分析计算非线性MEMS微梁问题。通过处理MEMS微梁的几何通过假设初始函数,将微梁非线性控制方程转换为线性化微分方程,建立逼近非线性微分方程的线性化迭代格式。采用重心有理插值配点法求解线性化微分方程,提出了数值分析MEMS微梁非线性弯曲问题的重心插值迭代配点法。给出了非线性微分方程的直接线性化和Newton线性化计算公式,详细讨论了非线性积分项的计算方法和公式。利用重心有理插值微分矩阵,建立了矩阵-向量化的重心插值迭代配点法的计算公式。数值算例结果表明,重心插值迭代配点法求解微梁非线性弯曲问题,具有计算公式简单、程序实施方便和计算精度高的特点。     

Dynamical behavior of nonlinear viscoelastic beams

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＧａｌｅｒｋｉｎｔｒｕｎｃａｔｉｏｎｉｓｗｉｄｅｌｙｕｓｅｄｔｏｓｔｕｄｙｔｈｅｄｙｎａｍｉｃａｌｂｅｈａｖｉｏｒ(ｅｓｐｅｃｉａｌｌｙｔｈｅｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃａｌｂｅｈａｖｉｏｒ)ｏｆｓｔｒｕｃｔｕｒｅｓ[1].Ｈｏｗｅｖｅｒ,ａｓｆａｒｔｈｅｒｅｉｓｎｏｄｉｒｅｃｔｅｖｉｄｅｎｃｅｔｏｐｒｏｖｅｔｈｅｐｌａｕｓｉｂｉｌｉｔｙｏｆｔｈｅｌｏｗｏｒｄｅｒＧａｌｅｒｋｉｎｔｒｕｎｃａｔｉｏｎ,ａｌｔｈｏｕｇｈｉｔｃａｎｂｅｉｎｆｅｒｒｅｄｆｒｏｍｃｅｒｔａｉｎｉ…     

Nonlinear forced vibration analysis of postbuckled beams

A numerical solution methodology is proposed herein to investigate the nonlinear forced vibrations of Euler–Bernoulli beams with different boundary conditions around the buckled configurations. By introducing a set of differential and integral matrix operators, the nonlinear integro-differential equation that governs the buckling of beams is discretized and then solved using the pseudo-arc-length method. The discretized governing equation of free vibration around the buckled configurations is also solved as an eigenvalue problem after imposing the boundary conditions and some complicated matrix manipulations. To study forced and nonlinear vibrations that take place around a buckled configuration, a Galerkin-based numerical method is applied to reduce the partial integro-differential equation into a time-varying ordinary differential equation of Duffing type. The Duffing equation is then discretized using time differential matrix operators, which are defined based on the derivatives of a periodic base function. Finally, for any given magnitude of axial load, the pseudo -arc-length method is used to obtain the nonlinear frequencies of buckled beams. The effects of axial load on the free vibration, nonlinear, and forced vibrations of beams in both prebuckling and postbuckling domains for the lowest three vibration modes are analyzed. This study shows that the nonlinear response of beams subjected to periodic excitation is complex in the postbuckling domain. For example, the type of boundary conditions significantly affects the nonlinear response of the postbuckled beams.     

Approximate solutions of axially moving viscoelastic beams subject to multi-frequency excitations

Approximate solutions of axially moving viscoelastic beams subject to multi-frequency excitations are presented. A non-linear partial-differential equation governing the transverse vibration of the beams is derived from Newton's second law, the Kelvin constitutive relationship, and the Lagrangian strain. Based on 1-term Galerkin's truncation, the governing equation is reduced to an ordinary differential equation. Three cases, including superharmonic resonance case, subharmonic resonance, and combination resonance are studied. The approximate solutions of the transverse vibration of the beams are obtained. Numerical results show that the approximate solutions are in good agreement with numerical results.     

Isospectral Euler-Bernoulli beams via factorization and the Lie method

We obtain isospectral Euler-Bernoulli beams by using factorization and Lie symmetry techniques. The canonical Euler-Bernoulli beam operator is factorized as the product of a second-order linear differential operator and its adjoint. The factors are then reversed to obtain isospectral beams. The factorization is possible provided the coefficients of the factors satisfy a system of non-linear ordinary differential equations. The uncoupling of this system yields a single non-linear third-order ordinary differential equation. This ordinary differential equation, called the principal equation, is analyzed, reduced or solved using Lie group methods. We show that the principal equation may admit a one-dimensional or three-dimensional symmetry Lie algebra. When the principal system admits a unique symmetry, the best we can do is to depress its order by one. We obtain a one-parameter family of invariant solutions in this case. The maximally symmetric case is shown to be isomorphic to a Chazy equation which is solved in closed form to derive the general solution of the principal equation. The transformations connecting isospectral pairs are obtained by numerically solving systems of ordinary differential equations using the fourth-order Runge-Kutta method.     

考虑恒载效应的拱形梁静力近似解

应用虚功原理,推导了考虑恒载效应影响时拱形梁在活载作用下的非线性微分方程,得到了方程的近似闭合解。根据方程的解,讨论了恒载大小及结构自身刚度(矢高、跨度、惯性矩及惯性半径等)不同因素在考虑恒载效应时对拱形梁静力特性的影响。通过与Takabatake得到的直梁解析解结果及作者在其他文献提出的有限元方法对拱形梁分析结果的比较,验证了本文非线性微分方程及其求解公式。结果表明,本文给出的非线性微分方程对于拱形梁和直线梁具有通用性,初始恒载的存在减小了拱形梁在活载作用下的静力反应,这种影响与恒载的大小及结构自身的刚度有关,对轻型结构的设计提出了一些建议。     

考虑收缩徐变影响的新旧混凝土叠合梁单元刚度矩阵改进

混凝土自身的收缩徐变会在新旧混凝土叠合梁中使应力重分布.为了计算重分布应力,首先推导以挠度表达的叠合梁非线性微分方程,然后通过求解该微分方程,引入位移形函数、刚度形函数和等效节点载荷形函数,最后得出新混凝土梁、旧混凝土梁和Goodman弹性夹层三合一的叠合梁改进型单元刚度矩阵和等效节点载荷,从而为收缩徐变影响下的混凝土的内力计算提供了一种有效的新方法.文中还进行了实例验证分析,并从中得出了一些有益的结论.     

Analytical and experimental studies on nonlinear characteristics of an L-shape beam structure

This paper focuses on theoretical and experimental investigations of planar nonlinear vibrations and chaotic dynamics of an L-shape beam structure subjected to fundamental harmonic excitation,which is composed of two beams with right-angled L-shape.The ordinary differential governing equation of motion for the L-shape beam structure with two-degree-of-freedom is firstly derived by applying the substructure synthesis method and the Lagrangian equation.Then,the method of multiple scales is utilized to obtain a four-dimensional averaged equation of the L-shape beam structure.Numerical simulations,based on the mathematical model,are presented to analyze the nonlinear responses and chaotic dynamics of the L-shape beam structure.The bifurcation diagram,phase portrait,amplitude spectrum and Poincare map are plotted to illustrate the periodic and chaotic motions of the L-shape beam structure.The existence of the Shilnikov type multi-pulse chaotic motion is also observed from the numerical results.Furthermore, experimental investigations of the L-shape beam structure are performed,and there is a qualitative agreement between the numerical and experimental results.It is also shown that out-of-plane motion may appear intuitively.     

Timoshenko-Euler楔形梁有限元

本文首先建立楔形梁包含轴力和剪切变形效应的平衡微分方程，由于该方法是二阶变系数微分方程，其争析解很难得到，本文通过将该方程中的变系数和方程的解用Chebyshev多项式逼近得到了Timoshenko-Euler楔形梁的单元刚度方程，最后通过算例检验了所得单元刚度方程的对称性，以及验证了计算悬臂梁挠度和悬臂柱弹性临界力的正确性及其收敛性，本文提出的方法可适用于任意变截面Timoshenko-Euler梁单元刚度方程的求解，运用此方法，除可以考虑轴力和剪切变形的影响外，还可以减少结构分析中的单元数和自由度，提高包含楔形构件的结构分析的精度和速度。