结构模态分析的有限元并行模拟

以子结构模态综合分析为基础,提出一种求解大型结构特征值问题的并行解法．采用子结构模态综合算法,结构特征模态采用子空间迭代方式并行求解．这种子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和质量阵求解综合系统的特征值问题．数值结果表明这种求解大型结构特征值问题的并行算法是可行有效的．     

基于模态缩减和模态实验数据的结构连接参数识别方法

利用在结构系统可测自由度上获得的不完备模态参数和子结构的有限元模型，根据模态缩减理论，建立了识别子结构间连接子结构参数的优化模型，采用逐次二次规划法求解，改善了测试噪声和模态截断误差的影响。该方法识别精度高、收敛速度快、计算量小，便于工程应用。     

大型结构特征值问题的Lanczos子结构并行算法

本文利用子结构和Lanczos方法,提出了大型结构固有频率与模态的并行解法。该方法在Lanczos方法的求解过程中,仅利用子结构刚度阵和质量阵并行进行凝聚,进而求得新的迭代矢量,最终求得三对角阵对应的特征值和特征向量。该算法在西安交通大学ELXSI-6400并行计算机上程序实现,计算结果表明能有效地节省计算时间和计算机的内存,为一种有效的大型工程结构动力问题的求解方法。     

基于高精度模态综合法的螺栓连接组合结构不确定模型更新

由于螺栓结构中存在复杂的微观特征,数值模型中无法避免的会引入不确定性。针对螺栓设计中参数不确定性量化问题,一种基于高精度模态综合法的螺栓连接组合结构随机模型更新方法被提出。首先考虑高阶剩余模态的影响,推导了适用于螺栓组合结构弹性交界面的高精度模态综合法;之后,在贝叶斯推理框架下,通过组合结构的自振频率和振型的概率分布,建立螺栓连接参数的后验概率密度函数,并借助DRAM(delayed rejection adaptive metropolis)抽样方法量化参数的不确定性。数值算例表明,针对各子结构之间自振频率差别较大的结构,与模态综合法相比,本研究方法在较少子结构模态参与的情形下,也可以保证良好的模型更新结果。     

基于键合图法的多个子结构自由界面模态综合法研究和应用

将键合图方法用于动态子结构研究,提出了一种求解多个子结构自由界面模态综合法的新方法。通过一系列物理和数学上的分析,详细推导出多个子结构自由界面模态综合法的计算过程。在本文给出的算例中,基于键合图方法的自由界面模态综合法,通过建立各个子结构的状态空间方程,计算子结构相应的特征值矩阵和振型矩阵;这些子结构在进行模态综合后,获得整体结构的特征值矩阵和振型矩阵,该结果和原来整体结构完全相同,这进一步说明了本文提出的新方法的正确性。运用本文提出的方法建立状态空间方程,在子结构的模态综合过程中,不会产生系统特征值的增失根问题,确保了子结构在综合后其整体结构模态信息的完整性。     

随机桁架结构的非平稳随机动力响应分析

本文研究了随机桁架结构在非平稳随机激励下的动力响应问题。在利用随机因子法分析随机结构动力特性的基础上，给出了一种分析随机结构非平稳随机响应的新方法。从结构非平稳随机响应的表达式出发，同时考虑桁架结构的物理参数、几何尺寸的随机性，利用求解随机变量函数矩的方法和求解随机变量数字特征的代数综合法，导出了随机桁架结构在非平稳随机激励下位移响应均方值和应力响应均方值的均值、方差和变异系数的计算表达式。通过算例，分析了结构物理参数和几何尺寸的随机性对结构位移响应均方值和应力响应均方值随机变量随机性的影响。本文方法具有对随机结构进行一次动力分析便可求得动力响应的数字特征，且可以考察结构任一参数的随机性对结构非平稳随机响应分析结果的影响之优点。     

大型结构特征值问题的混合粒度并行算法

本文提出一种求解大形结构特征值问题的粗细粒度混合并行算法：在子结构模态综合粗粒度并行算基础上,综合系统的特性值问题采用细粒度并行方式求解。细粒度并行包括子空间迭代法的子结构并行算法、雅可比分块并行计算的方法和一种Ｎｅｗｔｏｎ－Ｒａｐｈｏｎ迭代法在多处理器上任力均衡分配的有效策略。子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和国求综合系统的特征值问题。利用雅     

线性随机结构的非平稳随机响应变异性分析

对于具有随机参数的结构受到非平稳随机激励的问题，给出了结构随机响应变异系数的虚拟激励摄动算法。它应用虚拟激励法先将随机荷载转化为确定性荷载，以使随机问题精确地转化为仅结构参数具有随机性的问题，从而将问题归结为应用随机摄动法求解单随机问题。求解过程简单高效，且有较高的精度。     

柴油机曲轴整体结构的三维有限元分析

本文采用子结构方法对柴油机曲轴整体结构进行了三维有限元分析，对曲轴这种复杂的工程结构问题通过分成子结构进行求解可以建立更加灵活合理的结构模型，也可降低方程组的阶数，节省大量的重复计算工作，使我们能够在没有大机器的情况下，用很少的内存容量解算大自由度数的题目。     

柴油机曲轴整体结构的三维有限元分析

本文采用子结构方法对柴油机曲轴整体结构进行了三维有限元分析，对曲轴这种复杂的工程结构问题通过分成子结构进行求解可以建立更加灵活合理的结构模型，也可降低方程组的阶数，节省大量的重复计算工作，使我们能够在没有大机器的情况下，用很少的内存容量解算大自由度数的题目。     

Theoretical analysis of transient waves in a simply-supported Timoshenko beam by ray and normal mode methods

The dynamic transient responses of a simply-supported Timoshenko beam subjected to an impact force are investigated by two theoretical approaches – ray and normal mode methods. The mathematical methodology proposed in this study for the ray method enable us to construct the solution for the interior source problem and to extend to solve the complicated problem for the multi span of the Timoshenko beam. Numerical results based on these two approaches are compared. The comparison in this study indicates that the normal mode method is more computationally efficient than the ray method except for very short time after the impact. The long-time transient responses are easily calculated using the normal mode method. It is shown that the average long-time transient response converges to the corresponding static value. The Timoshenko beam theory is more accurate than the Bernoulli–Euler beam theory because it includes shear and rotary inertia. This study also provides the slender ratio for which the Bernoulli–Euler beam can be used for the transient-response analysis of the displacement. Moreover, the resonant frequencies obtained from finite element calculation based on the three-dimensional model are compared with the results calculated using the Timoshenko beam and Bernoulli–Euler beam theories. It is noted in this study that the resonant frequency can be accurately determined by the Timoshenko beam theory if the slender ratio is larger than 100, and by the Bernoulli–Euler beam theory if the slender ratio is larger than 400.     

Transient wave analysis of a cantilever Timoshenko beam subjected to impact loading by Laplace transform and normal mode methods

This study applies two analytical approaches, Laplace transform and normal mode methods, to investigate the dynamic transient response of a cantilever Timoshenko beam subjected to impact forces. Explicit solutions for the normal mode method and the Laplace transform method are presented. The Durbin method is used to perform the Laplace inverse transformation, and numerical results based on these two approaches are compared. The comparison indicates that the normal mode method is more efficient than the Laplace transform method in the transient response analysis of a cantilever Timoshenko beam, whereas the Laplace transform method is more appropriate than the normal mode method when analyzing the complicated multi-span Timoshenko beam. Furthermore, a three-dimensional finite element cantilever beam model is implemented. The results are compared with the transient responses for displacement, normal stress, shear stress, and the resonant frequencies of a Timoshenko beam and Bernoulli–Euler beam theories. The transient displacement response for a cantilever beam can be appropriately evaluated using the Timoshenko beam theory if the slender ratio is greater than 10 or using the Bernoulli–Euler beam theory if the slender ratio is greater than 100. Moreover, the resonant frequency of a cantilever beam can be accurately determined by the Timoshenko beam theory if the slender ratio is greater than 100 or by the Bernoulli–Euler beam theory if the slender ratio is greater than 400.     

复合材料叠层梁和金属梁的固有振动特性

对根据三种梁理论得到的金属梁和复合材料叠层梁的固有振动特性进行了对比性的研究对常用的三种梁理论在弹性碰撞分析中的应用进行了分析和比较     

Asymptotic frequencies of various damped nonlocal beams and plates

A striking difference between the conventional local and nonlocal dynamical systems is that the later possess finite asymptotic frequencies. The asymptotic frequencies of four kinds of nonlocal viscoelastic damped structures are derived, including an Euler–Bernoulli beam with rotary inertia, a Timoshenko beam, a Kirchhoff plate with rotary inertia and a Mindlin plate. For these undamped and damped nonlocal beam and plate models, the analytical expressions for the asymptotic frequencies, also called the maximum or escape frequencies, are obtained. For the damped nonlocal beams or plates, the asymptotic critical damping factors are also obtained. These quantities are independent of the boundary conditions and hence simply supported boundary conditions are used. Taking a carbon nanotube as a numerical example and using the Euler–Bernoulli beam model, the natural frequencies of the carbon nanotubes with typical boundary conditions are computed and the asymptotic characteristics of natural frequencies are shown.     

An efficient procedure to find shape functions and stiffness matrices of nonprismatic Euler–Bernoulli and Timoshenko beam elements

Nonprismatic beam modeling is an important issue in structural engineering, not only for versatile applicability the tapered beams do have in engineering structures, but also for their unique potential to simulate different kinds of material or geometrical variations such as crack appearing or spreading of plasticity along the beam. In this paper, a new procedure is proposed to find the exact shape functions and stiffness matrices of nonprismatic beam elements for the Euler–Bernoulli and Timoshenko formulations. The variations dealt with here include both tapering and abrupt jumps in section parameters along the beam element. The proposed procedure has found a simple structure, due to two special approaches: The separation of rigid body motions, which do not store strain energy, from other strain states, which store strain energy, and finding strain interpolating functions rather than the shape functions which suffer complex representation. Strain interpolating functions involve low-order polynomials and can suitably track the variations along the beam element. The proposed procedure is implemented to model nonprismatic Euler–Bernoulli and Timoshenko beam elements, and is verified by different numerical examples.     

Enhancement of elementary beam theories in order to obtain exact solutions for elastic rectangular beams

Boley's method is utilized in order to show that the elementary Bernoulli–Euler beam theory can be enhanced such that exact solutions of the plane-stress theory of linear elasticity are obtained for force loaded rectangular beams. An equivalent enhancement is derived for the elementary Timoshenko theory of beams. The enhancement terms act analogous to thermal loadings; they follow from the force loading of the rectangular beam in an explicit form. The resulting boundary value problem of fourth order can be efficiently solved by means of symbolic computer codes. As an illustrative example, a redundant beam is studied, which is simply supported at one end, and which is clamped at the other end. Outcomes for three alternative clamped end boundary conditions are compared.     

NODAL DISPLACEMENT ACCURACY OF ONE-DIMENSIONAL FINITE ELEMENTS

研究高次杆单元和梁单元的节点位移精度问题.首先求出一端固支均匀杆和悬臂梁在任意次多项式形式分布载荷作用下的位移精确解,然后用二次杆单元、五次欧拉梁单元和三次铁木辛柯梁单元求得了节点位移.通过比较有限元解与精确解以及利用静力凝聚方法,发现一次以上杆单元、三次以上欧拉梁单元以及三次以上铁木辛柯梁单元都可以给出精确的端点位移.     

从一个怪例看细长梁理论的基本假定

分别采用欧拉和铁木辛柯梁理论分析了均匀分布力偶作用下的两端固支等截面匀质细长 梁, 并通过ABAQUS有限元分析了一个实例, 验证了铁木辛柯梁理论分析的结果. 对比证明在 这种载荷及边界条件下即使细长梁, 也必须考虑剪切效应的影响.     

A micro scale Timoshenko beam model based on strain gradient elasticity theory

A micro scale Timoshenko beam model is developed based on strain gradient elasticity theory. Governing equations, initial conditions and boundary conditions are derived simultaneously by using Hamilton's principle. The new model incorporated with Poisson effect contains three material length scale parameters and can consequently capture the size effect. This model can degenerate into the modified couple stress Timoshenko beam model or even the classical Timoshenko beam model if two or all material length scale parameters are taken to be zero respectively. In addition, the newly developed model recovers the micro scale Bernoulli–Euler beam model when shear deformation is ignored. To illustrate the new model, the static bending and free vibration problems of a simply supported micro scale Timoshenko beam are solved respectively. Numerical results reveal that the differences in the deflection, rotation and natural frequency predicted by the present model and the other two reduced Timoshenko models are large as the beam thickness is comparable to the material length scale parameter. These differences, however, are decreasing or even diminishing with the increase of the beam thickness. In addition, Poisson effect on the beam deflection, rotation and natural frequency possesses an interesting “extreme point” phenomenon, which is quite different from that predicted by the classical Timoshenko beam model.     

经典理论与一阶理论之间简支梁特征值的解析关系

利用Euler-Bernoulli梁理论(EBT)和Timoshenko梁理论(一阶理论,TBT)之间,梁的特征值问题在数学上的相似性,研究了不同梁理论之间特征值的关系。将特征值问题的求解转化为一个代数方程的求解,并导出了不同梁理论之间梁的特征值之间的精确解析关系。因此,只要已知梁的经典结果(临界载荷和固有频率),便很容易从这些关系中获得一阶梁理论下的相应结果。这些解析结果清楚地显示了横向剪切变形对经典结果影响的本质特点。另外,从这些关系中获得的含有剪切变形影响的结果,可以用于检验一阶理论下梁特征值数值结果的有效性、收敛性以及精确性等问题。