柔性圆环壳在子午面内整体弯曲的复变量方程及细环壳的一般解

在З.Л.Аксельрад(E.L.Axelrad)非轴对称载荷下柔性旋转壳线性方程的基础上,导出了圆环壳在子午面内整体弯曲的复变量方程和相应的细环壳方程.该方程可与钱伟长给出了一般解的В.В.Новожилов(V.V.Novozhilov)轴对称环壳方程相类比.通过类比,给出了细环壳在子午面内整体弯曲的一般解.所给出的解可以用来计算波纹管整弯曲的应力和端面位移.     

弹性圆柱壳扭转屈曲研究

本文给出两端固支的弹性圆柱壳扭转屈曲实验与理论计算结果.实验发现,对于较长的壳,其屈曲后的变形并不占据整个壳体的长度.另外在计算中仅考虑壳体的法向边界条件,而不考虑其周向和轴向边界条件,结果和Yamaki精确解以及本文实验结果相符较好,说明周向和轴向边界条件对圆柱壳的扭转屈曲影响较小.     

旋转壳在子午面内整体弯曲的几何非线性摄动有限元肖及在波纹管计算中的应用（Ⅰ）

为切实有效地计算波纹管，建立了旋转壳在子午面内整体弯曲几何非线性问题的摄动有限元法。以结构环向应变的均方根为摄动小参数，将有限元节点位移列式和节点力列式直接展开。通过摄动小参数将非线性有限元的载荷分级和迭代过程有机地统一起来，即载荷的分级是有约束的，每一级载荷增量和所对应的位移增量之间的关系是已知的，每一级的计算一步到位。为叙述方便并具实用性，将旋转壳用截锥壳单元进行离散。位移分量和载荷分量沿环向按Fourier级数展开，沿子午线用多项式插值，端面弯矩和横向力化成载荷分量离散到节点上。采用Sanders中小转角非线性几何方程和各向同性广义Hooke定律，对多层材料叠合而成的旋转壳按各层薄膜应变、弯曲应变、扭转应变相等的原则进行处理，该方法能方便有效地计算单层和多层波纹管整体纯弯曲、横向弯曲的几何非线性问题。并为有限元处理非线性问题提供了一条新途径。     

圆柱壳开孔的应力集中──非圆孔问题的一般解

本文从Donnell型圆柱壳的基本方程出发,利用复变函数方法和保角映射技术,对圆柱壳开非圆形孔的问题进行了研究.首先给出了逼近具有非圆形孔的圆柱壳开孔问题一般解的完备函数序列,构造出了问题的一般解;其次利用有关圆柱壳开小孔的假设概念,给出了圆柱壳开非圆孔时边界条件的一般表达式.进而利用正交函数展开的方法,将待解的问题归结为一组无穷代数方程组的求解问题,并进行直接求解.在本文最后,对圆柱壳开圆孔.椭圆孔附近的应力集中问题进行了数值计算,给出了分析结果.     

旋转壳在子午面内整体弯曲的几何非线性摄动有限元法及在波纹管计算中的应用(Ⅰ)

为切实有效地计算波纹管,建立了旋转壳在子午面内整体弯曲几何非线性问题的摄动有限元法。以结构环向应变的均方根为摄动小参数,将有限元节点位移列式和节点力列式直接展开。通过摄动小参数将非线性有限元的载荷分级和迭代过程有机地统一起来,即载荷的分级是有约束的,每一级载荷增量和所对应的位移增量之间的关系是已知的,每一级的计算一步到位。为叙述方便并具实用性,将旋转壳用截锥壳单元进行离散。位移分量和载荷分量沿环向按Fourier级数展开,沿子午线用多项式插值,端面弯矩和横向力化成载荷分量离散到节点上。采用Sanders中小转角非线性几何方程和各向同性广义Hooke定律。对多层材料叠合而成的旋转壳按各层薄膜应变、弯曲应变、扭转应变相等的原则进行处理,该方法能方便有效地计算单层和多层波纹管整体纯弯曲、横向弯曲的几何非线性问题。并为有限元处理非线性问题提供了一条新途径。     

旋转薄壳自由振动中的转点问题及其奇异摄动解

在用渐近方法求解任意旋转薄壳(圆柱壳和球壳除外)自由振动的微分方程组时,在一定的参数范围内,存在转(向)点问题。其中,对于存在唯一简单转点的情况(此时,该转点具有枝点奇性),至今未获解决。本文解决了这一问题。文中仅讨论截顶旋转壳。此时,转点是唯一的奇点。最后,对于两端固定的任意旋转薄壳,给出了频率方程。     

具有不等厚表层的夹层旋转壳问题

本文探讨了具有各向同性夹心及不等厚表层组成的三层旋转壳的小挠度问题,采用广义变分法,在考虑表层抗弯刚度的情况下,导出了基本方程,其特例即为文献[3],[4]的结果. 根据所导出的基本方程,本文求解对称加载的夹层旋转壳问题,将其归结为只含一个广义位移的六阶微分方程. 本文的理论可应用于复合装甲及其它工程设计.     

弹性介质中受轴向冲击载荷作用的圆柱壳的屈曲临界载荷计算分析

建立并求解了弹性介质中圆柱壳的径向位移控制方程,考虑边界条件及相容条件,得到了应力波传播及反射过程中圆柱壳的动力屈曲分叉条件.通过计算得到了不同时间段屈曲临界载荷与应力波波阵面到达圆柱壳的位置、弹性介质的刚度、壳体未嵌入弹性介质部分的长度与总长之比的关系.数值计算结果表明,弹性介质中的圆柱壳发生轴对称屈曲和非轴对称屈曲趋势一致;嵌入弹性介质部分越深、弹性介质刚度越大圆柱壳越难屈曲;屈曲临界载荷随着弹性介质刚度的增大经历了增长缓慢、增长迅速以及增长较慢3个阶段;应力波反射前波阵面通过分界面后,屈曲仅发生在应力波传播区域,反射波波阵面通过分界面前,临界载荷较小时屈曲先发生在反射端部,随着轴向阶数增大,屈曲覆盖整个圆柱壳区域,反射波波阵面通过分界面后,壳体发生的屈曲始终覆盖整个圆柱壳区域.     

球壳和柱壳振动中一类方程组的求解

本文利用矩阵型式的Ｆｒｏｂｅｎｉｕｓ级数方法求解了球壳和柱壳振动中一类含正则奇点的常微分方程组，考虑了指标方程根的相互关系，从而全面地得到了不同情况下解的表达形式，为解析求解工作奠定了基础。     

自然弯扭梁广义翘曲坐标的求解

提出了自然弯扭梁受复杂载荷作用时静力分析的一种理论方法,重点在于对控制方程的求解,其中考虑了与扭转有关的翘曲变形和横向剪切变形的影响．在特殊的情况下,可以比较容易地得到这些方程的解答,包括各种内力、应力、应变和位移的计算．算例给出了平面曲梁受水平和垂直分布载荷作用时广义翘曲坐标的求解方法．计算结果表明,求得的应力和位移的理论值和三维有限元结果非常接近．此外,该理论不限于具有双对称横截面的自然弯扭梁,同样可推广至具有一般横截面形状的情况．     

正交各向异性体梁弯曲的弹性理论

本文由文献[1]横观各向同性板的弯曲弹性理论关于二维问题的特例,通过比拟,得到了正交各向异性梁弯曲的弹性理论,文中给出了求解正交各向异性梁弯曲问题的一种方法.提出了一种新的深梁理论,并指出了考虑横向剪切变形影响的Reissner理论对于应力分量的近似程度较差.     

新型空间薄壁梁单元

基于Timoshenko梁理论和Vlasov薄壁杆件约束扭转理论,建立了具有内部结点的新型空间薄壁截面梁单元．通过对弯曲转角和翘曲角采取独立插值的方法,考虑了横向剪切变形,扭转剪切变形及其耦合作用,弯曲变形和扭转变形的耦合以及二次剪应力等因素影响,由Hellinger-Reissner广义变分原理,推得单元刚度矩阵．算例表明所建模型具有良好的精度,可用于空间薄壁杆系结构的有限元分析．     

Multi-scale modeling of beam-like structures: A new boundary condition concept for the RVE

In this work a coupled two-scale beam model using Timoshenko beam elements [1] with finite displacements on the macro scale and fully non-linear 3D brick elements on the micro scale is proposed. The calculation is carried out with the so-called FE² concept. To achieve the coupling between the beam and the brick elements, the algorithm from [2] is adapted. Within the degenerated concept of the Timoshenko beam, the introduction of a pure shear deformation leads to significant problems concerning the equilibrium condition on the micro scale. Applying this deformation mode on the RVE with periodic boundary conditions results in a rigid body rotation. Using linear displacement boundary conditions instead, the wrapping deformation is suppressed on the boundary, leading to a length dependency in the torsional deformation mode. In addition, the shear forces introduce a bending moment, which depends on the length of the RVE and adds spurious normal stresses and a length dependency of the shear stiffness. To overcome these problems, periodic boundary conditions are applied and the displacement assumptions are modified such that the shear deformation is achieved with force pairs on both ends of the RVE. The resulting model leads to length independent results in tension, bending and torsion and a domain which is able to produce a pure shear stress state. Consequently, only this domain of the model should be homogenized which can be accomplished by modifying the variations in the algorithm [2]. The concept is validated by simple linear and non-linear test problems. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

考虑屋盖整体空间工作时的单层偏心厂房的平扭耦联地震反应分析

本文研究了单层偏心厂房的平动-扭转耦联地震反应.建立了考虑屋盖整体空间工作时的地震反应运动方程,并进行了分析,讨论.指出在包含有地震动转动分量作用时,地震动各运动分量独立作用的迭加原理和动力反应分析的振型分解法已不再适用.本文在不计地震转动分量的前提下,对单层偏心厂房的平扭耦联地震反应进行了大量的计算及对比分析,结论是:考虑屋盖变形后的计算结果能更好地解释地震破坏现象,并指出了不考虑屋盖变形的适用范围.     

Natural damped vibrations of anisotropic box beams of polymer composite materials. 2. Numerical experiments

The influence of the orientation of reinforcing fibers on the natural frequencies and mechanical loss coefficient of coupled vibrations of unsupported symmetric and asymmetric box beams, as evaluated in numerical experiments, is discussed. The calculations were performed under the assumption that the real parts of the complex moduli and mechanical loss coefficient are frequency-independent. Vibration modes were identified by their surface shapes. The boundaries of the regions of mutual transformation of interacting vibration modes were determined by the joint analysis of the dependences of the coupled and partial eigenfrequencies and the mechanical loss coefficients on the orientation angle of reinforcing fibers. It is established that vibrations of a symmetric box beam give rise to two primary interactions: bending–torsional and longitudinal–shear ones, which are united into a unique longitudinal–bending–torsional–shear interaction by the secondary interaction caused by transverse shear strains. Vibrations of an asymmetric box beam give rise to longitudinal–torsional and bending–bending (in two mutually orthogonal planes) interactions. It is shown that in a number of cases variation in the orientation angle of reinforcing fibers is accompanied with a mutual transformation of coupled vibration modes. If the differential equations for natural vibrations involve odd-order derivatives with respect to the spatial variable (a symmetric beam and the bending–bending interaction of an asymmetric beam), then, with variation in the orientation angle of reinforcing fibers, the mutual transformation of coupled vibration modes proceeds. If the differential equations for natural vibrations involve only even-order derivatives (the longitudinal–torsional interaction of an asymmetric beam), no mutual transformation of coupled vibration modes occurs.     

Selective finite element refinement in torsional problems based on the membrane analogy

This work presents a selective finite element refinement strategy based on the h-refinement type, in the context of a posteriori error estimates considerations (error computed after the application of the proposed refining scheme), based on a graphical procedure to determine progressively better estimates for the maximum shearing stress in prismatic torsional members. It is structured in an integrated FORTRAN code and DELPHI based environment to refine an initial arbitrary finite element mesh. The proposed procedure is founded on the membrane analogy that exists between membrane deflections and the torsion problem in the sense that the location of the membrane largest gradient drives the refining procedure. It is shown that multiple level application of the proposed method to two members with different cross sectional geometries with known analytic solutions leads to progressively more accurate estimates (< 1.0% error in most cases) for the maximum shearing stresses calculations. Finally, the proposed method is applied to the torsional analysis of an L section member, showing that for this practical case the procedure results in a very accurate calculation as well.     

热荷载作用下薄板的优化设计

板壳结构是一大类广泛使用的结构元件.在热荷载作用下,当热膨胀受到约束时,板壳结构产生内力及挠度,严重时影响结构的正常服役.由于热荷载的特殊性,简单地均匀加大板壳结构的厚度并不能有效地减少热变形和热应力,热结构设计因此特别困难.该文研究在给定材料体积的条件下,通过优化板壳结构的厚度分布来减少弹性薄板结构在热载荷下的变形.以结构的变形能为优化目标,在给定材料体积的条件下,建立了设计板壳结构厚度分布的优化问题列式,并采用变分法,推导出优化准则,给出了修改厚度的迭代公式.应用商用有限元软件的热结构分析功能,对程序进行二次开发,从而实现该优化算法.算例结果表明,采用该方法优化弹性薄板的厚度分布,可以大幅度地减小结构热变形,是一种有效的热结构设计方法.     

A two-variable first-order shear deformation theory considering in-plane rotation for bending,buckling and free vibration analyses of isotropic plates

This paper presents a two-variable first-order shear deformation theory considering in-plane rotation for bending, buckling and free vibration analyses of isotropic plates. In recent studies, a simple first-order shear deformation theory (S-FSDT) was developed and extended. It has only two variables by separating the deflection into bending and shear parts while the conventional first-order shear deformation theory (FSDT) has three variables. However, the S-FSDT provides incorrect predictions for the transverse shear forces on the insides and the twisting moments at the boundaries except simply supported plates since it does not consider in-plane rotation. The present theory also has two variables but considers in-plane rotation such that it is able to correctly predict the responses of plates with any boundary conditions. Analytical solutions are obtained for rectangular plates with two opposite edges that are simply supported, with the other edges having arbitrary boundary conditions. Numerical results of deflections, stress resultants, buckling loads and natural frequencies are presented with the FSDT, the S-FSDT and the present theory. Comparative studies demonstrate the effects of in-plane rotation and the accuracy of the present theory in predicting the bending, buckling and free vibration responses of isotropic plates.     

A mathematical analysis of the elasto-plastic plane stress problem of a power-law material

In this paper a generic study on the plane stress problem ofa power-law material undergoing infinitesimal deformations iscarried out, and a general solution for the stress and strainfields is derived using a stress function method and analyticfunction theory. Hencky's deformation theory and von Mises'yield criterion are used, and a differential transformationis invoked in the analysis. As an example, the closed-form solutionof the pure bending problem of a thin beam of power-law materialis obtained by applying the general solution directly.     

A homotopy analysis solution to large deformation of beams under static arbitrary distributed load

In this article, homotopy analysis method (HAM) is employed to investigate non-linear large deformation of Euler–Bernoulli beams subjected to an arbitrary distributed load. Constitutive equations of the problem are obtained. It is assumed that the length of the beam remains constant after applying external loads. Different auxiliary parameters and functions of the HAM and the extra auxiliary parameter, which is applied to initial guess of the solution, are employed to procure better convergence rate of the solution. The results of the solution are obtained for two different examples including constant cross sectional beam subjected to constant distributed load and periodic distributed load. Special base functions, orthogonal polynomials e.g. Chebyshev expansion, are employed as a tool to improve the convergence of the solution. The general solution, presented in this paper, can be used to attain the solution of the beam under arbitrary distributed load and flexural stiffness. Ultimately, it is shown that small deformation theory overestimates different quantities such as bending moment, shear force, etc. for large deflection of the beams in comparison with large deformation theory. Finally, it is concluded that solution of small deformation theory is far from reality for large deflection of straight Euler–Bernoulli beams.