完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>

|
问题描述:
为了实现《液氯运输泄漏事故扩散风险分析》
液氯运输泄漏事故扩散风险分析.pdf
(1.05 MB, 下载次数: 160
)
论文中的一个改进的高斯烟团模型,移动点源的高斯烟团模型。 主要公式: 主要公式中涉及到的一个式子: 绘图最后结果: 而我的结果却是不管车速v设的有多大,浓度值最大点的横坐标总是21附近,只受风速的影响不受车速的影响,为了便于运算我的积分时间仅为10s:。 我的程序如下:%% clc clear %% %根据论文中设置参数 t=10; % 10秒 H=1; % 泄露点距离地面的高度为1m u=2; % 平均风速u=2m/s v=12.5; % 车速v=45km/h=12.5m/s z=1.5; %% for x=1:130 %deta_x、deta_y、deta_z由大气稳定度决定 E级 deta_x=(0.06* x *((1+0.0001*x).^(-0.5))); deta_y=deta_x; deta_z=(0.03 * x * ((1+0.0003*x).^(-0.5))); for y=-100:100 s(x,y+101)=integral(@(tba) (((6500000.*(1-0.007.*(tba.^0.5)))./(((2.*pi).^1.5).*deta_x.*deta_y.*deta_z)).* (exp(-((x-u*t+u*tba-v*tba)).^2) ./(2.*(deta_x.^2))).*(exp(-((y).^2)./(2.*(deta_y.^2)))).*(exp(-((z-H).^2)./(2.*(deta_z.^2)))+exp(-((z+H).^2)./(2.*(deta_z.^2))))),0,10); end end x=1:130; y=-100:100; s=s';%转置 s=flipud(s);%上下翻转 save test_movecontinue1_5.mat load test_movecontinue1_5.mat %移动点源持续泄露 figure surf(x,y,s);%画三维图 [cs, h]=contour(x,y,s,'LevelList',[6 90 300 2500]); %1.5m时画出值为等值线 grid minor,xlabel( ' x-axis ' ),ylabel( ' y-axis '),zlabel( ' z-axis ' ); title( '1.5m高度的扩散情况' ) ; |
|
相关推荐
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
我用matlab的coder 封装了一个c语言的dll, 但是在用labview调用时,会出现识别不到库的问题,有大神遇到过吗
2673 浏览 0 评论
5246 浏览 0 评论
在matlab中如何计算含有第一类修正的贝塞尔函数的积分算不出的问题?
9043 浏览 0 评论
怎么利用matlab得到95%,80%和70%的置信区间,并生成不同区间下的功率误差贝塔分布?
10485 浏览 0 评论
请问simulink的s-function模块如何添加多输入输出接口
14660 浏览 2 评论
/9
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2025-12-3 13:22 , Processed in 0.645405 second(s), Total 73, Slave 54 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191

淘帖
6497