A distributed-mass approach for dynamic analysis of Bernoulli-Euler plane frames

The exact dynamic analysis of plane frames should consider the effect of mass distribution in beam elements, which can be achieved by using the dynamic stiffness method. Solving for the natural frequencies and mode shapes from the dynamic stiffness matrix is a nonlinear eigenproblem. The Wittrick-Williams algorithm is a reliable tool to identify the natural frequencies. A deflated matrix method to determine the mode shapes is presented. The dynamic stiffness matrix may create some null modes in which the joints of beam elements have null deformation. Adding an interior node at the middle of beam elements can eliminate the null modes of flexural vibration, but does not eliminate the null modes of axial vibration. A force equilibrium approach to solve for the null modes of axial vibration is presented. Orthogonal conditions of vibration modes in the Bernoulli-Euler plane frames, which are required in solving the transient response, are theoretically derived. The decoupling process for the vibration modes of the same natural frequency is also presented.     

Exact spectral elements for follower tension buckling by power series

The spectral dynamic stiffness method using exact solutions of the governing equations as shape functions has been popular for vibration and dynamic stability analyses of framed structures consisting of uniform members. Since non-uniform members do not generally have closed form solutions, special cases only have been considered. However, exact solutions are still possible for generally non-uniform members using power series. The paper studies the exact dynamic stability of columns with distributed axial force by power series. Both uniform and distributed, compression and tension, and conservative and non-conservative axial forces are considered. Interaction diagrams of various kinds of axial loads on the natural frequencies including different intensities of the distributed loads and degree of tangency are given. Follower tension buckling is reported for the first time. It is found that the power series outperforms the dynamic stiffness method in terms of versatility in applications and numerical stability at the very low and high ends of the frequency spectrum.     

Consideration of mechanical properties of single-walled carbon nanotubes under various loading conditions

In this article, mechanical properties of single-walled carbon nanotubes (SWCNTs) with various radiuses under tensile, compressive and lateral loads are considered. Stress–strain curve, elastic modulus, tensile, compressive and rotational stiffness, buckling behaviour, and critical axial compressive load and pressure of eight different zigzag and armchair SWCNTs are investigated to figure out the effect of radius and chirality on mechanical properties of nanotubes. Using molecular dynamic simulation (MDS) method, it can be explained that SWCNTs have higher Young’s modulus and tensile stiffness than compressive elastic modulus and compressive stiffness. Critical axial force of zigzag SWCNT is independent from the radius, but that of armchair type rises by increasing of radius, also these two types show different buckling modes.     

Coupled bending and torsional vibration of axially loaded Bernoulli-Euler beams including warping effects

The dynamic transfer matrix method for determining natural frequencies and mode shapes of the bending-torsion coupled vibration of axially loaded thin-walled beams with monosymmetrical cross sections is developed by using a general solution of the governing differential equations of motion based on Bernoulli-Euler beam theory. This method takes into account the effect of warping stiffness and gives allowance to the presence of axial force. The dynamic transfer matrix is derived in detail. Two illustrative examples on the application of the present theory are given for bending-torsion coupled beams with thin-walled open cross sections. The effects of axial load and warping stiffness on coupled bending-torsional frequencies are discussed. Compared with those available in the literature, numerical results demonstrate the accuracy and effectiveness of the proposed method.     

空间弹性细杆在扭矩作用下的屈曲

用有限元方法研究了三维弹性细杆在扭矩作用下的屈曲.利用自然坐标形式的细长空间曲杆的能量方程和2节点12个自由度的自然坐标形式的三维曲梁单元,采用特征值分析方法,研究分析了同时受有轴力和扭矩作用时的空间弹性细杆的屈曲问题.数值结果与存在的理论解极为吻合.具有一定曲率和挠率的空间细长曲杆,其临界扭矩值与扭矩的指向有着极其明显的关系.     

Solution to the problem of Nicolai

We consider the problem of Nicolai on dynamic stability of an elastic cantilever rod loaded by an axial compressive force P and a twisting tangential torque L in continuous formulation. The problem is to find the stability region for non-equal principal moments of inertia of the rod in the space of three parameters: P, L and the parameter α for the ratio of principal moments of inertia. New governing equations and boundary conditions, which form the basis for analytical and numerical studies, are derived. An important detail of this formulation is that the pre-twisting of the rod due to the torque L is taken into account. The singular point on the stability boundary at the critical Euler force P_E is recognized and investigated in detail. For an elliptic cross-section of a uniform rod the stability region is found numerically with the use of the Galerkin method and the exact numerical approach. The obtained numerical results are compared with the analytical formulas of the asymptotic analysis.     

Exact dynamic and static stiffness matrices of shear deformable thin-walled beam-columns

Total potential energy of non-symmetric thin-walled beam-columns in the general form is presented by introducing the displacement field based on semitangential rotations and deriving transformation equations between displacement and force parameters defined at the arbitrary axis and the centroid-shear center axis, respectively. Next, governing equations and force-deformation relations are derived from the total potential energy for a shear-deformable, uniform beam element and a system of linear eigenproblem with non-symmetric matrices is constructed based on 14 displacement parameters. And then explicit expressions for displacement parameters are derived and exact dynamic stiffness matrices are determined using force-deformatin relationships. In addition, the modified numerical method to eliminate multiple zero eigenvalues and to evaluate the exact static stiffness matrix is developed for spatial stability analysis. Finally, in order to demonstrate the validity and the accuracy of this study, the spatially coupled natural frequencies and buckling loads are evaluated and compared with analytical solutions or results analyzed by thin-walled beam elements and ABAQUS's shell elements.     

定跨度变截面梁的弯曲问题

一.引言 在处理梁的弯曲问题时,人们经常利用函数级数来表示有关的各量,并后而得到各该量的近似值。胡海昌曾经指出:在横向载荷和轴向力同时作用下,适宜于用梁的屈曲的本徵函数展开式来表示梁的挠度;其中φ_n是满足所给的梁的支座情况的屈曲本徵函数,a_n是常数系数。他求得一个相当简单的公式以已知的本徵函数和本徵值表示诸系数     

Comparison of moments and couple-generating forces near discontinuities in structural-acoustic systems

In modeling vibration isolators in structural-acoustic systems, the isolator's dynamic properties are often treated as acting only in the axial direction as moments are often neglected. Furthermore, the size, or scale, of the isolator is often neglected and the isolator is assumed to act at single points on the connected structures. Previous work has shown that concentrated moments can be particularly important when located near a fixed support or a structural discontinuity. This research extends that work to examine the importance of moment scale effects for a system containing a distributed structural discontinuity with its own scale. Moment scale effects are examined by determining the difference in radiated acoustic power for a simple system that is excited by a couple-generating distributed force and a concentrated moment. The distributed force produces a couple that is equivalent to the concentrated moment. As a result, only the scale is being examined. Particular interest occurs when the excitation is located near the structural discontinuity. Based on the cases studied here, moment scale is shown to be important at lower frequencies when the excitation is located near the edges of the discontinuity. At higher frequencies, any overlap of the excitation and discontinuity may warrant the need to consider moment scales.     

Study of a misaligned flexibly coupled shaft system having nonlinear bearings and cyclic coupling stiffness—Theoretical model and analysis

A model which enables dynamic analysis of flexibly coupled misaligned shafts is presented. The model is setup to account for both angular and parallel misalignment in the presence of mass unbalance and incorporates a coupling having angular, torsional and axial flexibility. Among the important features is the ability to simulate both nonlinear bearing stiffness and coupling angular-stiffness anisotropy. The equations of motion are derived for the linear system, extended to include nonlinear bearing effects and subsequently transformed into non-dimensional form for general application. A series of numerical analyses are performed and the influence of important system parameters assessed thereby providing insight to the resulting static and dynamic forces and motions. Angular and parallel misalignments are shown to produce fundamentally different system response. It is found that the static preload induced by both types of misalignment can play a key role in producing complex vibration resulting from it's interaction with rotating-element anisotropy and bearing nonlinear properties. Bearing static forces are altered and rotating elements are subjected to alternating forces which could affect fatigue life. Bearing forces can be further modified by the application of transmitted torque. The potential for great variability in system response is shown to exist due to the participation of numerous influential variables.     

Molecular dynamic investigation of mechanical properties of armchair and zigzag double-walled carbon nanotubes under various loading conditions

Using molecular dynamic simulation (MDS), effects of chirality and Van der Waals interaction on Young's modulus, elastic compressive modulus, bending, tensile, and compressive stiffness, and critical axial force of double-walled carbon nanotube (DWCNT) and its inner and outer tubes are considered. Achieving the highest safety factor, mechanical properties have been investigated under applied load on both inner and outer tubes simultaneously and on each one of them separately. Results indicate that as a compressive element, DWCNT is more beneficial than single-walled carbon nanotube (SWCNT) since it carries two times higher compression before buckling. Except critical axial pressure and tensile stiffness, in other parameters zigzag DWCNT shows higher amounts than armchair type. Outer tube has lower strength than inner tube; therefore, most reliable design of nanostructures can be attained if the mechanical properties of outer tube taken as the properties of DWCNT.     

Nonlocal shear deformable shell model for bending buckling of microtubules embedded in an elastic medium

A nonlocal shear deformable shell model is developed for buckling of microtubules embedded in an elastic matrix of cytoplasm under bending in thermal environments. The results reveal that the lateral constraint has a significant effect on the buckling moments of a microtubule when the foundation stiffness is sufficiently large.     

Natural vibration of pre-twisted shear deformable beam systems subject to multiple kinds of initial stresses

Free vibration and buckling of pre-twisted beams exhibit interesting coupling phenomena between compression, moments and torque and have been the subject of extensive research due to their importance as models of wind turbines and helicopter rotor blades. The paper investigates the influence of multiple kinds of initial stresses due to compression, shears, moments and torque on the natural vibration of pre-twisted straight beam based on the Timoshenko theory. The derivation begins with the three-dimensional Green strain tensor. The nonlinear part of the strain tensor is expressed as a product of displacement gradient to derive the strain energy due to initial stresses. The Frenet formulae in differential geometry are employed to treat the pre-twist. The strain energy due to elasticity and the linear kinetic energy are obtained in classical sense. From the variational principle, the governing equations and the associated natural boundary conditions are derived. It is noted that the first mode increases together with the pre-twisted angle but the second decreases seeming to close the first two modes together for natural frequencies and compressions. The gaps close monotonically as the angle of twist increases for natural frequencies and buckling compressions. However, unlike natural frequencies and compressions, the closeness is not monotonic for buckling shears, moments and torques.     

Synchronisation and general dynamic symmetry of a vibrating system with two exciters rotating in opposite directions

We derive the non-dimensional coupling equation of two exciters, including inertia coupling, stiffness coupling and load coupling. The concept of general dynamic symmetry is proposed to physically explain the synchronisation of the two exciters, which stems from the load coupling that produces the torque of general dynamic symmetry to force the phase difference between the two exciters close to the angle of general dynamic symmetry. The condition of implementing synchronisation is that the torque of general dynamic symmetry is greater than the asymmetric torque of the two motors. A general Lyapunov function is constructed to derive the stability condition of synchronisation that the non-dimensional inertia coupling matrix is positive definite and all its elements are positive. Numeric results show that the structure of the vibrating system can guarantee the stability of synchronisation of the two exciters, and that the greater the distances between the installation positions of the two exciters and the mass centre of the vibrating system are, the stronger the ability of general dynamic symmetry is.     

考虑介质膨胀速率的锂离子电池管状电极中扩散诱导应力及轴向支反力分析

硅作为锂离子电池阴极材料相对于传统负极材料具有高比容量,价格低廉等优势.本文针对充电过程中锂离子电池中电极建立力学模型和扩散模型,并在扩散模型引入考虑介质膨胀速率的影响.以硅空心柱形电极为例,分析了恒流充电下介质膨胀速率对电极中扩散诱导应力分布的影响,并研究了不同内外半径比、充电速率、材料参数以及锂化诱导软化系数(lithiation induced softening factor,LISF)对轴向的支反力达到临界欧拉屈曲力所需时间的影响.结果表明,随着电极中锂浓度上升,介质膨胀速率对应力分布的影响增大,对轴向的支反力影响较小.弹性模量和应力成正比,但其与轴向的支反力达到临界欧拉屈曲力所需时间无关;扩散系数与所需时间成反比;偏摩尔体积增大时,达到临界屈曲力所需时间减少;随着LISF绝对值增大,完全锂化时轴向力降低.     

轴向冲击载荷作用下双壁碳纳米管的动力屈曲

应用改进的有限元方法,建立考虑层间范德华力作用的壳-弹簧非线性有限元模型,基于B-R运动准则,系统地研究了双壁碳纳米管的动力屈曲问题,得到了轴向冲击载荷作用下双壁碳纳米管的临界动力屈曲载荷和临界动力失效载荷. 研究结果表明,在动力屈曲过程中,双壁碳纳米管层间距的变化非常小,各管的变形相互协调;碳纳米管中应力波的传播导致碳纳米管出现非对称屈曲模态,可明显观测到四个环向波瓣,沿着碳纳米管的轴线方向,四个波瓣的波峰和波谷交替变化. 对碳纳米管动力屈曲问题的研究表明,冲击载荷的大小和持续时间对碳纳米管的动力屈曲有 关键词： 碳纳米管 动力屈曲 冲击载荷     

STOCHASTIC DYNAMIC SENSITIVITY OF UNCERTAIN STRUCTURES SUBJECTED TO RANDOM EARTHQUAKE LOADING

The present study involves computation of stochastic sensitivity of structures with uncertain structural parameters subjected to random earthquake loading. The formulations are provided in frequency domain. A strong earthquake-induced ground motion is considered as a random process defined by respective power spectral density function. The uncertain structural parameters are modelled as homogeneous Gaussian stochastic field and discretized by the local averaging method. The discretized stochastic field is simulated by the Cholesky decomposition of respective co-variance matrix. By expanding the dynamic stiffness matrix about its reference value, the advantage of Neumann Expansion technique is explored within the framework of Monte Carlo simulation, to compute responses as well as sensitivity of response quantities. This approach involves only a single decomposition of the dynamic stiffness matrix for the entire simulated structure and the facility that several stochastic fields can be tackled simultaneously are basic advantages of the Neumann Expansion. The proposed algorithm is explained by an example problem.     

Dynamic stiffness method for inplane free vibration of rotating beams including Coriolis effects

The paper addresses the in-plane free vibration analysis of rotating beams using an exact dynamic stiffness method. The analysis includes the Coriolis effects in the free vibratory motion as well as the effects of an arbitrary hub radius and an outboard force. The investigation focuses on the formulation of the frequency dependent dynamic stiffness matrix to perform exact modal analysis of rotating beams or beam assemblies. The governing differential equations of motion, derived from Hamilton's principle, are solved using the Frobenius method. Natural boundary conditions resulting from the Hamiltonian formulation enable expressions for nodal forces to be obtained in terms of arbitrary constants. The dynamic stiffness matrix is developed by relating the amplitudes of the nodal forces to those of the corresponding responses, thereby eliminating the arbitrary constants. Then the natural frequencies and mode shapes follow from the application of the Wittrick–Williams algorithm. Numerical results for an individual rotating beam for cantilever boundary condition are given and some results are validated. The influences of Coriolis effects, rotational speed and hub radius on the natural frequencies and mode shapes are illustrated.     

Dynamic column buckling of multi-walled carbon nanotubes under axial impact load

This paper studies the dynamic column buckling of multi-walled carbon nanotubes (MWNTs) under axial impact load. The analysis is based on the continuum mechanics model and a simplified model for the van der Waals forces between adjacent layers. By introducing initial imperfections for MWNTs and applying the method of preferred mode, a buckling condition is derived for the buckling load and associated buckling mode. In particular, explicit expressions are obtained for double-walled carbon nanotubes (DWNTs). Finally, numerical calculations are worked out for a DWNT and a five-layer MWNT with different length-to-radius ratios.     

Estimating ensemble average power delivered by a piezoelectric patch actuator to a non-deterministic subsystem

Engineering systems such as aircraft, ships and automotive are considered built-up structures. Dynamically they are taught of as being fabricated from many components that are classified as ‘deterministic subsystems’ (DS) and ‘non-deterministic subsystems’ (Non-DS). Structures' response of the DS is deterministic in nature and analysed using deterministic modelling methods such as finite element (FE) method. The response of Non-DS is statistical in nature and estimated using statistical modelling technique such as statistical energy analysis (SEA). SEA method uses power balance equation, in which any external input to the subsystem must be represented in terms of power. Often, input force is taken as point force and ensemble average power delivered by point force is already well-established. However, the external input can also be applied in the form of moments exerted by a piezoelectric (PZT) patch actuator. In order to be able to apply SEA method for input moments, a mathematical representation for moment generated by PZT patch in the form of average power is needed, which is attempted in this paper. A simply-supported plate with attached PZT patch is taken as a benchmark model. Analytical solution to estimate average power is derived using mobility approach. Ensemble average of power given by the PZT patch actuator to the benchmark model when subjected to structural uncertainties is also simulated using Lagrangian method and FEA software. The analytical estimation is compared with the Lagrangian model and FE method for validation. The effects of size and location of the PZT actuators on the power delivered to the plate are later investigated.