Piecewise constrained optimization harmonic balance method for predicting the limit cycle oscillations of an airfoil with various nonlinear structures

A method is proposed to calculate the periodic solutions of piecewise nonlinear systems. The method is based on analytical derivation of nonlinear multi-harmonic equations of motion. Since periodic variations of nonlinear forces are characterized by different states, the vibration cycle is broken into sequential transition intervals according to the instant sets of state transitions. Analytical formulations of the harmonic coefficients of the nonlinear forces and its derivatives with respect to the harmonic coefficients of displacements are developed. Sensitivities of the harmonic coefficients of periodic solutions are determined for constructing explicit expressions for vibration amplitude levels as a function of structural parameters. Numerical investigations of the limit cycle oscillations and its sensitivities of an airfoil with different piecewise nonlinearities have been performed. The results show that the developed method is capable of determining the periodic solutions and its sensitivities with respect to the structural parameters. In order to guarantee time continuity of the nonlinear force, for the hysteresis model it is not right to track the periodic solutions by using the preload or freeplay as the continuation parameters.     

Mathematical modeling of brush seals

A computational fluid dynamics (CFD)-based model of brush seals has been developed and tested against other workers' experimental data. In the model, the brush is treated as an axisymmetric, anisotropic porous region with nonlinear resistance coefficients. The resistance coefficients are chosen through calibration against measurements. The CFD model gives predictions of flow rate, pressure distribution, velocity field, and bending forces on the bristles. The bristle forces are used in a separate calculation to estimate bristle bending and reaction forces on the shaft and backing plate. Bending in both the axial direction and the orthogonal plane are considered.     

Bending of Bars under a Follower Load

Exact solutions of the problem of nonlinear bending of a thin bar under a point follower load are given. The problem is studied for an arbitrary follower angle and the particular cases of axial and transverse follower forces are considered. The solutions are written in unified parametric form and expressed in terms of Jacobi elliptic functions.     

Refined geometrically nonlinear formulation of a thin-shell triangular finite element

A refined geometrically nonlinear formulation of a thin-shell finite element based on the Kirchhoff-Love hypotheses is considered. Strain relations, which adequately describe the deformation of the element with finite bending of its middle surface, are obtained by integrating the differential equation of a planar curve. For a triangular element with 15 degrees of freedom, a cost-effective algorithm is developed for calculating the coefficients of the first and second variations of the strain energy, which are used to formulate the conditions of equilibrium and stability of the discrete model of the shell. Accuracy and convergence of the finite-element solutions are studied using test problems of nonlinear deformation of elastic plates and shells. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 160–172, September–October, 2007.     

FINITE ELEMENT DISPLACEMENT PERTURBATION METHOD FOR GEOMETRIC NONLINEAR BEHAVIORS OF SHELLS OF REVOLUTION OVERALL BEDING IN A MERIDIONAL PLANE AND APPLICATION TO BELLOW (Ⅱ)

The finite-element-displacement-perturbation method (FEDPM)for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes (Ⅰ) was employed to calculate the stress distributions and the stiffness of the bellows. Firstly, by applying the first-order perturbation solution (the linear solution)of the FEDPM to the bellows, the obtained results were compared with those of the general solution and the initial parameter integration solution proposed by the present authors earlier, as well as of the experiments and the FEA by others.It is shown that the FEDPM is with good precision and reliability, and as it was pointed out in (Ⅰ) the abrupt changes of the meridian curvature of bellows would not affect the use of the usual straight element. Then the nonlinear behaviors of the bellows were discussed. As expected, the nonlinear effects mainly come from the bellows ring plate,and the wider the ring plate is, the stronger the nonlinear effects are. Contrarily, the vanishing of the ring plate, like the C-shaped bellows, the nonlinear effects almost vanish. In addition, when the pure bending moments act on the bellows, each convolution has the same stress distributions calculated by the linear solution and other linear theories, but by the present nonlinear solution they vary with respect to the convolutions of the bellows. Yet for most bellows, the linear solutions are valid in practice.     

Nonlinear bending of size-dependent circular microplates based on the modified couple stress theory

The present study proposes a nonclassical Kirchhoff plate model for the axisymmetrically nonlinear bending analysis of circular microplates under uniformly distributed transverse loads. The governing differential equations are derived from the principle of minimum total potential energy based on the modified couple stress theory and von Kármán geometrically nonlinear theory in terms of the deflection and radial membrane force, with only one material length scale parameter to capture the size-dependent behavior. The governing equations are firstly discretized to a set of nonlinear algebraic equations by the orthogonal collocation point method, and then solved numerically by the Newton–Raphson iteration method to obtain the size-dependent solutions for deflections and radial membrane forces. The influences of length scale parameter on the bending behaviors of microplates are discussed in detail for immovable clamped and simply supported edge conditions. The numerical results indicate that the microplates modeled by the modified couple stress theory causes more stiffness than modeled by the classical continuum plate theory, such that for plates with small thickness to material length scale ratio, the difference between the results of these two theories is significantly large, but it becomes decreasing or even diminishing with increasing thickness to length scale ratio.     

Nonlinear Bending of thin Elastic Rods

Exact solutions of the problem of nonlinear bending of thin rods under various fixing conditions and point dead loads are obtained. The solutions written in a unified parametric form and expressed in terms of the elliptic Jacobi functions are classified. These solutions depend on a single parameter — modulus of elliptic functions.     

弹性常数对薄板大挠度弯曲的影响

本文针对薄板大挠度弯曲问题，导出了不同弹性常数对应解之间的等价关系式，研究了弹性常数对横向荷载－挠度、弯曲应力－挠度、中面应力－挠度等非线性表达式系数的影响，还分析了文献［１］中近似换算公式的误差。     

Thermal bending analysis of laminated orthotropic plates by the generalized differential quadrature method

The interlaminar stresses in a thin laminated rectangular orthotropic plate with four sides simply supported edges under bending was determined by using the generalized differential quadrature (GDQ) method involving the effects of thermal expansion strain and transverse load. The approximate stress and displacement solutions are obtained under the effects of thermal expansion force and uniform pressure load for eight-layer unidirectional laminates, symmetric cross-ply laminates. Numerical results on the dominant interlaminar stresses and displacement of bending analysis are compared to the Navier solution. The thermal induced forces have significant effect on the bending of plates.     

Bending of thick plates with a concentrated load

In this paper, according to the simplified theory of [1], the bending of rectangular plates with two opposite edges simply supported and other two opposite edges being arbitrary under the action of a concentrated load is treated by means of properties of two-variable-function and the method of series^[2]. The effect of transverse shearing forces on the bending of plates is considered. When the thickness h of plates is small, and the term, whose orders are more than order of h² are neglected, then the results agree with the solutions corresponding to the problem of thin plates^[3]. At the end, the solutions of the bending problem of plates with arbitrary linear distributed load are also obtained.     

Stability of the in-plane bending mode of thin-walled framed structures

The paper outlines an approach to solving the stability problem for framed structures under arbitrary transverse loading. The available methods are limited by one law of variation in the bending moment responsible for loss of stability. The equilibrium equations for a thin-walled bar are integrated assuming that the bending moment is constant. The solution of the Cauchy problem is given in normal form. The arbitrary varying bending moment is approximated by a piecewise-constant function, which will be a little different from the original if the bar is partitioned into a great number of segments. The equations of the boundary-value problem for a discretized framed structure are derived using the boundary-element method. The critical forces and moments are determined from a transcendent equation. Numerical solutions are presented to demonstrate the high accuracy and efficiency of the approach. The solutions of test problems are in agreement with those obtained by Timoshenko     

Nonlinear dynamics of three-dimensional belt drives using the finite-element method

In this paper, new nonlinear dynamic formulations for belt drives based on the three-dimensional absolute nodal coordinate formulation are developed. Two large deformation three-dimensional finite elements are used to develop two different belt-drive models that have different numbers of degrees of freedom and different modes of deformation. Both three-dimensional finite elements are based on a nonlinear elasticity theory that accounts for geometric nonlinearities due to large deformation and rotations. The first element is a thin-plate element that is based on the Kirchhoff plate assumptions and captures both membrane and bending stiffness effects. The other three-dimensional element used in this investigation is a cable element obtained from a more general three-dimensional beam element by eliminating degrees of freedom which are not significant in some cable and belt applications. Both finite elements used in this investigation allow for systematic inclusion or exclusion of the bending stiffness, thereby enabling systematic examination of the effect of bending on the nonlinear dynamics of belt drives. The finite-element formulations developed in this paper are implemented in a general purpose three-dimensional flexible multibody algorithm that allows for developing more detailed models of mechanical systems that include belt drives subject to general loading conditions, nonlinear algebraic constraints, and arbitrary large displacements. The use of the formulations developed in this investigation is demonstrated using two-roller belt-drive system. The results obtained using the two finite-element formulations are compared and the convergence of the two finite-element solutions is examined.     

Stability of nonlinear periodic vibrations of 3D beams

The geometrically nonlinear periodic vibrations of beams with rectangular cross section under harmonic forces are investigated using a p-version finite element method. The beams vibrate in space; hence they experience longitudinal, torsional, and nonplanar bending deformations. The model is based on Timoshenko’s theory for bending and assumes that, under torsion, the cross section rotates as a rigid body and is free to warp in the longitudinal direction, as in Saint-Venant’s theory. The theory employed is valid for moderate rotations and displacements, and physical phenomena like internal resonances and change of the stability of the solutions can be investigated. Green’s nonlinear strain tensor and Hooke’s law are considered and isotropic and elastic beams are investigated. The equation of motion is derived by the principle of virtual work. The differential equations of motion are converted into a nonlinear algebraic form employing the harmonic balance method, and then solved by the arc-length continuation method. The variation of the amplitude of vibration in space with the excitation frequency of vibration is determined and presented in the form of response curves. The stability of the solution is investigated by Floquet’s theory.     

THE PROBLEMS OF THE NONLINEAR UNSYMMETRICAL BENDING FOR CYLINDRICALLY ORTHOTROPIC CIRCULA RPLATE（I）

ＴＨＥＰＲＯＢＬＥＭＳＯＦＴＨＥＮＯＮＬＩＮＥＡＲＵＮＳＹＭＭＥＴＲＩＣＡＬＢＥＮＤＩＮＧＦＯＲＣＹＬＩＮＤＲＩＣＡＬＬＹＯＲＴＨＯＴＲＯＰＩＣＣＩＲＣＵＬＡＲＰＬＡＴＥ（Ｉ）ＱｉｎＳｈｅｎｇ－Ｉｉ（秦圣立）ＨｕａｎｇＪｉａ－ｙｉｎ（黄家寅）（Ｑｕｆ...     

Nonlinear analysis of shear deformable beam-columns partially supported on tensionless three-parameter foundation

In this paper, a boundary element method is developed for the nonlinear analysis of shear deformable beam-columns of arbitrary doubly symmetric simply or multiply connected constant cross-section, partially supported on tensionless three-parameter foundation, undergoing moderate large deflections under general boundary conditions. The beam-column is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM-based method. Application of the boundary element technique yields a system of nonlinear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member as well as the shear forces along the span induced by the applied axial loading. Numerical examples are worked out to illustrate the efficiency, wherever possible, the accuracy and the range of applications of the developed method.     

Karman型正交异性矩形薄板非线性弯曲的统一求解方法

本文对Karman型四边支承正交异性薄板在5种不同边界条件下的几何非线性弯曲进行了统一分析。所设的位移函数均为梁振动函数。它们精确地满足边界条件，利用Galerkin方法和位移函数的正交属性，转换控制方程为非线性代数方程。用“稳定化双共轭梯度法”求解稀疏矩阵线性方程组以及“可调节参数的修正迭代法”求解非线性代数方程组，最后给出了相应的数值结果。     

Nonlinear dynamic analysis of Timoshenko beams by BEM. Part I: Theory and numerical implementation

In this two-part contribution, a boundary element method is developed for the nonlinear dynamic analysis of beams of arbitrary doubly symmetric simply or multiply connected constant cross section, undergoing moderate large displacements and small deformations under general boundary conditions, taking into account the effects of shear deformation and rotary inertia. Part I is devoted to the theoretical developments and their numerical implementation and Part II discusses analytical and numerical results obtained from both analytical or numerical research efforts from the literature and the proposed method. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique yields a nonlinear coupled system of equations of motion. The solution of this system is accomplished iteratively by employing the average acceleration method in combination with the modified Newton–Raphson method. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member, as well as the shear forces along the span induced by the applied axial loading.

     

Nonobvious features of dynamics of circular cylindrical shells

In the framework of the nonlinear theory of flexible shallow shells, we study free bending vibrations of a thin-walled circular cylindrical shell hinged at the end faces. The finite-dimensional shell model assumes that the excitation of large-amplitude bending vibrations inevitably results in the appearance of radial vibrations of the shell. The modal equations are obtained by the Bubnov-Galerkin method. The periodic solutions are found by the Krylov-Bogolyubov method. We show that if the tangential boundary conditions are satisfied “in the mean,” then, for a shell of finite length, significant errors arise in determining its nonlinear dynamic characteristics. We prove that small initial irregularities split the bending frequency spectrum, the basic frequency being smaller than in the case of an ideal shell.     

NONLINEAR THEORY OF MULTILAYER SANDWICH SHELLS AND ITS APPLICATION(Ⅲ)──LARGE DEFLECTION AND POSTBUCKLING OF SHALLOW SHELLS

In this paper, exact solutions of large deflection of multilayer sandwich shallow shells under transverse forces and different boundary conditions are presented. Exact results of postbuckling of multilayer sandwich plates, shallow cylindrical shells and nonlinear deflection of general shallow shells such as spherical shells under inplane edge forces are also obtained by the same procedure.     

转换矩阵在非线性材料杆系分析中的应用

证明了在杆系中,力的转换矩阵与位移的转换矩阵互为转置矩阵,当静不定非线性杆系静力平衡方程确定,而变形协调条件难以确定时,利用转置矩阵可以方便求得静不定非线性杆系的内力及有关节点位移。非线性材料杆系应力-应变关系σ=Bε^1/n中的幂n=2时,非线性材料静不定桁架有可能存在两个解;而采用常规方法求解静不定非线性杆系内力时有可能存在漏解现象,即出现仅能得到一个解的现象。非线性材料杆系应力-应变关系σ=Bε^1/n中的幂n=1时,假设非线性材料杆系各杆内力全部受拉力,或按各杆内力真实受力情况去求各杆内力及节点位移,求得结果的绝对值都是相同的,仅存在符号的差异;与按非线性材料杆系应力-应变关系σ=Bε^1/n中幂n=2时,求得的各杆内力及节点位移的其中一个解的绝对值是一致的。