转动盘一片系的循环对称模综

本文考察在任意恒定转速下叶片/轮盘系的固有模态分析,计及叶片预扭所引起的大变形效应.把群论算法用于子结构技术,通过联立迭代求出叶片/轮盘系在离心场中的平衡位置和应力分布,进而确定初应力刚度矩阵、大变形刚度矩阵,以及相应的动频模态.由于本方法考虑了叶片大变形的影响,且可消去子结构间的全体内连自由度。从而在确保精度的前提下,显著地压缩了计算的规模.实践证明本文提出的子结构联立迭代法有很好的收敛性,它使计算量有巨大的节省.     

Convergence predictions for aeroelastic calculations of tuned and mistuned bladed disks

Mistuning changes the dynamics of bladed disks significantly. Frequency domain methods for predicting the dynamics of mistuned bladed disks are typically based on iterative aeroelastic calculations. Converged aerodynamic stiffness matrices are required for accurate aeroelastic results of eigenvalue and forced response problems. The tremendous computation time needed for each aerodynamic iteration would greatly benefit from a fast method of predicting the number of iterations needed for converged results. A new hybrid technique is proposed to predict the convergence history based on several critical ratios and by approximating as linear the relation between the aerodynamic force and the complex frequencies (eigenvalues) of the system. The new technique is hybrid in that it uses a combined theoretical and stochastic/computational approach. The dynamics of an industrial bladed disk is investigated, and the predicted convergence histories are shown to match the actual results very well. Monte Carlo simulations using the new hybrid technique show that the aerodynamic ratio and the aerodynamic gradient ratio are the two most important factors affecting the convergence history.     

一种大型结构特征值问题的并行解法

本文提出了一种求解大型结构固有频率与模态的并行解法。该方法在子空间迭代过程中,利用子结构刚度阵和质量阵并行进行凝聚,求得下一次的迭代基矢量,直到收敛。该算法在西安交通大学ELXS1-6400并行机上程序实现,计算结果表明能大幅度节省计算时间,同时也有效地节省了内存。     

The worst response of mistuned bladed disk system using neural network and genetic algorithm

The objective of this paper is to develop an integrated approach using artificial neural networks (ANN) and genetic algorithms (GA) for predicting the worst response of mistuned bladed disk. ANN is used to predict the responses of bladed disk system which are used further in evaluation of fitness and constraint violation in GA process. A multilayer back-propagation neural network is trained with the results obtained from finite element model for different bladed disk configurations. Subsequently, GA is employed for arriving at optimum configuration of the bladed disk system by maximizing the blade responses. By integrating ANN with GA, the computational time required for obtaining optimal solution could be reduced substantially. The efficacy of this approach is demonstrated by carrying out studies on mistuned bladed disk systems for different sets of mistuning parameters, namely mistuning in modulus of elasticity and length of blades. Finally, the effect of adding shroud at the tip of blades in reducing the maximum response of the bladed disk system was investigated.     

Localization in Nonlinear Mistuned Systems with Cyclic Symmetry

In forced systems with cyclic symmetry localization can occur due toparameter uncertainties. Often, Monte-Carlo simulations are used to findregions, where the system response is sensitive to parameteruncertainties. These simulations require a large computation time.Therefore, an approximate method to calculate the envelopes of thefrequency response functions is developed in this paper. An example of anonlinear system with cyclic symmetry is a bladed disk assembly withfriction dampers. Friction dampers can be installed underneath the bladeplatforms of turbine blades. Due to dry friction and the relative motionbetween blades and dampers, energy is dissipated, which results in areduction of blade vibration amplitudes. By optimizing the mass of thefriction dampers, the best damping effects are obtained, which lead toan increase in the reliability of the turbine. In this paper, thecalculated response of a mistuned bladed disk assembly with frictiondampers is discussed. An approximate method is developed to calculatethe envelopes of the corresponding frequency response function forstatistically varying eigenfrequencies of the blades. Regions wherelocalization can occur with a high probability, are calculated by thismethod.     

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

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

基于模板的子结构多分辨率拓扑优化

多分辨率拓扑优化(multi-resolution topology optimization, MTOP)方法将有限元网格和密度网格解耦, 采用较粗的网格(超单元)进行有限元分析, 从而大大降低了拓扑优化过程中的结构分析成本. 但MTOP方法每次迭代都需要根据超单元内的平均密度计算有限元单刚, 不仅精度不够且在过滤半径较小的情况容易出现棋盘格现象和QR模式. 为解决相应问题, 本文将超单元视为子结构, 通过静态凝聚得到超单元刚度阵, 并进一步根据拓扑优化过程中子结构的密度分布特征组建了其模板库, 从而省去了超单元单刚的重复计算, 显著提高了MTOP方法的分析精度, 有效抑制了数值不稳定现象.      

A variational differential quadrature solution to finite deformation problems of hyperelastic shell-type structures: a two-point formulation in Cartesian coordinates

A new numerical approach is presented to compute the large deformations of shell-type structures made of the Saint Venant-Kirchhoff and Neo-Hookean materials based on the seven-parameter shell theory. A work conjugate pair of the first Piola Kirchhoff stress tensor and deformation gradient tensor is considered for the stress and strain measures in the paper. Through introducing the displacement vector, the deformation gradient, and the stress tensor in the Cartesian coordinate system and by means of the chain rule for taking derivative of tensors, the difficulties in using the curvilinear coordinate system are bypassed. The variational differential quadrature (VDQ) method as a pointwise numerical method is also used to discretize the weak form of the governing equations. Being locking-free, the simple implementation, computational efficiency, and fast convergence rate are the main features of the proposed numerical approach. Some well-known benchmark problems are solved to assess the approach. The results indicate that it is capable of addressing the large deformation problems of elastic and hyperelastic shell-type structures efficiently.     

界面连接刚度参数辨识的子结构分析法

以试验模态参数为基础,提出一种通过特征方程反问题辨识子结构界面连接刚度参数的子结构分析法。新方法以子结构动柔度矩阵特征方程为基础,建立求解界面结点内力和位移的方程,从而由子结构内部结点可测自由度上的位移用广义逆理论估计界面结点内力和位移。并通过迭代修正内部结点可测自由度上的试验值,以提高界面内力和位移的估计精度。最后通过连接子结构刚度矩阵建立的平衡方程求解相应的刚度参数。文中以太阳电池阵板间铰链副刚度参数辨识为例,将铰链副简化为两端结点各有6个自由度的弹簧连接元,考虑到自由度之间的耦合,推导了连接元的刚度矩阵。用上述方法辨识了铰链副6个自由度的刚度参数,得到满意的辨识结果。     

Fluid–structure interaction analysis of flexible turbomachinery

A method for the performance computation of an expandable-impeller pump is developed and validated. Large deformations of the highly flexible pump impellers result in a strong coupling between the impeller and fluid flow. The computational method therefore requires simultaneous solution of fluid flow and structural response. OpenFOAM provides the flow and mesh motion solvers and is coupled to an author-developed structural solver in a tightly coupled approach using a fixed-point iteration. The structural deformations are time-dependent because the material exhibits stress relaxation. The time-constant of the relaxation, however, is very large, thereby allowing quasi-steady simulations. A water-tunnel test of a viscoelastic hydrofoil is employed to validate the solver. Simulations of the test problem show good agreement with the experimental results and demonstrate the need for several sub-iterations of the solver even for the quasi-steady simulations.     

基于三阶剪切变形理论的悬臂输流管道自由振动

采用三阶剪切变形理论,结合有限元法研究了悬臂输流管道的自由振动问题．利用虚功原理建立了输流管系统的有限元方程,同时将悬臂端弹性支承以势能的形式引入到系统方程中,求解了系统前三阶的复频率．分别探讨了流体速度和弹簧刚度对系统复频率实部和虚部的影响,重点分析了弹簧刚度与前三阶固有频率间的关系．在弹性支承刚度为零的特例下,对比了本文结果与Timoshenko梁理论的结果,证明了本文方法的可靠性．研究发现系统固有频率的实部恒为负值,表明一端带有弹性支承的约束形式有利于提高悬臂输流管道自由振动的稳定性;流体的流动对管道振动起到了阻尼作用,在流动速度足够大的情况下,各阶振动固有频率均趋于零;当弹簧刚度为无穷大,且流体速度足够大时,输流管道将发生失稳．     

弹塑性大变形分析的一致性高阶无单元伽辽金法

采用无单元伽辽金法求解弹塑性大变形问题。充分利用无单元法易于建立高阶近似函数的优点，位移采用二阶移动最小二乘近似。在更新拉格朗日方法的框架下，通过对控制方程弱形式的线性化建立了内力率的表达式，并区分为材料和几何两部分。采用Hughes-Winget算法更新应力，建立了Newton-Raphson迭代求解所需的一致切线刚度阵。刚度阵的数值积分采用近来针对小变形分析建立的二阶一致三点积分格式QC3（Quadratically Consistent 3-point integration scheme）。数值结果证明了本文方法分析弹塑性大变形问题的有效性和优越性。     

非线性复合材料的一种杂交应力有限单元方法

建立了非线性复合材料模型的杂交应力有限元方法,并在材料主坐标系下提出直接方法计算单元非线性应力场,然后由此计算单元切线刚度矩阵和剩余载荷并转换到整体坐标系下,利用Newton-Raphson方法进行结构的位移迭代。在Hahn-Tsai非线性复合材料杂交元分析中,由位移和应力方程所导出求解单元非线性应力场的简单迭代法是条件收敛的,对较大载荷当迭代位移增加到一定程度以后无法得到应力收敛解。但是,利用本文提出的直接法由于完全避免了非线性应力场迭代,不仅很好地解决了这一问题,而且极大地提高了计算效率。数值算例说明该方法是确实有效的。     

无额外自由度广义有限元非线性分析

将无额外自由度的广义有限元法由线弹性分析扩展到弹塑性大变形分析.局部强化函数的构建依赖于已有节点,不引入额外自由度,避免了线性相关性问题.在更新拉格朗日框架下,通过控制方程弱形式的线性化推导得到了节点内力的率形式,并分为材料和几何两部分.考虑超弹性和亚弹-塑性两种材料模型,采用Newton-Raphson迭代求解,给出了相关的一致切线刚度阵.三个典型算例的数值结果表明,本文发展的非线性无额外自由度广义有限元方法不仅能够准确求解超弹性和弹塑性大变形问题,同时相比于传统的线性有限元方法具有更高的精度.本文工作进一步拓宽了无额外自由度广义有限元方法的应用领域.     

Hole-Drilling Residual Stress Measurement with Artifact Correction Using Full-Field DIC

A full-field, multi-axial computation technique is described for determining residual stresses using the hole-drilling method with DIC. The computational method takes advantage of the large quantity of data available from full-field images to ameliorate the effect of modest deformation sensitivity of DIC measurements. It also provides uniform residual stress sensitivity in all in-plane directions and accounts for artifacts that commonly occur within experimental measurements. These artifacts include image shift, stretch and shear. The calculation method uses a large fraction of the pixels available within the measured images and requires minimal human guidance in its operation. The method is demonstrated using measurements where residual stresses are made on a microscopic scale with hole drilling done using a Focused Ion Beam – Scanning Electron Microscope (FIB-SEM). This is a very challenging application because SEM images are subject to fluctuations that can introduce large artifacts when using DIC. Several series of measurements are described to illustrate the operation and effectiveness of the proposed residual stress computation technique.     

考虑塑性和刚度不匹配的焊接残余应力估算

现有残余应力计算方法未能考虑材料塑性变形和焊接接头刚度不匹配的影响,使得焊接残余应力计算结果和实际残余应力存在较大偏差.在2219-T87铝合金钨极氩弧焊焊接头残余应力测试基础上,提出一种基于非线性有限元和材料弹性模量分区的残余应力—释放应变曲线的残余应力计算方法,研究了材料塑性变形和接头刚度不匹配对焊接残余应力计算的影响.结果表明,焊接接头中非均质材料塑性不匹配可以引起对于残余应力计算的较大误差;材料塑性变形对残余应力的影响大于接头刚度不匹配对残余应力的影响.所提出方法修正了传统方法在焊接接头的残余应力计算中由于未考虑接头非均质材料塑性不匹配而引起的误差.     

自适应摄动法在失谐叶盘受迫振动分析中的应用

针对自适应摄动法在失谐叶盘受迫振动分析中,选择接近激振力频率的相应模态数目难于控制的问题,基于质量-弹簧模型,提出一种自动选择模态的方法。该方法可根据激振力频率与基系统固有频率集的关系,自动选择最接近激振力频率的相应模态个数。计算结果显示,对失谐叶盘系统,特别是在弱耦合情况下,该方法比人为选择模态更为可靠,可以使用较少机时,获得较高的求解精度。     

提高多体系统离散时间传递矩阵法计算精度的研究

深入探讨多体系统离散时间传递矩阵法对平面、空间刚体-光滑铰多体系统运动响应的研究。提出提高该方法计算精度和计算稳定性的方法，导出相应的多端刚体传递矩阵。设计了角坐标的迭代循环，该变量不必采用近似形式，因而提高了计算精度和计算稳定性。其它方法有：1)增加泰勒展开式的高阶项；2)合理选择将速度、加速度表示为位移的线性函数的方法(本文简称线性化方法)；3)合理确定迭代循环的初值等。对4种多体系统进行了计算机仿真研究，表明本文提出的方法有效。     

Analytical formulation of the Jacobian matrix for non-linear calculation of the forced response of turbine blade assemblies with wedge friction dampers

A fundamental issue in turbomachinery design is the dynamical stress assessment of turbine blades. In order to reduce stress peaks in the turbine blades at engine orders corresponding to blade natural frequencies, friction dampers are employed. Blade response calculation requires the solution of a set of non-linear equations originated by the introduction of friction damping.

Such a set of non-linear equations is solved using the iterative numerical Newton–Raphson method. However, calculation of the Jacobian matrix of the system using classical numerical finite difference schemes makes frequency domain solver prohibitively expensive for structures with many contact points. Large computation time results from the evaluation of partial derivatives of the non-linear equations with respect to the displacements.

In this work a methodology to compute efficiently the Jacobian matrix of a dynamic system having wedge dampers is presented. It is exact and completely analytical.

The proposed methods have been successfully applied to a real intermediate pressure turbine (IPT) blade under cyclic symmetry boundary conditions with underplatform wedge dampers. Its implementation showed to be very effective, and allowed to achieve relevant time savings without loss of precision.