对于变上限积分的求解
求裂纹附加柔度矩阵时，遇到的变上限积分，不知道大家都是怎么解决的？？？ 原帖由 malong 于 2007-5-29 10:00 发表 http://www.chinavib.com/forum/images/common/back.gif求裂纹附加柔度矩阵时，遇到的变上限积分，不知道大家都是怎么解决的？？？
我是不懂的了，xjzuo 版主、rocwoods 兄 在积分方程求解问题上都颇有心得，让他们尽情发挥好了:loveliness: 把方程贴出来吧,让大家一起分析. 大家可以先看看，不知道数值积分工具箱好用不，还有不少人编程序来解这个方程！ 其中的F2是什么? 不好意思，F2忘掉了 现在补充 请自己先动手写一下表达式, 或自己计算的代码.
另:这个问题应该可以参考我的一个帖子:"关于带参数的积分问题". 必须指出的是,你这不是所谓的变上限积分,而就是普通的二重积分.上下限可以是确实的数,也可以是外圈变量的表达式,你这个就是这种情况.一般地,上下限是确实的数时的情况比较容易求出,matlab 中就有dblquad来求.上下限是表达式时就必须自己编了,你可以去查查.
给定的R,a,E,v,alpha(中间那个还像是alpha吧),应该积出来的就是一个确实的值.我用maple只能积出第二项,第一项积不出来,需要特殊的处理,找找matlab这方面的书,应该会有的. 仔细看了一下word文档,发现果然不是变上限积分，而是一般情形下的二重积分．
这个问题可以用＂数值积分工具箱＂函数quad2dggen求解，工具箱可从Mathworks网页上下载．
当然，自己编程也可解决，只是要花点时间． 好的 我再试试 谢谢大家！
[ 本帖最后由 eight 于 2007-6-1 09:36 编辑 ]
页:
[1]