复合材料三角形网格加筋圆锥壳体位移型稳定性方程及其总体稳定性分析

本文采用Donnell型扁壳理论,首先利用最小势能原理和广义平均筋条刚度法推导出用位移分量表示的复合材料三角形网格加筋叠层圆锥壳体的稳定性方程,考虑了蒙皮最一般的拉弯与拉扭耦合关系和加筋筋条的偏心效应,并讨论了该方程的基本性质.根据外压实验观察结果,通过选取适当的位移分量表达式,并运用Galerkin法分析了在均布外压作用下复合材料三角形网格加筋叠层圆锥壳体总体稳定性,得到了临界载荷的解析表达式,并对某一类C/E复合材料三角形内网格加筋圆锥壳体的临界外压进行了计算,所得理论值与实验结果很好地吻合.最后,讨论了有关参数对临界载荷的影响.本文所建立的新方程和所得结果对于航空航天结构非常有用.     

阶梯柱屈曲的改进Fourier级数分析

该文对阶梯柱的弹性屈曲问题进行了研究。首先基于改进Fourier级数法采用局部坐标逐段建立阶梯柱的位移函数表达式,然后由带约束的势能变分原理得到含屈曲荷载的线性方程组,利用线性方程组有非零解的条件把问题转化为矩阵特征值问题得到临界载荷,最后讨论方法中的参数取值,并把结果与已有文献和有限元的结果比较,从而验证方法的精度。所提模型在阶梯柱的两端和变截面处引入横向弹簧和旋转弹簧,通过改变弹簧的刚度值模拟不同的边界。所提方法在工程设计中能比较精确地确定各种弹性边界条件下阶梯柱的临界载荷。     

薄球壳在均布外压与温度耦合作用下的热屈曲研究

从张量方法推导出的轴对称薄球壳屈曲方程出发,推导出在均布外压与温度耦合作用下用位移表示的薄球壳热屈曲方程;应用虚功原理建立薄球壳屈曲最小势能泛函;进一步用Ritz(里兹)法分析了周边简支的半球壳的3种热屈曲问题.得到了: 1) 温度不超过屈曲临界温度值时,均布外压的临界载荷;2) 均布外压载荷为0时,屈曲临界温度值;3) 均布外压载荷不超过临界载荷时,屈曲临界温度值.     

论非线性稳定分析中的临界载荷并建议一组求解临界载荷的变分原理

本文从变分原理的概念出发,给出了一组求解非线性稳定临界载荷的变分公式.从本文的变分公式中可非常方便地得到失稳临界载荷的上限值.     

随机最优控制的二阶必要条件综述

文章介绍作者(含合作)近期在随机最优控制的二阶必要条件方面的工作.首先,在凸控制约束情况下给出经典意义下随机奇异最优控制的逐点型二阶必要条件.其次,利用针状变分获得具有非凸控制约束的Pontryagin最大值原理意义下随机奇异最优控制的逐点型二阶必要条件.最后利用变分分析的工具进一步改进非凸控制约束情形的结果,并将其推广到具有状态约束的情形.     

均衡约束为二阶锥约束广义方程的数学规划问题的二阶充分条件(英)

本文考虑一类均衡约束为二阶锥约束广义方程的数学规划问题. 我们通过一个非光滑映射的方向导数, 给出了临界锥的定义, 并建立它在可行点处的等价形式. 基于此临界锥, 我们提出了均衡约束为二阶锥约束广义方程的数学规划问题的二阶充分性条件, 并且验证了在适当的条件下, M-稳定点处的二阶充分性条件是二阶增长条件成立的充分条件.     

粘弹性方程一种新的分裂正定混合元法

本文用分裂正定混合有限元方法研究二阶粘弹性方程. 首先构造一种新的分裂正定混合变分形式和基于这种分裂正定混合变分形式关于时间的半离散格式, 然后绕开关于空间变量的半离散化格式, 直接从时间半离散出发构造出全离散化的分裂正定混合有限元格式, 并给出这种分裂正定混合有限元解的误差估计. 这种研究思路使得理论论证变得更简单,这是处理二阶粘弹性方程的一种新的尝试.     

拟抛物方程的两类分裂正定混合元方法

本文对拟抛物方程构造两种分裂对称正定混合元方法.通过适当选取变分形式,格式分裂成两个独立对称正定子格式,并且方法不需要验证LBB条件.收敛性分析表明方法关于变量u和引进的变量σ分别具有L 2(Ω)和H(div;Ω)范数意义下的最优收敛阶.最后,通过数值实验验证了方法的有效性.     

非保守力作用下杆的塑性动态稳定性

本文用Liapunov第二方法分析非保守力作用下直杆的塑性动态稳定性.杆处于粘性阻尼介质中,并受到切向均布的随动载荷作用.分析中在应力-应变关系中引入了应变率效应.导出了一个稳定性条件,并求出了临界失稳载荷,讨论了应变率效应对杆的稳定性的影响.     

变截面二维功能梯度微梁的振动和屈曲特性北大核心CSCD

基于修正的偶应力理论和Timoshenko梁理论,应用变分原理建立了变截面二维功能梯度微梁的自由振动和屈曲力学模型.模型中包含金属组分和陶瓷组分的材料内禀特征尺度参数,可以预测微梁力学行为的尺度效应.采用Ritz法给出了任意边界条件下微梁振动频率和临界屈曲载荷的数值解.数值算例表明:微梁厚度减小时,无量纲一阶频率和无量纲临界屈曲载荷增大,尺度效应增强.锥度比对微梁一阶频率的影响与边界条件密切相关,同时,对应厚度和对应宽度锥度比的影响也有明显差异.变截面微尺度梁无量纲一阶频率随着陶瓷和金属的材料内禀特征尺度参数比的增加而增大,且不同边界条件时增大程度不同.厚度方向和轴向功能梯度指数对微梁的一阶频率和屈曲载荷也有显著的影响.     

The Stability of a Flexible Vertical Rod on a Vibrating Support

The classical Kapitsa problem of the inverted flexible pendulum is generalized. We consider a thin homogeneous vertical rod with a free top end and pivoted or rigid attached lower end under the weight of the pendulum’s action and vertical harmonic vibrations of the support. In both cases of attachment, we have stability conditions for the vertical rod position. We take the influence of axial and bending rod vibrations and describe the bending vibrations using the Bernoulli–Euler beam model. The solution is built as a Fourier expansion by eigenfunctions of auxiliary boundary-value problems. As a result, the problem is reduced to the set of ordinary differential equations with periodic coefficients and a small parameter. The asymptotic method of two-scale expansions is used for its solution and to determine the critical level of vibration. The influence of longitudinal waves in the rod essentially decreases the critical load. The single-mode approximation has an acceptable accuracy. With pivoting support at the lower end of the rod, we find the explicit approximate solution. For the rigid attachment, we conduct numerical analysis of the critical level of vibrations depending on the problem parameters.     

Dynamic analysis of short and long journal bearings in laminar and turbulent regimes,application in critical shaft stiffness determination

Linear and non-linear stability of a flexible rotor-bearing system supported on short and long journal bearings is studied for both laminar and turbulent operating conditions. The turbulent pressure distribution and forces are calculated analytically from the modified Reynolds equation based on two turbulent models; Constantinescu's and Ng–Pan–Elrod. Hopf bifurcation theory was utilized to estimate the local stability of periodic solutions near bifurcating operating points. The shaft stiffness was found to play an important role in bifurcating regions on the stable boundaries. It was found that for shafts supported on short journal bearings with shaft stiffness above a critical value, the dangerous subcritical region can be eliminated from a range of operating conditions with high static load. The results presented have been verified by published results in the open literature.     

Stability of a square plate with an elastic border

We study the problem of stability of a square plate subject to uniform compression in both directions. The elastic border of the plate is characterized by the four stiffness coefficients (corresponding to the number of sides) in relation to the angle of rotation. We obtain an approximate formula for computing the critical load. For the cases when the plate has hinge support along the entire border or is rigidly clamped the results determined by this formula practically coincide with the exact solutions. One table. Bibliography: 2 titles. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 22, pp. 46–50, 1991.     

A cylindrical plate with a weakly fixed curvilinear edge made of a transversally isotropic material

Free oscillations and stability under an axial compression of a thin cylindrical plate with a weakly fixed rectilinear edge made of a transversally isotropic material with low stiffness with respect to transverse displacements are considered. The curvilinear edges of the plate are assumed to be hingedly supported. The oscillation frequencies and the critical load for a plate with a free or weakly fixed edge are smaller than those for a shell closed in the circumferential direction. The shapes of oscillations and the forms of stability loss localized near the weakly fixed edge and damped at a distance from it are considered. The Timoshenko-Reissner model is used. Localized forms are analyzed by using a system of equations for Timoshenko-Reissner shallow shells, which is derived for this purpose. The main special feature of this system is that it contains a separate equation describing a solution with large variability. For the example of the stability problem under consideration, the error involved in the system of equations for Timoshenko-Reissner shallow shells is studied. The critical load values obtained with the use of the Kirchhoff-Love and Timoshenko-Reissner models are compared.     

Convergence of a non-stiff boundary integral method for interfacial flows with surface tension

Boundary integral methods to simulate interfacial flows are very sensitive to numerical instabilities. In addition, surface tension introduces nonlinear terms with high order spatial derivatives into the interface dynamics. This makes the spatial discretization even more difficult and, at the same time, imposes a severe time step constraint for stable explicit time integration methods.

A proof of the convergence of a reformulated boundary integral method for two-density fluid interfaces with surface tension is presented. The method is based on a scheme introduced by Hou, Lowengrub and Shelley [ J. Comp. Phys. 114 (1994), pp. 312-338] to remove the high order stability constraint or stiffness. Some numerical filtering is applied carefully at certain places in the discretization to guarantee stability. The key of the proof is to identify the most singular terms of the method and to show, through energy estimates, that these terms balance one another.

The analysis is at a time continuous-space discrete level but a fully discrete case for a simple Hele-Shaw interface is also studied. The time discrete analysis shows that the high order stiffness is removed and also provides an estimate of how the CFL constraint depends on the curvature and regularity of the solution.

The robustness of the method is illustrated with several numerical examples. A numerical simulation of an unstably stratified two-density interfacial flow shows the roll-up of the interface; the computations proceed up to a time where the interface is about to pinch off and trapped bubbles of fluid are formed. The method remains stable even in the full nonlinear regime of motion. Another application of the method shows the process of drop formation in a falling single fluid.

     

One-dimensional dynamic thermoviscoelastic contact with damage

The problem of thermoviscoelastic dynamic contact between a rod and a rigid obstacle, when the material damage is taken into account, is modeled and analyzed. The contact is modeled by the normal compliance condition and the stress-strain constitutive equation is of Kelvin-Voigt type. The damage, which describes the reduction of the load carrying capacity of the rod, evolves because of the opening of microcracks as a result of tension or compression. When the damage reaches a critical value at a point on the rod the material cannot carry any load and the system breaks down. Mathematically, this is expressed by the quenching of the solution. The existence of a local weak solution is established using penalization and a priori estimates.     

Instability of a rotating fluid

The paper is concerned with the problem of stability of a uniformly rotating viscous incompressible self-gravitating liquid bounded by a rotationally symmetric free surface. It is proved that it is unstable if the second variation of the energy functional is not positive. Bibliography: 11 titles.     

Robust stabilizing compensators for flexible structures with collocated controls

For flexible structures with collocated rate and attitude sensors/actuators, we characterize compensator transfer functions which guarantee modal stability even when stiffness/inertia parameters are uncertain. While the compensators are finite-dimensional, the structure models are allowed to be infinite-dimensional (continuum models), with attendant complexity of the notion of stability; thus exponential stability is not possible and the best we can obtain is strong stability. Robustness is interpreted essentially as maintaining stability in the worst case. The conditions require that the compensator transfer functions be positive real and use is made of the Kalman-Yakubovic lemma to characterize them further. The concept of positive realness is shown to be equivalent to dissipativity in infinite dimensions. In particular we show that for a subclass of compensators it is possible to make the system strongly stable as well as dissipative in an appropriate energy norm.This research was supported in part under NASA Grant No. NCC 2-374.     

The stability of the shape of an unattached elastic rod with stiff flanges on its ends

The stability of an unattached column consisting of an elastic rod with stiff flanges on its ends under longitudinal compression is investigated. The load under which the plane of the flange surface is tilted from the plane of the support surface is found. This tilting is accompanied by considerable rotation (reversing) of the flanges and corresponding bending of the rod axis. Abrupt replacement of the rectilinear or bent equilibrium shape by an equilibrium shape that is non-contiguous to it occurs. It is established that the columns behave differently when this a change in the equilibrium shape occurs, depending on the ratio of the length of the rod to the length of the flanges.     

Static and dynamic beam forms of the loss of stability of a long orthotropic cylindrical shell under external pressure

A cylindrical shell with end sections which are closed and supported by hinges, in accordance with the concepts of the rod theory, is considered to be under the action of an omnidirectional external pressure which remains normal to the lateral surface during the deformation process. It is shown that, for such shells, the previously constructed consistent equations of the momentless theory, reduced using the Timoshenko shear model to the one-dimensional equations of the rod theory, describe three forms of loss of stability: (1) static loss of stability, which occures through a bending mode from the action of the total end axial compression force since, under the clamping conditions considered, its non-conservative part cannot perform work on deflections of the axial line; (2) also a static loss of stability but one which occurs through a purely shear mode with the conversion of a cylinder with normal sections into a cylinder with parallel sloping sections and a corresponding critical load which is independent of the length of the shell; (3) dynamic loss of stability which occurs through a bending-shear form and can only be revealed by a dynamic method using an improved shear model.