一类生物流体力学模型的稳定性

在研究颅内压集中参数模型的基础上，改进自适应的龙格-库塔法对一类生物流体力学模型进行数值模拟。通过合理控制计算量，得到了微分方程近似解的局部截断误差的估计。使用矩阵特征值分析微分方程的稳定性，在实际生理范围内变化模型参数，得到特征值随参数变化的关系，指出模型中存在实部为正的特征值。文章还讨论了控制矩阵特征值的变化对数值方法稳定步长的影响，并得到步长的取值范围。通过理论分析。得到控制方程的解随时间的发展和模型中生理参数的选取相关。分析了特征值变化的稳定性和参数的关系，并对模型进行改进，讨论其稳定性的情况。     

复弯矩表示的非圆截面杆平衡的Schrǒdinger方程及其半解析解

弹性细杆的平衡和稳定性问题的研究在工程和分子生物学中有重要的应用背景。利用文中提出的复柔度概念，建立了用复弯矩表示的非圆截面杆平衡的Schrǒdinger方程。借助复曲率概念，导出以杆的曲率、挠率和截面相对Frenet坐标系的扭角为未知变量的2阶常微分方程，此方程与传统使用的Kirchhoff方程等价。文献中仅适用于圆截面杆平衡问题的Schrǒdinger方程为本文导出方程的特例。对于准对称截面杆，用小参数法分别建立了零次和一次近似方程，其中零次近似方程存在解析解。对于截面的主轴坐标轴与中心线的Frenet坐标轴重合的无扭转杆特殊情形，Schrǒdinger方程转化为Duffing方程，应用数值方法作出了Duffing杆变形后的三维几何图形。     

强流长距离加速过程中的束包络数值计算

直线感应加速器束包络的计算编码（ABC）能给出束流从注入器经加速器直到最后聚焦的包络图，可以方便地研究各种磁场位形及隙压参数对束流输运的影响。作为算例，本文给出了束流为3kA，在长达40余米直线感应加速器的多种磁场位形下的束包络。束包络的计算对调试加速器磁场位形保证束流的正常输运及抑制束流不稳定性可提供帮助。     

Bifurcation of flow and mass transport in a curved blood vessel

A numerical analysis of flow and concentration fields of macromolecules in a, slightly curved blood vessel was carried out. Based on these results, the effect of the bifurcation of a flow on the mass transport in a curved blood vessel was discussed. The macromolecules turned out to be easier to deposit in the inner part of the curved blood vessel near the critical Dean number. Once the Dean number is higher than the critical number, the bifurcation of the flow appears. This bifurcation can prevent macromolecules from concentrating in the inner part of the curved blood vessel. This result is helpful for understanding the possible correlations between the blood dynamics and atherosclerosis. The project supported by National Natural Science Foundation of China (10002003), JSPS Postdoctoral Fellowship for Foreign Researcher and Foundation for University Teachers, the Ministry of Education     

促进其线性频散特征另一种形式的Bousinesq方程

Ｂｏｕｓｉｎｅｓｑ方程能够用于模拟表面重力波传播过程中的折射、绕射、反射以及浅化，非线性作用等现象．用不同垂直积分方法所得到的二维Ｂｏｕｓｓｉｎｅｓｑ方程形式具有不同的线性频散特征．采用两个不同的水深层的水平速度变量组合，推导出一个新形式的Ｂｏｕｓｉｎｅｓｑ方程．通过对其参数的设置可得到精确的线性频散解Ｐａｄｅ近似４阶精度．其适用范围已由原来的浅水，向深水拓进．相速误差小于２％，其拓展适用范围可达到０８个波长水深．应用所得到的新型Ｂｏｕｓｉｎｅｓｑ方程，采用有限差分法，对经典工况进行了数值模拟，其计算结果表明，计算值与物模实验值吻合较好．这说明本文新形式的Ｂｏｕｓｓｉｎｅｓｑ方程对变水深非线性效应所产生的能量频散有着较为精确的描述     

球形封闭容器内一个简单的煤粉燃烧爆炸模型

在分析了大量球形封闭容器内煤粉燃烧爆炸实验数据基础上，考虑了煤粉燃烧爆炸机理所涉及的湍流燃烧、相变、各种化学反应动力学过程等复杂因素，并且对球形封闭容器内由于煤粉混合不均匀造成的燃烧不充分给予了考虑，得到了球形封闭容器内煤粉燃烧爆炸特征的数值计算结果，计算的压力－时间曲线与实验结果符合较好。     

Melnikov方法在输流管混沌运动研究中的应用

对基础简谐运动激励下两端固定输流管道的混沌运动进行了研究，推导出了系统的运动方程，确定了系统存在的平衡点及其稳定性，计算出了未扰系统的同宿轨道，并利用Melnikov方法得到了系统发生混沌运动时参数需满足的临界条件，同时还利用相平面图和：Poincare映射等方法对管道的混沌运动进行了数值模拟，通过比较发现，由Melnikov方法确定的临界参数值要稍小于数值模拟中首次观察到混沌运动时对应的临界值。     

A quasi-flow constitutive model with strain-rate dependence

In this paper, the proposed is a quasi-flow constitutive model with strain-rate sensitivity for elastic plastic large deformation. The model is based on the Quasi-flow Corner theory, and is suitable for the sheet metal forming process simulation with a variable punch machine velocity. Uniaxial tensile tests and deep-drawing tests of a circular blank with square punch are carried out and numerically simulated. The consistency between the experimental and the numerically simulated results shows the validity of the present new constitutive model. The project supported by the Scientific Foundation of National Outstanding Youth of China (10125208), the National Natural Science Foundation of China (19832020), and the National Education Committee of China     

冲击消振器的概周期碰振运动分析

建立了冲击消振器对称周期运动的Poincar啨映射方程 ,讨论了对称周期运动的稳定性与局部分岔。通过数值仿真研究了冲击消振器在非共振、弱共振和强共振条件下的概周期碰振运动及其向混沌的转迁过程。     

Theoretical statistical solution and numerical simulation of heterogeneous brittle materials

The analytical stress-strain relation with heterogeneous parameters is derived for theheterogeneous brittle materials under a uniaxial extensional load, in which the distributions of theelastic modulus and the failure strength are assumed to be statistically independent. This theoreticalsolution gives an approximate estimate of the equivalent stress-strain relations for 3-D heterogeneousmaterials. In one-dimensional cases it may provide comparatively accurate results. The theoreticalsolution can help us to explain how the heterogeneity influences the mechanical behaviors, Further, anumerical approach is developed to model the non-linear behavior of three-dimensional heterogeneousbrittle materials. The lattice approach and statistical techniques are applied to simulate the initialheterogeneity of heterogeneous materials. The load increment in each loading stage is adaptivelydetermined so that the better approximation of the failure process can be realized. When the maximumtensile principal strain exceeds the failure strain, the elements are considered to be broken, which canbe carried out by replacing its Young‘s modulus with a very small value. A 3-D heterogeneous brittlematerial specimen is simulated during a full failure process. The numerical results are in good agreementwith the analytical solutions and experimental data.