时域边界元法分析撞水响应

基于势流理论，考虑流场的可压缩性，首先利用积分变换导得了势流问题的一个动力学倒易定理，在此基础上，进而求得问题对应的时空边界积分方程，然后通过对边界和时间轴同时离散，建立了一组有递推形式的时间边界元方程最后结合液面条件和物体运动方程耦全求解得到了刚体的撞水响应。     

非稳态弹流润滑状态下滚动体-滚道摩擦副的动力特性研究

针对滚动体-滚道摩擦副,建立了点接触非稳态弹流润滑数学模型,利用FFT技术和半解析算法数值求解了接触体在自由振动过程中油膜压力和膜厚的变化,同时结合有阻尼系统的自由振动模型,给出了预测点接触摩擦副动力特性的方法,在较宽的载荷和速度范围内分析了接触副的等效刚度系数和阻尼系数的变化.结果表明：根据接触副的实际工作载荷和速度所确定的无量纲自然频率来进行非稳态弹流的计算所得到的膜厚结果更接近实验值;在接触体的振动过程中油膜的压力和厚度在平衡位置附近上下波动,且由于润滑油膜的作用接触体的振动幅值逐渐减小;刚度系数随载荷参数的增加而增加,随速度参数的增加而减小,而阻尼系数的变化规律较复杂,在不同的载荷和速度范围内呈现出不同的变化趋势.     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。     

砂土类材料离散元细观参数的自动识别

准确的细观参数是采用离散元进行合理模拟的前提,基于最优化理论,本文提出了一种砂土类材料离散元细观参数的自动标定方法:采用PFC离散元软件模拟砂土的双轴压缩试验,以使PFC数值模拟的应力–应变曲线与某标准砂土室内双轴压缩实验结果的误差最小为优化目标,采用模式搜索法实现砂土类材料离散元细观参数的自动标定。本文提出的离散元细观参数自动标定方法具有识别效率快,精度高的优点。     

分块法研究圆柱绕流升阻力

 使用新的分块耦合方法，分别对单圆柱和串列双圆柱绕流进行了数值 模拟. 对于单圆柱绕流，低$Re$下计算所得到的定常涡尺寸与实验非常接近. 对于 串列双圆柱绕流，研究分析了改变双圆柱中心间距对上下游圆柱的升阻力系数和脉动频率所 产生的影响，计算结果与实验非常吻合，为进一步研究涡致振动提供了依据.     

滚滑工况弹流反常温度场的研究

研究了弹流反常温度场的形成机理及影响因素，指出入口温升是压缩功发热和逆流剪切热所致，而出口局部低温是负压缩功吸热的结果，出口温度的再次微幅上升则是压缩功消失后剪切热作用的结果．研究结果表明，入口温升随载荷的增加而增大，随卷吸速度的增加显著升高而几乎与滑滚比无关；在高速小滑滚比工况下，接触区的最高温度有可能出现在入口位置；入口温升增加了材料在工作中经受高温的次数，对其接触疲劳寿命有不利影响；在保证润滑性能的前提下，适当减少供油量可以减小逆流，从而降低入口温升。     

纳米二氧化硅水溶胶成膜特性研究

研究了纳米SiO2水溶胶的固液二相流成膜特性．结果表明，水溶胶中的Si02固体颗粒的粒径及质量分数对粘度的影响较小，对成膜能力和承载能力的影响较大；在薄膜润滑阶段，固体颗粒的粒径和质量分数对润滑状态的稳定性具有一定的影响，随着速度的增加，液体膜厚度逐渐增大，固体颗粒对液体润滑膜稳定性及性能的影响逐渐减弱；在弹流润滑阶段，固体颗粒的影响很小。     

土的工程力学性质的颗粒流模拟

基于颗粒流理论,引入不同的颗粒接触连接本构模型,分别建立了砂土和粘性土的颗粒流模型．通过颗粒流数值模型试验,对砂土和粘性土的室内平面应变试验及其剪切带形成和发展进行了数值模拟,分别对比了不同围压下颗粒流试样与室内试验的应力应变关系曲线,基本再现了砂土和粘性土试样应力．应变关系．通过砂土和粘性土PFC试样剪切带模拟表明,当围压较小时试样内部颗粒位移量小而且分布范围较广,当围压增大时,试样内部颗粒位移量也增大,而且发生较大位移颗粒的分布范围趋于集中,同时随着围压的增大试样内部形成明显的剪切带．无论砂土还是粘性土的PFC试样,随着围压的增加剪切带的形状趋于集中,而且剪切带宽度在减小．在围压很小时,试样内形成大的破坏区域,在围压较大时出现明显的线破坏区．这些规律基本与室内试验结果相似。     

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

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

高压下受限齿轮油成膜特性研究

本文揭示了在实际应用中极端工况下齿轮油的润滑特性.为了测量在高接触压力和不同温度条件下齿轮油的成膜能力,制作了高精度的膜厚测量仪,采用相对光强法实现纳米级膜厚测量.试验中共采用了5种性质不同的齿轮油,在纯滚动条件下测量接触区的油膜厚度,结果表明：随着接触压力的升高油膜厚度明显降低,但压力对润滑状态影响不大;温度的改变不但能影响油膜厚度,对润滑状态的影响也很明显;GL-5 85W/190和GL-5 85W/90齿轮油在接触压力达3 GPa温度提高到120℃时仍能形成很厚的油膜,但其余3种润滑剂在极端工况下成膜能力不足.本文最后根据试验结果还提出了含有丰富添加剂的齿轮油润滑模型.