1. 引言
自然界中,拥有纳米孔洞结构的生物材料包括羽毛、鹿角、皮质骨和硅藻等,这些结构对它们的功能起到了至关重要的作用 [1]。自20世纪20年代,美国工程师M. Raney发现了雷尼镍之后,多孔金属就引起了全世界科学家的关注。近些年来,纳米多孔金属的发展更是突飞猛进,作为一类新型纳米材料,纳米多孔金属孔径在0.1~100 nm之间,其纳米级的双连续孔洞结构和巨大的比表面积,使得它们兼具泡沫材料的结构和纳米材料的功能特性:特别的物理、化学和机械性能、优越的催化性能、表面增强拉曼散射、应变可逆性和高屈服强度,在传感、催化、燃料电池等领域具有广阔的应用前景 [2]。
随着工业发展,人类社会对资源的需求不断增加。资源的可持续利用变得越来越重要。资源短缺问题,引起了人们对绿色化工的重视,尤其是对生产条件更环保且原料成本更低、更易获取的化工产品的生产与加工工艺的重视。相比于传统的催化剂,纳米多孔金具有更高的比表面积以及更多的化学活性位,同时金元素无毒且具有优异的耐腐蚀性,低温下具有良好的催化效果,是一种理想催化系统,对于现代工业技术的革新来说具有重要意义。
能量储存是现代社会面临的一项重大问题。为了解决全球变暖和不可再生化石燃料日益短缺等问题,开发低成本、高效率和更环保的能源转换、储存设备迫在眉睫。设备的性能与材料的性能密切相关,因此材料技术的突破是开发新一代储能和转换设备的关键。近年来纳米技术的发展表明,纳米材料用作锂电池、燃料电池和超级电容器的电极或电解液,表现出独特的性能 [3]。纳米多孔材料用作超级电容器电极材料的关键问题,是在比表面积(确保高电容)和孔径分布(方便电解液流通)之间找到平衡。因此,具有高表面积、高孔隙率和良好导电性的纳米多孔金属,是理想的超级电容器电极材料。
与其他材料一样,在将纳米多孔金属推向大规模工业应用之前,对其力学性能和力学可靠性进行研究是必不可少的环节。目前,人们已经从实验现象上对纳米多孔金属的力学性能有了一定的知识积累,并且以泡沫金属的研究为基础,对其破坏机理和力学行为等领域开展了一些理论研究。
纳米多孔金属的性能表征主要通过试验和理论分析两种手段进行研究,前者需要花费大量人力、物力和财力;后者由于成分、结构的可设计性和制备工艺、加工方法的差异,使得纳米多孔金属的性能数据呈现出极大的分散性,其性能是多种因素互相影响和作用的结果,而材料力学性能与影响因素之间通常是典型的非线性关系,难于用数学模型予以描述。并且,纳米多孔金属的设计与应用主要取决于材料的性能,材料的性能又受到许多复杂因素的影响,这给材料的性能测试带来很大的困难。
随着人工神经网络技术不断地成熟与发展,该技术在纳米多孔金属力学性能建模方面越来越显示出明显的优越性。利用已得的试验数据,将人工神经网络技术应用于纳米多孔金属的力学性能测试之中,力求以较少的试验量获得充分可靠的信息,从而为降低纳米多孔金属的试验和设计成本,减少投资,缩短研制周期,为扩大其应用提供一种新的思路。
2. 纳米多孔金属的结构与力学性质
纳米多孔金属材料骨架的微观尺寸处于纳米尺度,大多数制得的纳米材料呈现出由尺寸大小不一的开孔构成的双连续微观结构,如图1所示,其具有体积比高、导电性好、化学稳定性好、生物相容性好、渗透性好等优点。2001年,Erlebacher等 [5] 通过实验与计算机模拟结合的方法,发现纳米多孔结构的形成是一个内在动力学过程。在纳米多孔结构形成过程中,电位差驱动的活泼原子溶解,惰性原子在界面上的扩散,材料表面积的不断增加,三个因素同时作用。制备纳米多孔金属的主要方法有:模板法、电化学法以及脱合金法等。除了通过实验方法对制备的多孔金属进行力学分析测试外,Newman等人 [6] 通过模拟由A和B两种组分组成的二元流体混合物的旋态分解得到的模型与纳米多孔金属微观结构非常相似,可以在分子动力学(molecular dynamics, MD)模拟用于深度研究纳米多孔金属变形机制和力学性能。
![](//html.hanspub.org/file/8-1280912x9_hanspub.png)
Figure 1. SEM micrographs of dealloyed 25% Au foam by (a) Free corrosion and (b) Potentiostatically driven [4]
图1. (a) 自由腐蚀和(b) 恒电势驱动下脱合金25%纳米多孔金的SEM照片 [4]
为了探索纳米多孔材料的性能,人们进行了大量的实验和理论研究。特别是纳米孔金(nanoporous gold, np-Au, NPG),具有制备简单、抗腐蚀、抗氧化稳定性好 [7] 等优点,是研究纳米孔金属力学行为的理想选择。Np-Au是一种孔洞介于2 nm和10 μm之间的双连续孔洞可调材料,大多数的研究现状表明 [8] [9],其弹性模量在5~13 GPa之间,相对密度和屈服强度的取值范围分别是0.25~0.40和16~150 MPa。Biener等人 [8] 的实验表明,即使相对质量密度低至0.25,np-Au的承载能力(如屈服强度)仍然远远大于大块Au。众所周知,若材料的尺寸小于10 μm,其力学性能与宏观大尺寸情况下相比有明显改变。其主要原因是在微观情况下,位错、孪晶、晶界滑移等塑性变形 [10] 受到长度尺度如晶粒尺寸、试样大小、大应变梯度、钝化层和位错源的缺失等因素 [11] 的直接影响。到目前为止,对纳米多孔材料力学性能的理论研究大多基于Gibson-Ashby模型,然而Gibson-Ashby模型并不能解释微观结构的尺寸效应,因此,一些研究通过实验和计算对方程进行了修正 [9] [12]。此外,Feng等人 [13] 提出了一种考虑表面能和残余表面应力影响的单细胞微力学模型来测试纳米孔材料的尺寸依赖弹性性能,Xia等 [14] 从理论上研究了实验中观察到的具有分层结构的纳米多孔材料的力学性能。
3. 基于神经网络的力学性能分析
实际的工程应用问题往往是非线性的,Gibson-Ashby及其修正的模型都不能很准确的刻画出np-Au的力学性能,反向传播(backpropagation, BP)神经网络研究方法已经广泛应用于材料科学研究的各个领域,然而用于np-Au力学性能定量研究的神经网络模型却未见有报道。究其原因在于影响np-Au力学性能的因素较多且定量化较为困难,而BP神经网络通过自身的自学习自适应性,不需要任何的先验函数或者既定公式,它直接作用于已经获得的样本数据,用训练好的网络模型进行np-Au力学性能的推理测试,从而对np-Au的力学性能进行准确的预报,减少了大量繁琐的步骤和工作量,缩短了研究时间,为np-Au的研究提供了新的思路和方法。
3.1. 人工神经网络的特点
人工神经网络对不确定的非线性关系可以不断学习和适应,来不断的完成自身网络,进而可以学习处理和识别相似的图像和数据,呈现出非常强的自主学习能力。
BP神经网络即目前神经网络模型中,运用的最多的、最具有代表性的误差反向传播多层前馈式网络。BP神经网络一般由输入层(Input layer)、隐藏层(Hidden layer)和输出层(Output layer)三层结构组成,隐层可以有一层或多层。BP算法包括信号的前向传播和误差的反向传播两个部分。正向传播即为输入信号通过隐藏层作用于输出神经元,经过非线性函数变换,产生输出信号。通过计算输出与期望输出的误差,则进入了反向传播过程:输出误差经由隐藏层向输入层逐层反向传播,逐层修改权重值,使得误差沿梯度方向下降,经过反复学习训练,最终获得一个使得误差尽量小的网络参数(权值和阈值),训练结束。BP神经网络具有非常优秀的非线性映射能力,可以通过最简单的三层网络以极高的精度逼近非线性连续函数,对于工程上面临绝大多数的高度非线性映射和不确定性的实际问题有着良好的应用前景。然而,BP神经网络也有如下一些缺点,如:算法可能陷入局部极值,收敛速度慢,学习复杂性高等。
3.2. BP神经网络参数的选择
在实现NPG力学模型的过程中,需要自行设置的超参数一共有:
1) 隐藏层层数,设置为3层。一般的BP神经网络构输入层–隐藏层–输出层。隐藏层数的多少决定了神经网络的精度和计算量大小,然而实际问题的绝大多数均可以由三层BP神经网络来完成,任意的连续函数都可以使用其进行收敛运算。因此在建立NPG力学模型时采用三层神经网络;
2) 隐藏藏神经元的数量,设置为8。若隐含神经元数太少,会使神经网络误差过大,失去了模型建立的意义。隐含神经元数过多,会使得计算量过大训练时间太长甚至不能收敛从而形成死循环。一些经验公式可用于确定隐含层神经元m:
(1)
(2)
(3)
(4)
公式中n, t, α分别代表输入量个数、输出量个数和小于10的非负常数。
确定BP神经网络隐含层的方法一般除了上述的经验公式外还可以根据计算神经网络的均方根误差J来确定,具体的计算公式为:
(5)
公式中m为样本个数,
为期望输出值,y为BP网络实际输出值。经过多次计算得出当隐含层数神经元为8时,BP神经网络模型的均方根误差最小。因此,最终确定神经网络结构为2-8-1三层全连接BP神经网络结构如图2所示。
![](//html.hanspub.org/file/8-1280912x16_hanspub.png)
Figure 2. Prediction BP network model for Nanoporous gold
图2. NPG屈服强度BP网络模型
3) 训练样本和检测样本的划分比例,本文将总数据集60%的样本作为训练集;
4) 训练方法的选择,可以选择的训练方法有梯度下降法、RMSProp、Momentum和Adam优化算法等,本文采用梯度下降法对网络模型进行训练;
5) 训练次数,训练次数设定为5000次。
3.3. 样本数据
本文所搜集的数据均来自于实验,均为毫米尺寸NPG样品的微拉伸实验和纳米压痕实验所得。NPG的宏观弹性模量和流动应力依赖于韧带尺寸、相对密度和微观结构。本文总共搜集了34组数据,包括实际和模拟的NPG的相对密度、韧带尺寸、弹性模量和屈服强度。
对数据的统计分析情况如下:
1) 相对密度:指的是NPG的密度与密实块体金的密度之比,考虑了体积收缩效应。其中,通过脱合金制备和相场法模拟所得的NPG的相对密度都在20%~40%之间。
2) 韧带尺寸:如图3,某一纳米多孔金样品的韧带尺寸并不是统一的,常用韧带直径的平均值来表示韧带尺寸。不同的脱合金方法可以制造不同韧带尺寸的NPG,其中常用HNO3、HCl/NaOH对Au-Ag、Al-Au进行化学脱合金,其中Au-Ag脱合金所制得的NPG韧带大小平均值为15 nm,Al-Au所制得的韧带大小平均值在10~80 nm之间。
![](//html.hanspub.org/file/8-1280912x17_hanspub.png)
Figure 3. Distribution of ligament diameter [9]
图3. 韧带直径分布示意 [9]
表1是对收集到的NPG样本数据的统计表,其中相对密度用小数表示,韧带尺寸用韧带直径表示,单位为nm,屈服强度的单位为MPa。
![](Images/Table_Tmp.jpg)
Table 1. Statistical table of NPG sample data
表1. NPG样本数据的统计表
3) 屈服强度,这里主要讨论实验所得的屈服强度数据。实验所得的屈服强度包括微拉伸所得和纳米压痕所得。纳米压痕中直接获得的数据是np-Au的硬度而非屈服强度,由于应力集中等,其与屈服强度之间的差异是可以预料的,所以需要对纳米压痕实验所得硬度数据进行处理。
在纳米压痕试验过程中,通常会对纳米多孔金属的行为进行假设:纳米多孔金属的压痕可以认为相当于多孔材料的正常压缩。对于相对密度小于0.3的泡沫材料,Wilsea等人 [15] 确定材料在压头穿透过程中不会横向膨胀,并会容纳施加的位移,因为孔隙只是在压头下方塌陷。这一行为对应于近似于零的塑性泊松比(塑性变形时压头下的横向应变与轴向应变之比)的假设,这将代表一种直至韧带塌陷和致密化完全可压缩的材料。这是一个典型的假设,即在塑性变形过程中,多孔结构不受周围环境的约束,所以,屈服强度通常被假定为与测量的硬度相等。一些计算研究发现,多孔材料的塑性泊松比接近于0,但这些研究中的相对密度范围一般在0.1以下。同时,随着塑性泊松比从0增加到0.5,硬度与屈服强度比从1迅速上升到3 (塑性流动过程中不可压缩的致密材料的典型特征:硬度与屈服强度之比接3)。综上所述,如果纳米多孔金属的纳米压痕直接等效于传统多孔材料压缩试验,那么对于单轴压缩的np-Au试样,应测量非常近似为零的塑性泊松比。
然而,这一假设并不成立。Volkert等人 [16] 报道,在压缩np-Au微柱时,测定了非零的泊松比:在测试前后分别测量微柱的直径,并假设在压缩过程中(试样变形主要为塑性)保持恒定,他们计算出平均塑性泊松比为0.20。此外,Nicolas等人的两项 [17] [18] 研究表明,在毫米尺度的拉伸和压缩试验中,连续测量到了弹性状态下的非零泊松比。在单晶np-Au的研究中,弹性下的平均弹性泊松比为0.22,这个值非常类似于由Volkert等人确定的0.20塑性泊松比。由于np-Au宏观极脆,无法直接测量毫米试样在拉伸作用下的塑性变形,因此无法通过拉伸试验确定塑性泊松比。根据Shaw等 [19] 报道的塑性泊松比与硬度-屈服强度比的演化关系,平均泊松比为0.22,表明np-Au的硬度-屈服强度比为2.65。
将该值作为校正因子,类似于文献 [20] 中较早的校正,将该因子应用于本研究中获得的纳米压痕结果,可以得到一组新的屈服强度值。这些新值可以与毫米级拉伸试验结果进行比较。
3.4. BP神经网络训练结果
用上述条件所形成的训练方案,建立相应的纳米多孔金的屈服强度(Yield Strength, YS)的BP神经网络模型,进行训练得到以下结果。
实验样本数据为微拉伸实验和纳米压痕实验所得,已经归一化处理,处理公式为:
(6)
式中:xscaled为归一化后的数据,x为样本数据,xmin和xmax分别样本最小值和最大值。
在神经元层数为3层(1层输入层,1层隐藏层,1层输出层),隐藏层神经元数为8,训练集占总数据集60%,学习率为0.2。训练模型表示通过有标签样本来学习(确定)所有权重和偏差的理想值。在监督式学习中,机器学习算法通过检查多个样本并尝试找出可最大限度地减少损失的模型,损失是对糟糕预测的惩罚。也就是说,损失是一个数值,表示对于单个样本而言模型预测的准确程度。如果模型的预测完全准确,则损失为零,否则损失会较大,训练模型的目标是从所有样本中找到一组平均损失“较小”的权重和偏差。本模型采用的是公式(5),也就是均方差(或者L2)作为损失函数,在训练次数为5000次的条件下,训练得到训练集均方误差为0.0013101,损失曲线如图4,测试集均方误差为0.00290034。
(a) (b)
Figure 4. Loss function value of testing set
图4. 训练集损失值变化
训练集损失值变化曲线如图4,其中(a)是损失曲线的整体下降趋势,(b)为损失值下降细节,可以见到,神经网络的损失值在训练初期下降很快,在训练不到100次的情况下,就从初始的~0.14下降到4 × 10−3以下,训练200次后损失值保持稳定并略有上升,但很快又降低,训练500次后,下降趋于平缓,收敛速度变慢。
图5是训练集数据的真实值和测试值的比较,其中蓝色圆点为真实值,红色折线为测试值的连线。图6是将训练好的神经网络参数用于测试集所得到的测试值和真实值的比较,其中蓝色线条为测试值构成的折线,橙色线条表示真实值。从图5可以看出,神经网络对训练集的拟合效果并不好,但相对来说,图6所示的测试集的拟合效果反而表现得不错,训练集和测试集的损失都在10−3数量级,误差都大于一般所要求的小于10−5。阅读表2,最小误差为3.0%,预测值的平均值非常接近真实值,但同时也看到误差绝对值最大达到了43.04%,回看训练数据,依旧可以发现第2、3、10、14组等数据在数据集中为孤例,若去除误差值在10%以上的孤例,可得到9.12%的误差,误差可降低6.24个百分点。从另一个方面来看,BP神经网络作为数据驱动方法的一个分支,其做出精准预测的基础是较大数量的样本数据,本文收集到的数据有限,可以从某种程度上对本文的方法进行验证,然而样本数据量偏少对本文的预测的准确性也有一定的影响。
![](//html.hanspub.org/file/8-1280912x20_hanspub.png)
Figure 5. Real value and test value of training set
图5. 训练集的真实值和测试值
![](//html.hanspub.org/file/8-1280912x21_hanspub.png)
Figure 6. Real value and test value of testing set
图6. 测试集的真实值和预测值
![](Images/Table_Tmp.jpg)
Table 2. Comparison between predict value based on BP neural network and real value of yield strength for NPG
表2. BP神经网络模型NPG屈服强度预测值和真实值的比较
Xia [14] 和Nicolas [20] 的研究表明,屈服强度随韧带直径的变化呈幂律关系,与相对密度之间呈近似线性关系,Nicolas还给出了用微拉伸实验和纳米压痕实验数据推导所得的屈服强度公式:
(7)
式中*表示泡沫性能,s表示固体性能,ρ*/ρs为相对质量密度,l为韧带直径。
本文和文献 [17] 数据来源都是微拉伸实验和纳米压痕实验数据,同样对纳米压痕的数据都采用2.65的校正因子,表3中列出了通过公式(7)计算的屈服强度值,同时也列出了相同参数下利用BP神经网络预测的数据,通过表中的数据可以看出,本文BP神经网络的平均误差值为7.56%,而式(7)所预测的屈服强度误差平均值高达65.17%,也从一个侧面反映了BP神经网络的有效性。
![](Images/Table_Tmp.jpg)
Table 3. Comparison of result from equation and BP neural network
表3. BP神经网络和公式预测值对比
4. 结论
本文选用BP神经网络模型分析纳米多孔金的屈服强度,通过搜集整理34组力学性能数据样本,建立纳米多孔金以相对密度和韧带直径为输入,屈服强度为输出的三层BP网络预测模型,测试集均方根误差数量级为0.145830 (模拟数据)、0.00290034 (实验数据),屈服强度平均误差控制在10%以内。
通过对纳米多孔金力学性能的模拟和研究分析,可以得出以下结论:
1) 纳米多孔金的力学性能是由其相对密度、韧带尺寸甚至是微观结构决定的,在相同结构条件下,结合纳米多孔金的相对密度和韧带直径典型数据可以用BP神经网络拟合其屈服强度的经验公式。
2) BP神经网络进行纳米多孔金力学性能的预测,与公式推导相比,不需要研究其力学性能和本构参数之间线性关系或者幂律关系,为研究纳米多孔金属的力学性能提供了一个新的思路。
3) 同时可以看到,对纳米多孔金力学性能的研究缺乏有效的数据库,所搜集到的数据比较少,通过BP神经网络作基础性研究,其可靠性和精度还需要大量的数据作为支撑。
基金项目
国家自然科学基金(51401148)。