计算机磁头／磁盘超薄气膜润滑压强的算子分裂算法

以任意拉森数的超薄气体润滑方程为基础，采用算子分裂法和非结构三角网格的有限元法计算Ω型磁头的空气轴承气垫面(ABS)气膜压强分布、气浮力和纵倾力矩；在分析比较流量系数的各种算法的基础上，确定采用多项式拟合数据库计算流量系数．计算结果表明：采用算子分裂法可以有效克服在高轴承数时的数值不稳定性，消除数值振荡；在小轴承数时，气浮力随轴承数增加而增大，当轴承数增大到某一数值后，气浮力趋近某一稳定值，此时气膜压强分布与磁头造型基本一致；气膜的纵倾力矩在轴承数的某一临界值附近出现最大值．     

海岸碎波拍的计算

研究了一维碎波拍的计算方法．波浪控制方程由欧拉方程在短波周期上平均而得到．为了考虑低频波浪在海岸上的破碎，采用了ＷＡＦ（ＷｅｉｇｈｔｅｄＡｖｅｒａｇｅＦｌｕｘ）方法求解长波方程．对由辐射应力和海底摩擦力引起的方程非齐次项采用了时间算子分裂法进行了处理，并对由此产生的误差进行了修正．计算了单色波和双色波垂直入射海岸的情况，并与实验结果进行了对比．     

基于新型解耦算法的激波诱导燃烧过程数值模拟

对一种模拟化学非平衡流动的时间和空间二阶精度新型解耦算法进行两方面改进,流动算子采用基于Runge-Kutta方法的时间格式以后, 可以推广到更多的空间差分格式,化学反应源项求解算子可以采用梯形公式、拟稳态逼近法和变系数常微分方程求解器. 对H化学非平衡流动; 解耦算法; 计算方法对一种模拟化学非平衡流动的时间和空间二阶精度新型解耦算法进行两方面改进,流动算子采用基于Runge-Kutta方法的时间格式以后,可以推广到更多的空间差分格式,化学反应源项求解算子可以采用梯形公式、拟稳态逼近法和变系数常微分方程求解器.对H_2/Air预混气体中激波诱导振荡燃烧的Lehr试验进行模拟,考察了化学动力学模型、网格尺寸和差分格式耗散大小对计算结果的影响,同时对不同的化学反应源项算子求解算法的计算效率进行了比较.     

七方程可压缩多相流模型的HLLC格式及应用

针对Saurel和Abgrall提出的两速度两压力的七方程可压缩多相流模型,改进了其数值解法并应用于模拟可压缩多介质流动问题.在Saurel等的算子分裂法基础上,根据Abgrall的多相流系统应满足速度和压力的均匀性不随时间改变的思想,推导了与HLLC格式一致的非守恒项离散格式以及体积分数发展方程的迎风格式.进一步,通过改变分裂步顺序,构造了稳健的结合算子分裂的三阶TVD龙格-库塔方法.最后通过几个一维和二维高密度比高压力比气液两相流算例,显示了该方法在计算精度和稳健性上的改进效果.      

AN HLLC SCHEME FOR THE SEVEN-EQUATION MULTIPHASE MODEL AND ITS APPLICATION TO COMPRESSIBLE MULTICOMPONENT

针对Saurel和Abgrall提出的两速度两压力的七方程可压缩多相流模型,改进了其数值解法并应用于模拟可压缩多介质流动问题．在Saurel等的算子分裂法基础上,根据Abgrall的多相流系统应满足速度和压力的均匀性不随时间改变的思想,推导了与HLLC格式一致的非守恒项离散格式以及体积分数发展方程的迎风格式．进一步,通过改变分裂步顺序,构造了稳健的结合算子分裂的三阶TVD龙格一库塔方法．最后通过几个一维和二维高密度比高压力比气液两相流算例,显示了该方法在计算精度和稳健性上的改进效果．     

广义Stokes方程的完全边界积分表示式及其在求解N-S方程中的应用

为了分解N-S方程组各变量相互偶合,本文采用Peaceman-Rachford算子分裂法,将时间相依的N-S方程组分解成不存在上述偶合特性的线性和非线性的子问题。线性子问题具有广义Stokes方程类型。本文采用多重互易法,即采用多阶拉普拉斯算子基本解逐步变换,将其解表示成完全边界积分形式,从而使问题的计算维数降低一维。广义Stokes方程的算例以及二维圆柱在剪切流中的Stokes绕流解,都表明多重互易算法具有高效特点,而且后者与文[3]解析解吻合得非常好。     

不同设计参数下螺旋槽径向液体润滑轴承油膜稳定性的比较计算

采用Ｇａｌｅｒｋｉｎ有限元法求解油膜压力的雷诺方程和扰动压力的摄动方程,并分析了螺旋槽径向液体润滑轴承的动态性能。采用逐步修正的迭代法实现正压力区和油膜破裂区交界处的雷诺边界条件。在不同参数下比较了螺旋槽轴承油膜不产生涡动的转轴临界质量Ｍｃ,并进行了参数的比较计算,发现合理地选择参数可以得到良好的油膜稳定性,甚至可达到绝对稳定。     

二维磁流体力学ALE数值模拟

介绍采用ALE方法进行的电磁内爆二维磁流体（MHD）力学模拟。二维MHD力学模拟的ALE方法分成拉格朗日、网格重分、对流输运三步。拉格朗日步采用文献[4]的时间分裂法,分成辐射冷却、热扩散、磁扩散、Lagrangian流动四步求解,并对Lagrangian流动步的方法进行了改进,消除了其拉氏计算的质量损失等。编制了二维计算程序并通过验证,获得的二维电磁内爆不稳定性发展与文献[2,5,11]吻合。     

单节流孔静压球面气体轴承动态特性的有限元分析

采用分布参数和小参数摄动方法简化动态雷诺润滑方程，推导了用于计算小孔节流方式的单节流孔静压球面气体轴承气膜动态特性的一系列计算公式，并通过动态压力分布方程的有限元求解得到了气膜的动态刚度、动态阻尼曲线以及轴承的稳定工作范围图谱，进而给出了静压气体轴承的稳定性判据．结果表明，小孔节流方式静压气体轴承的稳定性主要取决于节流气腔的容积、供气压力和承载质量．     

航空发动机轴承腔润滑的气液两相均匀流研究

基于轴承腔中润滑油的两相均匀流动模型,采用湍流模式和有限差分数值方法计算轴承腔内三维定常N-S方程,对腔内润滑油的气液两相均匀流动特性进行研究,以获得气液两相均匀流条件下润滑油流场、压力场和速度场在轴承腔内的分布情况,分析转子转速、含气率和润滑油进口速度对润滑油出口压力以及润滑油与壁面之间剪切力的影响,同时对单相流和两相均匀流润滑性能差异进行比较.结果表明,转子转速、含气率和润滑油进口速度对润滑油出口压力和腔内壁面与润滑油间的平均剪切力具有不同影响,而采用2种流动模型计算出的轴承腔润滑油出口压力的差异较大,同时支持了开展航空发动机轴承腔润滑两相流动分析的必要性.     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

两自由度振动系统的斜碰撞分析

研究斜碰撞振动系统动力学的一个关键问题是对系统在碰撞前后的状态进行合理描述和正确计算．针对两弹性体斜碰撞问题，基于瞬间碰撞假设，提出了采用步进冲量来分析和求解斜碰撞前后的状态关系；并以弹簧摆和振子组成的两自由度斜碰撞振动系统为例，具体介绍了该算法如何实现．用解析方法讨论了该系统在斜碰撞过程中可能出现的各种力学现象，将冲量步进算法得到的数值解与解析结果进行对比，取得了完全一致的结果．该数值方法能适应多种斜碰撞问题的计算．     

Application of the theory of branching of solutions of nonlinear equations in the estimation of pole perturbance of linear closed-loop systems

When the state and input matrices of a multivariable linear time-invariant system are perturbed, the problem of the estimation of pole perturbance of closed-loop system is considered by making use of the theory of branching of solutions for nonlinear equations.     

The distribution of zeroes of solutions of neutral equations

1.IntroductionThestudyofoscillator}'behaviorofneutraldifferentialequationsisarelativealnewfieldandisveryinterestinginapplications.Inrecentyears.theoscillatorytheoryfortheseequationshasbeenextensivelydeveloped.WerefertotherecentbookbyGyoriandLadaSI'I,Baino…     

CHARACTERISTICS OF SUBDIFFERENTIALS OF FUNCTIONS

ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ（郭兴明）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ￥ＧｕｏＸｉｎｇｍｉｎｇ（ＲｅｃｅｉｖｅｄＪｕｎｅ１６，１９９５...     

SWT-120风洞稳定段的性能测量

在普通超音速风洞中,由于受到噪声干扰,很难进行有效的边界层转捩特性试验研究,针对超声速流动特点发展较低噪声风洞十分必要,而稳定段设计的好坏直接影响到下游试验段噪声水平。本文介绍一座低噪声风洞稳定段的结构和性能测量结果。先对稳定段的结构设计做了简单介绍,然后对实验结果进行分析,实验结果表明在大角度扩散段内装置孔锥,稳定段安装消音夹层和阻尼网组等部件后,气流的速度脉动和压力脉动明显降低,其中压力脉动降低一个量级,速度脉动为1%。进一步优化设计和改进工艺,速度脉动还可进一步降低。测量结果表明SWT-120稳定段的设计是成功的,对我们以后发展更高性能的静风洞有借鉴和参考价值。     

综合管廊早龄期混凝土温度应力仿真模拟

以汉中市综合管廊为工程背景,通过有限元软件,首先模拟分析综合管廊早龄期温度场,随后将温度场结果附加到管廊结构中求解温度应力,从而分析综合管廊早龄期混凝土温度应力与浇筑长度的关系。模拟分析发现:综合管廊在混凝土浇筑早期,温度场沿长度方向变化不显著,且随着浇筑长度的增加,温度场不变;在墙体温度上升阶段,墙体内部沿长度方向的应力主要为压应力,当墙体混凝土内部温度达到最高,混凝土开始降温,结构内变形由膨胀逐渐变为收缩,内力由压应力逐渐变为拉应力,随着时间的变化,墙体内外温差趋于稳定,墙体温度应力也逐渐平稳,此时温度应力主要受外界环境变化影响;当综合管廊浇筑长度在60m范围内,浇筑长度越大,温度应力越大,但当浇筑长度超过60m,温度应力不再增加。     

Evolution of viscoelastic properties of polyurethane in the course of curing

Simultaneous measurements of the components of dynamic modulus at the set of frequencies were made by Fourier-Transform Mechanical Spectroscopy (FTMS) to investigate the curing process of a segmented polyurethane.The whole process of curing is treated as consisting of three stages. In the first stage, the material can flow and this stage is complete at the gel-point. In the second stage, the curing proceeds in the network and ends with phase separation of aggregating segments (blocks) of polymeric chains. In the third stage, the curing continues in the two-phase system. The object of investigation was the transition from the second to the third stage. It was shown that the position of the maximum of loss tangent depends on frequency and does not coincide with different characteristic points marking the time dependencies of G and G. This is due to the formation of micro-particles of the newly formed phase. The FTMS method is used for separation of the modulus changes into parts related to formation of chemical and physical networks.Special means of treating calorimetric measurements identified a weak maximum heat output at the moment corresponding to the transition point.     

Complete Asymptotic Expansions of Solutions of the System of Elastostatics in the Presence of an Inclusion of Small Diameter and Detection of an Inclusion

We consider the system of elastostatics for an elastic medium consisting of an imperfection of small diameter, embedded in a homogeneous reference medium. The Lamé constants of the imperfection are different from those of the background medium. We establish a complete asymptotic formula for the displacement vector in terms of the reference Lamé constants, the location of the imperfection and its geometry. Our derivation is rigorous, and based on layer potential techniques. The asymptotic expansions in this paper are valid for an elastic imperfection with Lipschitz boundaries. In the course of derivation of the asymptotic formula, we introduce the concept of (generalized) elastic moment tensors (Pólya–Szegö tensor) and prove that the first order elastic moment tensor is symmetric and positive (negative)-definite. We also obtain estimation of its eigenvalue. We then apply these asymptotic formulas for the purpose of identifying with high precision the order of magnitude of the diameter of the elastic inclusion, its location, and its elastic moment tensors.