Symbolic generation of large multibody system dynamic equations using a new semi-explicit Newton/Euler recursive scheme

Summary The aim of this paper is to show that multibody systems with a large number of degrees of freedom can be efficiently modelled, taking conjointly advantage of a recursive formulation of the equations of motion and of the symbolic generation capabilities.Recursive schemes are widely used in the field of multibody dynamics since they avoid the explosion of the number of arithmetical operations in case of large multibody models. Within the context of our field of applications (railway dynamics simulation), explicit integration schemes are still prefered and thus oblige us to compute the generalized accelerations at each time step. To achieve this, we propose a new formulation of the well-known Newton/Euler recursive method, whose efficiency will be compared with a so-called O(N) formulation.A regards the symbolic generation, often decried due to the size of the equations in case of large systems, we have recently implemented recursive multibody formalisms in the symbolic programme ROBOTRAN [1]. As we shall explain, the recursive nature of these formalisms is particularly well-suited to symbolic manipulation.All these developments have been successfully applied in the field of railway dynamics, and in particular allowed us to analyse the dynamic behaviour of several railway vehicles. Some typical results related to a completely non-conventional bogie will be presented before concluding.     

船舶密闭舱室火灾过程及其模拟计算

密闭舱室火灾是一个三维非定常过程，压力升高是这一过程的重要特征。本文采用质量守恒的方法进行压力校正，建立了计算密闭舱室火灾过程的数学物理模型，数值计算采用非均匀交错网格，结果与实验和理论分析相符较好。     

随机选取法计算化学反应流爆轰过程引起的参数涨落

利用ＲＣＭ（随机选取法）解反应流爆轰问题时，掌握影响流动参数随机涨落的主要因素是选择适当计算参数、保证得到正确结果的关键。各种流动参数的随机涨落相互关联，只需要搞清楚炸药未反应质量分数的随机涨落并做出分析，其它流动参数涨落的情况也就清楚了。给出了影响未反应质量分数涨落的因数，同时还给出了具体算例，检验所得结论的正确性。     

采动诱发岩体移动破坏过程数值模拟研究

本文基于岩石介质的宏观非线性主要是由非均质性和各向异性造成的, 应用新的数值计算软件RFPA(2D), 对采动引起岩体失稳破坏的全过程进了数值模拟研究。     

L-leap-accelerating the stochastic simulation of chemically reacting systems

Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems,in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap.The L-leap method can better satisfy the leap condition.Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.     

八节点等参元滑动接触算法的改进

分析认为采用八节点等参元计算滑动接触问题时,往往存在较大的接触应力计算误差的主要原因在于这种单元具有两种不同的节点。为此,本文提出了双模拟接触边界法,即分别采用两种不同的节点,形成两个模拟接触边界,并限制“从接触边界”上的节点只能和“主接触边界”上其同类节点形成的模拟接触边界形成接触关系,从而有效的减小了采用八节点等参元计算滑动接触问题时的误差。     

NUMERICAL SIMULATION OF OIL MIGRATION-ACCUMULATION OF MULTILAYER AND ITS APPLICATION

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｏｎｅｏｆｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ.Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃ…     

单相旋转型驻波超声电机的数学模型及仿真

尽管单相驱动的旋转型驻波超声电机的原理样机早巳研制成功，但对这类电机的运动机理的研究却一直局限在定性分析上。本文将从能量角度，根据Hamilton变分原理，定、转子摩擦界面采用纯滑动模型建立了该种电机的数学模型；并根据此模型，用MATLAB语言作了计算机仿真，得到了该种电机的开、关过程动态响应。     

CAD软件在工程地质三维建模中的应用

如何快速、准确地建立地质体的三维模型一直是众多岩土工程数值模拟工作者所面临的难题。虽然三维地学模拟软件具有很好的三维地质建模能力,但是由于数据结构的差异,采用他们现行三维地学模拟软件建立的地质模型难以导入数值模拟分析软件中,以为相应工程问题的数值模拟服务。目前,随着各种CAD、CAM软件行业的的飞速发展,涌现出了许多优秀的三维建模软件,而且这些软件大都与现行数值分析软件有着良好的数据接口功能。据此,本文提出了采用现行CAD软件来建立工程地质体的三维模型,使得建立的模型达到既"可视"又"可算"的目的。将其应用于云南某高速公路边坡的三维建模中,证明了该法具有方便、快捷和合理等优点。     

The use of renormalization for calculating effective permeability

There is a need in the numerical simulation of reservoir performance to use average permeability values for the grid blocks. The permeability distributions to be averaged over are based on samples taken from cores and from logs using correlations between permeabilities and porosities and from other sources. It is necessary to use a suitable effective value determined from this sample. The effective value is a single value for an equivalent homogeneous block. Conventionally, this effective value has been determined from a simple estimate such as the geometric mean or a detailed numerical solution of the single phase flow equation.If the permeability fluctuations are small then perturbation theory or effective medium theory (EMT) give reliable estimates of the effective permeability. However, for systems with a more severe permeability variation or for those with a finite fraction of nonreservoir rock all the simple estimates are invalid as well as EMT and perturbation theory.This paper describes a real-space renormalization technique which leads to better estimates than the simpler methods and is able to resolve details on a much finer scale than conventional numerical solution. Conventional simulation here refers to finite difference (or element) techniques for solving the single phase pressure equation. This requires the pressure and permeability at every grid point to be stored. Hence, these methods are limited in their resolution by the amount of data that can be stored in core. Although virtual memory techniques may be used they increase computer time. The renormalization method involves averaging over small regions of the reservoir first to form a new averaged permeability distribution with a lower variance than the original. This pre-averaging may be repeated until a stable estimate is found. Examples are given to show that this is in excellent agreement with computationally more expensive numerical solution but significantly different from simple estimates such as the geometric mean.