Solve Finite Element Analysis equation

2020 4 17
1 个回答
回答已采纳
更新时间：2020 4 20
13 次查看（30 天）

0 个投票

Hi,
I'm trying to solve the following FEM equation for 1D unsteady heat conduction using MATLAB. I'm new to FEM and to the mathematical approach in MATLAB. So, I'd appreciate any help/suggestion.
What I need help with is to get the converged solution.
At i=0, T=10 and at i=5, T=50. When n=0, T_i^(n+1)=10, 73, 100, 85, 50 respectively at i=1, 2, 3, 4, 5.
dt=0.1 and (alpha/dx)=2
Again, I'd appreciate any and every suggestion here.
Thank you in advance!

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Jake
Jake 2020-4-19
Sorry about the late response.
Yes, I figured that the equation can be simplified. I will try approaching it that way. Thanks :)
Jake
Jake 2020-4-20
The grid representation? I'm assuming you're referring to the FEM part here and not the MATLAB part.

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

darova
darova 2020-4-20

1 个投票

Here are some of my thoughts

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Jake
Jake 2020-4-20
Yes! I figured the same way but my concern is what to have on left side of the equation (or the right side, because as I mentioned in my original post, having two equations on either sides gives an error in MATLAB.) For instance, T1 and T2 are both unknowns.
Pardon my lack of knowledge in both FEM and MATLAB at this point but I'm learning as I go.
This, however, is exactly the confirmation I needed. Thank you very much, Darova.
darova
darova 2020-4-20
编辑：darova 2020-4-20
As you see we have 3 equations with 3 uknowns(T1 T2 and T3)
All you have to do is to solve this system:
example
a = value
dt = value
dx = value
syms T1 T2 T3 % uknowns
eq1 = 1/6*(10-10)/dt + 2/3*(T1-73)/dt ...
eq2 = ...
eq3 = ...
res = vpasolve([eq1 eq2 eq3],[T1 T2 T3],[73 100 85])
I think it's can be automatized for many rows( n )
Jake
Jake 2020-4-20
编辑：Jake 2020-4-20
Thank you! I have yet to learn vpasolve and now I have a reason to!
This is great. Thank you, again.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Structural Mechanics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by