Non-linear flexural-flexural-torsional-extensional dynamics of beams—I. Formulation

The non-linear differential equations of motion, and boundary conditions, for Euler-Bernoulli beams able to experience flexure along two principal directions (and, thus, flexure in any direction in space), torsion and extension are formulated. The beam's material is assumed to be Hookean but its properties may vary along its span. The nonlinearities present in the differential equations include contributions from the curvature expression and from inertia terms. A set of differential equations with polynomial nonlinearities to cubic order, suitable for a perturbation analysis of the motion, is also developed and the validity of the inextensional approximation is assessed. The equations developed here reduce to those for an inextensional beam. In Part II of this paper, a specific example of application is analyzed and the results obtained are compared with those available in the literature where several non-linear terms have been neglected a priori.     

THE STABILITY OF AN AXIALLY ACCELERATING BEAM ON SIMPLE SUPPORTS WITH TORSION SPRINGS

The axially moving beams on simple supports with torsion springs are studied. The general modal functions of the axially moving beam with constant speed have been obtained from the supporting conditions. The contribution of the spring stiffness to the natural frequencies has been numerically investigated. Transverse stability is also studied for axially moving beams on simple supports with torsion springs. The method of multiple scales is applied to the partialdifferential equation governing the transverse parametric vibration. The stability boundary is derived from the solvability condition. Instability occurs if the axial speed fluctuation frequency is close to the sum of any two natural frequencies or is two fold natural frequency of the unperturbed system. It can be concluded that the spring stiffness makes both the natural frequencies and the instability regions smaller in the axial speed fluctuation frequency-amplitude plane for given mean axial speed and bending stiffness of the beam.     

GEOMETRICALLY NONLINEAR FINITE ELEMENT MODEL OF SPATIAL THIN-WALLED BEAMS WITH GENERAL OPEN CROSS SECTION

Based on the theory of Timoshenko and thin-walled beams, a new finite element model of spatial thin-walled beams with general open cross sections is presented in the paper, in which several factors are included such as lateral shear deformation, warp generated by nonuni- form torsion and second-order shear stress, coupling of flexure and torsion, and large displacement with small strain. With an additional internal node in the element, the element stiffness matrix is deduced by incremental virtual work in updated Lagrangian （UL） formulation. Numerical examples demonstrate that the presented model well describes the geometrically nonlinear property of spatial thin-walled beams.     

Free torsion of thin-walled structural members of open- and closed-sections

Free torsion of thin-walled structures of open- and closed-sections is a classical elastic mechanics problem, which, in literature, is often solved by the method of membrane analogy. The method of membrane analogy, however, can be only applied to structures of a single material. If the structure consists of both open- and closed-sections, the method of membrane analogy is difficult to be applied. In this paper, a new method is presented for solving the free torsion of thin-walled structures of open- and/or closed- sections with multiple materials. By utilizing a simple statically indeterminate concept, torsional equations are derived based on the equilibrium and compatibility conditions. The method presented here not only is very simple and easy to understand but also can be applied to thin-walled structures of combined open- and closed-sections with multiple materials.     

轴对称弹性体扭转问题的无网格自然邻接点Petrov-Galerkin法

本文将无网格自然邻接点Petrov-Galerkin 法应用于轴对称弹性体扭转问题的求解．无网格自然邻接点Petrov-Galerkin 法采用自然邻接点插值构造试函数，并且采用三角形线性单元的形函数作为加权残值法的加权函数．自然邻接点插值构造的试函数满足Kronecker delta 函数性质，因此本质边界条件的施加十分方便．由于几何形状和边界条件的轴对称特点，原来的空间问题简化为二维问题求解，因此计算时只需要横截面上离散节点的信息．数值算例结果表明，所提出的方法对求解轴对称弹性体扭转问题是行之有效的．     

Almansi-type boundary conditions for electric potential inducing flexure in linear piezoelectric beams

Using the recent results found in [1, 2] we prove that it is possible to induce flexure in linear piezoelectric beams by means of quadratic Almansi type boundary conditions for the electric potential. Beams constituted by transversely isotropic piezoelectric materials whose symmetry axis is parallel to the axis of the beam are considered. Our choice of boundary conditions for the electric potential has been suggested by the results found in [1, 3]. An explicit expression of material parameters that influence flexure is given in terms of piezoelectric moduli. Received: October 21, 1996     

Some characteristic points of the cross-sections of anisotropic beams

Summary The axis of twist in pure torsion and the centre of flexure in bending by a terminal load for beams with a general rectilinear anisotropy, as well as the axis of flexure for orthotropic beams acted on by a uniform load are investigated. The notions of the axes of pure rotation and of pure translation are introduced. For definiteness, some particular problems treated in more detail involve a restricted type of orthotropy.This work has been sponsored by the United States Army under contract No. DA-11-022-ORD-2059.     

A relook at Reissner’s theory of plates in bending

Shear deformation and higher order theories of plates in bending are (generally) based on plate element equilibrium equations derived either through variational principles or other methods. They involve coupling of flexure with torsion (torsion-type) problem and if applied vertical load is along one face of the plate, coupling even with extension problem. These coupled problems with reference to vertical deflection of plate in flexure result in artificial deflection due to torsion and increased deflection of faces of the plate due to extension. Coupling in the former case is eliminated earlier using an iterative method for analysis of thick plates in bending. The method is extended here for the analysis of associated stretching problem in flexure.     

A structural gradient theory of torsion,the effects of pretwist,and the tension of pre-twisted DNA

A structural gradient theory of torsion of thin-walled beams is developed. A non-local estimate of the mean value of the angle of twist of the beam leads to a shear gradient that is energetically consistent with a bi-moment, in the spirit of the averaging theory of Vardoulakis and Giannakopoulos (2006). The geometric details of the cross section play the role of the microstructure of the beam, introducing a size effect in the torsion problem. The appropriate boundary conditions are derived from the variational formulation of the problem. The proposed gradient elasticity theory is identical to Vlasov’s torsion theory of thin walled elastic beams. The tension of pre-twisted DNA is analyzed at high axial loads, where enthalpic elasticity prevails. A size effect is naturally introduced, indicating that shorter DNA lengths lead to stiffer response in torsion. It is shown also that the complete unwinding of DNA triggers the debonding of its strands.     

Shear stresses in elastic beams: an intrinsic approach

The theory of non-uniform flexure and torsion of Saint-Venant's beam with arbitrary multiply connected cross section is revisited in a coordinate-free form to provide a computationally convenient context. Numerical implementations, by Matlab, are performed to evaluate the maximum elastic shear stresses in beams with rectangular cross sections for different Poisson's ratios. The deviations between the maximum and mean stresses are then diagrammed to adjust the results provided by Jourawski's method.     

Nonlinear dynamic analysis of Timoshenko beams by BEM. Part II: applications and validation

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. In Part I the governing equations of the aforementioned problem have been derived, leading to the formulation of five boundary value problems with respect to the transverse displacements, to the axial displacement and to two stress functions. These problems are numerically solved using the Analog Equation Method, a BEM based method. In this Part II, numerical examples are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. Thus, the results obtained from the proposed method are presented as compared with those from both analytical and numerical research efforts from the literature. More specifically, the shear deformation effect in nonlinear free vibration analysis, the influence of geometric nonlinearities in forced vibration analysis, the shear deformation effect in nonlinear forced vibration analysis, the nonlinear dynamic analysis of Timoshenko beams subjected to arbitrary axial and transverse in both directions loading, the free vibration analysis of Timoshenko beams with very flexible boundary conditions and the stability under axial loading (Mathieu problem) are presented and discussed through examples of practical interest.

     

Parametric stress-concentration factors in greenstick bone fracture

The reported work is a part of an ongoing research program concerned with structural analysis of fractured long bone and methods of internal fixation. The stress-concentration factors for equine metacarpus bones containing greenstick fractures and “through” fractures (surgically repaired) were determined for the compression, flexural and torsional modes of loading based on whole bone (unfractured) strengths. The greenstick type of fracture was simulated with saw cuts at the mid-span of the bone, and the parameters varied were depth of fracture and orientation of fracture. All specimens consisted of fresh dead bone which had been placed in a freezer within 4 hr after expiration. The maximum stress-concentration factors for the simulated greenstick fractures studied were about 3.4 for compression, 4.3 for torsion and 16 for flexure. The stressconcentration factors for fractured bones surgically repaired with commercial plates were about 3.0 for compression, 2.7 for torsion and 6.1 for flexure.     

用样条边界元与奇异积分方程法联合求解多裂纹柱体的扭转

通过间解的分离，本文将径向多裂纹柱体的导曲函两个调和函数表示，使问题归为解一组混混合型积分方程。针对方程的特点，本文联合使用三次样条边界法与奇异积分方程的数值方法对所得方程建立了数值法，并对裂纹相交情形作了特殊处理。最后对工程中感兴趣的一些典型的多裂纹柱体的扭转作了例题计算，结果表明，本文方法具有收敛快，精度高的特点。     

The ellipse of torsional elasticity

Summary In perfect analogy with Culmann's theory of the elasticity ellipse a flexo-torsional ellipse is devised for rapid calculation of the displacements, hyperstatic reactions and influence lines in beams with a straight axis and variable cross section even in conditions of non uniform torsion.

---

Sommario In perfetta analogia con la classica teoria di Culmann dell'ellisse di elasticità, si instaura una ellisse di carattere flesso-torsionale che permette il calcolo rapido degli spostamenti, delle iperstatiche, delle linee di influenza nelle travi ad asse rettilineo a sezione comunque variabile, anche in regime di torsione non uniforme.

     

弹性板边界元中计算边界高阶导数的渐近收敛积分方程

本文针对板弯曲边界元方法中计算边界曲率等高阶导数项时边界积分方程中出现的高阶奇异积分项,通过对未知挠曲函数作渐近展开并加以适当摄动,获得了渐近收敛的边界积分方程。采用这一方法计算板边界上的曲率分布,获得了满意的数值结果。     

“Negative roughness” and polymer drag reduction

Based on an analogy to the Colebrook-White equation, a technique has been developed to allow polymer-solution extrapolation or “scaling” from one pipe size to another at constant values of ΔB. Each experimental data point can be transferred to a new pipe size by a simple, pocket-calculator method which preserves the experimental value of ΔB exactly. Thus scaling can be easily accomplished, without resorting to iteration or graphical techniques. The “negative-roughness” idea can also explain the loss of ΔB or drag reduction with increasing flow velocity.     

采用有限单元法计算梁任意形状截面特性

采用将梁截面离散化的方式,用数值积分计算截面的几何特性,并根据梁剪切变形和扭转理论,利用变分原理建立截面的有限元法方程,求解任意形状截面的扭转常数、剪切中心以及剪切面积修正系数等特性.本方法适用于各种形式的截面,具有计算精度高及适应性强的特点.根据上述理论编制了相应程序,按照不同的单元划分方式,分别计算出矩形截面截面特性,与理论解进行比较;又对舟山市定海长峙至岙山预应力混凝土连续箱梁截面进行了计算,并与Ansys结果进行比较,均证明采用本文的计算方法能得到满意的结果,且该方法适用于各种形状的截面形式.     

A new finite element of spatial thin-walled beams

Based on the theories of Timoshenko＇s beams and Vlasov＇s thin-walled members, a new spatial thin-walled beam element with an interior node is developed. By independently interpolating bending angles and warp, factors such as transverse shear deformation, torsional shear deformation and their Coupling, coupling of flexure and torsion, and second shear stress are considered. According to the generalized variational theory of Hellinger-Reissner, the element stiffness matrix is derived. Examples show that the developed model is accurate and can be applied in the finite element analysis of thinwalled structures.     

Vibration and stability of an axially moving viscoelastic beam with hybrid supports

Vibration and stability are investigated for an axially moving beam constrained by simple supports with torsion springs. A scheme is proposed to derive natural frequencies and modal functions from given boundary conditions of an elastic beam moving at a constant speed. For a beam constituted by the Kelvin model, effects of viscoelasticity on the free vibration are analyzed via the method of multiple scales and demonstrated via numerical simulations. When the axial speed is characterized as a simple harmonic variation about the constant mean speed, the instability conditions are presented for axially accelerating viscoelastic beams in parametric resonance. Numerical examples show the effects of the constraint stiffness, the mean axial speed, and the viscoelasticity.