| 查看: 1270 | 回复: 2 | ||
herozhu木虫 (正式写手)
|
[求助]
关于MATLAB有限差分模拟激光加工过程的边界条件
|
|
在用matlab编写程序模拟激光加工过程，由于涉及相变，故选用了热焓法（enthalpy method），即某一方格（element）的能量达到某一值，该单元格发生相应变化。例如图中所示，用方格中心点的温度来代表整个单元格，如单元格（1，n）的能量超过了汽化对应能量，则认为该单元被去除掉。暂未考虑液态飞溅、等离子体等等因素。 关于控制方程，是从单元格能量守恒推导而得 关于边界条件，简单起见，左右下三边都设为绝缘，上表面考虑对流与辐射。激光为纳秒脉冲，高斯分布，从上表面入射。 一直没有考虑明白的是上表面上的边界条件问题。我一直将激光光源作为“体热源”处理，即在每一个单元格内，根据激光能量密度，在x,z两方向及时间t上做积分，带入单元格能量守恒公式中，从而推导下一时刻单元格的enthalpy，并得到相应温度T。但对于辐射及对流计算，必须使用材料表层的温度进行计算。则具体到这一层上，边界条件应为 之前的做法是暂不考虑表层，对所有单元格进行求解，得到下一时刻的单元格温度，再由第一层单元格温度按照对流及辐射（不考虑激光）写出方程fsolve解出该时刻的表层温度（相当费时），但观察结果有问题，会出现单元格能量条跳跃或是温度低于初始值，考虑是边界条件的问题。希望能得到大家的意见及建议，此过程的任意方面都欢迎讨论。  20131215_153418.jpg |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求解微分方程
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 该如何列全偏微分方程组的边界条件？
已经有14人回复
- 有限差分离散二阶偏微分方程后如何编译matlab程序解方程组？
已经有3人回复
- 有限元第二类边界条件 matlab仿真怎么编
已经有3人回复
- 求Matlab仿真模拟一级、二级倒立摆模型方程及matlab仿真模拟程序
已经有6人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- [交流]分享Matlab写的PML吸收边界条件。代码很精简
已经有11人回复
- 怎么用matlab求出满足以下条件的所有情况
已经有4人回复
- 【求助】同样一个问题用matlab和fortran编程计算，结果相差比较大，不知道为什么
已经有4人回复
- 【求助】偏微分方程边界条件的处理
已经有4人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】请教虫友关于matlab差分方程的求解和画图程序
已经有4人回复
- 【原创】关于MATLAB 描述分子布朗运动等行为的过程
已经有59人回复
- 【求助】用matlab编过有限差分法的请进
已经有9人回复
- 【求助】关于MATLAB中元胞自动机的边界条件和均匀分布
已经有4人回复
herozhu
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2408.6
- 散金: 31
- 红花: 1
- 帖子: 514
- 在线: 118.9小时
- 虫号: 1772719
- 注册: 2012-04-23
- 性别: GG
- 专业: 零件加工制造
2楼2014-03-26 12:38:11
3楼2015-03-15 20:41:33
10