| 查看: 1922 | 回复: 10 | |||
| 【悬赏金币】回答本帖问题，作者kevin123581将赠送您 10 个金币 | |||
kevin123581金虫 (正式写手)
|
[求助]
激光加工温升计算，一维傅里叶传导或者双温模型，福柯普朗克模型
|
||
小弟本科学的机械，现在跟的导师是做激光加工的，从0开始在学。但是导师叫我做一个脉冲激光的照射温升计算，弄了3个月了，还是各种不对。12ps，10^14 w/m^2的单脉冲强度我算下来居然只升高14度，材料是碳化硅。 双温模型的图也是怪怪的。所以想请问下各位前辈有没做过这方面研究的能指导一二，最好能借下程序，因为是半路出家，matlab也不是很熟练，但是感觉又没什么问题，结果就是不对。 望大家相助啊，实在无脸见导师了，谢谢 方便的话 请加我qq 584566830 ，大谢了  一维傅里叶传导，图形呈线性，最高点太低，脉冲结束了还没有温降  和上面同样的程序，改了脉冲长度，强度和光吸收系数。图形形状不错，但是没温升几乎。  双温模型的程序图，4h-sic材料，晶格温度也是基本没变 [ Last edited by kevin123581 on 2012-3-13 at 17:06 ] |
回复此楼» 猜你喜欢
难道我真的要放弃吗？
已经有44人回复
-
【45岁以上博士】柔性人才引进项目，有补贴
已经有18人回复
-
Pattern Recognition期刊，二审审稿邀请发出一个月了，但有一个审稿人一直未接受审稿
已经有12人回复
-
神奇的中医
已经有10人回复
-
怎么带研究生？
已经有39人回复
-
脱酯基合成步骤
已经有7人回复
-
省基金对国自然申请是有利还是有弊
已经有5人回复
-
博后面上今天下午会公布吗？有无消息？
已经有41人回复
-
评审专家会不会很在意申请人的单位啊
已经有13人回复
-
老博士选高校
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 济南大学材料学院招聘信息（微纳粉体专业）
已经有22人回复
2楼2012-03-13 21:32:20
kevin123581
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 761.3
- 散金: 81
- 红花: 6
- 帖子: 539
- 在线: 238.7小时
- 虫号: 1658947
- 注册: 2012-03-02
- 性别: MM
- 专业: 光学
3楼2012-03-14 11:44:14
龙山狮崖
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 5877.1
- 红花: 6
- 帖子: 391
- 在线: 249.4小时
- 虫号: 1683187
- 注册: 2012-03-11
- 专业: 凝聚态物性 II ：电子结构
4楼2012-03-14 14:30:00
kevin123581
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 761.3
- 散金: 81
- 红花: 6
- 帖子: 539
- 在线: 238.7小时
- 虫号: 1658947
- 注册: 2012-03-02
- 性别: MM
- 专业: 光学
5楼2012-03-14 16:01:31
龙山狮崖
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 5877.1
- 红花: 6
- 帖子: 391
- 在线: 249.4小时
- 虫号: 1683187
- 注册: 2012-03-11
- 专业: 凝聚态物性 II ：电子结构
6楼2012-03-14 18:36:22
kevin123581
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 761.3
- 散金: 81
- 红花: 6
- 帖子: 539
- 在线: 238.7小时
- 虫号: 1658947
- 注册: 2012-03-02
- 性别: MM
- 专业: 光学
7楼2012-03-14 19:03:24
kevin123581
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 761.3
- 散金: 81
- 红花: 6
- 帖子: 539
- 在线: 238.7小时
- 虫号: 1658947
- 注册: 2012-03-02
- 性别: MM
- 专业: 光学
|
%热源Q项当成偏微分方程里面的项 close all L=3e-4; L=3e-4; dx=L/200; dx=9.829e-7 % 空间步长 x=0:dx:L; x=x'; Nx=length(x); tend=50e-12; dt=0.5e-12; %时间步长 t_1=12e-12; %12ps去掉加热 t=0:dt:tend; % D=1.81e-4; %热扩散系数 cp=690; %比热 rho=3200; %密度 I0=1.69e14; %激光能量密度 w/m2 R=0.27; R1=1-R; %反射率 alpha=21000; %光吸收系数 beta=D*dt/dx/dx gamma=I0*R1*alpha*exp(-alpha*x)/rho/cp*dt; T0=300; T=T0*ones(size(x)); %一行一个温度 A=diag(-beta*ones(Nx-1,1),-1)+diag((1+2*beta)*ones(Nx,1))+diag(-beta*ones(Nx-1,1),1); A(1,1)=1; A(1,2)=-1; A(end,end)=1; A(end,end-1)=-1; for i=1:length(t) B=T(:,i)+gamma; B(1)=0; B(end)=0; T(:,i+1)=A\B; if (t(i)>=t_1) gamma=zeros(size(x)); end end plot(t*1e12,T(1,2:end)) xlabel('t(ps)') ylabel('t_s_u_r_f') =============================== 原方程 dT/dt=D*(d^2T/dt^2)+I0*(1-R)*alpha*exp(-alpha*x)/(cp*rho) IC: T(x,t=0)=300K BC: dT/dx(x=0)=dT/dx(x=L)=0 |
8楼2012-03-14 19:11:55
龙山狮崖
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 5877.1
- 红花: 6
- 帖子: 391
- 在线: 249.4小时
- 虫号: 1683187
- 注册: 2012-03-11
- 专业: 凝聚态物性 II ：电子结构
9楼2012-03-14 20:01:25
kevin123581
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 761.3
- 散金: 81
- 红花: 6
- 帖子: 539
- 在线: 238.7小时
- 虫号: 1658947
- 注册: 2012-03-02
- 性别: MM
- 专业: 光学
10楼2012-03-15 11:34:29
5