材料科学与工程前沿——特种陶瓷材料内容•特种陶瓷简介•特种陶瓷的应用与课题•空天陶瓷•陶瓷制备过程进展2016年3月科技史、人类史上的重大事件计算机(人工智能)计算机不能做什么?材料人应该用计算机做什么?现状如何?数学在材料研究中的作用有没有用?用没用?为什么?计算在材料研究中的作用实验数据处理材料计算计算机在材料研究中的作用1.数据处理软件(excel,origin)2.计算机编程(C、VB)3.科学计算专业平台(MATLAB)4.材料计算专业软件(MaterialStudio,Ansys)5.材料基因组计划材料基因组计划MaterialsGenomeInstitute(MGI)材料创新基础数字化信息平台材料特性——组成、结构、工艺根据应用需求,通过计算选材、配方、选择制备工艺及其条件实验技术、计算技术和数据库数学模型及算法材料基因组3大系统材料大通量计算系统材料性能大通量测试技术系统材料设计性能数据库与信息平台系统材料基因与性能的关系原子参数(半径、价态、电负性等)分子参数(拓扑参数、构象等)凝聚态结构(晶型、径向分布函数等)温度、压力、外场等复杂数学关系(多因子、非线性、高噪声、非高斯分布、非均匀分布等)光、电、磁特性热力学特性力学特性化学特性生物活性。。。。。材料基因组计划的本质依据大量的理论计算及实验数据建立影响因素(基因)与材料性能的数学模型给该数学模型输入各因素的值预测材料性能依据用户对材料性能的需求,反算(设计)材料体系及制备工艺、条件等数学模拟运用数学再现一个系统、过程或现象研究其原理、规律及控制方法数学模型对于一个真实过程,为了一个特定的目的根据内在规律,做出必要的简化假设运用适当的数学工具得到一个数学结构数学模型的内容F(x,y,β)=0y:因变量,x:自变量,β:参数;F:数学结构(代数方程、微分方程)x,y往往是可观测量,数学建模的主要工作是建立合适的数学结构,根据x,y的大量观测值,获得合适的该模型的参数数学模型的分类按经验成分机理模型:基本定律推到而得。偏微分方程半经验模型:主要依据物理定律,同时包括一定经验假设经验模型:以具体系统的考察结果为基础,不反应过程的内部本质和特征扩散过程机理模型)(xCDxtC半经验模型模型22CxCDt经验模型模型taawt102D是具有物理意义的参数;而a0,a1不具物理意义数学模型的分类按表现特征确定性与随机性模型:是否有随机因素的影响静态与动态模型:是否与时间有关线性与非线性模型:微分方程是非线性离散与连续模型:变量是否连续本质上讲大多数问题是随机、动态、非线性数学模型的分类确定性:)(xfy)(xfy随机性:静态:xCDJ22CxCDt动态非线性2210xaxaay线性xaay10经验模型的建立•经验公式bxeyy02210xaxaay模型参数的获取?模型观测:5432154321yyyyyxxxxx拟合参数二元模型F(x)x1x2y线性:2121110xaxaay二次非线性:模型参数的回归分析变量之间的相关关系•自变量x,应变量y•y=f(x,β)•函数f:数学形式和参数β–y=ax3+bx2+cx+d–y=aebx•函数形式表示一般性,具体的参数表示特殊性•如果已知函数形式及其参数,即可获得任意自变量时的因变量值函数的确立•形式的确定:理论或经验•参数的确定:回归分析参数的回归分析•真实函数:y=ax+c•估计函数:ŷ=âx+ĉ•做实验获取一定的x时的y的观测值•x可以是随机的,也可以的确定的;y往往是随机的•通过获得大量的(x,y)数据,估计出函数参数的估计值â,ĉ多项式模型的参数回归•y=f(x)=a0xn+a1xn-1+…+an-1x+an•实验数据:x=[x1,x2,……xm];y=[y1,y2,……ym];mna=polyfit(x,y,n);a=[a0,a1,……,an]非线性数学模型的参数获取非线性函数y=F(x,β);模型的观测值xdata=(x1,x2,….,xn),ydata=(y1,y2,….,yn)搜索合适的参数β使得模型预测因变量y与其观测值ydata的误差平方和最小。(最小二乘)β=lsqcurvefit(fun,β0,xdata,ydata);β0:初始参数不知数学形式的复杂模型•模型涉及多因素、多状态函数(性能)。并且交互作用强。•涉及多种复杂的物理、化学过程•不确定性因素很多•只能实验获得大量的各因素及其性能的观测值,中间作用机制不明。•黑箱模型多元模型F(x)x1x2y1x3x4y2y3数学表达式?复杂;黑箱F11(x)x1x2y1x3x4y2y3F12(x)F13(x)F23(x)F22(x)F21(x)人工神经网络人工神经网络•人工神经网络是由大量处理单元互联组成的非线性、自适应信息处理系统。•人工神经网络是并行分布式系统,采用了与传统人工智能和信息处理技术完全不同的机理,克服了传统的基于逻辑符号的人工智能在处理直觉、非结构化信息方面的缺陷,具有自适应、自组织和实时学习的特点。人工神经网络结构组成温度压力湿度强度韧性化学稳定性神经元权重wi隐单元输入单元输出单元神经元工作原理XiWibiYi)(iiiibxwfyxi:输入向量wi:权重系数bi:阀值yi:输出值f:激发函数可以是线性的,也可以是非线性的激发函数•输入量•线性:y=k*net+c•非线性斜面函数•阀值阶跃函数(k=∞)•S型Sigmoid函数)(netiiibxwnetnetnetknetnetf||*)(]1,1[e11)(]1,0[e11)(22xxxenetfnetf神经网络拓扑结构•多层结构•输入层(0层)•输出层(n层)•隐藏层:隐藏层数多,能够描述更复杂的交互作用,结构越复杂,网络训练数据越多,精度不一定越高。一般1层•各层神经元的个数与连接:增加隐层节点数可提高精度神经网络的训练•用大量的模型的观测值(输入与输出值)训练网络,获得预测误差较小的网络参数(权重)•算法:反向传播算法(BackProagation)(BP网络)人工神经网络的建立与应用•设计网络类型与结构•搜集大量输入、输出数据(x,y)•训练网络,确定网络权重wi•运用网络预测性能(输出)数学模型在间接测试中的应用16:24:07测试表征可观测量x与待测参数β–可观测量直接测得β=x:空间尺度、质量、温度……–可由可观测量用显函数算出β=f(x):强度、密度……–不可由可观测量用显函数算出f(x,β)=0:扩算系数、热导率……–反问题反问题:扩散xCDJ22CxCDt微分方程的求解•边界条件,初始条件•解析解?扩散实验多孔材料PVA溶液tw实验结果多项式拟合结果扩散系数?扩散数学模型1222299)/exp(161nnntqrtDqMM22CxCDt该实验条件下的解233)tan(nnnqqq已知:r、p、V、Co实测:Mt、t待定:D非线性拟合可获得反问题f(x,β)=0的求解方法建立隐函数F=f(x,β),用数值分析求解F=0时的β0β是n维,观测值x为m个m=n:恰定方程,唯一解β0mn:超定方程,最小二乘法,求最优β0MATLABβF0β0热导率稳态、均质:xqxTkxtTrccTAQ/问题:–高隔热材料温度梯度很大–高温辐射传热的贡献高温传热有限体积法模型将隔热材料划分为n层i层的能量方程导热热流本身辐射热流本身辐射热流得到的辐射热流得到的辐射热流导热热流4'''kikiTqqxTTqkikiki11I''xTTqkikiki11I'层间的辐射换热其中ɛ为辐射发射率,ɣ为反射率,θ为透过率方程组联立4k1-n1-4k2-n2-4k334k224k114kn4k1-n1-4k444k334k22n1-n432-1-n-2-n-3-2-11-n1-n2-n2-n33221n1-n1-n443322)(T)(T......)(T)(T)(T)(T)(T......)(T)(T)(TGG......GGGGG......GGG1--1--............1--1--1--1--1............--1--1--1nnnnkkkkkkkkkk求解数学模型f(x,β)=0观测值x:表背温T、热流量Q、试样尺寸传热参数β:热导率a+bT、发射率ɛ、透过率θ4个待定参数,4套以上观测值x,最小二乘法求出其他数值分析的应用微分方程的求解(有限元方法)插值数值微分。。。。。。结论数据挖掘是高水平研究的关键除了会用excel,origin,更要会用MATLAB