基于修正偶应力理论的Timoshenko微梁模型和尺寸效应研究

基于修正偶应力理论,将Timoshenko微梁的应力、偶应力、应变、曲率等基本变量,描述为位移分量偏导数的表达式.根据最小势能原理,推导了决定Timoshenko微梁位移场的位移场控微分方程.利用级数法求解了任意载荷作用下Timoshenko简支微梁的位移场控微分方程,得到了反映尺寸效应的挠度、转角及应力的偶应力理论解.通过对承受余弦分布载荷Timoshenko简支微梁的数值计算,研究了Timoshenko微梁的挠度、转角和应力的尺寸效应,分析了Poisson比对Timoshenko微梁力学行为及其尺寸效应的影响.结果表明:当截面高度与材料特征长度的比值小于5时,Timoshenko微梁的刚度和强度均随着截面高度的减小而显著提高,表现出明显的尺寸效应;当截面高度与材料特征长度的比值大于10时,Timoshenko微梁的刚度与强度均趋于稳定,尺寸效应可以忽略;材料Poisson比是影响Timoshenko微梁力学行为及尺寸效应的重要因素,Poisson比越大Timoshenko微梁刚度和强度的尺寸效应越显著.该文建立的Timoshenko微梁模型,能有效描述Timoshenko微梁的力学行为及尺寸效应,可为微电子机械系统（MEMS）中的微结构设计与分析提供理论基础和技术参考.     

Formulation and evaluation of a finite element model for the linear analysis of stiffened composite cylindrical panels

A finite element model for linear static and free vibration analysis of composite cylindrical panels with composite stiffeners is presented. The proposed model is based on a cylindrical shell finite element, which uses a first-roder shear deformation theory. The stiffeners are curved beam elements based on Timoshenko and Saint-Venant assumptions for bending and torsion respectively. The two elements are developed in a cylindrical coordinate system and their stiffness matrices result from a hybrid-mixed formulation where the element assumed stress field is such that exact equilibrium equations are satisfied. The elements are free of membrane and shear locking with correct satisfaction of rigid body motions. Several examples dealing with stiffened isotropic and laminated plates and shells with eccentric as well as concentric stiffeners are analyzed showing the validity of the models.     

Stress and free vibration analysis of functionally graded beams using static Green's functions

In this paper static Green's functions for functionally graded Euler-Bernoulli and Timoshenko beams are presented. All material properties are arbitrary functions along the beam thickness direction. The closed-form solutions of static Green's functions are derived from a fourth-order partial differential equation presented in [2]. In combination with Betti's reciprocal theorem the Green's functions are applied to calculate internal forces and stress analysis of functionally graded beams (FGBs) under static loadings. For symmetrical material properties along the beam thickness direction and symmetric cross-sections, the resulting stress distributions are also symmetric. For unsymmetrical material properties the neutral axis and the center of gravity axis are located at different positions. Free vibrations of functionally graded Timoshenko beams are also analyzed [3]. Analytical solutions of eigenfunctions and eigenfrequencies in closed-forms are obtained based on reference [2]. Alternatively it is also possible to use static Green's functions and Fredholm's integral equations to obtain approximate eigenfunctions and eigenfrequencies by an iterative procedure as shown in [1]. Applying the Sensitivity Analysis with Green's Functions (SAGF) [1] to derive closed-form analytical solutions of functionally graded beams, it is possible to modify the derived static Green's functions and include terms taking cracks into account, which are modeled by translational or rotational springs. Furthermore the SAGF approach in combination with the superposition principle can be used to take stiffness jumps into account and to extend static Green's functions of simple beams to that of discontinuous beams by adding new supports. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于非约束模态的中心刚体-Timoshenko梁动力学建模与分析北大核心CSCD

梁的横向变形会导致梁纵向缩短,建模过程中考虑梁横纵变形二次耦合项则存在动力刚化现象,这说明梁的纵向变形会对模型的广义刚度造成影响.对于做旋转运动的梁结构,旋转运动时还会受到离心力的作用而产生轴向拉力,轴向拉力同样也会引起梁的轴向变形,这种影响对粗短梁更加明显.以大范围运动中心刚体-Timoshenko梁模型为研究对象:首先,运用Timoshenko梁理论以及Hamilton原理建立含离心力的动力学模型;其次,引入非约束模态概念,采用Frobenius方法求解非约束模态振型函数以及固有频率;最后,通过数值仿真探究不同恒定转速时非约束模态与约束模态广义刚度的差异和非约束模态条件下离心力对模型的影响.     

Dynamic modeling for rotating composite Timoshenko beam and analysis on its bending-torsion coupled vibration

A formulation is presented for steady-state dynamic responses of rotating bending-torsion coupled composite Timoshenko beams (CTBs) subjected to distributed and/or concentrated harmonic loadings. The separation of cross section's mass center from its shear center and the introduced coupled rigidity of composite material lead to the bending-torsion coupled vibration of the beams. Considering those two coupling factors and based on Hamilton's principle, three partial differential non-homogeneous governing equations of vibration with arbitrary boundary conditions are formulated in terms of the flexural translation, torsional rotation and angle rotation of cross section of the beams. The parameters for the damping, axial load, shear deformation, rotation speed, hub radius and so forth are incorporated into those equations of motion. Subsequently, the Green's function element method (GFEM) is developed to solve these equations in matrix form, and the analytical Green's functions of the beams are given in terms of piecewise functions. Using the superposition principle, the explicit expressions of dynamic responses of the beams under various harmonic loadings are obtained. The present solving procedure for Timoshenko beams can be degenerated to deal with for Rayleigh and Euler beams by specifying the values of shear rigidity and rotational inertia. Cantilevers with bending-torsion coupled vibration are given as examples to verify the present theory and to illustrate the use of the present formulation. The influences of rotation speed, bending-torsion couplings and damping on the natural frequencies and/or shape functions of the beams are performed. The steady-state responses of the beam subjected to external harmonic excitation are given through numerical simulations. Remarkably, the symmetric property of the Green's functions is maintained for rotating bending-torsion coupled CTBs, but there will be a slight deviation in the numerical calculations.     

Green's functions for forced vibration analysis of bending-torsion coupled Timoshenko beam

A method based on Green's functions is proposed for the analysis of the steady-state dynamic response of bending-torsion coupled Timoshenko beam subjected to distributed and/or concentrated loadings. Damping effects on the bending and torsional directions are taken into account in the vibration equations. The elastic boundary conditions with bending-torsion coupling and damping effects are derived and the classical boundary conditions can be obtained by setting the values of specific stiffness parameters of the artificial springs. The Laplace transform technology is employed to work out the Green's functions for the beam with arbitrary boundary conditions. The Green's functions are obtained for the beam subject to external lateral force and external torque, respectively. Coupling effects between bending and torsional vibrations of the beam can be studied conveniently through these analytical Green's functions. The direct expressions of the steady-state responses with various loadings are obtained by using the superposition principle. The present Green's functions for the Timoshenko beam can be reduced to those for Euler–Bernoulli beam by setting the values of shear rigidity and rotational inertia. In order to demonstrate the validity of the Green's functions proposed, results obtained for special cases are given for a comparison with those given in the literature and they agree with each other exactly. The influences of external loading frequency and eccentricity on Green's functions of bending-torsion coupled Timoshenko beam are investigated in terms of the numerical results for both simply supported and cantilever beams. Moreover, the symmetric property of the Green's functions and the damping effects on the amplitude of Green's functions of the beam are discussed particularly.     

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)     

A micromechanics-based Cosserat-type model for dense particulate solids

A two-dimensional continuum theory is presented for cohesionless granular media consisting of identical rigid disks. While the normal deformation of contacting particles is constrained, the tangential frictional contact is modelled by a line spring with a constant stiffness. To describe the static frictional system transmitting couples at contacts, a Cosserat-type continuum including rotational degrees of freedom is appropriate. Contrary to the classical elastic medium, movement of particles within a granular system in response to applied loads can give rise to localisations of force chains and large voids. In addition to relative displacement and rotation, a director governing the direction of interparticle forces and a phase field delineating density variation, are therefore introduced. Total work done involving these two order parameters for a particle is attained on an orientation average. Based on the formulation of free energy, a concentration- and anisotropy-dependent formulation for static quantities (stress and couple stress) in the rate form is derived in light of the principles of thermodynamics. It is consistent with the requirement of observer independence and material symmetry. The governing equations for two order parameters are derived, in which void concentration and stress anisotropy are related to relative displacement and rotation. As an example, the proposed model is applied to the hardening regime of deformation of a dense particle assembly with initial perfect lattice under simple shear. It is demonstrated that in the presence of dilatancy and director variation, there exists a linear relation between the shear stress and strain, in coincidence with experimental observations. Received: February 24, 2005     

A micromechanics-based Cosserat-type model for dense particulate solids

A two-dimensional continuum theory is presented for cohesionless granular media consisting of identical rigid disks. While the normal deformation of contacting particles is constrained, the tangential frictional contact is modelled by a line spring with a constant stiffness. To describe the static frictional system transmitting couples at contacts, a Cosserat-type continuum including rotational degrees of freedom is appropriate. Contrary to the classical elastic medium, movement of particles within a granular system in response to applied loads can give rise to localisations of force chains and large voids. In addition to relative displacement and rotation, a director governing the direction of interparticle forces and a phase field delineating density variation, are therefore introduced. Total work done involving these two order parameters for a particle is attained on an orientation average. Based on the formulation of free energy, a concentration- and anisotropy-dependent formulation for static quantities (stress and couple stress) in the rate form is derived in light of the principles of thermodynamics. It is consistent with the requirement of observer independence and material symmetry. The governing equations for two order parameters are derived, in which void concentration and stress anisotropy are related to relative displacement and rotation. As an example, the proposed model is applied to the hardening regime of deformation of a dense particle assembly with initial perfect lattice under simple shear. It is demonstrated that in the presence of dilatancy and director variation, there exists a linear relation between the shear stress and strain, in coincidence with experimental observations.     

考虑非局部剪切效应的碳纳米管弯曲特性研究

基于Hamilton(哈密顿)变分原理和非局部连续介质弹性理论,建立了新型非局部Timoshenko(铁木辛柯)梁模型（ANT）,推导了碳纳米管（CNT）的ANT弯曲平衡方程以及两端简支梁、悬臂梁和简支 固定梁的边界条件表达式,分析了剪切变形效应和非局部微观尺度效应对碳纳米管弯曲特性的影响.数值计算结果显示,碳纳米管的弯曲刚度随着小尺度效应的增强而升高.其次,这种小尺度效应对自由端受集中力的悬臂梁碳纳米管有明显作用,其刚度变化规律和其它约束条件的碳纳米管一样,这一点是ANT模型区别于普通非局部纳米梁模型的主要特点.经分子动力学模拟验证,ANT模型是合理分析碳纳米管力学特性的有效方法.     

具有分数导数本构关系的粘弹性Timoshenko梁的静动力学行为分析

利用粘弹性材料的三维分数导数型本构关系,建立粘弹性Timoshenko梁的静、动力学行为研究的数学模型;分析Timoshenko梁在阶跃载荷作用下的准静态力学行为,得出了问题的解析解,考察了一些材料参数对梁的挠度的影响。基于模态函数讨论了粘弹性Timoshenko梁在横向简谐激励作用下的动力响应,并考察了剪切和转动惯性对梁振动响应的影响。     

Vibrations of Timoshenko beams with isogeometric approach

A study on the free vibration analysis of Timoshenko beams is presented here. In order to determine natural frequencies of beams, a thick beam element is developed by using isogeometric approach based on Timoshenko beam theory which allows the transverse shear deformation and rotatory inertia effect. Three refinement schemes such as h-, p- and k-refinement are used in the analysis and the identification of shear locking is also conducted by using numerical examples. From numerical results, the present element can produce very accurate values of natural frequencies and the mode shapes due to exact definition of the geometry. With higher order basis functions, there is no shear locking phenomenon in very thin beam situations. Finally, the benchmark tests described in this study are provided as future reference solutions for Timoshenko beam vibration problem.     

Static and dynamic analysis of cracked or weakened structures

Stiffness modifications in engineering structures, for example due to damage and cracking, will inevitably also lead to changes in deformations, internal forces, natural frequencies and mode shapes of the structures. In this paper, an efficient and simple method for sensitivity analysis of cracked or weakened structures under time-harmonic loading is presented. The method is based on a comparison between the strain energy and the kinetic energy of an uncracked structure and that of a cracked structure in conjunction with the application of exact or approximate Green's functions as described in [3] for the static case. The present analysis enables the prediction of any changes in the displacements and stresses and has a lower computational effort as compared to available classical methods, because only the damaged region has to be re-considered in the method. Green's functions are taken as a basis of the approach, which have the ability to weight the influence of the stiffness modifications in a region of a structure and show how sensitive other regions respond to the stiffness modifications. Based on linear elastic fracture mechanics, cracked or damaged regions are approximated by spring models in the analytical solution of some simple beam problems, while cracked finite elements are used for complicated cases where analytical solutions cannot be obtained. Sensitivity analysis with Green's functions (SAGF) approach is applied to static and dynamic analysis of cracked and weakened structures, which consist of homogeneous materials or fiber reinforced composites like reinforced concretes. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于剪切效应纤维梁单元的结构非线性有限元数值模拟

基于Euler-Bernoulli梁理论的经典纤维模型忽略了剪切变形给截面带来的影响,为了得到更加精确的梁单元模型,该文基于考虑剪切效应的纤维梁单元,根据Timoshenko梁理论,推导了该纤维梁单元的刚度矩阵,并结合弹塑性增量理论,同时考虑了几何非线性和材料非线性的双重影响,建立了压弯剪复杂应力状态下结构非线性有限元...     

准三维功能梯度微梁的尺度效应模型及微分求积有限元

基于修正的偶应力理论与四参数高阶剪切-法向伸缩变形理论,提出了一种具有尺度依赖性的准三维功能梯度微梁模型,并应用于小尺度功能梯度梁的静力弯曲和自由振动分析中.采用第二类Lagrange方程,推导了微梁的运动微分方程及边界条件.针对一般边值问题,构造了一种融合Gauss-Lobatto求积准则与微分求积准则的2节点16自由度微分求积有限元.通过对比性研究,验证了理论模型以及求解方法的有效性.最后,探究了梯度指数、内禀特征长度、几何参数及边界条件对微梁静态响应与振动特性的影响.结果表明,该文所发展的梁模型及微分求积有限元适用于研究各种长细比的功能梯度微梁的静/动力学问题,引入尺度效应会显著地改变微梁的力学特性.     

中心刚体-外Timoshenko梁系统的建模与分岔特性研究

对于中心刚体固结悬臂梁系统,当不考虑梁剪应力(即Euler-Bernoulli梁)影响时,匀速转动梁的平凡解是稳定的。而对于深梁,有必要考虑剪应力(即Timoshenko梁)的影响,此时其匀速转动平凡解将出现拉伸屈曲。为此采用广义Hamilton变分原理建立了中心刚体固结Timoshenko梁这类刚-柔耦合系统的非线性动力学模型,应用数值方法研究了匀速转动Timoshenko梁非线性系统的分岔特性,以及失稳的临界转速。     

Elasticity solution of multi-span beams with variable thickness under static loads

This paper studies the stress and displacement distributions of continuously varying thickness multi-span beams simply supported at two ends and under static loads. The intermediate supports of the beam may be elastic and/or rigid in one or two directions. On the basis of the two-dimensional plane elasticity theory, the general solution of stress function, which exactly satisfies the governing differential equations and the simply supported boundary conditions, is deduced. In the present analysis, the reaction forces of the intermediate supports are regarded as the unknown external forces acting on the lower surface of the beam under consideration. The unknown coefficients in the solutions are determined by using the Fourier sinusoidal series expansions to the boundary conditions on the upper and lower surfaces of the beam and using the linear relations between reaction forces and displacements of the beam at intermediate supports. The solution obtained is exact and excellent convergence has been confirmed. Comparing the numerical results obtained from the proposed method to those obtained from the Euler beam theory, the Timoshenko beam theory and those obtained from the commercial finite element software ANSYS, high accuracy of the present method is demonstrated.     

Exact analytic solutions to the problem on plane buckling modes of rectangular orthotropic plates with free edges under biaxial loading

A two-dimensional linearized problem on plane buckling modes (BMs) of a rectangular plate with free edges, made of an elastic orthotropic material, underbiaxial tension-compression is considered. With the use of double trigonometric basis functions, displacement functions exactly satisfying all static boundary condition on plate edges are constructed. It is shown that the exact analytic solutions found describe only the pure shear BMs, and if the normal stress in one direction is assumed equal to zero, an analog of the solution given by the kinematic Timoshenko model can be obtained. Upon performing the limit passage to the zero harmonic in the displacement functions of one of the directions, the solution to the problem of biaxial compression can be obtained by equating the Poisson ratio to zero; in the case of uniaxial compression, this solution exactly agrees with that following from the classical Bernoulli-Euler model. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 2, pp. 149–170, March–April, 2007.     

A smart displacement based (SDB) beam element with distributed plasticity

The standard displacement based inelastic beam element suffers of approximations related to the inability of the cubic polynomial interpolation functions to properly describe the displacement response of the beam when exhibiting inelastic behaviour. The increase of the number of finite elements, or the use of higher order functions with additional internal degrees of freedom, are common remedies suggested to improve the approximation leading to an unavoidable reduction of the computational efficiency. Alternatively, it has been shown that the development of force based finite elements, based on the adoption of exact force shape functions, lead to more accurate results, although requiring different and more complicated iterative solution strategies. Within this scenario, this paper proposes a new inelastic beam element, within the context of the displacement based approach, based on variable displacement shape functions, whose analytic expressions are related to the plastic deformation evolution in the beam element. The adaptive generalised displacement shape are obtained by identifying, at each step, an equivalent tangent beam, characterised by abrupt variations of flexural stiffness, as a suitable representation of the current inelastic state of the beam. The presented approach leads to the formulation of a Smart Displacement Based (SDB) beam element whose accuracy appears to be comparable to those obtained through a force based approach but requiring a reduced implementation effort and a more straightforward approach. The term ‘smart’ aims at emphasizing the ability of the element to upgrade the displacement field according to the current inelastic state.     

轴向运动导电导磁梁的磁弹性振动方程

针对磁场环境中轴向运动导电导磁梁磁弹性耦合振动的理论建模问题进行研究.基于Timoshenko(铁木辛柯)梁理论并考虑几何非线性因素,给出轴向运动弹性梁在横向双向振动下的形变势能、动能计算式以及电磁力和机械力的虚功表达式.应用Hamilton(哈密顿)变分原理,推得磁场中轴向运动Timoshenko梁的非线性磁弹性耦合振动方程,并给出了简化形式的Euler-Bernoulli(欧拉 伯努利)梁磁弹性振动方程.根据电磁理论和相应的电磁本构关系,得到载流导电弹性梁所受电磁力的表达式,基于磁偶极子-电流环路模型给出铁磁弹性梁所受磁体力和磁体力偶的表述形式.通过算例,分析了轴向运动导电弹性梁的奇点分布及其稳定性问题.