最近在作论文的时候,碰到一个积分,由于积分式很长,好像还是不可积的,用了好几种方法也没有积出来,希望大家能帮我看看,赐教一下,小弟先谢谢大家了!程序在下面
-----------------------------
Q=10; % 冷却率
T=980; %初始温度
Ae3=870;
time=1; % 时间步
T=T-Q*time; %温度的改变
% 此处开始if语句
if T<=Ae3
Tc=Ae3-T; % 过冷度
syms T; %声明符号变量
% 计算Dr=Dr0*exp(-Qr/(R*T))的参数
Dr0=1.75/10^5;
Qr=143.32;
R=8.3144;
%K1,K2
K1=2.07/10^5;
K2=6.33/10^15;
%k是Boltzmann常数(玻尔兹曼)
k=1.38/10^23;
%超元素S的自由能中的相关参数值
Tmsi=0;
Tnsi=-3;
Tmmn=-39.5;
Tnmn=-37.5;
%Xsi,Xmn为置换型合金元素si,mn的摩尔分数
Xsi=0.2;
Xmn=0.3;
%超元素S的活度系数的相关参数值
%Wr为奥氏体中碳原子间的交换能
Wr=1250;
%Xc是奥氏体中碳元素的初始摩尔分数
Xc=0.5;
%计算:超元素S的活度系数as
Zr=14-12*exp(Wr/(R*T));
%奥氏体到铁素体相变中超元素S的自由能
Gfe=20853.06-466.35*T-0.046304*T^2+71.147*T*log(T);
Gs=141*(Xsi*(Tmsi-Tnsi)+Xmn*(Tmmn-Tnmn))+Gfe;
%铁素体形核的驱动力
Gn=Gs-R*T*(log((1-Zr*Xc)/(1-Xc)))/(Zr-1);
%奥氏体中碳溶质的扩散率
Dr=Dr0*exp(-Qr/(R*T));
% I是在某一温度T时奥氏体晶粒表面单位面积的铁素体形核率
I=K1*(1/sqrt(k*T))*Dr*exp(-K2/(k*T*(Gn)^2));
%f为被铁素体占据的奥氏体晶粒表面分数
f=0;
%积分函数
ft=(I*(1-f))/Q;
%此处是一个关于ft的积分函数,积分限是(0,Tc),积分变量是(Ae3-T),该如何对 ft 函数进行积分呢?????
%Sr 为该元胞的表面积
%Lca,两个近邻元胞的距离。
Lca=20;
%此处元胞类型取为四边形元胞(正方体)
L=Lca;
Sr=Lca*Lca;
nuc_pf=Np*Sr;
end
------------------------------------------
麻烦大家帮我看一下了。谢谢了