An Implicit Multistage Integration Method Including Projection for the Numerical Simulation of Constrained Multibody Systems

To be efficient, the simulation of multibody system dynamics requires fast and robust numerical algorithms for the time integration of the motion equations usually described by Differential Algebraic Equations (DAEs). Firstly, multistep schemes especially built up for second-order differential equations are developed. Some of them exhibit superior accuracy and stability properties than standard schemes for first-order equations. However, if unconditional stability is required, one must be satisfied with second-order accurate methods, like one-step schemes from the Newmark family.Multistage methods for which high accuracy is not contradictory with stringent stability requirements are then addressed. More precisely, a two-stage, third-order accurate Implicit Runge–Kutta (IRK) method which possesses the desirable properties of unconditional stability combined with high-frequency dissipation is proposed.Projection methods which correct the integrated estimates of positions, velocities and accelerations are suggested to keep the constraint equations satisfied during the numerical integration. The resulting time integration algorithm can be easily implemented in existing incremental/iterative codes. Numerical results indicate that this approach compares favourably with classical methods.     

注射成型中充填流动对短纤维取向的影响研究

在建立聚合物熔体在型腔中流动以及纤维取向的数学模型的基础上，提出一种数值方法，同时进行充填流动的模拟和纤维取向的预测，构造包含结点的体积微元，求解结点处的速度梯度。并且给出两个例子，说明流动对纤维取向的影响，预测的纤维取向分布与Ｌｏｃｋｅｔ的定性预测相符合。     

混合电路表面温度分布的一种新计算方法

研究混合电路的矩形三层有限模型，通过计算机解含有特殊函数的数值积分来计算混合电路上表面温度分布，并将这一方法编制成计算机软件，计算结果和实测结果基本一致。     

基于规则矩阵的数值化推理算法

何华灿等于１９９４年首次提出了数值化推理技术的概念，给出了一种基于规则矩阵的数值化推理算法，它同时完成规则的激活和可信度传播计算。本文提出另一种基于规则矩阵的数值化推理算法，它先判定规则是否被激活，再进行被激活规则的可信度传播计算。这两种方法各有优缺点，适用于不同场合，可根据应用需要选择     

SF_6断路器灭弧室结构对电场的影响

运用四边形八结点等参元有限元法对两种不同结构的SF_6断路器灭弧室内的静电场进行了数值计算。通过对100％开距时数值结果的分析,讨论了SF_6断路器灭弧室结构对电场分布的影响,分析了两种结构的优劣并提出了对其中一种结构的改进方案。     

全天棚送风相对两侧墙回风式洁净室特性研究

对全天棚送风相对两侧墙下部回风式洁净室的空态工况下的气流流型以及污染物分布特性进行了模拟，并将模拟结果与已有实验，实测结果进行了比较，分析，提出一些关于这种洁净室的新观点。     

流线型散流器顶送洁净室气流组织的研究

采用Ｋ－ε两方程紊流数学模型，对流线型散流器顶送洁净室的流场和污染物的分布进行了数值模拟。通过对计算结果的分析，认清了各种影响因素对洁净室气流组织的作用，根据流场特点对洁净室进行分类，确定了每一类的应用范围，并提出了优化设计参数。     

淬火工件/介质边界换热条件研究

在试验的基础上，开发了淬火工件／介质边界换热条件程序，该程序通过反算法计算边界换热系数．利用该程序与ＡＢＡＱＵＳ非线性有限元程序对淬火过程温度场进行数值模拟，模拟结果与实测值比较，误差在１０％以内．     

TNT 粉尘爆炸压力发展过程研究

在２０Ｌ球形粉尘爆炸装置中，对ＴＮＴ粉尘爆炸压力发展过程进行了研究，主要讨论了浓度、粒度对其最大爆炸压力（ｐｍ）及最大压力上升速率（（ｄｐ／ｄｔ）ｍ）的影响；同时，理论建模对不同粒度、浓度、湍流度等情况下压力发展进行了研究，发现与实验基本吻合．     

硫化天然橡胶溶胀曲线的一种拟合方法

探讨采用１／Ｙ＝Ａ＋Ｂ／Ｘ数学模型拟合硫化天然橡胶溶胀曲线。结果表明，该模型能好地描述硫化天然橡胶溶胀的动力学过程，对三种不同交联密度的样品拟合后，其相关系数均大于０．９７４８。