离散体高速绕球时球体表面受力分析.docVIP

离散体高速绕球时球体表面受力分析 　　【摘 要】 基于离散群动力学模型和对流体绕球流动的认识，本文通过设计一套风洞试验，用应力应变片采集球体表面的压力曲线图，并辅以高速摄像系统观测拍照试验中的粉体的运动情况。对拍摄的图片和应力应变系统采集的信息数据与绕球流动数学模型进行对比分析。结果表明：离散体高速绕球体运动时，球体表面可以划分三个应力区域，分别为正应力区域、剪切力区域和雷诺应力区域。在球体表面0～45°、315°～360°是正应力分布区域，90°～135°、225°～270°是剪切力分布区域，45°～90°、270°～315°是正应力与剪切力混合区域。球体表面的正应力越大则剪切力越小，反之剪切力越大则正应力越小。 　　【关键词】 离散体 绕球运动 风洞试验 压力分布 　　在离散体高速绕球运动时离散体与球体发生着撞击、剪切、挤压、摩擦等各种运动[1-2]，这些运动都是非线性且相对独立的，这些运动时而接触时而分离，其运动过程在力学上面及其复杂，造成用传统的分析方法很难以解决。之前比较主流研究的人也最多的观点是把离散体颗粒当作一个理想的质点，对单个质点进行各种动力学研究并对其进行数学描述[3-6]。区别于以往的用单个离散体颗粒的力学研究分析球体表面的受力情况，本文通过具体的风洞试验，借鉴流体力学绕球流动的方法建立离散体群动力学模型，对整个离散体群绕运动进行分析从而研究得出球体表面的受力情况。 　　1 研究方法 　　建立风洞模型进行离散体高速绕球运动的风洞实验。将球体固定在风洞中，同时风洞中充满离散体颗粒，这些离散体颗粒不停围绕球体进行各种运动。为了得到球体表面的压力分布，在球体表面以球体的半径为半径的剖面圆上布置8个，再在这个剖面的两边以26为半径的剖面圆上分别布置8个应变片，这样总共24个应变片。通过应变片测得球体在风洞中表面的动应变数值。然后通过一个固定力产生的应变和力的关系的标值，将所测得的球体表面的动应变转化为相应的力，然后与理论分析的压力情况相比较，得出球体表面大概的压力分布情况。 　　1.1 风洞实验 　　选择的风洞尺寸为风洞尺寸为Φ290mm、长2000mm。参考???验型振动磨机参数：滚筒直径：126mm；滚筒长度：100mm；偏心距：20mm；钢锻直径分别为28.0mm、26.0mm、19.0mm、15.0mm，长度均为15mm。球体直径为：60mm。以球体的半径为半径的剖面圆上布置8个应变片，这个剖面的两边以26为半径的剖面圆上分别布置8个应变片，总共24个应变片。 　　根据试验方案拟采用200目石膏粉、方解石粉和石英粉作为风洞试验的离散体材料。在本风洞试验中，在球体的24个点上通过安装的应变片用三种材料共测得了72组数据。试验所选用的粉体材料为200目方解石粉、200目石膏粉、200目石英粉，每种粉体材料测得24组数据。试验测得24个预定点的动应变数据。测得数据有几百个，但是在某一应变大小范围内波动，求得这些数据的平均值，然后再加以处理分析。采集所得数据是应变值，但是我们想要的是力，所以必须把应变转换为力。取与风洞试验相同的球体材料，加载一个数值固定的力，测得应变，即可得出力与应变关系。 　　试验测得是球体表面的动应变数值，要把它们转化为对应力的数值，必须测得相同材料的标值。所谓标值即是施加一固定不变的力在相同材料上，测出材料的应变数值，得出力与应变的一一对应关系。选择质量为200克的一个砝码， 　　其应变为175.345.根据湖南湘潭市本地的重力加速度为 9.794m/s2由F=mg得200克标值对应的力为1.9588牛。处理所测得的72组数据得到三种粉体材料分别在所测点的应变平均值，再通过前面所做的力和应变的标值通过一一对应，将应变转化为相应的力（如表1）。 　　1.2 球体表面理论无因此压力分布曲线 　　工程上常用无因次压力来表示球体上任一点处的压力。这样可以使得无因次压力分布规律不受速度及压力p的绝对值的影响，使用时更方便。 　　无因次压力的定义 　　分析可得离散体绕旋转球体运动时，球体体表面上无因次压力分布规律，表达式为 　　本试验所选取的振动磨机参数为：滚筒直径Φ=126mm，长为100mm；试验用有机玻璃管直径Φ=290mm，长为2000mm；= 22.8m/s。可得球体体表面无因此压力分布曲线图如下图1： 　　2 实验结果及分析 　　将表1中的数据，将同一剖面上的8个点如上表所示的所测得的力用软件绘图，经处理后可得球体表面试验无因此压力分布曲线如下图2所示： 　　由图可以看出，在θ=±45°的区域中，理论曲线与多种速度下的离散体绕球体情况比较接近。试验表明，当极角θ比较小时，不同速度的离散体其压应力分布几乎合成为一条曲线。随着极角θ的增大，不同速度的试验曲