屋檐水槽构建模型.doc555VIP免费

组别：第六小组姓名：梁伟仙、叶佳盛、冼深发屋檐水槽模型摘要对于倾斜屋顶采用何种结构进行有组织排水，是民用建筑一个普通而有实际的问题。在民用建筑中一般是在房顶的边缘安装一个檐槽和一个竖立的排水管来排水，这种排水方式的可行性探讨实际上是解决水槽的容量在单位时间内能否足以排除雨水的问题。构建这种排水方式的微分方程模型，运用Maple对设计方案进行了可行性论证，提出了屋檐水槽模型的优化方案。关键词:速度平衡原理；微分方程；数值解；maple一、问题的重述为了雨天出入方便，房屋管理部门想在房顶边缘安装一个檐槽。现在有一个公司想承担这项业务，他们承诺：提供一种新型的可持久的檐槽，这种檐槽不管天气如何都能排掉房顶的雨水。房屋部门希望检验公司的承诺能否实现。简单来说,从屋脊到屋檐的房顶可以看成是一个12米长,6米宽的矩形平面.房顶与水平方向的倾斜角度要视具体的房屋而定,一般来说,这个角度通常在20°~50°之间。现在有一个公司想承接这项业务,他们允诺:提供一种新型的可持久的檐槽,它包括一个横截面为半圆形(半径为7.5厘米)的水槽和一个竖直的排水管(直径为10厘米),并且不管天气情况如何,这种檐槽都能排掉房顶的雨水。如图1所示。二、模型的假设2.1、合理假设(1)、降雨分布均匀并以垂直下落，并且直接落在房顶上；(2)、所有落在房顶上的雨水迅速流入水槽中；(3)、雨水不断从水槽中溅出；(4)、排水管道顺畅，没有任何障碍或阻塞；(5)、假设雨开始下时槽内有雨水，且水深高度为0.01m。2.2、符号说明：如下图表所示各因素与符号说明有关的因素因素类型符号单位降水强调输入变量rm/s时间变量ts房顶的倾斜度输入参数弧度房顶长度输入参数dm房顶的宽度输入参数bm水槽的半径输入参数am水槽中水的深度输入参数hm水槽中水的容量变量Vm流入水槽的流速变量Qm/s流出水槽的流速变量Qm/s排水管横截面积参数Am重力加速度常数gm/s表1三、模型的建立参照图1,根据速度平衡的原理,对于房顶排水系统有:水槽中水的流量的变化率=雨水的流入流量-排水流出的流量。即Vc(t)=Q-Q,这里Q、Q分别是单位时间流入水槽和从水槽流出的雨水流量。房顶雨水的流动情况,如图2所示。房顶的面积是bd,由于房顶是倾斜的,根据合理假设(1),实际受雨的水平面积应为bdcos,房顶上雨水的流量就是r(t)bdcos,雨水的流动是沿倾斜的房顶向下的,从而流入水槽的流量应该是它在铅垂方向的分量,即Q=r(t)cossin(1)而根据能量守恒定律得出，（2）所以，Vc(t)=Q-Q水槽中水的深度hV(t):=a^2*d*(arcos((a-h(t))/a)-(a-h(t))*sqrt(2*a*h(t)-h(t)^2)/a^2)定义v(t);V(t):=adarccos->V1：=diff(v(t),t);#求v’（t）>v2:=a^2*d(diff(h(t),t)/(2*a*h(t)-h(t)^2)^(1/2)+diff(h(t),t)*(2*a*h(t)-h(t)^2)^(1/2)/a^2-1/2*(a-h(t))/(2*a*h(t)-h(t)^2)^(1/2)/a^2*(2*a*diff(h(t),t)-2*h(t)*diff(h(t),t)));>v(t):=simplify(%);#对上式化简v(t):=-（8）所以根据（1）、（2）、（8）这样就得到模型：（9）即（10）：四、模型的求解与分析设定一组数据：a=0.075m,b=6m,d=12m,g=9.8m/s,A=0.0025m,,h(0)=0.01m（假设槽内有些积水）。下面，用Maple求（10）数值解。>restart;eq:=diff(h(t),t)=(r*b*d*sin(alpha)*cos(alpha)-A*sqrt(2*g*h(t)))/(2*d*sqrt(2*a*h(t)-(h(t))^2));#定义（10）式；>eq:=subs(a=0.075,b=6,d=12,g=9.8,A=0.0025*Pi,alpha=Pi/6,%);把这组数据代到上式eq:=>eq:=evalf(eq);化简针对房屋管理部门的要求，所以我们可能遇到两种情况：1、r(t)为常数；2、r(t)为周期函数：（1）当r(t)为常数时。考虑水槽的深度趋于一个低于0.075m的稳定值，即h’(t)=0>eq1:=0=.4166666667e-1*(31.17691454*r-.3477105892e-1*h(t)^(1/2))/(.150*h(t)-1.*h(t)^2)^(1/2);eq1:=0=>solve(eq1,h(t));8.03952576210r>eq2:=803952.5762*r^2>=0.075;eq2:=0.075<=8.03952576210r>solve(eq2,{r});{r<=-0.0003054326190},{0.0003054...