数学建模(5)---煤矸石堆积问题
煤矸石堆积问题<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" />
摘 要
我们对煤矸石的堆积储存问题进行了研究。根据煤矸石的堆积要求,建设一段与地面角度约为β=25°的直线型上升轨道(角度过大,运矸车无法装满),用在轨道上行使的运矸车将矸石运到轨道的顶端后向两侧倾倒,待矸石堆高后,再借助矸石堆延长轨道,这样逐渐堆起如下图1所示的一座矸石山来。
首先,对堆积煤矸石最终所形成的规则的几何体分析,所形成的几何题是数学中较常见的几何体,因此可以根据几何关系找出这几个体中的相关因数。在整个运输煤矸石的工程中,随着对放煤矸石山的坡面长度的增加,整个运输煤矸石的费用也会相应的增加。整个费用包括土地征用费用和运输煤矸石的费用。土地征地费用与堆积的煤矸石所占地面积有关,而占地面积又与煤矸山坡面长度的增加而增加;同时,运输煤矸石所用的电费也是随着坡面长度的增加而增加的。将费用问题转化成为研究煤矸石山坡面长度,而整个坡面长度是时间的函数,于是总费用又与煤矸石的开采时间建立了联系,最终我们通过给定使用时间,得出坡面长度,确定处理矸石的总经费,在处理征地费时,由于地价的年涨幅大于了银行的贷款利息,所以我们用开始由银行贷款购足使用年限内所需的所有土地者正方式,结果表明这种操作后,设计中的年处理经费除部分用来缴纳电费外,其余全部偿还银行的贷款,按处理经费为100万元/年,适用年限为20年出矸率为10%计算,该设计经费是过用的且有结余。而且可以在第14年时还清银行贷款,在20年后,可以盈利4203637.83元。
关键词
使用年限 坡面长度 机械能 几何体体积 处理总费用 地价涨幅 贷款利率
一、 问题的提出
煤矿采矿时,会产出废料美感使在平原地区,煤矿只得征用土地堆放煤矸石。堆放煤矸石时,需架设一端与地面角度为一定值β的轨道(角度过大煤车无法装满),矸石的自然安息角а也为定值,运送煤矸石的矸石车机械效率也随轨道的增加而下降,土地的征用费存在年涨幅,在众多因素的影响下,我们需解决如下问题:
1. 判断设计处理经费是否够用;
2. 根据有关数据制定合理征地计划;
3. 对不同出矸率预测处理矸石的最低费用;
二、 问题分析
根据题意,煤矸石堆的几何形状是一定的(如图1),我们要解决的就是建立煤矸石堆积总费用与几何形状的体积的关系,由于地价存在年涨幅,征用土地需向银行贷款,又由于地价涨幅和银行贷款利率有所不同,还需考虑贷款方案,根据煤矿使用数年后的矸石体积,计算出坡面长度x确定占地面得出处理总费用,列出在不同出矸率时的最小费用。
<?xml:namespace prefix = v ns = "urn:schemas-microsoft-com:vml" />
图1
三、 模型假设
1. 电费在当年付清、不拖欠;
2. 银行存、贷款利率和地价涨幅都固定不同;
3. 运矸所消耗的电能出机械损失外,全部转化成煤矸石的重力势能(运矸车匀速运动,且矸石车与坡面轨道见摩擦不计);
4. 征地为每年年初至征地一次;
5. 每年的处理经费按年拨对,且除交付当年电费外,其余全部用来征地或偿还贷款,严格做到专款专用,
四、 名词及符号约定
1. x ——煤矸石山坡面长度;
2. ρ ——容重;
3. а ——矸石自然堆放安息角;
4. k1 ——用矸车所需的电费;
5. η0 ——运矸车机械效率初始值;
6. a ——轨道每延长<?xml:namespace prefix = st1 ns = "urn:schemas-microsoft-com:office:smarttags" />10米,机械效率在原有的基础上下降的百分数;
7. k2 ——土地征用费的现值;
8. b ——地价年涨幅;
9. c1 ——银行的存款利率;
10. c2 ——银行的贷款利率;
11. M ——煤矿设计的原煤年产量;
12. N ——煤矿的设计寿命;
13. d ——煤矿的出矿率;
14. S(x) ——坡面长度为x时,矸石堆实际占地面积;
15. V(x) ——坡面长度为x时,矸石堆的实际体积;
16. η(x) ——坡面长度为x时,运矸车的机械效率;
17. w(x) ——坡面长度为x时,运矸车所要消耗的电能;
18. G(n) ——煤矿开采n年时,矸石的总体积;
19. H(n) ——煤矿开采n年时,所用的总的电费;
20. S(x(N)) ——在使用年限范围内,所用的土地总面积;
五、 模型的建立
1) 对几何体的分析
由图1知,A-SBOD为一个棱锥部分,A-DCOB为圆锥部分,SB与圆锥地面相切,图1中各个关系如下:
AO=sin(β)·x CO=sin(β)·cot(а) ·x SO=cos(β)·x
∠BOS=arc cos(tan(β)·cot(а))
2) 矸石山的底面积S(x)
S(x)=S四边形SDOB+S扇形BODC
S四边形SDOB=x2·cos(β)·sin(β)·cot(а)·sin[arc cos(tan(β)·cot(а))]
S扇形BODC= x2·[pi-arc cos(tan(β)·cot(а))]·sin2(β)·cot2(а)
S(x)= x2·{cos(β)·sin(β)·cot(а)·sin[arc cos(tan(β)·cot(а))]+[pi-arc cos(tan(β)·cot(а))]·sin2(β)·cot2(а)}=K(а, β) ·x2
(其中:K(а, β)= cos(β)·sin(β)·cot(а)·sin[arc cos(tan(β)·cot(а))]+[pi-arc cos(tan(β)·cot(а))]·sin2(β)·cot2(а)
3) 矸石山体积V(x)
V(x)=1/3·S(x) ·AO=1/3·K(а, β)·sin(β)·x3
4) 机械效率
由题意知,在轨道每延长10米,效率在原有基础上下降2%,所以η(x)=30%·(1-2%)x/10
5) 机械能在坡面长度为x时,对微元增量dx存在体积微增量dv(x)
dw(x)=dv(x) ·ρ·g·x·sin(β)/ η(x)= K(а, β) ·sin2(β) · ρ·g·x3·dx/(30%·(1-2%)x/10)
w(x)=∫0x (K(а, β) ·sin2(β) · ρ·g·x3 /(30%·(1-2%)x/10))dx
6) 开采n年时出矸石的总体积;
G(n)=n·M·d/ρ=V(x(n))
X(n)=(3·n·M·d/ρ·K(а, β) ·sin(β))1/3
7) n年所用总电费
H(n)=w(x(n)) ·0.5/(3.6·106)
第n年需交电费为H(n)-H(n-1)
8) 使用年限N内,所用土地总面积S(x(N))= K(а, β) ·x2(N)
附录一:
1.
function [rst]=t1(nYear)
if nargin~=1
error('Error input arguments!');
else
b=25/180*pi;a=55/180*pi;
p=2000;g=9.8;
sinb=sin(b);
cosb=cos(b);
cota=cot(a);
arccos=acos(tan(b)*cot(a));
kab1=cosb*sinb*cota*sin(arccos);
kab2=(pi-arccos)*sinb^2*cota^2;
kab=kab1+kab2;
xf=(4500000*nYear*0.1/(kab*sinb))^(1/3);
kw=kab*sinb^2*p*g/0.3;
sym x;
wIntegral=int('x^3/(0.98^(x/10))','x',0,xf);
W=kw*wIntegral;
monW=0.5/(3.6*10^6);
Area=kab*xf^2;
monA=1.1*80000/666.6667;
AllMoney=W*monW+Area*monA;
% disp('The total area is (Unit:MU);');disp(Area/666.6667);
% disp('Area money (Unit:RMB):');disp(Area*monA);
% disp('Eletricty money (Unit:RMB) ');disp(W*monW);
% disp('Need money (Unit:RMB):');disp(AllMoney);
rst=W*monW;
% rst=AllMoney;
end;
2.
function t2()
for iCount=1:19
rstW(iCount,1)=t1(iCount+1)-t1(iCount);
end;
abc=7595900*1.05-(1000000-73391);
for nCount=1:19
disp('The year is; ');disp(nCount+1);
% disp('All money have returned.');break;
% else
abc=abc*1.05-(1000000-rstW(nCount,1))
% end;
end;
转载于:https://www.cnblogs.com/AloneSword/archive/2004/09/12/2237721.html
数学建模(5)---煤矸石堆积问题相关推荐
- 清风数学建模笔记——Excel画图保姆级教学
备注:本篇内容是在Office环境下作图,WPS不支持直方图和箱线图好像.免费的Office软件可以在微信搜索中搜"Office"然后挑一个下载即可 本文借鉴了数学建模清风老师的课 ...
- 2023五一数学建模该选哪道题?选题建议及初步思路来啦
大家好呀,2023五一数学建模开赛了,在这里给大家带来初步的选题建议及思路,之后我会一直更新更加详细的完整思路以及后续也会有完整成品的讲解视频,手把手教大家如何去做哈,目前的初步思路讲解视频: 202 ...
- 2021 年高教社杯全国大学生数学建模竞赛A题分析
2021 年高教社杯全国大学生数学建模竞赛A题分析 题目 赛题分析 前言 问题一分析 问题二分析 问题三分析 题目 A 题 "FAST"主动反射面的形状调节 中国天眼--500 米 ...
- 2021 年第十一届 MathorCup 高校数学建模挑战赛A题分析
2021 年第十一届 MathorCup 高校数学建模挑战赛A题分析 题目:自动驾驶中的车辆调头问题 问题分析 前言 综述 问题一分析 问题二分析 问题三分析 问题四分析 问题五分析 问题六分析 题目 ...
- 全国计算机建模三等奖,2009年全国数学建模真题(论文获国家 三等奖).doc
2009年全国数学建模真题(论文获国家 三等奖) 制动器试验台的控制方法分析 摘要 制动器的作用在于使行驶减速或者停止,其设计是车辆设计中最重要的环节之一,直接影响着人身和车辆的安全.本文分析了制动器 ...
- 《数学建模:基于R》一一2.2 方差分析
本节书摘来自华章计算机<数学建模:基于R>一书中的第2章,第2.2节,作者:薛 毅 更多章节内容可以访问云栖社区"华章计算机"公众号查看. 2.2 方差分析 方差分析是 ...
- 数学建模优化模型简单例题_数学建模之优化模型:存储模型
点击上方「蓝字」关注我们 最近,为申报市级精品课程,我为我校"数学建模与科学计算"课程录制了讲课视频,下面是3.1节优化模型的第一个例子:存储模型.敬请大家批评指正! 优化模型是数 ...
- 资料分享:数学建模资料分享 -- 神经网络部分
背景 周日的时候,为数学建模俱乐部的同学们进行了一场有关人工神经网络方面的分享.虽然在这个方面有一些积累,但过于零散,所以拿了一堆文件拼凑成整体的内容. 幸亏自己有分享的习惯,学会一些知识就写下来,这 ...
- 资料分享:数学建模资料分享 -- 图论部分
背景 今天上午,在教六第一阶梯教室为数学建模俱乐部的同学们分享了有关图论的基本知识和应用. 课后,为同学们留了一个算法实现的小练习,大家可以先做一下.在本图文的末尾处,我把上课的资料以及代码分享出来, ...
最新文章
- java 异常 理解_java中的异常理解
- SBULL块设备驱动程序分析
- decimal 类型数据怎么引用_Java基本数据类型和引用类型
- Python中的这几种报错你遇到过吗?
- 数据库的事务,隔离级别和3大范式
- web前端网页开发离不开的三要素:HTML、CSS、JavaScript
- 一朝不知IP事,一世妄为测试人...
- php xampp 安装yaf,在apache中,yaf访问yaf,apache会异常退出是为何?
- Java中this关键字的使用
- easyui tree复选框是否打钩状态_原创深度:如何利用OpenVINO工具套件监控机器操作员情绪状态(一)...
- http之httpClient工具类
- 游戏音效常见格式以及优缺点
- 重磅:使用UOS微信桌面版协议登录,wechaty免费版web协议重放荣光
- 【Jersey】Jersey框架介绍与使用
- VAM语料库--学习笔记
- linux课程--实验二 Linux 基本命令操作2
- 【JS】json导出到excel,自定义文件名和后缀名
- 谈谈结果导向的做事方式
- 2.4G与5GWiFi频段的区别
- C语言---集合---合取-析取-异或-条件语句-双条件语句
热门文章
- Windows Server 2016补丁更新机制
- linux /dev/null,Shell中 /dev/null和 /dev/null 21
- 50道编程小题目之【完数】
- swift 拖动按钮_Swift 简单控件示例:滑块(UISlider)
- What's the difference between Unicode and UTF-8?
- 各执一词,民用安防市场现状看法PK
- eclipse 基础设置
- bzoj2423[HAOI2010]最长公共子序列
- Redis 的持久化
- 在Linux中实现https访问站点