弹性连接旋转柔性梁动力学分析

采用Chebyshev 谱方法对考虑根部连接弹性的平面内旋转柔性梁动力学特性进行研究. 基于Gauss-Lobatto 节点与Chebyshev 多项式方法对柔性梁变形场进行离散,通过投影矩阵法施加固定及弹性连接边界条件. 利用Chebyshev 谱方法获得了系统固有频率和模态振型数值解,通过与有限元方法及加权残余法的比较,验证了方法的有效性. 分析了弹性连接刚度、角速度比率、系统径长比及梁的长细比等参数对系统固有频率及模态振型的影响. 研究发现:由于系统弯曲模态、拉伸模态的频率随各参数的变化规律不一致,将出现频率转向与振型转换现象;随着弹性连接刚度、角速度比率及系统径长比的增大,低阶弯曲模态频率增大并超过高阶拉伸模态频率,随着梁的长细比的增大,低阶拉伸模态频率增大并超过高阶弯曲模态频率.      

Geometrically nonlinear finite-element models for thin shells with geometric imperfections

A finite-element method to analyze the stress–strain state and stability of thin shells with geometric imperfections is proposed. An arbitrary curvilinear finite element with vector approximation of the displacement function is used. To solve the systems of nonlinear algebraic equations by iteration methods, linearized stiffness matrices of finite elements and residual and load vectors are formed. The stress–strain state of a thin-walled shell with real geometric imperfections under surface pressure and axial compression is analyzed. The effect of geometric imperfections on the critical combination of loads is evaluated     

3D numerical simulation of drilling residual stresses

Drilling can affect the integrity of the surface of a mechanical component and reduce its fatigue life. Thus, drilling parameters such as lubrication or drilling velocity must be optimized to ensure a satisfactory residual mechanical state of the hole surfaces. Unfortunately, experimental tests are time consuming and it is not easy to observe the cutting process because of the confinement of the drill zone. The literature does not exhibit any numerical simulation capable of simulating 3D thermomechanical phenomena in the drill zone for large depth holes. Therefore, residual stresses cannot be easily simulated by means of the sole drilling parameters. The aim of this article is to propose a new numerical approach to compute drilling residual stresses for large-depth holes. A first simulation is developed to simulate heat transfer by means of a 3D thermoviscoplastic simulation in a new Rigid-ALE framework allowing the use of large calculation time steps. Then, a time interpolation and a spatial projection are implemented to rebuild the Lagrangian thermal history of the machined component. Finally, a thermo-elastoplastic simulation is carried out to compute residual stresses in the final workpiece. In this paper, the method is applied to a 316L austenitic stainless steel in the case of an unlubricated hole. The computed residual stresses are compared to experimental measurements.     

Application of the random eigenvalue problem in forced response analysis of a linear stochastic structure

The random eigenvalue problem arises in frequency and mode shape determination for a linear system with uncertainties in structural properties. Among several methods of characterizing this random eigenvalue problem, one computationally fast method that gives good accuracy is a weak formulation using polynomial chaos expansion (PCE). In this method, the eigenvalues and eigenvectors are expanded in PCE, and the residual is minimized by a Galerkin projection. The goals of the current work are (i) to implement this PCE-characterized random eigenvalue problem in the dynamic response calculation under random loading and (ii) to explore the computational advantages and challenges. In the proposed method, the response quantities are also expressed in PCE followed by a Galerkin projection. A numerical comparison with a perturbation method and the Monte Carlo simulation shows that when the loading has a random amplitude but deterministic frequency content, the proposed method gives more accurate results than a first-order perturbation method and a comparable accuracy as the Monte Carlo simulation in a lower computational time. However, as the frequency content of the loading becomes random, or for general random process loadings, the method loses its accuracy and computational efficiency. Issues in implementation, limitations, and further challenges are also addressed.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

三维结构安定分析的直接算法

基于线性随动强化理论和Von. Mises屈服准则,对蒙板结构直接安定分析法进行了扩展,建立了结构的三维安定直接分析法。根据投射原理,推导出结构发生塑性安定的存在条件,便于调整控制加载步长和载荷历程。采用逐次增量加载方式,确定出背应力的偏移范围,克服了原始直接分析法不能获得安定极限的缺陷,并得到安定极限条件下结构中残余应力与应变的分布状况。该数值方法将弹塑性问题分解为弹性问题和特征应变决定的残余问题,节约计算时间,提高计算效率,将该算法应用于相关算例,并与有关数值结果相比较,验证了该算法的有效性。     

Application of conjugate-gradient-like methods to a hyperbolic problem in porous-media flow

This paper presents the application of a preconditioned conjugate-gradient-like method to a non-self-adjoint problem of interest in underground flow simulation. The method furnishes a reliable iterative solution scheme for the non-symmetric matrices arising at each iteration of the non-linear time-stepping scheme. The method employs a generalized conjugate residual scheme with nested factorization as a preconditioner. Model runs demonstrate significant computational savings over direct sparse matrix solvers.     

移动荷载作用下地基动力分析的有限元方法

通过对地基动力问题的基本方程进行变换,把基本方程变换到随荷载移动的运动坐标系中,通过加权残数法推导了相应的单元刚度矩阵,从而建立了移动问题的有限元格式,并发现移动荷载问题的单元刚度矩阵是对相应静力问题单元刚度矩阵的修正,在静力单元刚度矩阵的主对角元素上增加与移动速度有关的项,即可得到移动问题有限元的单元刚度矩阵,这样就将动力学问题转化为“拟静力”问题处理。文中用移动问题有限元方法计算了地基的动力响应,并与解析解进行了对比,以说明本方法具有较好的精度。     

Iterative solution techniques for the stokes and Navier-Stokes equations

The linear system arising from a Lagrange-Galerkin mixed finite element approximation of the Navier–Stokes and continuity equations is symmetric indefinite and has the same block structure as a system arising from a mixed finite element discretization of a Stokes problem. This paper considers the iterative solution of such a system, comparing the performance of the one-level preconditioned conjugate residual method for indefinite matrices with that of a more traditional two-level pressure correction approach. Asymptotic estimates for the amount of work involved in each method are given together with the results of related numerical experiments.     

Multi-fidelity surrogate reduced-order modeling of steady flow estimation

A multi-fidelity reduced-order model (ROM), which incorporates low-fidelity data to improve the prediction of high-fidelity results, is proposed for the reconstruction of steady flow field at different conditions. The spatial basis functions of low-fidelity and high-fidelity data, which are generated for all training sets are extracted separately by proper orthogonal decomposition. Then a surrogate model is used to construct mappings between the mode coefficients obtained from low-fidelity and the high-fidelity data. In the online stage, both the low-fidelity flow at the predicted state and the surrogate model are needed to predict the mode coefficients of the high-fidelity flow, and the high-fidelity flow field is subsequently reconstructed. This method differs from existing surrogate-based reduced-order modeling method because it allows the use of partial physical information for flow estimation, which is coming from the low-fidelity data instead of adopting a black-box mapping between system state and the projection coefficients. Numerical studies are presented for a lid-driven cavity problem and transonic flow past a NACA0012 airfoil. Two low-fidelity models, with either a coarse mesh or a lower numerical order, are considered respectively. Results show that the proposed multi-fidelity ROM predicts the flow field accurately and outperforms the traditional methods in both interpolated and extrapolated conditions.     

Nuclei of strain for bi-material elastic media with sliding interface

The nuclei of strain method, which was previously developed for a perfectly bonded bi-material elastic medium, is modified to obtain the fundamental solutions for the elastic state of a bi-material, isotropic elastic medium with sliding interface. The solutions are conveniently expressed in terms of strength vectors at object and image points which are related by the transmission and reflection matrices given in the paper.     

THE MATHEMATICAL MODELLING OF NEAR COAST SHALLOW WATER CIRCULATION

This paper presents a finite element method to solve the shallow water circulation problem numerically. Considering the Coriolis effect, bottom friction and eddy viscosity, the continuity equation and momentum equation are integrated vertically. Using Galerkin weighted residual method, the weak variational formulation is derived for the finite element analysis. The split-time method is applied for the numerical integration instead of iteration for nonlinear terms. Moreover, an artificial smooth approach is proposed to suppress the short wavelength noise. In order to save Computer storage units, a densed storage scheme is set up, where all the zero elements in large scaled and sparse matrices are excluded.     

Discrete time transfer matrix method for dynamics of multibody system with flexible beams moving in space

In this paper,by defining new state vectors and developing new transfer matrices of various elements moving in space,the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space.Formulations and numerical example of a rigidflexible-body three pendulums system moving in space are given to validate the method.Using the new method to study the dynamics of multi-rigid-flexible-body system moving in space,the global dynamics equations of system are not needed,the orders of involved matrices of the system are very low and the computational speed is high,irrespective of the size of the system.The new method is simple,straightforward,practical,and provides a powerful tool for multi-rigid-flexible-body system dynamics.     

FAST ALGORITHMS FOR DETERMINING POSITIVE DEFINITENESS OF LARGE SCALE MATRICES IN SOLID STRUCTURE ANALYSIS

对于结构稳定性分析中超大规模矩阵正定性判定,必须采用并行计算方法,传统方法如计算特征值、主子式行列式及LDLT等直接方法难以实现.本文提出了一些可适用于并行的迭代判定算法.借鉴力学系统中能量下降的思想,发展了一种判定矩阵正定性的新思路,即将矩阵的正定性判定问题转化为一个优化问题,并基于优化算法来判定矩阵的正定性.提出了基于最速下降法和共轭梯度法来进行矩阵正定性判定的算法.然后考虑到力学系统刚度矩阵的稀疏性和结构刚失稳状态的弱非正定性,提出可以先截超平面后解方程求驻值点的方法来判定弱非正定矩阵的正定性.为了保证对强非正定矩阵判定的准确性,本文提出可以高效混杂使用截平面法和共轭梯度法.数值实验结果表明,本文提出的算法具有准确性和高效性.对于强非正定矩阵而言,共轭梯度算法更加高效;而对于弱非正定矩阵,则是截平面法和混杂算法更加高效.这些算法都容易在机群上实现并行计算,能够快速判定大规模矩阵的正定性.     

Finite time chaos control for a class of chaotic systems with input nonlinearities via TSM scheme

This paper investigates nonsingular terminal sliding mode control for a class of uncertain systems with nonlinear inputs and its application in chaos control. When some of the system states are finite-time stable, the nonlinear items that coupled with these states may come into zeros in other subsystems. This will simplify the stability analysis of the whole system greatly. Compared with the traditional finite-time stabilization design method, the introduction of the terminal sliding mode can reduce the input dimensions. Only one control input is requested to realize chaos control of the Liu system when unmatched uncertainties and input nonlinearity coexist. The parameter matrices in the TSM can be determined through the solution of LMIS. Simulation results are given to demonstrate the effectiveness of the proposed method.     

Solution of free-boundary problems using finite-element/Newton methods and locally refined grids: Application to analysis of solidification microstructure

A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.     

阶梯压电层合梁的波动动力学特性

采用行波理论系统地研究了压电阶梯梁的自由振动分析以及强迫响应的分析方法. 基于分布 参数理论研究了压电阶梯梁的波传播特性，忽略柔性梁横向剪切和转动惯量的影响，给出了 梁的轴向和横向的简谐波解. 将压电阶梯梁离散化为单元，考虑压电片的刚度和质量的影响， 建立了节点散射模型. 应用位移连续和力平衡条件，推导了节点的波反射和波传递矩阵，在 此基础上，引入波循环矩阵的概念，给出波循环矩阵、波传递系数矩阵的确定方法. 应用波 循环矩阵可以有效地计算结构的固有频率. 另外，应用波传递系数研究了压电陶瓷作动器位 置对其驱动能力的影响. 得出两个主要结论：1)作动器靠近悬臂梁固定端将有较强的驱动 能力，悬臂梁边界反射行波产生弯曲消失波有利于增大压电波的模态传递系数；2)模态传 递系数与固有频率的灵敏度密切相关，波传递系数越大, 对应该处固有频率变化灵敏度越大. 另外，数值算例表明了行波方法比有限元方法具有更高的计算精度.     

Optimization of multibody systems based on the generalized-α projection method for DAEs

Efficient optimization strategy of multibody systems is developed in this paper. Augmented Lagrange method is used to transform constrained optimal problem into unconstrained form firstly. Then methods based on second order sensitivity are used to solve the unconstrained problem, where the sensitivity is solved by hybrid method. Generalized-α method and generalized-α projection method for the differential-algebraic equation, which shows more efficient properties with the lager time step, are presented to get state variables and adjoint variables during the optimization procedure. Numerical results validate the accuracy and efficiency of the methods is presented.     

约束层阻尼板动力学问题的半解析解

利用条形传递函数方法(SDTFM)得到了约束层阻尼(CLD)板动力学问题的半解析解.首先对CLD板沿纵向离散成多个条形单元,基于Hamilton原理推导出条形单元的刚度矩阵和质量矩阵,仿照有限元法组集得到系统的总刚度矩阵和总质量矩阵.经Laplace变换后引入状态向量,采用分布参数传递函数方法在状态空间内建立CLD板的控制方程并进行求解.最后以对边固支和悬臂CLD板为例,得到了板的动力学特性和频响曲线,并与NASTRAN或相关文献结果进行了比较,吻合良好,验证了该方法的有效性.从推导过程和算例可以看出,该方法所需的单元数目少,获得的是半解析解,计算效率高且准确可靠.     

Continuous-Time Bilinear System Identification

The objective of this paper is to describe a new method for identification of a continuous-time multi-input and multi-output bilinear system. The approach is to make judicious use of the linear-model properties of the bilinear system when subjected to a constant input. Two steps are required in the identification process. The first step is to use a set of pulse responses resulting from a constant input of one sample period to identify the state matrix, the output matrix, and the direct transmission matrix. The second step is to use another set of pulse responses with the same constant input over multiple sample periods to identify the input matrix and the coefficient matrices associated with the coupling terms between the state and the inputs. Numerical examples are given to illustrate the concept and the computational algorithm for the identification method.