SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

二自由度耦合非线性随机系统的最大Lyapunov指数和稳定性

利用摄动方法讨论了一类耦合二自由度非线性系统，在小强度白噪声参数激励下系统运动模态的稳定性，获得了系统扩散过程的稳态概率密度的渐近表达式，由此获得了系统运动模态几乎必然稳定的充分必要条件。     

一种岩石损伤本构模型在地下强爆炸中的应用

为了描述地下爆炸波的传播,建立了一种考虑塑性硬化、剪胀和损伤软化效应的岩石本构模型,并应用于地下强爆炸自由场的数值计算,获得的速度和位移波形以及峰值速度衰减曲线等与国外地下试验测量数据和同类计算相比十分接近,从而验证了本构模型的有效性。     

偏压电场对压电板中Lamb波相速度的影响

本文研究了偏压电场作用下,Lamb波在压电板中的传播行为，首先给出了偏压电场作用时压电板中的应力场及电位移场，然后通过求解含初应力及初电位移的小幅波动问题的耦合方程，分别给出了Lamb波的对称模态和反对称模态的相速度方程，以典型的PZT－5H压电陶瓷板为例进行了数值计算,并讨论了偏压电场对Lamb波相速度及频散曲线的影响，结果表明，偏压电场可以显著地改变Lamb波的传播速度，借此可使声波器件获得延时效果。     

一种确定均匀动脉壁面切应力的非线性方法

从Ling和Atabek提出的``局部流'理论出发，提出一种利用测量血液黏度、管轴上 的血流速度、压力和管径波形计算均匀动脉管壁切应力的非线性方法. 将这种方法与柳兆荣 等提出的利用测量血液黏度、管轴上的血流速度和平均管径计算切应力的线性方法比较，结 果表明，当管壁脉动幅度较小时，两种方法计算的压力梯度、流速剖面和管壁切应力差别较 小；而当管壁脉动幅度增大时，两种方法计算的压力梯度、流速剖面和管壁切应力差别增大. 对于小幅脉动均匀动脉，用线性方法计算管壁切应力有较高的精度；而对于大变形 均匀动脉，则需要考虑非线性因素对管壁切应力的影响. 由于作为输入量的血液黏度、轴心 血流速度、压力波形和管径波形可在活体上通过无损伤或微损伤的检测方法得到， 所提出的计算切应力的方法为在体或离体研究切应力与动脉重建的关系提供了方法学基础.     

ANALYTICAL EVALUATION OF PERMANENT DEFLECTION OF A THIN CIRCULAR PLATE STRUCK NORMALLY AT ITS CENTER BY A PROJECTILE

The permanent deflection of a thin circular plate struck normally at its center by a projectile is studied by an approximate theoretical analysis,FEM simulation and experiment.The plate made of rate sensitive and strain-hardening material undergoes serious local deformation but is not perforated during the impact.The theoretical analysis is based on an energy approach, in which the Cowper-Symonds equation is used for the consideration of strain rate sensitive effects and the parameters involved are determined with the aid of experimental data.The maximum permanent deflections predicted by the theoretical model are compared with those of FEM sim- ulation and published papers obtained both by theory and experiment,and good agreement is achieved for a wide range of thickness of the plates and initial impact velocities.     

磁探针速度传感器研究

应用Hopkinson杆装置,对磁探针速度计进行了实验测定。发现对待测自由面速度阶跃,磁探针有两种类型稳定的输出信号,分别对应于不同的速度范围。对某个特定磁探针来说,其初始信号峰值与速度成正比的条件是速度大于某个阈值。实验还测定了间距大小对磁探针输出信号的影响,给出了磁探针对阶跃实验的响应波形。基于实验条件,本文探讨了磁探针原理、测量电路特性以及靶板涡流的穿透深度等理论问题,得到的结果可以定性解释磁探针的两种特征信号以及它们间的转化机理。本文工作表明,磁探针可以用来定量测量自由面速度的初始阶跃幅度,准确测定全时域的速度阶跃时刻,但不足以确定整个速度历史。为使这种传感器更加成熟,其结构和线路有待改进。     

半浮区液桥热毛细对流速度场特征

液桥流场的速度分布是液桥表面张力梯度驱动对流的一个重要研究内容。通过地面模拟微重力环境的液桥实验以及数值计算,研究了外加温差、液桥腰径等对流场速度分布的影响,实验结果与数值计算结果相当符合。结果还表明,当液桥腰径减小时,流场的结构发生了变化,涡心由一个转化成两个。     

麦克斯韦速率分布的特色推导方法

大学物理教学的一个基本点是培养学生的空间想象力.麦克斯韦速率分布律是气体分子运动论的中心内容,是大学物理气体运动理论中讲授的一个难点,其公式抽象、繁难,学生不易理解.本文根据速度空间概念,给出速度球的表面积相当于气体分子微观状态数的观点,利用拉郎格日函数,推导理想气体平衡态下气体分子的速率分布函数.这种推导方法相对比较简单、易学,便于学生理解和培养学生的空间想象力.