薄板弯曲理论中一个新的广义变分原理及其交替max,min性质

本文在薄板弯曲的经典理论中提出一个以挠度、横向剪应变、曲率、弯矩、横向剪力、边界上的未知反力为自变函数的新的变分原理,在对自变函数预加不严厉的约束后,新泛函具有一连串交替max,min的性质。     

用广义变分原理和重调函数级数求解圆板的混合边值问题

1.前言弹性力学中的广义变分原理是一般性的变分原理.在这一原理中,自变函数可以任意选取,而自变函数问的相互关系(几何、物理、平衡三个方面)和边界条件由泛函的驻值来保证.这一变分原理广泛应用于有限元方法和弹性力学数值分析等问题中.利用广义变分原理求解弹性薄板弯曲问题的开创性工作可见文献[3].然而,在广义变分原理的具体应用方面,仍然存在着许多问题.例如,在弹性力学空间问题中,有位移,应变和应力等15个自变函数,人们还不太清楚怎样具体选择这些自变函数为好.又如,若选择的自变函数和受力物体的真实变形状态不适应时,此时广义变分原理不能导致近似解,有时甚至会得到错误的解答.     

荷载横向变位下变截面连续箱梁剪滞分析的有限梁段法

为研究荷载横向作用位置变化对箱梁剪滞效应的影响,对箱梁顶、底板、悬臂板分别设置了不同的剪滞纵向位移差函数;假定纵向翘曲位移沿横向分布为k次抛物线,并考虑剪滞和剪切双重效应的影响,通过能量变分法推导出了荷载横向变位时梁段单元的平衡控制微分方程组及其闭合解;提出了能对工程中常见的变截面连续箱梁剪滞效应进行分析的有限梁段法。该方法计算结果与有限元模型、已有模型试验结果的最大误差在5.95%~9.74%之间,两种工况下计算结果的叠加与有限元结果相对误差在0.07%~19.18%之间,均吻合良好,说明将基于有限梁段法的剪滞效应变分解和叠加原理用于求解复杂力状态下的剪滞效应是可行的。剪滞翘曲位移横向分布函数精度选择的研究结果表明:均布荷载分别作用于腹板顶部、顶板中心时,翘曲位移横向分布函数宜分别选用三次、二次抛物线。     

略论Hellinger-Reissner和胡海昌-鹫津久一郎两种广义变分原理的联系

Hellinger和Reissner先后于1914年和1950年提出了弹性力学中的一种广义变分原理,其中位移和应力看作是二类独立的自变函数.后来这种变分原理常叫做Hellinger-Reissner变分原理.本人和鹫津久一郎先后于1954年和1955年提出了另一种广义变分原理,其中位移、应变和应力三类变量都看作是独立的自变函数.后来这种变分原     

广义变分原理与无条件变分原理

在讨论变分原理和有限元素法的许多文献中,常常出现“广义变分原理”和“无条件变分原理”这样两个名词。在不同的文献中,例如在专著[1]～[3]中,这两个名词有不同的含义,给读者带来许多不便。本文对此进行讨论,并希望能求得一致的理解和用法。有无条件是一个古老的数学概念,指的是变分式中的自变函数事前要不要满足某种条件(不包括函数的连续性、可导性、可积性等一系列定性的要求,有时甚至不包括边界条件,见文献[3]第21页)。如果在某变分式中,自变函数事前不用满足什么条件,那末     

粘弹层合板的稳态振动和层间应力

应用混合分层理论和Ressiner混合变分原理，在板厚方向取二次位移插值函数和三次、四次横向应力插值函数推导出粘弹层合板的动力学方程，得出简支粘弹层合板稳态振动的解。不仅得出与三层弹性板精确的自振频率吻合良好的解，而且对于粘弹层合板，所计算的自振频率和结构损耗因子也与三维结果吻合较好。计算了自由阻尼层合板对应的低阶法向位移响应幅值和层问横向应力的幅值。结果表明，较高的层间横向正应力是低频稳态振动中引起粘弹层合板分层破坏的主要因素，采用适当模量和厚度的粘弹性材料将有效地降低粘弹层合板的层间横向正应力的幅值。     

面向变厚度柔性轧制工艺的帽型梁横向冲击吸能优化设计

作为汽车主要吸能构件的帽型梁的吸能提升设计是备受关注的问题.研究表明,通过优化薄壁结构的厚度可有效提升吸能性能,但复杂的厚度分布造成制造困难.针对可实现厚度调控的工艺,发展易制造的结构设计方法极为必要.本文基于变厚度柔性轧制工艺(variable gauge rolling, VGR)可实现厚度调控的特点,发展建立帽型梁横向冲击吸能优化设计方法.基于变厚度柔性轧制工艺生产的柔性轧制板(tailor rolled blanks,TRB)的特点,将受横向冲击的帽型薄壁梁设计成沿轴线分段变厚度、分段间设梯度过渡段的结构形式,通过调整各段厚度、分段位置和过渡层梯度变化规律,实现性能的优化.以应变能密度分布均匀为优化准则、基于混合元胞自动机(hybird cellular automata, HCA)方法构建优化模型和求解方法,并在迭代过程中施加满足轧制约束的过滤函数,使结构满足轧制工艺要求.其中,轧制约束的过滤函数由粒子群算法自动寻找.基于本文方法,具体设计了柔性轧制帽型梁横向冲击吸能最优的分段位置、各段厚度及过渡段厚度的梯度过渡方式,设计结果验证了方法的有效性.     

面向变厚度柔性轧制工艺的帽型梁横向冲击吸能优化设计

作为汽车主要吸能构件的帽型梁的吸能提升设计是备受关注的问题.研究表明,通过优化薄壁结构的厚度可有效提升吸能性能,但复杂的厚度分布造成制造困难.针对可实现厚度调控的工艺,发展易制造的结构设计方法极为必要.本文基于变厚度柔性轧制工艺(variable gauge rolling, VGR)可实现厚度调控的特点,发展建立帽型梁横向冲击吸能优化设计方法.基于变厚度柔性轧制工艺生产的柔性轧制板(tailor rolled blanks, TRB)的特点,将受横向冲击的帽型薄壁梁设计成沿轴线分段变厚度、分段间设梯度过渡段的结构形式,通过调整各段厚度、分段位置和过渡层梯度变化规律,实现性能的优化.以应变能密度分布均匀为优化准则、基于混合元胞自动机(hybird cellular automata, HCA)方法构建优化模型和求解方法,并在迭代过程中施加满足轧制约束的过滤函数,使结构满足轧制工艺要求.其中,轧制约束的过滤函数由粒子群算法自动寻找.基于本文方法,具体设计了柔性轧制帽型梁横向冲击吸能最优的分段位置、各段厚度及过渡段厚度的梯度过渡方式,设计结果验证了方法的有效性.      

有限变形弹性杆中的非线性波及周期解

利用Ham ilton变分原理,导出了计及有限变形和横向Possion效应的弹性杆中非线性纵向波动方程.利用Jacob i椭圆正弦函数展开和第三类Jacob i椭圆函数展开法,对该方程和截断的非线性方程进行求解,得到了非线性波动方程的两类准确周期解及相应的孤波解和冲击波解,讨论了这些解存在的必要条件.     

功能梯度变曲率曲梁的几何非线性模型及其数值解

基于弹性曲梁平面问题的精确几何非线性理论,建立了功能梯度变曲率曲梁在机械和热载荷共同作用下的无量纲控制方程和边界条件,其中基本未知量均被表示为变形前的轴线坐标的函数。以椭圆弧曲梁为例,采用打靶法求解非线性常微分方程的两点边值问题,获得了两端固定功能梯度椭圆弧曲梁在横向非均匀升温下的热弯曲变形数值解,分析了材料梯度指数、温度参数、结构几何参数等对曲梁受力及变形的影响。     

Representation theorem for linear,isotropic tensor functions of a skew argument

The present paper offers two proofs of the representationtheorem for linear,isotropic tensor functions of a skew ar-gument.The first proof is new.The second one is basicallyalong the lines of reasoning exploited in[1].but more con-cise,and it corrects some errors committed in[1].     

考虑横向剪切变形的扁球壳在集中力作用下的分析

采用适于夹层壳的直线假设扁壳理论，应用三角级数法，导出了扁球壳齐次方程的解析解。进而分析了在顶点作用法向集中力和在偏心集中力作用下的解。计算了在偏心集中力作用下带孔球壳的位移和应力，并将结果与经典理论的结果进行了比较分析，结果表明，在集中力作用处和孔边处两种理论结果明显不同。     

A two variable refined plate theory for orthotropic plate analysis

A new theory, which involves only two unknown functions and yet takes into account shear deformations, is presented for orthotropic plate analysis. Unlike any other theory, the theory presented gives rise to only two governing equations, which are completely uncoupled for static analysis, and are only inertially coupled (i.e., no elastic coupling at all) for dynamic analysis. Number of unknown functions involved is only two, as against three in case of simple shear deformation theories of Mindlin and Reissner. The theory presented is variationally consistent, has strong similarity with classical plate theory in many aspects, does not require shear correction factor, gives rise to transverse shear stress variation such that the transverse shear stresses vary parabolically across the thickness satisfying shear stress free surface conditions. Well studied examples, available in literature, are solved to validate the theory. The results obtained for plate with various thickness ratios using the theory are not only substantially more accurate than those obtained using the classical plate theory, but are almost comparable to those obtained using higher order theories having more number of unknown functions.     

Reduced-order models for microelectromechanical rectangular and circular plates incorporating the Casimir force

We consider the von Kármán nonlinearity and the Casimir force to develop reduced-order models for prestressed clamped rectangular and circular electrostatically actuated microplates. Reduced-order models are derived by taking flexural vibration mode shapes as basis functions for the transverse displacement. The in-plane displacement vector is decomposed as the sum of displacements for irrotational and isochoric waves in a two-dimensional medium. Each of these two displacement vector fields satisfies an eigenvalue problem analogous to that of transverse vibrations of a linear elastic membrane. Basis functions for the transverse and the in-plane displacements are related by using the nonlinear equation governing the plate in-plane motion. The reduced-order model is derived from the equation yielding the transverse deflection of a point. For static deformations of a plate, the pull-in parameters are found by using the displacement iteration pull-in extraction method. Reduced-order models are also used to study linear vibrations about a predeformed configuration. It is found that 9 basis functions for a rectangular plate give a converged solution, while 3 basis functions give pull-in parameters with an error of at most 4%. For a circular plate, 3 basis functions give a converged solution while the pull-in parameters computed with 2 basis functions have an error of at most 3%. The value of the Casimir force at the onset of pull-in instability is used to compute device size that can be safely fabricated.     

Representations of linear,isotropic tensor functions of an asymmetric argument

This paper offers a mathematical demonstration of the representation theorems for linear, isotropic scalar- and tensor-valued tensor functions of an asymmetric argument.     

Vibration and buckling of laminated sandwich plates having interfacial imperfections

The vibration and buckling characteristics of sandwich plates having laminated stiff layers are studied for different degrees of imperfections at the layer interfaces using a refined plate theory. With this plate theory, the through thickness variation of transverse shear stresses is represented by piece-wise parabolic functions where the continuity of these stresses is satisfied at the layer interfaces by taking jumps in the transverse shear strains at the interfaces. The transverse shear stresses free condition at the plate top and bottom surfaces is also satisfied. The inter-laminar imperfections are represented by in-plane displacement jumps at the layer interfaces and characterized by a linear spring layer model. It is quite interesting to note that this plate model having all these refined features requires unknowns only at the reference plane. To have generality in the analysis, finite element technique is adopted and it is carried out with a new triangular element developed for this purpose, as any existing element cannot model this plate model. As there is no published result on imperfect sandwich plates, the problems of perfect sandwich plates and imperfect ordinary laminates are used for validation.     

Shear correction factors and an energy-consistent beam theory

Presented here is a new derivation of shear correction factors for isotropic beams by matching the exact shear stress resultants and shear strain energy with those of the equivalent first-order shear deformation theory. Moreover, a new method of deriving in-plane and shear warping functions from available elasticity solutions is shown. The derived exact warping functions can be used to check the accuracy of a two-dimensional sectional finite-element analysis of central solutions. The physical meaning of a shear correction factor is shown to be the ratio of the geometric average to the energy average of the transverse shear strain on a cross section. Examples are shown for circular and rectangular cross sections, and the obtained shear correction factors are compared with those of Cowper (1966) . The energy-averaged shear representative is also used to derive Timoshenkos beam theory.     

Elastodynamic analysis of low tension cables using a new curved beam element

In this paper, we address and overcome the difficulties associated with the use of the classic cable theory to treat low tension cables by developing a new three-noded locking-free nonlinear curved beam element. Based upon nonlinear generalized curved beam theory, large deformations and rotations in the new element are formulated in terms of Updated Lagrangian framework. Consistently coupled polynomial displacement fields are used to satisfy the membrane locking-free condition and the requirement of being able to recover the inextensible bending modes. Quintic transverse displacement interpolation functions are used to represent the bending deformation of the beam, while the axial and torsional displacement fields are derived by integration of the presumably linear membrane and torsional shear strain fields, which are coupled with the transverse displacement fields. Numerical results are presented to demonstrate the superior accuracy and the high convergence rate of the newly developed curved beam element. The stability and accuracy of the new element are further validated by experiments of an instrumented free-swinging steel cable experiencing slack and low tension. Good agreements in cable position and tension are observed between the experimental results and the finite element predictions.