没有合适的资源?快使用搜索试试~ 我知道了~
利用相场法模拟金属晶粒生长的计算代码
需积分: 48 141 下载量 78 浏览量
2016-03-20
16:49:37
上传
评论 15
收藏 4KB TXT 举报
温馨提示
试读
3页
利用相场法模拟金属晶粒生长的计算代码
资源推荐
资源详情
资源评论
%模拟镍的相场计算(无任何耦合、粗网格、变边界条件)-----------------------------
clear;clc;clf;tic
%参数----------------------------------------------------------------------
m=0.05; eps=0.005; r=0.05; %动力学系数、界面相关参数、各向异性系数
Tm=1728; L=2350; DT=0.155; Cp=5.42; %熔点、潜热、热扩散率、比热
Delta=-0.6; R=0.05; alfa=400; %无量纲过冷度、形核半径、耦合相关参数alfa
EPS=eps^2;A=m/EPS;B=(m/EPS)*EPS;C=eps*alfa*(-Delta);D=30/(-Delta);
%网格划分-------------------------------------------------------------------
N=100; NTimeSteps=20000; Dx=0.02; Dy=Dx; Dt=1e-5; E=(m*Dt)/eps^2; DXX=Dx^2;
[x,y]=meshgrid(0:Dx:2); t=20; T=0; h=Dt/Dx/Dx;
%初始条件------------------------------------------------------------------
nn = length(x);
for i=1:nn
for j=1:nn
%phy_n(i,j)=0.5*(tanh((sqrt(x(i,j).^2+y(i,j)^2)-R)/(2*sqrt(2)*eps))+1);
phy_n(i,j)=0.5*(tanh((sqrt((x(i,j)-x(N/2,N/2)).^2+(y(i,j)-y(N/2,N/2))^2)-R)/(2*sqrt(2)*eps))+1);
U_n(i,j)=Delta* phy_n(i,j);
end
end
%主程序---------------------------------------------------------------------
for n=1:NTimeSteps
%边界条件------------------------------------------------------------------
phy_n(:,1)=phy_n(:,2); phy_n(:,N+1)=phy_n(:,N);
U_n(:,1)=U_n(:,2); U_n(:,N+1)=U_n(:,N);
phy_n(1,:)=phy_n(2,:); phy_n(N+1,:)=phy_n(N,:);
U_n(1,:)=U_n(2,:); U_n(N+1,:)=U_n(N,:);
%--------------------------------------------------------------------------
for i=2:nn-1
for j=2:nn-1
phy_x=(phy_n(i+1,j)-phy_n(i-1,j))./2./Dx;
clear;clc;clf;tic
%参数----------------------------------------------------------------------
m=0.05; eps=0.005; r=0.05; %动力学系数、界面相关参数、各向异性系数
Tm=1728; L=2350; DT=0.155; Cp=5.42; %熔点、潜热、热扩散率、比热
Delta=-0.6; R=0.05; alfa=400; %无量纲过冷度、形核半径、耦合相关参数alfa
EPS=eps^2;A=m/EPS;B=(m/EPS)*EPS;C=eps*alfa*(-Delta);D=30/(-Delta);
%网格划分-------------------------------------------------------------------
N=100; NTimeSteps=20000; Dx=0.02; Dy=Dx; Dt=1e-5; E=(m*Dt)/eps^2; DXX=Dx^2;
[x,y]=meshgrid(0:Dx:2); t=20; T=0; h=Dt/Dx/Dx;
%初始条件------------------------------------------------------------------
nn = length(x);
for i=1:nn
for j=1:nn
%phy_n(i,j)=0.5*(tanh((sqrt(x(i,j).^2+y(i,j)^2)-R)/(2*sqrt(2)*eps))+1);
phy_n(i,j)=0.5*(tanh((sqrt((x(i,j)-x(N/2,N/2)).^2+(y(i,j)-y(N/2,N/2))^2)-R)/(2*sqrt(2)*eps))+1);
U_n(i,j)=Delta* phy_n(i,j);
end
end
%主程序---------------------------------------------------------------------
for n=1:NTimeSteps
%边界条件------------------------------------------------------------------
phy_n(:,1)=phy_n(:,2); phy_n(:,N+1)=phy_n(:,N);
U_n(:,1)=U_n(:,2); U_n(:,N+1)=U_n(:,N);
phy_n(1,:)=phy_n(2,:); phy_n(N+1,:)=phy_n(N,:);
U_n(1,:)=U_n(2,:); U_n(N+1,:)=U_n(N,:);
%--------------------------------------------------------------------------
for i=2:nn-1
for j=2:nn-1
phy_x=(phy_n(i+1,j)-phy_n(i-1,j))./2./Dx;
资源评论
dayou2016
- 粉丝: 7
- 资源: 2
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 2SJ343-VB一款P-Channel沟道SOT23的MOSFET晶体管参数介绍与应用说明
- 一款基于Matlab的像素级图像标记工具,支持FITS、JPG、PNG、TIF等图像格式,可导出所有标记点坐标+源代码+文档说明
- 基于javaweb的在线答题平台.zip
- 三菱PLC例程源码丹佛斯变频器与FX通讯实例(适用于所有丹佛斯变频器)
- 微信每日早安问候(用微信传递爱情).zip
- 基于matlab-模板匹配-角度修正的铭牌识别+源代码+文档说明
- 三菱PLC例程源码丹佛斯VT2800与FX2N用RS485通讯示例
- 微信 macOS 客户端增强 Tweak 动态库.zip
- UAM概念飞机 UAM概念飞机
- 一个使用 Java 开发的、类似 WordPress 的产品,支持多站点、多语种自动切换等.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功