matlab小球下降时间最短,基于Matlab的小球在粘滞液体中下落运动的研究

粘滞系数是表征液体粘滞性强弱的一个重要物理量,其准确测量被广泛的应用于机械、化学、材料、医学、工程、卫生保健和其他工业领域[1-2]。用落球法测液体的粘滞系数,是大学物理实验的一个基本而又经典的实验,通过实验可以测量液体的粘滞系数[3]。落球在粘性液体运动过程中,可以看成由加速下落和匀速下落两个过程组成,在测量时间时,因只取小球匀速运动的那段时间,所以要在小球下落一段距离后再计时。有文献研究表明,直径为3mm的钢球,从无初速度进入液面到匀速运动,用时很短,经过距离也很短,不需要预留一段距离[4]。本文利用不同直径的小球无初速度地在甘油中下落,测量出了小球通过中间AB段所需的时间,并通过MATLAB软件对实验数据进行了拟合和计算,数据表明,小球直径很小时,可以忽略量筒的直径对小球的影响,小球的直径较大时,要考虑量筒的直径对小球的影响,同时,还要考虑小球下落产生的湍流对小球的影响,要用不同的粘滞系数公式进行计算。1实验原理1.1小球在液体中下落的加速度已知密度为,直径为d的小球在高度为h密度为0的液体中所受的力包括:重力G=mg,浮力F=gV,以及黏粘力f=6rv;由牛顿第二定律知ma=mg-F-f,整理得ma=(-0)gV-6rv,代入小球体积公式V=43r3=16d3得:ma=(-0)g6d3-3dv;在小球匀速v0(收尾速度)下落时刻解得粘滞系数为=(-0)gd218v0,一般情况下,当小球在容器直径为D、液体高度为H的液体内下落时,液体各方向上都是无限广阔的这一假设条件是不能成立的。因此,考虑到容器壁的影响f=6rv(1+2.4dD)(1+3.3rH)(1+613Re-191080Re2+…)在直径为D的容器中,如果在钢球下落速度很小的情况下,不会有湍流,可忽略了修正项1+613Re-191080Re2+……于是修正后公式得=(-0)gd218v0(1+2.4dD)将其带入加速度公式得:a=(-0)g6md3-3(-0)gd318mv0(1+2.4dD)v令M=(-0)g,K=1(1+2.4dD)v0故得到简化后的加速度公式a=M(1-Kv)。1.2小球在液体中的运动速度对上述加速度公式a=M(1-Kv)进行积分,可得到小球运动的速度公式v=1K(1-e-KMt)。1.3小球在液体中下落的深度对上述速度公式v=1K(1-e-KMt)进行积分得小球下落深度公式h=1Kt+1K2M(e-KMt-1)。以上分析得出了小球在液体中的运动规律以及各种运动参数与时间的函数关系:速度公式为v=1K(1-e-KMt);加速度公式为a=M(1-Kv)=Me-KMt;深度公式为h=1Kt+1K2M(e-KMt-1)。1.4加速段的位移可由上述公式求得加速度公式a=Me-KMt,可以取加速末端时刻为=5KM,代入其中得出a=0.00674M=0.0543ms-2,M的大小为8.055m/s2,所以这时的小球的运动可以近似看作匀速运动,故可以求出小球加速段的位移为h=1K2M(e-5+4),可以近似取为h=4K2M;从而求得h与小球直径的关系式为h=(-0)g812d4。同时,这时的速度称为收尾速度v=1K,把AB间运动看成匀速运动时,则期间运动的时间为tAB=HABv。2实验过程及数据处理2.1实验过程实验装置如图1所示。图1实验装置图将直径分别为0.003241m,0.003992m,0.004780m,0.005683m的四种钢球,钢球的密度为=7.6103kg/m3,无初速放入直径D为0.062287m的量筒中,量筒中放入密度为0=1.35103kg/m3的