A novel semi-inverse solution method for elastoplastic torsion of heat treated rods

Torsion rods are a primary component of many power transmission and other mechanical systems. The behavior of these rods under elastoplastic torsion is of major concern for designers. Different methods have so far been proposed which deal with the elastoplastic torsion of rods, most of which assume constant yield stress. This assumption produces rough and inaccurate results when the rods are heat treated, since in the process of heat treatment the form of yield stress distribution across the rod cross section changes. We propose a new method for calculating elastoplastic torsion of rods of simply connected cross section which is based on heat treatment observations. In our method the full plastic stress function is obtained by using the semi-inverse method. Elastoplastic stress function is obtained by generalizing the idea of the membrane analogy and using a piecewise continuous stress function. Since the proposed form of yield stress distribution can not be handled by the current Finite Element packages, we produce a computer package with a 3D graphical interface capable of calculating and displaying the 3D elastoplastic stress function, shear stress contours, and torque-angle of rotation per unit length. We show that our method produces excellent agreement for several known cross sections in comparison to methods which assume constant yield stress.     

Meshless methods for the inverse problem related to the determination of elastoplastic properties from the torsional experiment

The problem of determining the elastoplastic properties of a prismatic bar from the given experimental relation between the torsional moment M and the angle of twist per unit length of the rod’s length θ is investigated as an inverse problem. The proposed method to solve the inverse problem is based on the solution of some sequences of the direct problem by applying the Levenberg-Marquardt iteration method. In the direct problem, these properties are known, and the torsional moment is calculated as a function of the angle of twist from the solution of a non-linear boundary value problem. This non-linear problem results from the Saint-Venant displacement assumption, the Ramberg–Osgood constitutive equation, and the deformation theory of plasticity for the stress–strain relation. To solve the direct problem in each iteration step, the Kansa method is used for the circular cross section of the rod, or the method of fundamental solutions (MFS) and the method of particular solutions (MPS) are used for the prismatic cross section of the rod. The non-linear torsion problem in the plastic region is solved using the Picard iteration.     

Compensation of torsion in rods by piezoelectric actuation

This paper considers the compensation of torsional deformations in rods with the help of thin integrated piezoelectric actuator layers. A laminated orthotropic rod is considered, for which the material properties of each layer are assumed to be homogenous. For the sake of a generalization, the piezoelectric actuation is expressed in terms of eigenstrains. The main scope is the derivation of a distribution of eigenstrains that is able to completely compensate the angle of twist caused by external torsional moments. Saint Venant’s theory of torsion for laminated orthotropic rods is extended for the presence of eigenstrains, which is performed by introducing an additional warping function. It is shown that the actuating torsional moment is a function of the eigenstrains and the additional warping function. For the example of a rectangular cross section, an analytic solution for the actuating moment and the additional warping function is presented. The results are verified by three-dimensional finite-element computations showing a very good accordance with the theoretical results over a large parameter range.     

Stability and vibration of a helical rod constrained by a cylinder

The stability and vibration of an elastic rod with a circular cross section under the constraint of a cylinder is discussed. The differential equations of dynamics of the constrained rod are established with Euler’s angles as variables describing the attitude of the cross section. The existence conditions of helical equilibrium under constraint are discussed as a special configuration of the rod. The stability of the helical equilibrium is discussed in the realms of statics and dynamics, respectively. Necessary conditions for the stability of helical rod are derived in space domain and time domain, and the difference and relationship between Lyapunov’s and Euler’s stability concepts are discussed. The free frequency of flexural vibration of the helical rod with cylinder constraint is obtained in analytical form. The project supported by the National Natural Science Foundation of China (10472067). The English text was polished by Yunming Chen.     

Equations of oscillations of rods of variable composition

We obtain differential equations for the general case of longitudinal, torsional, and transverse oscillations of rods to some parts of which masses are being added or detached. We solve certain special problems concerning the oscillations of such rods of variable composition. In deriving generalized equations of oscillations of rods of variable composition we employ the assumption of planar sections, the assumption of small deformations, and other customary simplifications. We also employ the simplifying assumption of close action; i.e., we assume that the masses being detached and added interact with the rod only at the instant of direct contact. Forces of internal nonelastic resistance are not taken into account. We assume also that in the undeformed state the elastic axis is rectilinear and that the centers of gravity of cross sections are not displaced from their initial positions relative to the cross sections. There may be a change of mass per unit length of the rod both on account of a change in density as well as on account of a change in area of a cross section, the latter being understood to be the union of the initial area of the cross section and the areas of the parts being added and detached. In addition, with the rod there may be associated particles of variable mass distributed continuously or discretely along the length of the rod. We assume that these particles do not interact among themselves but only with the rod.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 103–108, January–February, 1972.     

A new beam element for analyzing geometrical and physical nonlinearity

Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.     

Kovalevskaya弹性杆

应用Kirchhoff比拟讨论Kovalevskaya情况弹性细杆的平衡稳定性问题．导出Kirchhoff方程的解析积分．对于杆截面的主轴与Frenet坐标轴重合的无扭转杆的特殊情形作定性分析，讨论其平衡状态的稳定性与分岔．证明了判断受拉扭作用的圆截面直杆平衡稳定性的Greenhill公式也适用于Kovalevskaya情形的非圆截面杆．     

Post-critical behavior of damped beam columns with variable cross section subjected to distributed follower forces

In this paper the post-critical behavior of beam columns with variable mass and stiffness properties subjected to follower forces arbitrarily distributed along their length in the presence of damping (both internal and external) is investigated using a complete nonlinear dynamic analysis. Although the static nonlinear analysis is more economical in computational cost, it is associated only with the loss of local stability via flutter or divergence. Thus, the nonlinear dynamic analysis is adopted in order to examine the global stability of the system. The governing equations of hyperbolic type are derived in terms of the displacements by considering (a) nonlinear response including the axial deformation, (b) nonlinear response excluding the axial deformation and (c) linear response. Moreover, as the cross-sectional properties of the beam vary along its axis, the resulting coupled nonlinear differential equations have variable coefficients. Their solution is achieved using the analog equation method (AEM) of Katsikadelis. Besides its accuracy and effectiveness, this method overcomes the shortcoming of a possible FEM solution which may experience a lack of convergence. The problems treated in this investigation include beam columns with various load distributions, such as constant, linear and parabolic. Some of the conclusions detected in studying the nonlinear dynamic stability of Beck’s column with variable cross section (Katsikadelis and Tsiatas, Nonlinear dynamic stability of damped Beck’s column with variable cross section. Int. J. Non-linear Mech. 42, 164–171, 2007), are also valid for the case of distributed loads. The important, however, finding is that the post-critical response under distributed loads depends on the law of distribution of mass and stiffness properties, which may lead also to explosive flutter (unbounded amplitude), in contrast to Beck’s column (end-tip load) where the motion is always bounded.     

Bending of a fiber-reinforced viscoelastic composite plate resting on elastic foundations

Composite structures on an elastic foundation are being widely used in engineering applications. Bending response of inhomogeneous viscoelastic plate as a composite structure on a two-parameter (Pasternak’s type) elastic foundation is investigated. The formulations are based on sinusoidal shear deformation plate theory. Trigonometric terms are used in the present theory for the displacements in addition to the initial terms of a power series through the thickness. The transverse shear correction factors are not needed because a correct representation of the transverse shear strain is given. The interaction between the plate and the foundation is included in the formulation with a two-parameter Pasternak’s model. The effective moduli and Illyushin’s approximation methods are used to derive the viscoelastic solution. The effects played by foundation stiffness, plate aspect ratio, and other parameters are presented.     

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.     

Stability analysis of helical rod based on exact Cosserat model

Helical equilibrium of a thin elastic rod has practical backgrounds, such as DNA, fiber, sub-ocean cable, and oil-well drill string. Kirchhoff's kinetic analogy is an effective approach to the stability analysis of equilibrium of a thin elastic rod. The main hypotheses of Kirchhoff's theory without the extension of the centerline and the shear deformation of the cross section are not adoptable to real soft materials of biological fibers. In this paper, the dynamic equations of a rod with a circular cross section are established on the basis of the exact Cosserat model by considering the tension and the shear deformations. Euler's angles are applied as the attitude representation of the cross section. The deviation of the normal axis of the cross section from the tangent of the centerline is considered as the result of the shear deformation. Lyapunov's stability of the helical equilibrium is discussed in static category. Euler's critical values of axial force and torque are obtained. Lyapunov's and Euler's stability conditions in the space domain are the necessary conditions of Lyapunov's stability of the helical rod in the time domain.     

Specific features of the nonlinearly elastic behavior of cylindrical compressible bodies in torsion

The torsion problem of a cylinder with a circular transverse cross section twisted by end moments that are equal in magnitude and opposite in direction is considered for various models of nonlinearly elastic compressible media. The problem is solved by the semi-inverse method of elasticity theory. The Poynting effect, which consists of variation in the length of a shaft in torsion, is treated qualitatively and quantitatively. The results of the numerical and asymptotic (only terms that are quadratic relative to the displacement gradient are conserved) solutions for various models of the nonlinearly elastic behavior of materials are compared. An analysis of the results shows that in some cases, the quasilinear model is not applicable for studying the behavior of nonlinearly elastic compressible media. Rostov State Construction University, Rostov-on-Don 344022. Rostov State University, Rostov-on-Don 344090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 2, pp. 188–193, March–April, 2000.     

On the cross section of minimum stress concentration in the Saint-Venant theory of torsion

A lower bound is derived for the maximum stress in the torsion of cylindrical solids of simply-connected cross section. This bound, which is expressed in terms of the applied torque and the cross-sectional area, is isoperimetric in that it coincides with the maximum stress when the cross section is circular. It confirms the notion that given the applied torque and the area of the cross section, the least maximum stress occurs when the section is circular. Related isoperimetric upper bounds are derived for the minimum value of the stress at boundary points.     

Dynamic deformation of a thermoviscoelastic rod of triangular cross section in a coupled formulation

For the coupled model of a thermoviscoelastic rod of equilateral triangular cross section, two exact solutions are obtained for the cases where a normal displacement and a shear stress or a tangential displacement and a normal stress are specified on the lateral surface of the rod. A dimensionless parameter R₀ is introduced to judge the appropriateness of taking into account the coupling in the formulation of the problem. Formulas are given for the velocities and lengths of the temperature, shear, and longitudinal waves, which can be used in experiments to determine the physical properties of thermoviscoelastic materials. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 4, pp. 128–143, July–August, 2007.     

圆截面弹性细杆的平面振动

基于Kirchhoff理论讨论圆截面弹性细杆的平面振动．以杆中心线的Frenet坐标系为参考系建立动力学方程．杆作平面运动时,其扭转振动与弯曲振动解耦．讨论任意形状杆的扭转振动和轴向受压直杆在无扭转条件下的弯曲振动,证明直杆平衡的静态Lyapunov稳定性与欧拉稳定性条件为动态稳定性的必要条件．考虑轴向力和截面转动惯性效应的影响,导出弯曲振动的固有频率．     

TORSIONAL CENTER OF BEAMS OF ARBITRARY COMPLEX CROSS-SECTION 1)

为了计算任意复杂非圆截面梁横截面扭转中心的位置,用节线法将其约束受扭后所有横截面面外变形的形状用一族包含节线未知函数的曲面表示,建立梁约束受扭时的控制方程后,再用常微分方程求解器分别求出单纯扭矩与横向载荷单独作用时节线未知函数的数值解,最后用刚度等效原理导出复杂截面梁横截面扭转中心的位置。算例计算结果表明：该方法是合理的、有效的,是计算任意复杂非圆截面梁横截面扭转中心位置的可靠方法。     

开口薄壁梁的扭转理论与应用

以Vlasov薄壁构件理论为基础, 推导了开口薄壁构件一阶扭转理论. 该理论考虑了翘曲剪应力对截面转角的影响, 截面的转角分为自由翘曲转角和约束剪切转角, 在约束扭转中, St.Venant扭矩仅仅与自由翘曲转角有关, 而翘曲扭矩仅与约束剪切转角有关. 利用半逆解方法求出了约束扭转中薄壁构件的St.Venant扭矩表达公式; 依据能量方法, 建立了约束剪切转角和翘曲扭矩之间的关系, 并提出了翘曲剪切系数概念, 给出了一阶扭转理论的微分方程. 为了有效求解微分方程, 给出了求解微分方程的初参数法方程和相应的影响函数矩阵; 当St.Venant扭矩可以忽略时, 得到与一阶弯曲理论(Timoshenko梁理论)相似的一阶扭转理论简化形式. 最后利用算例证明了一阶扭转理论和简化理论的有效性.      

Shear stress analysis in engineering beams using deplanation field of special 2-D finite elements

This paper deals with the 2-D finite element shear stress analysis in beams, loaded by bending with shear and St. Venant’s torsion. The properties of these finite elements, like stiffness matrices as well as load vectors, are derived on the basis of their axial nodal displacements, e.g. by warping field. Proposed finite elements enable stress analysis independently of both cross-sectional member shape and material properties. Stiffness matrices and load vectors are derived for several finite element types. Material is assumed to be isotropic and linear elastic. For justification of the proposed stress analysis procedure, some examples are presented.     

Analytical solution of restrained torsional stresses and displacement for rectangular-section box bar with hontycomb core

Differential equation of restrained torsion for rectangular-section box bar with honeycomb core was established and solved by using the method of undetermined function. Non-dimension normal stress, shear stress acting in the faceplate and shear stress acting in the honeycomb-core and warping displacement were deduced. Numerical analysis shows the normal stress attenuates quickly along x-axis. Normal stress acting on the cross section at a distance of 20 h from the fixed end is only one per cent of that acting on the fixed end.     

A Methodology for Obtaining Material Properties of Polymeric Foam at Elevated Temperatures

A new methodology to characterise the elastic properties of polymeric foam core materials at elevated temperatures is proposed. The focus is to determine reliable values of the tensile and compressive moduli and Poisson’s ratio based on strain data obtained using digital image correlation (DIC). In the paper a detailed coverage of the source of uncertainties in the experimental procedure is provided. The uncertainties include those associated with the load introduction, the measurement and the data processing. The design of the specimens and loading jigs are developed and assessed in terms of the introduction of uniform strain. It is shown that due to the mismatch in stiffness between the jig material and the foam the introduction of a uniform strain through the cross section of the specimens is difficult to obtain. A means for correcting for the non-uniform strain across the specimen cross section is developed. To validate the methodology, tests are firstly conducted at room temperature on Divinycell PVC H100 foam. It is shown that the material is highly anisotropic with a stiffness of 50% less in the plane of the foam sheet compared to the through-thickness direction. It is also shown that because of the compliance of the foam, jig misalignment causes large errors in the measurement, and a means for correcting for this is defined. Tests are then conducted in a temperature controlled chamber at elevated temperatures ranging from 20°C to 90°C. A nonlinear reduction in Young’s modulus is obtained with significant degradation occurring after 70°C. The Poisson’s ratio remains fairly stable at different temperatures. A strong theme in the paper is the accuracy and precision of the DIC data and the factors which introduce scatter in the data, along with the uncertainties that this introduces. Particular attention is paid to the affect of the correlation parameters on the derived strain data.