计算瑞利积分
29日视图(30天)
显示旧的评论
目前我想计算压力的速度分布v(复数)在盘子里。一个好的方法应该使用瑞利积分方程。我已经实现了它,它很好地适合点远离板。
但是我有麻烦评估板上的瑞利积分因为它变得奇异。这是因为表面之间的距离r的评估点压力应该得到零计算。r站在分母上我无法评估。
澄清我修改图片公式和离散的版本我用。整合发生在整个表面板的年代。这将是伟大的听到一个方法来评估板表面上的公式。也许甚至有已经实现的积分方程。我也添加了我现有的代码。
提前感谢任何努力解决这个问题。
q = Element-surface *元素的速度
vorRay =((ω1我* * rhoL * Elsize ^ 2) /(2 *π));
祈祷= 0(长度(点),1);
[xx, yy] = meshgrid (x, y);
岩藏= 1 *(ω/ c);
为o = 1:长度(分)
距离=√(点(1,o) - xx)。^ 2 + (2 o) yy(点)。^ 2 +点(o) ^ 2)”;
祈祷(o) = vorRay *总和(和((v /距离)。* exp(岩藏*距离)));
结束
接受的答案
杰里米Kemmerer
2014年10月1日
你好马提亚,
提到的方法“特殊技术需要使用“计算板本身的压力(看在第五节)。他们给一个参考,可能是有用的。
同样,你看着用数值求积集成而不是明确的金额?一些MATLAB函数可能会对你有用: