简支梁大挠度弹性后屈曲载荷与变形的优化算法

针对简支梁结构大挠度后屈曲载荷与变形的计算问题,本文提出了一种直接求解其后屈曲载荷和变形的优化算法。在简支梁处于大挠度屈曲平衡状态下,将梁结构划分为有限子段,以待求后屈曲载荷为设计变量,根据起点的边界条件和每个子段满足的弯矩变形公式,累积计算出其他各个节点的坐标,以得到的终点坐标满足的边界条件构建目标函数模型。在此基础上,通过MATLAB编制优化程序分析了两个典型算例,并将理论结果与相关软件的计算结果进行对比,从而证明了本文算法的正确性。本文算法求解过程简单、快速,具有一定的实用性,为变截面结构大挠度弹性屈曲稳定性问题的研究提供了参考。     

小波加权残值法在斜板后屈曲上的应用

将斜板在大挠度理论下的平衡方程和变形协调方程转化为斜坐标系下的表达式,并无量纲化,对于四边固支斜板,选用尺度小于1的二维bior3.1重构尺度函数与小波作试函数,满足板在平面外的几何与自然边界条件,利用小波-Galerkin法将板的无量纲平衡控制方程与变形协调方程转化为两个非线性方程组,从而将后屈曲路径的求解转化为两非线性方程组的求解问题.以载荷作为迭代步长,采用Newton-Raphson迭代算法求得承受单向压缩四边固支斜板在不同边长比、不同斜角下的后屈曲平衡路径及四级渐近解.     

变截面压杆稳定非线性微分方程边值问题的最优化算法研究

针对任意约束类型的变截面受压杆件的稳定临界载荷计算问题,结合非线性微分方程数值算法和最优化方法,以起点边界的初始条件、待求临界荷载和附加约束力为设计变量;以终点边值条件满足的函数关系与位型条件构建目标函数,提出变截面压杆临界载荷和稳定位型的优化求解算法。应用VB编制通用的优化计算程序,分析了典型算例;通过对比发现,本文以较少设计变量实现了临界载荷的高精度计算,为工程应用提供参考。     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

初曲矩形薄板的非线性动力屈曲研究

对两参数冲击载荷面内压缩作用下初曲矩形薄板的非线性动力屈曲问题进行了理论研究。首先采用双重余弦函数的组合确定了面内冲击矩形薄板的艾雷应力函数和中面力的分布；其次根据伽辽金法求得了初曲矩形薄板非线性动力屈曲问题的控制方程，基于巴拿赫压缩映象原理，采用逐次逼近方法求解了该控制方程。最后，应用本文发展的理论，给出了面内两参数冲击载荷作用下初曲矩形薄板动力屈曲响应的计算实例，计算结果与已有的实验结果较吻合     

层合析参固支边界条件下的不稳定后屈曲分析

为了深入地研究复合材料层板所独具的后屈曲特性,利用能量变分原理和非线性几何方程建立了具有弹性约束的复合材料层板在面内载荷作用下的非线性稳定性控制方程组,并运用广义傅立叶级数法对其进行求解。重点分析了非对称层板在固支边界条件下的稳定性问题,发现层板在此条件下有可能存在非对称的失稳临界点和不稳定的后屈曲路径,进而构造了简化的物理模型进行解释,指出后屈曲的非对称性是由于结构关于Z轴不对称,而不稳定性是由于固支边界条件阻碍了前屈曲的发生。     

旋转导向系统底部钻具组合力学计算新方法

底部钻具组合(Bottom Hole Assembly, BHA)的力学计算是进行旋转导向系统钻具组合设计和井眼轨迹控制的基础．本文在纵横弯曲法、加权余量法和有限元法的基础上,提出了一种基于管柱单元组合的BHA力学求解方法．该方法以保持微分方程特解相同的条件作为管柱单元划分依据,对底部钻具组合进行单元划分;将管柱单元的挠曲线通解替代传统加权余量法的试函数,结合单元节点类型和组合关系,构建出包含边界条件和通解系数的线性方程组．对于边界条件未知问题,以支点反力符号条件为约束条件,以挠曲线越界程度、切点处弯矩平衡条件为目标函数,建立了以边界条件未知量为变量的多目标优化数学模型．应用多目标优化算法求解出边界变量,进而得出通解系数;最后,得出管柱单元内各点位移和载荷．以推靠式旋转导向系统BHA力学计算为案例,对比本文方法与加权余量法的计算结果,发现两者计算结果相差较小,表明本文方法可用于旋转导向系统底部钻具组合的力学计算．     

结构后屈曲分析中的混合Newton-Lanczos算法

本文提出一种适于结构非线性后屈曲分析的混合Newton-Lanczos算法。与当前流行的弧长法不同,本文提出的算法采用传统的载荷增量法进行逐步求解,可求出给定载荷下的结构变形且适于任意外加载荷。对于临界载荷附近的迭代应用了Lanczos法求解方程及相应变载技巧。文中给出的若干数值计算结果表明了该算法在结构非线性后屈曲分析中的适用性。     

热环境中功能梯度圆柱薄壳分岔屈曲的边界约束效应

基于前屈曲一致理论,研究了热环境中受轴压功能梯度圆柱薄壳分岔屈曲的边界约束效应问题.导出前屈曲变形的解析解,结合分离变量法与有限差分法求解分岔屈曲控制方程,由此导出确定临界轴压的非线性特征值问题.考虑材料热物性质与温度的相关性,分别就两端简支和两端固支边界条件,分析了温度梯度、初始几何缺陷、组分材料体积分数对分岔屈曲临...     

受火作用铝合金梁的临界温度及其参数分析

基于ANSYS数值模拟结果,首次对铝合金梁在火荷载作用下的响应行为进行模拟。系统比较了简支、铰支、固支、一端固支另一端简支等几种典型边界约束条件下铝合金梁的抗火特性,着重讨论了各自的临界温度并与已有规范进行了对比,在此基础上考察了端部约束条件、载荷比、载荷类型等因素对铝合金梁弯曲屈曲失效行为的影响。结果表明:设计规范公式计算出的临界温度偏于保守,比有限元结果低2.5%~16.7%,而有限元更能反映真实情况;对于不同的约束情形,铰支梁的抗火能力最强,其次是固支梁,简支梁的抗火能力最差;载荷比越大,铝合金梁的临界温度越低,抗火能力越弱,载荷比由0.3增加到0.7时,简支梁和一端固支另一端简支梁的临界温度减小约43℃,固支和铰支梁的临界温度分别减小37℃和30℃;载荷类型对铝合金梁的影响较小。本文研究为铝合金结构抗火设计提供了理论依据。     

Buckling and post-buckling analyses of size-dependent piezoelectric nanoplates

This paper attempts to investigate the buckling and post-buckling behaviors of piezoelectric nanoplate based on the nonlocal Mindlin plate model and von Karman geometric nonlinearity. An external electric voltage and a uniform temperature rise are applied on the piezoelectric nanoplate. Both the uniaxial and biaxial mechanical compression forces will be considered in the buckling and post-buckling analysis. By substituting the energy functions into the equation of the minimum total potential energy principle,the governing equations are derived directly, and then discretized through the differential quadrature(DQ) method. The buckling and post-buckling responses of piezoelectric nanoplates are calculated by employing a direct iterative method under different boundary conditions. The numerical results are presented to show the influences of different factors including the nonlocal parameter, electric voltage,and temperature rise on the buckling and post-buckling responses.     

Bifurcation and post-buckling analysis of bimodal optimum columns

A mathematical formulation of column optimization problems allowing for bimodal optimum buckling loads is developed in this paper. The columns are continuous and linearly elastic, and assumed to have no geometrical imperfections. It is first shown that bimodal solutions exist for columns that rest on a linearly elastic (Winkler) foundation and have clamped-clamped and clamped-simply supported ends. The equilibrium equation for a non-extensible, geometrically nonlinear elastic column is then derived, and the initial post-buckling behaviour of a bimodal optimum column near the bifurcation point is studied using a perturbation method. It is shown that in the general case the post-buckling behaviour is governed by a fourth order polynomial equation, i.e., near the bifurcation point there may be up to four post-buckling equilibrium states emanating from the trivial equilibrium state. Each of these equilibrium states may be either supercritical or subcritical in the vicinity of the bifurcation point. The conditions for stability of these non-trivial post-buckling states are established based on verification of positive semi-definiteness of a two-by-two matrix whose coefficients are integrals of the buckling modes and their derivatives. In the end of the paper we present and discuss numerical results for the post-buckling behaviour of several columns with bimodal optimum buckling loads.     

Initial thermo-mechanical post-buckling of beams with temperature-dependent physical properties

The objective of this work is to analyze the elastic buckling and initial post-buckling behavior of slender beams subjected to uniform heating. The beams are assumed to be double-hinged with fixed ends, preventing thermal expansion. Consequently, destabilizing compressive forces arise that may lead to beam buckling. When the temperature is further increased, the beam experiences finite displacements, with the result that the analysis is geometrically non-linear. The modulus of elasticity and the thermal induced strain, key material properties for this problem, are temperature-dependent. Thus, the coefficients of the governing equations are not constant. This suggests the physical non-linearity of the mathematical model. Hence, the analysis is geometrically and physically non-linear. The analysis is sensitive to the beam initial temperature, as the thermal strain is a function of the initial and final temperatures. The material is considered to be linear elastic, and consequently viscoelastic and plastic effects are not taken into account. Furthermore, the beam cross-section properties are assumed to be constant, which is consistent with the small strain formulation. A perturbation method is applied to the governing non-linear differential equations so that the initial post-buckling behavior may be analytically determined when temperatures above the critical temperature are applied to the beam. To illustrate the application of the formulation we present a case study for the aluminum 7075-T6 alloy, a material commonly used in aerospace and naval industries. Nonetheless, it is expected similar behavior for other metallic materials. The curves that define the variation of the modulus of elasticity, the thermal strain and the yield stress with temperature are considered in our analysis. The change in length, reaction forces at the supports and geometric configurations are obtained as a function of temperature and the beam slenderness ratio. The critical buckling loads and temperatures and the initial post-buckling analysis are also calculated in the context of the temperature-independent physical properties. Our results emphasize the importance of modeling the material's non-linearity if accuracy is required. However, from a practical application point of view results are acceptable if temperature-independent physical properties are employed, especially for large slenderness ratios.     

Experimental Validation of Reduction Methods for Nonlinear Vibrations of Distributed-Parameter Systems: Analysis of a Buckled Beam

An experimental validation of the suitability of reduction methods for studying nonlinear vibrations of distributed-parameter systems is attempted. Nonlinear planar vibrations of a clamped-clamped buckled beam about its first post-buckling configuration are analyzed. The case of primary resonance of the nth mode of the beam, when no internal resonances involving this mode are active, is investigated. Approximate solutions are obtained by applying the method of multiple scales to a single-mode model discretized via the Galerkin procedure and by directly attacking the governing integro-partial-differential equation and boundary conditions with the method of multiple scales. Frequency-response curves for the case of primary resonance of the first mode are generated using both approaches for several buckling levels and are contrasted with experimentally obtained frequency-response curves for two test beams. For high buckling levels above the first crossover point of the beam, the computed frequency-response curves are qualitatively as well as quantitatively different. The experimentally obtained frequency-response curves for the directly excited first mode are in agreement with those obtained with the direct approach and in disagreement with those obtained with the single-mode discretization approach.     

Creep buckling and post-buckling analysis of the laminated piezoelectric viscoelastic functionally graded plates

The creep buckling and post-buckling of the laminated piezoelectric viscoelastic functionally graded material (FGM) plates are studied in this research. Considering the transverse shear deformation and geometric nonlinearity, the Von Karman geometric relation of the laminated piezoelectric viscoelastic FGM plates with initial deflection is established. And then nonlinear creep governing equations of the laminated piezoelectric viscoelastic FGM plates subjected to an in-plane compressive load are derived on the basis of the elastic piezoelectric theory and Boltzmann superposition principle. Applying the finite difference method and the Newmark scheme, the whole problem is solved by the iterative method. In numerical examples, the effects of geometric nonlinearity, transverse shear deformation, the applied electric load, the volume fraction and the geometric parameters on the creep buckling and post-buckling of laminated piezoelectric viscoelastic FGM plates with initial deflection are investigated.     

Nonlinear free vibration and post-buckling analysis of functionally graded beams on nonlinear elastic foundation

In this study, simple analytical expressions are presented for large amplitude free vibration and post-buckling analysis of functionally graded beams rest on nonlinear elastic foundation subjected to axial force. Euler–Bernoulli assumptions together with Von Karman’s strain–displacement relation are employed to derive the governing partial differential equation of motion. Furthermore, the elastic foundation contains shearing layer and cubic nonlinearity. He’s variational method is employed to obtain the approximate closed form solution of the nonlinear governing equation. Comparison between results of the present work and those available in literature shows the accuracy of this method. Some new results for the nonlinear natural frequencies and buckling load of the FG beams such as the effect of vibration amplitude, elastic coefficients of foundation, axial force, and material inhomogenity are presented for future references.     

基于非线性分析的加肋板肋条位置无网格优化

在加肋板无网格模型中,肋条的位置对各种工况下加肋板受力性能的影响至关重要.文章基于一阶剪切变形和移动最小二乘法理论提出一种考虑非线性影响的加肋板无网格模型,并利用遗传算法优化肋条位置.首先,采用离散节点分别对平板和肋条进行离散,得到加肋板的无网格离散模型;其次,通过冯·卡门大挠度理论得到非矩形板几何非线性问题的弯曲控制方程;再次,通过哈密顿原理得到加肋非矩形板自由振动问题的控制方程;最后引入遗传算法,以肋条的位置为设计变量、非矩形加肋板中心点挠度最小或自振频率最大为目标函数,对肋条位置进行优化.在考虑了几何非线性影响的肋条位置优化过程中,肋条位置改变时只需重新计算位移转换矩阵,避免了网格重构.本文以全局荷载下单肋条菱形板为例与理论解进行对比,进行有效性验证.再以板的中点挠度最小和自振频率最大为优化目标,对局部荷载作用下不同形状、不同肋条布置方式的加肋板进行优化,分析方法的收敛性及稳定性.