《高大平房仓20仓的曲面拟合粮温分析研究》由会员分享，可在线阅读，更多相关《高大平房仓20仓的曲面拟合粮温分析研究（9页珍藏版）》请在金锄头文库上搜索。
1、高大平房仓20仓的曲面拟合粮温措施小草摘要：应用曲面拟合粮温的措施分析研究高大平房仓20仓的粮温状况，分析研究了粮温与平面几何空间,与空间、时间,与空间、时间和气温的曲面拟合，得出了曲面拟合粮温函数,作了等温线图和细分各层粮温图,并分析研究。核心词:曲面拟合粮温 粮温与空间、时间、气温的曲面拟合函数 等温线图 细分各层粮温图 分析研究前 言应用曲面拟合粮温的措施分析研究粮温状况是一种新措施。我应用该措施对高大平房仓20仓这个实例进行分析研究，以便打开分析研究粮温的新措施做点有益的摸索。重要对高大平房仓20仓的宽度方向上的垂直面、长度上的垂直面和水平面的第3层的粮温T(,Y)进行了曲面拟合，并作
2、了其等温线图,此外对测温线8号的粮温（中央点）进行了曲面拟合粮温T(X,t时间)和T(X,t气温，t时间)，并作了其等温线图和细分各层粮温变化图,而后进行了分析讨论。解决数据重要措施：曲面拟合法；重要工具：电子表格中的函数LES()和曲面俯视图功能。一、 计算成果1.1、20仓的宽度方向上的垂直面的粮温T（,Y）曲面拟合函数和其等温线图：其基函数为双3次多项式：i (i=,， j0,2,3)曲面拟合粮温: i jT(X，Y)=j XiYj i=0 =0表1：高大平房仓仓第6列垂直剖面拟合粮温系数Aij13030.42570.-.-016-0.-0.20-0.0.030.0.0.复有关系数R2=
3、0.9985，剩余原则差e=0127,自由度f1=5,f2=,F值1417724。-1135791113151719-0.40.20.81.422.63.23.84.455.60510152025303540北墙南墙图1：高大平房仓20仓第6列垂直剖面粮温俯视图（7月31日横峰县吴家）35-4030-3525-3020-2515-2010-155-100-51.2、20仓的长度方向上的垂直面的粮温（X,Y)曲面拟合函数和其等温线图：其基函数为双3次多项式：Xij （0，1,3 j=0,1,2,3)曲面拟合粮温: 3T(X，）=A XiY =0 j=0表2：高大平房仓2仓第3行垂直剖面拟合粮温系
4、数Aij0123i32.80.0.0-.-.0.-02-0.00.0.30.0.-0.复有关系数2=.9756，剩余原则差Se=1.496，自由度f1=1,f2=28，F值=74.12。1.3、0仓的水平面的第3层粮温T（X，)曲面拟合函数和其等温线图:其基函数为对称园锥面二次曲面函数： Uk（X,)Q（,Y，j ,Yj)=（(X-Xj)2+(Y-)2)05 式中：Xj ,Yj为基函数的中心点坐标。曲面拟合粮温: 了 mT（X，Y)=Aj Q（X，Y, ,j)+A0= Aj（(X-)+(-Y)2)0.5 +A0 j=1 j=1表:高大平房仓20仓第3层水平剖面拟合粮温的中心点方位和系数AjA0
5、表中心点坐标（0，0）(0,9.08)(,18.6)(0，27.2）（0,3.）(0,45.4)相应系数2749-1.40-04562-.0686-0.84030.514中心点坐标(6,)（4.6,.54)（4，13.2)(4.6,22.)(4.6，378)（4.6,4086)相应系数2.391.9870.723312111.205.475中心点坐标(9.2,)(9.2，98)(9.,1.16)（9.，27.24）相应系数01625.1040.4539-0.42常数项7.15080.719e2.717自由度f16自由度f2F76214、0仓中央测温线号的粮温(X，t)曲面拟合函数和其等温线图:
6、前面项都是粮温在平面几何空间的曲面拟合粮温,目前讨论一下粮温与深度（X)和时间（t）的曲面拟合问题。将.1中Y换成t就可以。X值取粮面处为0,时间从6月20日至月8日止，以6月0日为0。其基函数为双3次多项式:Xij （=0,1, j=0,1，2,3)曲面拟合粮温: 3 3T(X，Y）=Aij itj i=0 j=0表4:高大平房仓0仓第2号测温线垂直拟合粮温系数Aij023i32.6330.0.-0.1-14.-0.-0.0.0.0.-0.-0.复有关系数R.972,剩余原则差e=0.41，自由度1=5，2=8,F值=4021.485。1、20仓中央测温线28号的粮温T(X，t气温,）曲面拟
7、合函数和其等温线图：目前讨论一下粮温与深度(X）、气温和时间()的曲面拟合粮温T（X,t气温，t）问题。值取粮面处为0，时间t从6月7日至月8日止，以6月0日为，气温为气象平均气温，为粮温前一天的气温。其基函数为：Xtjt气温 (i0,1,2,3 j0,1 k=0,)曲面拟合粮温: 3 1 1T(X,t气温，t)Aijk ittk气温 = j=0 k0表5:ijk系数表(j，)（，0)(，1)（1,0)（,)i02.6484.0.-.-7.0.00.1.-0-0.3-.0.0.0复有关系数2=91,剩余原则差Se=0.958,自由度15,f2=6，F值=60.6。二、 分析讨论1、到7月日时，
8、从上图中可以得出粮温的重要状况:粮面及表层粮温很高，下降梯度很大:0以上粮温从粮面进一步到0.至.米，5以上粮温从粮面进一步到0.8至.0米，东边进一步到.8米,北边进一步到.2米。垂向面上:2以上粮温的等温线基本是水平线，0等温线是弧线，两头向下(底层向上)，中间平缓,15等温线是一种长椭圆。中央区的等温线近为水平线,也就是说中央区各层粮温是粮堆高度的一维函数。 水平面上：中间有一种近长椭圆状的15等温线,四周粮温较高,东边、北边粮温更高,四角粮温较高。冷心（粮温15）的大体位置：距粮面28米,距地面0.6-1米，距东墙4.8-7.米，距西墙1.-8米，距北墙0.3-43米，距南墙0-21米
9、，其体积大体在200立方米，占粮堆的35%。此外,东墙、北墙边粮温较高,墙与墙、墙与地面相交处粮温较高，距地面2米处墙边的粮温较低，冷心区在这里向外延伸达到最大边沿。2、粮温的影响因素：从上面各图可以得出粮温的影响有气温、仓温、墙壁(墙温)、地温、阳光、仓库构造和外部环境7大因素,以及它们共同作用因素。阳光因素：如图3的东墙和北墙边的粮温明显高于其她墙边，这是由于东墙和北墙晒太阳;仓库构造因素：上面各图中，墙与墙相交处和墙与地坪相交处的粮温相对较高,由于是平房仓,故有墙与墙相交处，园仓没有；外部因素:仓与21仓之间做了一种仓棚,南墙基本晒不着太阳,故南墙边的粮温较低；共同作用因素:如图中2等温
10、两头向下弯曲就是仓温和墙温共同作用的成果。、一种有趣的现象:在上面图2、中，X=-0.4米处(粮温)有一层薄薄的5-40区，当时仓温只有285,故当时粮面粮温不也许达这样高温，但高大平房在月份中曾经达至37（早上8时,粮面粮温)和42（晚上10时,仓温)，故曲面拟合粮温有历史恢复功能？后来进一步摸索;、测温线28号粮温拟合状况，有三点需要指出： 一是粮温30达到0.8米深处后,基本上停止进一步了,二是20线从6月20日在深.米处到8月8日在1.8米处,平均每天进一步.02米,进一步很慢,粮面到深1米以内的粮温受气温直接影响并随着深度而削弱随着其波动而波动随着其达到最大气温而先后达到粮温的最大值
11、。5、从粮温经验公式中懂得，在中央区中第层、第2层和第3层粮温与导温系数有如下关系： (T|n+1-Tn1)*Z( Tn-Tn+1)a=3600*4*K(1000*C*P）=0.135 式中:Tn为第+1层发展1天后的粮温，Z为第n+小层到n层的距离（米）。由拟合粮温懂得，(|n+-n+)/(T-n+k）(在计算时深度按0.2米分层，故有许多层)在深度.米到.6米时，在图上基本上是水平直线,也就是基本上是常数,其平均值为.3，与粮温经验公式一致,其上粮温受气温直接影响，其下粮温受地温影响,故无此常数。两个措施相印证，也为粮温的经验公式提供了根据和运用范畴。6、不完美之处: 由于是以某行或某列或某层实测来进行曲面拟合粮温函数的，故在拟定冷心位置时,有的数据相印证要差一点，图2的左下角数据明显偏小,再有用基点中心法所求的拟合粮温的剩余原则差偏大等。三、 结论通过以上分析研究，曲面拟合粮温是可以精确反映粮温（复有关系数2=0.97,剩余原则差S=4)状况的，是一种有用的新措施。它不仅把粮温与空间联系起来，还可以把粮温与空间、时间和气温联系起来，扩大了分析研究领域和运用领域,加深对粮温变化规律的结识。7月18日晚
