| 查看: 2822 | 回复: 7 | ||
[求助]
偏微分方程Matlab求解 已有1人参与
|
各位路过大神求帮忙，Matlab求下面的一个模型，方程和模型见图。我用matlab编程求解，结果有错误，下图2中的红色框内的数据按程序逻辑应该为0，但是计算出来不是，导致后面数据都有问题，而且数据每迭代几步以后就会增加一个数量级，问题在哪我也查不出来！希望会的帮运行下看看，是数据类型不对导致的计算误差，还是因为差分格式有问题！万分感谢，如果能解决问题金币不是好说，可以追加！！！！ 方程.jpg  图片2.jpg[ Last edited by 小学男生 on 2014-2-21 at 11:38 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : kyisunjian0220.m
2014-02-21 11:38:15, 2.62 K
» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
变系数偏微分方程matlab求解
已经有6人回复
- 一个偏微分方程组的求解问题
已经有5人回复
- matlab 非线性微分方程求解
已经有3人回复
- 偏微分方程求解！
已经有4人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- matlab解偏微分方程求助
已经有12人回复
- 求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~
已经有4人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 抛物型偏微分方程的数值解求解问题，谢谢！
已经有14人回复
- Matlab求解二阶偏微分方程组，希望能给出相应的求解代码
已经有10人回复
- 二阶非线性微分方程求解
已经有14人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- 有限差分离散二阶偏微分方程后如何编译matlab程序解方程组？
已经有3人回复
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- 三元偏微分方程求解
已经有4人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- 偏微分方程数值解有基本理论和一些例子（化工类）MATLAB源程序
已经有296人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】偏微分方程的基本解
已经有5人回复
- 【求助】用MATLAB 求微分方程
已经有18人回复
2楼2014-02-21 22:53:10
3楼2014-02-22 13:52:43
4楼2014-02-23 02:59:36
|
首先谢谢你如此细致的看！ y1代表的是图模型中横向放置的大矩形（主体）的位移，它是时间t的函数，初始条件y1(0)=0,dy1(0)=-3.1305(即自由跌落触地前瞬间的速度：负根号下2gH). y2(程序中使用的是U，y2在程序中为主体的速度)代表的是图模型中纵向放置的小矩形(易损件)的位移，它是时间t和x（以易损件和主体结合部位中点为坐标原点，竖直向上为正方建立x坐标轴）的函数。y2对时间的偏导代表易损件的速度，对x的偏导代表易损件的内应力分布。 初始条件y2(x,0)=0,Py2(x,0)(y2对时间的偏导)=-3.1305(即自由跌落触地前瞬间的速度：负根号下2gH). 边界条件y2(0,t)=y1(t),即x=0处的位移始终与主体的位移相同。Py2(x=L,t)=0,即x=L处（易损件自由端）的内应力始终为0. 系统中其他参数如下：m=3kg，易损件密度p=500kg/m3,弹性模量E=10^7Pa,易损件长度L=0.03，杆的横截面积A0=5*10^-6m2,a1=280，a2=8.94*10^4,a3=1.91,A=0.01m2,厚度h=0.035m，跌落高度H=0.5m，上述参数在y1的求解中我基本都带到方程中化简了，求解y1时不考虑红色方框中部分。二阶双曲偏微分方程中的alpha^2=E/P(弹性模量/密度)。 上述描述可能比较混乱，因为都是基于图中方程描述的，而在程序中，可能字母有变化，程序中部分有注释，如果还有看不明白的再告诉我。谢谢 |
5楼2014-02-23 10:38:27
|
楼主你好，你也说的很详细 不过我这里有几个疑问 首先，你把坐标设在结合部位，这一点是运动的还是一直在初始时刻的位置？如果是前者，那易损件实际是处于一个非惯性系，运动方程既要考虑其相对主体的运动，也要考虑主体本身的运动规律。 其次，即便你说了关于易损件的两个边界条件，我还得说，你少了一个关于易损件位移对时间一次导数在某个边界处的边界条件，也就是易损件在各个时刻某个边界处的速度要满足的条件。我没有看你的程序，但我看你的这个帖子，似乎依然没有提到 我还是觉得你的坐标系的设置有点不清晰，至少没有很明确的说明白。 方程本身来说，真心不复杂，那个双曲方程用有限差分法也是很容易算的，至于主体的常微分方程，也是很简单的 再好好审视你模型本身，尤其是一些基本的假定和坐标设置，一定要整理清楚，然后问题可能就解决了。 |
6楼2014-02-24 12:59:04
|
你分析的很对，尤其是坐标的设置。这个问题是别人发表过的一篇论文中（论文中只有结果与方程，导师要求我将过程以及程序写出），我也有疑问，文章中系统里是将y2作为易损件的绝对位移处理的（即易损件相对于跌落地面的位移，而不是主题），但是绝对位移对x的一阶导却又表示为易损件的内应力（不应该是相对位移吗），这个疑问想向你请教。 对于你提到的少了一个“关于易损件位移对时间一次导数在某个边界处的边界条件”，是我没有说清楚，易损件结合部位“y2(0,t)=y1(t),即x=0处的位移始终与主体的位移相同”，同时速度也相同即Py2(0,t)=dy1(t)。 y1以及速度dy1，我已经通过龙格库塔法数值求出，针对于易损件的双曲偏微分方程，只求位移的话直接利用差分法我也求出来了，但是要求易损件速度(y2对时间的一阶导)和内应力(位移对x的一阶导)，需要将双曲方程降阶为一阶耦合的双曲偏微分方程组，在编程求解的时候，原帖图中红色方框内的V2（内应力）数据，程序逻辑应该是0，但是结果不是，后面数据就开始波动了，而且每迭代2-3数量级就会增加一个，到后面数据就大的离谱，我纠结的在这个地方！ |
7楼2014-02-24 21:41:21
8楼2014-02-24 21:45:25
