气体动理学统一算法的隐式方法研究

目前的气体动理学统一算法(unified gas kinetic scheme, 简称UGKS) 在求解高速流动问题时的计算效率,难以满足求解复杂工程问题的需求. 为了提高该算法的计算效率, 本文对模型方程的对流项和碰撞项进行了隐式处理, 并针对UGKS 界面通量与演化时间相关的特点, 引入了演化时间平均界面通量, 通过对控制方程矩阵进行近似LU 分解(lower-upper decomposition), 实现了隐式UGKS. 不同来流马赫数的圆柱绕流算例测试表明, 只要演化时间选取得当, 隐式方法可以得到与显式方法完全相同的结果, 且计算效率可以提高1~2 个量级.      

Dynamic modeling and computational efficiency analysis for a spatial 6-DOF parallel motion system

A common approach for simplification analysis of complex dynamic model is presented, and the simplified dynamic model of a spatial 6-DOF parallel motion system with high computational efficiency is proposed for high real time control. By using Kane method, the full dynamic model of a spatial 6-DOF parallel motion system viewed as 13 rigid bodies is built. With rigid body decomposition, the full dynamic model is separated into several parts firstly, and then some separated parts are further divided into many dynamic components in terms of motion separation and the relationship with acceleration or velocity. The contribution of each dynamic term is analyzed for a specified spatial 6-DOF parallel motion system, and the simplified model is derived. Comparing with full dynamic model, the simplified error is analyzed, and the computational efficiency of the simplified model is discussed in a real-time industrial computer. The simplified strategy is confirmed in simulation. The simplified error is less than 8%, the simplified model can improve the computational efficiency by more than 70%, and the execution time is less than 0.1 ms, which can achieve the requirements of high real time control. The numerical results illustrate that the proposed approach is feasible and effective for simplification analysis of dynamics and the derived simplified dynamic model can be used in real-time control system with small simplified error.     

Three-dimensional unsteady flow simulations: Alternative strategies for a volume-averaged calculation

This work builds on a SIMPLE-type code to produce two numerical codes of greatly improved speed and accuracy for solution of the Navier–Stokes equations. Both implicit and explicit codes employ an improved QUICK (quadratic upstream interpolation for convective kinematics) scheme to finite difference convective terms for non-uniform grids. The PRIME (update pressure implicit, momentum explicit) algorithm is used as the computational procedure for the implicit code. Use of both the ICCG (incomplete Cholesky decomposition, conjugate gradient) method and the MG (multigrid) technique to enhance solution execution speed is illustrated. While the implicit code is first-order in time, the explicit is second-order accurate. Two- and three-dimensional forced convection and sidewall-heated natural convection flows in a cavity are chosen as test cases. Predictions with the new schemes show substantial computational savings and very good agreement when compared to previous simulations and experimental data.     

An explicit time domain solution for ground stratum response to harmonic moving load

Based on the thin layer method originally proposed in frequency domain, an explicit time domain semi–analytical solution has been developed for simulating three-dimensional layered ground responses to harmonic moving loads. The Fourier–Laplace transforms were applied to derive the transformed solution that satisfied the boundary conditions of horizontal infinities. The eigenvalue decomposition was performed with respect to Laplace parameter to express the ground motion corresponding to the eigenmodes. The formulation for each eigenmode incorporating the moving load expression was transformed back into time domain analytically, and the global system responses were given by means of the general mode superposition method. The proposed explicit time domain solution is suitable for studying various types of moving load acting on or inside the ground. In this paper a moving harmonic load with rectangular distribution was adopted to demonstrate the ground response simulation. Two illustrative examples for moving load with speeds below or above the ground Rayleigh wave velocity were presented to test the computational accuracy and efficiency of the proposed approach. A parametric study was also performed to investigate the influences of soil properties on the ground responses.The project is partially supported by the National Natural Science Foundation of China (50538010). The English text was polished by Yunming Chen.     

Use of Cholesky Coordinates and the Absolute Nodal Coordinate Formulation in the Computer Simulation of Flexible Multibody Systems

In a previous publication, procedures that can be used with the absolute nodal coordinate formulation to solve the dynamic problems of flexible multibody systems were proposed. One of these procedures is based on the Cholesky decomposition. By utilizing the fact that the absolute nodal coordinate formulation leads to a constant mass matrix, a Cholesky decomposition is used to obtain a constant velocity transformation matrix. This velocity transformation is used to express the absolute nodal coordinates in terms of the generalized Cholesky coordinates. The inertia matrix associated with the Cholesky coordinates is the identity matrix, and therefore, an optimum sparse matrix structure can be obtained for the augmented multibody equations of motion. The implementation of a computer procedure based on the absolute nodal coordinate formulation and Cholesky coordinates is discussed in this paper. Numerical examples are presented in order to demonstrate the use of Cholesky coordinates in the simulation of the large deformations in flexible multibody applications.     

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

圆形地下连续墙场地地震反应的简化计算方法

首先将水平地震荷载作用沿环向展开为傅立叶级数,利用沿圆周法向和切向两个方向的正交性,再将水平地震荷载转化为一致横向荷载组合,使三维求解问题转变为一系列旋转子午面上二维问题的叠加,进而提出圆形地下连续墙场地地震反应的简化计算方法。之后,采用等效线性化方法来考虑土体非线性性质,进一步将非线性地震反应问题转化为线性问题求解。...     

大跨度桥梁空间脉动风场的计算机模拟

针对现有Deodatis方法模拟大跨度桥梁空间脉动风场中存在的计算量问题,通过对谱分解矩阵引入插值近似,减少谱分解的次数,从而提高该谐波合成法的计算效率,并节省内存花费,实现了对三维空间脉动风场的有效模拟。改进方法模拟的脉动风速样本仍保持各态历经性,且逐渐收敛到目标功率谱。用改进的Deodatis方法模拟了润扬长江悬索桥桥面主梁上作用的纵向脉动风速。结果表明,该改进措施对Deodatis方法的应用效果非常明显,改进的Deodatis方法模拟脉动风速样本的相关函数与目标相关函数均吻合良好。尽管改进后的Deo-datis方法对谱分解矩阵采用了插值近似,但模拟的随机风速样本仍具有很好的精度。     

Computer Implementation of the Absolute Nodal Coordinate Formulation for Flexible Multibody Dynamics

Deformable components in multibody systems are subject to kinematic constraints that represent mechanical joints and specified motion trajectories. These constraints can, in general, be described using a set of nonlinear algebraic equations that depend on the system generalized coordinates and time. When the kinematic constraints are augmented to the differential equations of motion of the system, it is desirable to have a formulation that leads to a minimum number of non-zero coefficients for the unknown accelerations and constraint forces in order to be able to exploit efficient sparse matrix algorithms. This paper describes procedures for the computer implementation of the absolute nodal coordinate formulation' for flexible multibody applications. In the absolute nodal coordinate formulation, no infinitesimal or finite rotations are used as nodal coordinates. The configuration of the finite element is defined using global displacement coordinates and slopes. By using this mixed set of coordinates, beam and plate elements can be treated as isoparametric elements. As a consequence, the dynamic formulation of these widely used elements using the absolute nodal coordinate formulation leads to a constant mass matrix. It is the objective of this study to develop computational procedures that exploit this feature. In one of these procedures, an optimum sparse matrix structure is obtained for the deformable bodies using the QR decomposition. Using the fact that the element mass matrix is constant, a QR decomposition of a modified constant connectivity Jacobian matrix is obtained for the deformable body. A constant velocity transformation is used to obtain an identity generalized inertia matrix associated with the second derivatives of the generalized coordinates, thereby minimizing the number of non-zero entries of the coefficient matrix that appears in the augmented Lagrangian formulation of the equations of motion of the flexible multibody systems. An alternate computational procedure based on Cholesky decomposition is also presented in this paper. This alternate procedure, which has the same computational advantages as the one based on the QR decomposition, leads to a square velocity transformation matrix. The computational procedures proposed in this investigation can be used for the treatment of large deformation problems in flexible multibody systems. They have also the advantages of the algorithms based on the floating frame of reference formulations since they allow for easy addition of general nonlinear constraint and force functions.     

用虚拟激励法求解非比例阻尼线性体系的非平稳随机地震响应

应用复振型分解方法,将非比例阻尼线性体系在地震作用下的动力方程求解问题转化为若干个广义复振子的求解与叠加问题。通过假定地震地面运动为一零均值的非平稳随机激励,应用虚拟激励法原理,推导得到了广义复振子动力坐标计算的一般公式,进而得到了非比例阻尼线性体系非平稳随机地震响应计算的一般解答。由于可以选择少量共轭复振型的影响进行计算,对于大型复杂非比例阻尼结构,其随机地震响应计算工作量可以大幅度减小。算例证实了这种方法的可靠性及可行性。     

改进的泡泡布点方法研究

泡泡布点方法能够在复杂区域内生成高质量的节点集，但是其计算效率仍有待提高，为此本文做了两方面的改进。一是种子填充算法思想应用于节点的初始布置中，从而生成数目合适的节点集，省去了原布点方法中的节点增删过程，有效地节约模拟时间；二是简化了泡泡运动模拟的控制方程，使得运动模拟更加简单。数值算例表明，改进的泡泡布点方法生成的网格平均质量均高于0.94,且计算效率相比原布点方法提高90%以上。因此，改进的泡泡布点方法是一种高质量、高效率且适宜求解大规模问题的节点布置方法。     

非比例阻尼线性体系平稳随机地震响应计算的虚拟激励法

应用复振型分解方法,将非比例阻尼线性体系在地震作用下的动力方程求解问题转化为若干个广义复振子的求解与叠加问题。通过假定地震地面运动为一零均值的平稳随机激励,应用虚拟激励法原理,推导得到了广义复振子动力坐标的解析计算公式,进而得到了以复振型为基础的非比例阻尼线性体系随机地震响应计算的一般实数解析解答。算例证实了这种方法的可靠性及高效率。     

Comparative study on multibody vehicle dynamics models based on subsystem synthesis method using Cartesian and joint coordinates

The subsystem synthesis method has been developed in order to improve computational efficiency for a multibody vehicle dynamics model. Using the subsystem synthesis method, equations of motion of the base body and each subsystem can be solved separately. In the subsystem synthesis method, various coordinate systems can be used and various integration methods can be applied in each subsystem, as long as the effective mass matrix and the effective force vector are properly produced. In this paper, comparative study has been carried out for the subsystem synthesis method with Cartesian coordinates and with joint relative coordinates. Two different integration methods such as an explicit integrator and an explicit implicit integrator are employed. In order to see the accuracy and computational efficiency from the different models based on the different coordinate systems and different integration methods, a rough terrain run simulations has been carried out with a 6 × 6 off-road multibody vehicle model.     

Computation of frequencies and damping ratios for adaptive structures

In this paper, a new technique of state vector approach to solving damping ratios of adaptive structures is proposed, by which the state matrix can be analytically conducted the Cholesky decomposition, providing that the mass matrix is derived by lumped method. In this case, not only the computional accuracy is raised, but also the numerical operations are reduced. Finally, some numerical results are presented to show that the solution method is simple and efficient. Foundation item: the National 863 Program Foundation of China(863-2-4-4-4)     

建筑结构显式拟静力推覆分析方法研究

基于结构的精细化非线性有限元模型和显式拟静力求解方法,提出了结构拟静力推覆分析方法(EQPA);为控制结构惯性效应,给出了加载模式和基于能量的评价方法。在自主研发的CPU+GPU并行有限元软件中开发实现EQPA,可有效克服大规模复杂结构隐式Push-Over分析中计算量大和收敛性差的问题。使用EQPA对某超高层剪力墙结构的抗震性能进行分析,并与三组人工地震动的简化动力增量分析(sIDA)对比。结果表明,EQPA的结构惯性效应得到有效抑制;EQPA耗时仅为sIDA耗时的4.97%;两种分析方法的结果有一定差异,但两者的结构能力谱曲线、层间变形规律、薄弱楼层和构件损伤规律较为接近。     

大跨度空间网格结构多点输入反应谱计算方法的研究

根据大跨度空间网格结构的特点,忽略拟静力反应项与动力反应项的耦合项,提出了多点输入反应谱计算的简化方法.该方法所表示的结构地震动反应主要由拟静力响应和动力响应两部分组成,其中拟静力响应是由地面各点输入位移不一致引起的;动力响应部分主要是由结构在地面加速度作用下引起的,它在形式上与传统的一致输入反应谱保持一致,但考虑了多点输入的影响.结果表明,简化方法是可行的,它一方面减少了计算工作量,另一方面又保证了计算精度.     

基于虚拟激励法的大跨度桥梁非平稳随机振动分析闭合解

针对大跨度结构考虑地面空间运动的非一致地震响应分析问题,结合虚拟激励法,应用傅里叶分析建立了结构非平稳随机振动响应演变功率谱分析的频域方法。建立的方法完全基于频域执行,给出了响应演变功率谱的闭合解表达式。由于实现了确定性调制过程与随机过程的有效分离,应用离散傅里叶变换进行计算不需要较高的采样分析频率就可以获得较好精度的数值结果。数值算例研究了某斜拉桥的考虑地面运动空间效应非平稳随机地震响应,与通常时频分析方法进行对比,验证了本文频域方法的正确性和有效性。     

设防烈度下非比例阻尼结构地震随机响应峰值区间估计方法

建立与抗震规范设防水准相一致的地震动随机模型,给出了非比例阻尼结构体系地震反应运动方程,推导了非比例阻尼结构体系分析的复振型时域闭合解和地震随机响应的闭合解.通过随机响应的峰值因子和变异系数获得了非比例阻尼结构体系的地震随机响应峰值均值和方差,提出了设防烈度下非比例阻尼结构体系随机响应峰值区间估计方法.计算31条实际地震记录的复振型时域闭合解,并与地震随机响应的闭合解进行了比较,算例结果表明随机响应峰值区间估计方法的合理性和有效性.地震响应区间分析方法将是工程抗震设计的一个重要方法.     

带约束多体系统动力学方程的隐式算法

研究了带约束多体系统隐式算法,用子矩阵的形式推导出了多体系统正则方程的Ｊａｃｏｂｉ矩阵,它适用于多种隐式算法并给出了隐式Ｒｕｎｇｅ－Ｋｕｔｔａ算法,最后用一算例表明了隐式算法的计算效率和精度明显优于算法。     

一种结构动力响应分析的快速高阶算法

为了提高基于高阶格式的结构动力响应微分求积分析方法的计算效率,发展了一种求解动力方程的快速算法．利用微分求积原理将结构动力方程转化为标准Sylvester方程的形式,通过对系数矩阵进行矩阵分解,进而将动力响应Sylvester方程化为一系列标准线性方程组,采用相关成熟算法求解这些线性方程组后即可获得结构动力时程响应的全部解答．结构动力响应微分求积分析方法为高阶数值方法,一步计算可以获得多个时点处的动力响应．基于本文快速算法,不必直接对矩阵方程进行求解．数值算例表明,本文快速算法能够准确地计算出结构动力响应,具有数值精度高、收敛性好的优点．