纤维增强聚合物加固木柱的非线性稳定性分析

考虑加固层中纤维增强聚合物布（FRP布）拉伸与压缩时的不同弹性模量,基于梁大挠度变形假定,首先建立了FRP加固细长木梁大挠度弯曲的一般数学模型,给出了考虑梁弯曲二阶效应的非线性控制方程．其次,研究了FRP布加固细长简支木柱的非线性稳定性问题,得到了FRP加固简支木柱的临界载荷公式．理论证明了其过屈曲解的存在性,并利用摄动法,得到了临界载荷附近过屈曲状态的渐近解析解．进行了参数分析,结果表明:FRP加固层对临界载荷有显著的影响,而对其无量纲过屈曲状态影响较小．     

双层网格圆底扁球壳的非线性稳定问题

从基于等效夹层壳思想的双层网格扁壳,非线性弯曲理论的变分方程出发,利用坐标变换方法和驻值余能原理,导出双层网格圆底扁球壳,在均布压力作用下的轴对称大挠度方程和边界条件．采用修正迭代法,求得了两类边界条件下双层网格圆底扁球壳的非线性载荷-位移关系式和临界屈曲载荷的解析表达式,并讨论了几何参数对临界屈曲载荷的影响．     

复合材料多层板壳大挠度非线性问题的迭代解法

在王震鸣等人提出的各向异性多层扁壳的大挠度方程的基础上,提出了复合材料多层板壳大挠度非线性问题的迭代解法。分析了四边简支的复合材料多层矩形扁壳,与小挠度线性理论解析解及有限元非线性解进行了对比。结果表明,载荷较小并发生小挠度时,所得的大挠度解和小挠度解析解非常接近,载荷较大时,所得解和有限元非线性解非常接近。     

变厚度夹层截顶扁锥壳的非线性稳定性分析

对具有变厚度夹层截顶扁锥壳的非线性稳定问题进行了研究。利用变分原理导出表层为等厚度而夹心为变厚度的夹层截顶扁锥壳的非线性稳定问题的控制方程和边界条件,采用修正迭代法求得了具有双曲型变厚度夹层截顶扁锥壳的非线性稳定性问题的解析解,得到了内边缘与一刚性中心固结而外边缘为可移夹紧固支的变厚度夹层截顶扁锥壳临界屈曲载荷的解析表达式,讨论了几何参数和物理参数对壳体屈曲行为的影响。     

非Lévy型正交各向异性开口圆柱壳屈曲问题的辛叠加解析解

该文基于笔者提出的辛叠加方法得到了经典解法难以直接获得的典型非Lévy型正交各向异性开口圆柱壳屈曲问题的解析解.首先，基于Donnell薄壳理论建立了正交各向异性开口圆柱壳屈曲问题的Hamilton体系控制方程，然后将非Lévy型边界下的原问题拆分为两个子问题，在Hamilton体系下利用分离变量和辛本征展开等数学手段对子问题进行求解，最后基于原问题边界条件，通过子问题解的叠加求得原问题的解析解.数值算例表明，辛叠加解析解与有限元数值解结果吻合良好.同时，定量研究了长度和厚度等参数对屈曲载荷的影响.相比于半逆解法等传统解析方法，辛叠加方法基于严格的数学推导，无需假定解的形式，可以获得更多类似问题的解析解.     

正交各向异性椭圆板的弹性失稳

本文以von Kármán型方程为基础并利用一般分支理论讨论了正交各向异性椭圆板在面内边缘均布压力作用下的弹性失稳.利用Liapunov-Schmidt过程证明了单特征值处分支解的存在性并利用小摄动展开得到了分支解的渐近表达式.最后利用有限单元法计算了正交各向异性椭圆板的临界载荷并进行了板的过屈曲分析,还考察了材料和几何参数对稳定性的影响.     

关于钱氏摄动法的高阶解的计算机求解和收敛性的研究

本文借助于中心受集中载荷圆板小挠度问题的积分方程,获得了摄动参数为中心挠度的任意n阶摄动解的解析式.于是,任意次摄动解的所有待定系数能用计算机求解.因此,获得了相当高阶的摄动解.在此基础上,讨论了钱氏摄动法的渐近性和适用区.     

一维纳米准晶层合梁的非局部振动、屈曲与弯曲研究

基于非局部理论,建立了一维纳米准晶层合简支深梁模型,研究了其自由振动、屈曲行为及其弯曲变形问题.采用伪Stroh型公式,导出了纳米梁的控制方程,并通过传递矩阵法获得简支边界条件下纳米准晶层合梁固有频率、临界屈曲载荷及弯曲变形广义位移和广义应力的精确解.通过数值算例,分析了高跨比、层厚比、叠层顺序及非局部效应对一维纳米准晶层合简支梁固有频率、临界屈曲载荷和弯曲变形的影响.结果表明：固有频率和临界屈曲载荷随着非局部参数增大而减小;外层准晶弹性常数更高时,固有频率和临界屈曲载荷更大;叠层顺序对纳米准晶梁的力学行为有较大影响.所得的精确解可为纳米尺度下梁结构的各种数值方法和实验结果提供参考.     

环形板的屈曲状态

本文研究了内边界固定、外边界可移夹紧并受面内径向均匀压缩推力p作用的环形板的轴对称屈曲状态。首先论述了问题的合理提法,给出了未屈曲板的平凡解的解析公式;得到了在平凡解处线性化问题的特征方程,并对某些参数值给出了前两个特征值(即临界载荷)。然后,在一个适度的假设下,证明了所有的临界载荷都是分支点,并给出了临界载荷附近屈曲解的渐近公式。     

基于非局部应变梯度理论功能梯度纳米板的弯曲和屈曲研究

以纳米机器人等智能器件中的功能梯度纳米板结构为研究对象,基于非局部应变梯度理论,研究了其弯曲和屈曲问题.推导了一般情况下的功能梯度纳米板运动方程,弯曲和屈曲作为其特例可简化而成.分析了非局部尺度参数、材料特征尺度参数、梯度指数、纳米板尺寸等对弯曲挠度和临界屈曲载荷的影响.结果表明:不同高阶连续介质力学理论下的最大挠度都随梯度指数的增大而增大,正方形纳米板挠度较小,且板厚越大,弯曲挠度越小;最大挠度随非局部尺度参数的增大而增大,随材料特征尺度参数的增大而减小.临界屈曲载荷随梯度指数的增大而减小,随板厚、长宽比的增大而增大,随非局部尺度的增大而减小,随材料特征尺度的增大而增大.非局部应变梯度高阶弯曲和屈曲中存在结构软化与硬化机制,两个内特征参数之间具有耦合效应,当非局部尺度大于材料特征尺度时,非局部效应在功能梯度纳米板力学性能中占主导作用;当材料特征尺度大于非局部尺度时,应变梯度效应占主导作用.解析结果还证明了当非局部尺度等于材料特征尺度时,非局部应变梯度理论结果退化为经典结果.     

Fast solution of elliptic control problems

Elliptic control problems with a quadratic cost functional require the solution of a system of two elliptic boundary-value problems. We propose a fast iterative process for the numerical solution of this problem. The method can be applied to very special problems (for example, Poisson equation for a rectangle) as well as to general equations (arbitrary dimensions, general region). Also, nonlinear problems can be treated. The work required is proportional to the work taken by the numerical solution of a single elliptic equation.     

An analytical spectral stiffness method for buckling of rectangular plates on Winkler foundation subject to general boundary conditions

An analytical spectral stiffness method is proposed for the efficient and accurate buckling analysis of rectangular plates on Winkler foundation subject to general boundary conditions (BCs). The method combines the advantages of superposition method, stiffness-based method and the Wittrick–Williams algorithm. First, exact general solutions of the governing differential equation (GDE) of plate buckling considering both elastic foundation and biaxial loading is derived by using a modified Fourier series. The superposition of such general solutions satisfy the GDE exactly and BCs approximately, which guarantees the rapid convergence and high accuracy. Then, based on the exact general solution, the spectral stiffness matrix which relates the coefficients of plate generalized displacement BCs and force BCs is symbolically developed. As a result, arbitrary BCs can be prescribed straightforwardly in the stiffness-based model. As an efficient and reliable solution technique, the Wittrick–Williams algorithm with the J₀ problem resolved is applied to obtain the critical buckling solutions. The accuracy and efficiency of the method are verified by comparing with other methods. Benchmark buckling solutions are provided for plates with all possible boundary conditions. Also, dependence of various factors such as foundation stiffness, load combinations and aspect ratio on the buckling behaviors are investigated.     

弹性直杆中一类动力屈曲问题的渐近解

对于弹性杆受刚性块轴向撞击的动力屈曲问题而言,由于轴向载荷形式较为复杂,问题将归结为关于非线性偏微分方程组解的讨论,至今仍未能得到一个理论上的解析解,为此,讨论了有限长理想弹性直杆的此类动力屈曲问题,采用小参数的摄动展开和变分法,成功地得到了这一问题的一个理论上的近似解,并给出了相应的算例,从中得到了一些有益的结论.     

Buckling of a nonlinear elastic rod

The buckling of a pin-ended slender rod subjected to a horizontal end load is formulated as a nonlinear boundary value problem. The rod material is taken to be governed by constitutive laws which are nonlinear with respect to both bending and compression. The nonlinear boundary value problem is converted to a suitable integral equation to allow the application of bounded operator methods. By treating the integral equation as a bifurcation problem, the branch points (critical values of load) are determined and the existence and form of nontrivial solutions (buckled states) in the neighborhood of the branch points is established. The integral equation also affords a direct attack upon the question of uniqueness of the trivial solution (unbuckled state). It is shown that, under certain conditions on the material properties, only the trivial solution is possible for restricted values of the load. One set of conditions gives uniqueness up to the first branch point.     

On the integral equation method for buckling loads

An initial value method for the integral equation of the column is presented for determining the buckling load of columns. The differential equation of the column is reduced to a Fredholm integral equation. An initial value problem is derived for this integral equation, which is reduced to a set of ordinary differential equations with prescribed initial conditions in order to find the Fredholm resolvent. The singularities of the resolvent occur at the eigenvalues. Integration of the equations proceeds until the integrals become excessively large, indicating that a critical load has been reached. To check this method, numerical results are given for two examples, for which the critical load is well known. One is the Euler load of a simply supported beam, and the other case is the buckling load of a cantilever beam under its own weight. The advantage of this initial value method is that it can be applied easily to solve other nonlinear problems for which the critical loads are unknown. This approach will be illustrated in future papers.     

New results for aggregating integer-valued equations

A variety of results have been given for aggregating integer-valued (diophantine) equations whose variables are restricted to nonnegative integers. In each, integer weights are identified for the equations so that their linear combination yields a single equation with the same solution set of the original system of equations. Because the coefficients of the aggregated equation tend to achieve unwieldy sizes as the number of original equations increases, the goal is to identify weights so these coefficients will lie in a range as limited as possible. We give theorems which separately and in combination provide new methods for aggregating general integer-valued equations. Our results include formulations that do not require linearity of the original system, or nonnegativity of component variables. We also demonstrate that our theorems yield as special cases earlier results (analytical formulae) conjectured to yield the smallest possible weights for less general domains. As another application, the presented results were used to develop a highly efficient approach for the integer knapsack problem. Empirical outcomes show that the developed solution procedure is significantly superior to advanced branch and bound methods (previously established to be the most efficient knapsack solution procedures).     

Creep stability of polymer rods

The author examines the problem of the buckling of a hinger rod of rigid homogeneous polymer material under constant load. The results of an experimental investigation are presented. The theoretical calculations are based on the nonlinear generalized Maxwell equation. A numerical solution has been obtained on a computer. The results of this solution are compared with the experimental data. An analytic solution that makes it possible to estimate certain limiting values is obtained for the linearized equation.Mekhanika Polimerov, Vol. 4, No. 1, pp. 145–150, 1968     

对称迭层矩形板的平面应力分析

常用的对称迭层板为各向异性板．根据平面应力问题的基本方程精确地用应力函数解法求得了各向异性板的一般解析解．推导出平面内应力和位移的一般公式,其中积分常数由边界条件来决定．一般解包括三角函数和双曲函数组成的解,它能满足4个边为任意边界条件的问题．还有代数多项式解,它能满足4个角的边界条件．因此一般解可用以求解任意边界条件下的平面应力问题．以4边承受均匀法向和切向载荷以及非均匀法向载荷的对称迭层方板为例,进行了计算和分析．     

Problem of an elastic half-plane with a crack emerging orthogonally at the boundary

The problem of the half-plane, in which a finite crack emerges orthogonally at the boundary, is studied. On the edges of the crack a self-balancing load is applied. A detailed investigation is carried out for an integral equation with respect to the unknown derivative of the displacement jump, to which the problem can be reduced. The exact solution of the integral equation is constructed with the aid of the Mellin transform and the Riemann boundary value problem for the halfplane. The asymptotic behavior of the solution at both ends of the crack is elucidated. First the asymptotic behavior of the solution at the point of emergence of the crack is obtained and the dependence of this asymptotic behavior on the type of the load is established. For a special form of the load one obtains a simple expression of the stress intensity coefficient. In the case of a general load, the asymptotic behavior is used for the construction of an effective approximate solution on the basis of the method of orthogonal polynomials. As a result, the problem reduces to an infinite algebraic system, solvable by the reduction method.Translated from Dinamicheskie Sistemy, No. 4, pp. 45–51, 1985.     

非均匀变截面梁动力响应的一般解

本文利用精确解析法[1]给出非均匀变截面梁在任意谐振荷载和边界条件下的动力响应的一般解.问题最后归结为求解一个非正定微分方程.对于这一问题用一般变分法求解失败,利用本文的方法,这个问题的一般解可以写成解析的形式.因此对优化特别方便.本文给出收敛性证明.文末给出的算例表明本文的方法可获得满意的结果.