简支梁大变形的优化算法

提出一种求解几何非线性问题的优化算法,并研究了简支梁的几何非线性大变形问题。首先取简支梁大变形后的平衡状态为研究对象,分别创建它的变分模型和微分模型;然后基于微分模型,通过动坐标的迭代关系式求得微段端点坐标,构建微段端点未知坐标的目标函数;最后确定简支梁几何非线性大变形的最优化问题,并编制相应优化程序进行求解。通过分析典型算例,并同有限元方法的计算结果相比较,表明提出的优化算法在求解强几何非线性大变形问题中的正确性,为处理几何非线性大变形问题提供了一种新方法。     

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

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

逐段变形效应叠加法在简支梁的应用

逐段变形效应叠加法可应用于简支梁.分别计算了受对称载荷、非对称载荷的阶梯简支梁中截面的挠度及受对称载荷作用的阶梯简支梁端截面的转角.逐段变形效应叠加法可计算受任意载荷的阶梯简支梁中截面的挠度和端截面的转角.     

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

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

基于改进粒子群的薄壁变截面刚架临界载荷优化算法

针对大型变截面薄壁结构的稳定问题,以一类任意约束对称结构受非对称载荷的单跨刚架为研究对象,结构拆分为相关铁木辛柯（Timoshenko）梁,结合差分原理和最优化方法,以每段刚架的每个离散点挠度、临界载荷、轴力、剪力和梁端弯矩为设计变量,建立求解满足边界条件的非线性差分方程模型,提出基于优胜劣汰粒子更新的粒子群（IPSO）临界载荷优化算法。运用JAVA编程语言编制对应优化程序,分析典型算例并核实ABAQUS仿真结果。研究表明,本文提出的优化算法获得了有效的变形位型和高精度的临界载荷计算,能更好地描述刚架受力下位型和载荷的力学关系,进一步为工程设计与分析提供支持。     

基于改进粒子群的薄壁变截面刚架临界载荷优化算法

针对大型变截面薄壁结构的稳定问题,以一类任意约束对称结构受非对称载荷的单跨刚架为研究对象,结构拆分为相关铁木辛柯(Timoshenko)梁,结合差分原理和最优化方法,以每段刚架的每个离散点挠度、临界载荷、轴力、剪力和梁端弯矩为设计变量,建立求解满足边界条件的非线性差分方程模型,提出基于优胜劣汰粒子更新的粒子群(IPSO)临界载荷优化算法。运用JAVA编程语言编制对应优化程序,分析典型算例并核实ABAQUS仿真结果。研究表明,本文提出的优化算法获得了有效的变形位型和高精度的临界载荷计算,能更好地描述刚架受力下位型和载荷的力学关系,进一步为工程设计与分析提供支持。     

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

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

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

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

经典梁热过屈曲问题的解析解

基于非线性经典梁理论,建立了控制轴向和横向变形的基本方程,将两个非线性方程化简为一个关于横向挠度的四阶非线性积分-微分方程。对于本文所考虑的三类边界条件,该方程与相应的边界条件构成了微分特征值问题;直接求解该问题,得到热过屈曲构形的解析解,该解是外加热载荷的函数。为考察热载荷以及边界条件的影响,根据得到的解析解给出了一些数值算例,讨论了梁过屈曲行为的性质。本文得到的解析解可用于验证或改进各类近似理论和数值方法。     

考虑材料非线性的两端固支柱后屈曲载荷优化算法

针对满足Ramberg-Osgood本构关系的6082-T6铝合金固支柱的后屈曲问题进行了研究。利用泰勒展开式得到了用于计算横截面上弯矩的应力显式表达式，并通过弯矩-曲率关系导出了用剪力表示的控制方程。由于同时考虑材料非线性和几何非线性，此时控制方程为二阶非线性微分方程，本文提出了求解包含屈曲载荷的控制方程的优化算法。以屈曲载荷和固支端处的曲率为设计变量，以固支柱中点处的转角和挠度形成目标函数。利用进退法和黄金分割法改变设计变量的值，通过程序中包含的R-K法输出不同的结果，然后将输出的结果代入到目标函数中进行比较，获得包含目标函数极小值的最小区间，最终实现了对满足计算精度的设计变量值的确定。相较于打靶法复杂的分析过程，该优化算法优化过程简单，计算速度较快。为了验证本文算法的正确性，与两端固支6082-T6铝合金柱后屈曲时的数值解进行了对比。     

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.     

Development mechanism of local plastic buckling in bars subjected to axial impact

In order to clarify the developmental mechanism of the local plastic buckling and the interaction between axial wave and buckling deformation in an axially impacted slender-bar, the non-linear dynamic equations in the incremental form are derived and solved by use of the finite difference method, with the axial wave front treated as a moving boundary. The initial local-buckling deflection given by the characteristic-value analysis is used as the initial condition of the solution of the equations, instead of the initial imperfection that is assumed in literatures. It is found that the initial buckling deflection with one half-wave, occurring near the impacted end, develops into the higher post-buckling mode with several half-waves, as the axial compression waves propagate forward. The numerical results show that no strain reversal occurs at the early stage of post-buckling process, and the solution corresponding to the tangent-modulus theory is valid for the dynamic plastic post-buckling response of the bar at this stage. The theoretical results are in good agreement with the experimental results reported in the literature.     

Secondary buckling of a flat plate under uniaxial compression : Part 1: theoretical analysis of simply supported flat plate

While studies of post-buckling behavior and load-carrying capacities of thin plates subjected to uniaxial compression have been limited to stable conditions, further post-buckling loading generates an unstable condition. The secondary buckling which occurs with snap-through to higher-order deflections under such unstable conditions has not been analyzed in detail as yet. In the first part of this paper, a thin square plate under uniaxial compression, which is simply supported along four edges, is considered. A method based on the second variation of the total potential energy is then proposed for evaluating the stability of the post-buckling equilibrium state and inevitable secondary buckling is derived analytically. The effects of various factors, such as initial imperfections, assumed virtual displacement pattern, post-buckling deflection pattern and in-plane boundary conditions, on the secondary buckling values are discussed. In part 2, secondary buckling of clamped plates is analyzed by use of the finite element method and the resultant numerical results are compared with experimental results.     

考虑弯扭形变的薄壁箱形梁的非线性后屈曲

在大位移和扭转的前提下,通过一中等弯曲扭转的位移场描述了薄壁箱形梁在偏心载荷作用下的静稳定性问题.该非线性公式可用于分析简支薄壁箱形梁在不同载荷作用下的屈曲和后屈曲行为.采用伽辽金方法将非线性微分系统离散,并通过牛顿-拉普森增量迭代法求解得代数方程组.数值计算结果表明,当前屈曲位移不可忽略时,经典的横向屈曲预测是保守的...     

Buckling and post-buckling analysis for magneto-elastic–plastic ferromagnetic beam-plates with unmovable simple supports

This paper presents an analysis of buckling/snapping, bending and post-buckling/snapping behaviors of magneto-elastic–plastic interaction and coupling for soft ferromagnetic beam-plates with geometrically nonlinear deformation and unmovable simple supports at the ends of the plates. Based on the expression of magnetic force from the variational principle of ferromagnetic plates, the theory of thin plates with the nonlinear deformation of van Karman’s type, and the Mises yield criterion and the increment theory for plastic deformation, here, we establish a numerical code to quantitatively simulate the behaviors of the nonlinearly multi-coupling problems by the finite element method. Along with that the phenomena of buckling, bending, and post-buckling/snapping, or the characteristic curve of deflection versus magnitude of applied magnetic fields are numerically displayed, the critical values of buckling/snapping and yield magnetic fields, and the expansibility of plastic region after the plates undergo plastic deformation with increasing of the applied magnetic fields, as well as the evolvement of deflection configuration of the plate are numerically obtained in a case study.     

Buckling and post-buckling of annular plates under non-axisymmetric plane edge forces

On the basis of both the general theory [1,2] and the finite element method [4] of perforated thin plates with large deflection,the buckling and post-buckling of annular plates under non-axisymmetric plane edge forces are studied.The Project supported by the State Education Commission of the People's Republic of China     

细长柔韧压杆弹性失稳后挠曲线形状的计算机仿真

采用Maple编程对细长柔韧压杆弹性失稳后挠曲线形状进行了计算机仿真,进行了细长柔韧压杆弹性失稳后最大挠度和挠曲线封闭两种情况下的挠曲线形状仿真和详细的解答.分析计算了失稳后屈曲的力学特征,给出了解析表达式;分析计算了失稳后屈曲的平衡状态曲线的几何特征,绘出了计算机仿真曲线.结果表明:失稳后最大挠度和挠曲线封闭是属于两个完全不同的屈曲状态.     

功能梯度压杆后屈曲形态的数值计算

基于轴线可伸长杆的过屈曲精确数学模型,采用打靶法对两端简支功能梯度压杆后屈曲行为进行了数值分析.其中假设功能梯度梁的材料性质沿厚度方向按照幂函数连续变化.给出了不同梯度指标下FGM杆的后屈曲特征曲线,并与金属和陶瓷两种单相材料杆的相应特性进行了比较,分析和讨论了载荷、材料的梯度指数、长细比对杆后屈曲平衡路径的影响.     

Nonlinear thermo-elastic buckling characteristics of cross-ply laminated joined conical–cylindrical shells

Here, the nonlinear thermo-elastic buckling/post-buckling characteristics of laminated circular conical–cylindrical/conical–cylindrical–conical joined shells subjected to uniform temperature rise are studied employing semi-analytical finite element approach. The nonlinear governing equations, considering geometric nonlinearity based on von Karman’s assumption for moderately large deformation, are solved using Newton–Raphson iteration procedure coupled with displacement control method to trace the pre-buckling/post-buckling equilibrium path. The presence of asymmetric perturbation in the form of small magnitude load spatially proportional to the linear buckling mode shape is assumed to initiate the bifurcation of the shell deformation. The study is carried out to highlight the influences of semi-cone angle, material properties and number of circumferential waves on the nonlinear thermo-elastic response of the different joined shell systems.