摘 要：低速冲击后压缩强度是复合材料层合板的重要性能指标.在分析CCF300/5428复合材料层合板遭受不同能量等级低速冲击后压缩强度试验数据的基础上,提出一种基于广义回归神经网络技术的低速冲击后压缩强度估算模型.该模型以冲击能量、凹痕深度以及损伤面积为输入参数，以高斯函数为隐含层激励函数，使用部分试验数据训练，寻找最优光滑因子.此外，以最优光滑因子对所提模型进行重构并采用部分试验数据对该模型进行验证.结果表明：基于广义回归神经网络的模型具有较好的试验数据泛化能力，可用于估算低速冲击后复合材料的压缩强度.
关键词：神经网络；复合材料；低速冲击；压缩强度；估算
复合材料具有质量轻、强度高、可设计性好、成本低以及耐蚀耐候等优点，已在航空航天等领域得到了广泛应用[1-8].复合材料层合板是由2层或2层以上单层交错的材料粘合组成的结构板，每个单层又由平面铺设的增强体纤维和基体材料构成.研究表明[9-12]，复合材料层合板受到低速冲击后，其性能尤其剩余压缩强度会急剧下降，这在很大程度上会削弱构件的承载能力，并破坏材料结构的完整性[7,13].因此，研究低速冲击后复合材料层合板损伤程度对剩余压缩强度的影响以及剩余压缩强度的估算方法对复合材料结构设计具有重要的意义.
再来看看《现代汉语词典》中这类词的释义情况。[注]如没有特别说明,本文所引用的释义均依据中国社会科学院语言研究所词典编辑室:《现代汉语词典》,北京:商务印书馆,2012年6月第6版,2016年9月第7版。
复合材料层合板低速冲击后压缩强度的研究方法主要有4种.① 损伤演化法：利用动态有限元计算层合板的冲击损伤，以其对应的退化后刚度作为层合板的初始损伤，再用损伤演化模拟层合板的压缩破坏过程并计算剩余压缩强度[14-15]；② 软化夹杂法：将低速冲击后损伤视为规则形状的夹杂软化,再运用各种失效判据确定层合板的失效强度[16]；③ 子层屈曲法：将低速冲击后损伤看作多个面积不同、形状规则的分层,假设压缩破坏过程是各分层不断产生屈曲失效的过程,并认为所有分层屈曲破坏时材料结构发生破坏，以确定失效强度[3]；④ 开孔等效法：使用圆孔或椭圆孔等效替代低速冲击后损伤,再使用材料断裂韧性的判据来确定层合板破坏时的强度，作为失效强度[17].
本文以沪深两市高技术上市公司为样本，研究了管理层能力对创新投入的影响，并对比考察了不同产权性质下市场竞争和高管激励机制对二者关系的调节作用。结果发现：（1）高能力管理者抑制企业的研发投入，在国有背景企业中这种负向作用更明显；（2）充分的市场竞争可以有效缓解民营企业高能力管理者对研发投入决策的风险规避特征，而在国有企业表现为抑制性； （3）工资薪酬仍然是缓解代理冲突、促进高能力管理者开展创新活动的有效途径，而股权激励的促进作用尚不明显。
本文分析不同能量等级的低速冲击后CCF300/5428材料体系的复合材料层合板压缩强度数据，提出一种基于广义回归神经网络技术的低速冲击后复合材料层合板压缩强度的新型预测方法.考虑到样本试验数据的特殊性及样本容量的局限性，对影响复合材料层合板低速冲击后压缩强度的材料体系力学属性、铺层顺序、铺层数量、增韧技术、几何尺寸、冲头形状、冲击能量、环境等因素进行简化.在此基础上，利用广义回归神经网络建立估算模型，运用试验数据训练并结合轮询遍历交叉检验法计算出高斯函数的最佳光滑因子.重构了神经网络估算模型并进行层合板压缩强度估算验证.经过对比分析模型估算数据与试验数据，认为估算结果较为准确、可靠，且估算过程快速而简单.
党的群众路线，强调党员自身的认识和对党的认识，更重要的一点是更为科学、理性地认识群众的意愿。联系群众，团结群众，毛泽东曾经明确指出：“要团结一切可以团结的力量。”这是毛泽东在革命抗战时期所体现的伟大的马克思主义胸怀。这句话在当代也是我党需要坚定的信念，广大党员干部必须再次认清我党不仅仅是中国的执政党，更是中国社会的建设党和社会主义事业发展的领路人。邓小平曾明确提出：“把党建设成为有战斗力的马克思主义政党，成为领导全国人民进行社会主义物质文明和精神文明建设的坚强核心。”［4］9
1 试验
转借图书参与抽奖是奖转借系统与微信抽奖相结合，即只要读者手中有尚在借阅状态的图书发布出去，未有借阅的图书，则可以发布求书，一旦发布转借或是求书的信息后，读者的微信公众号即可转入微信抽奖平台的抽奖池，获得参与抽奖的机会。
按照ASTM D7136落锤冲击试验标准进行复合材料层合板低速冲击试验.试件尺寸为100 mm×150 mm×8 mm，试验区域的尺寸为75 mm×125 mm×8 mm.试验环境条件为干态室温，试验组数为24，冲头选用直径为 12.7 mm的半球形冲头.冲击能量(E)、冲击试验后试件表面凹痕深度(D)、冲击损伤面积(S)以及根据ASTM D7137标准测得试件冲击后剩余压缩强度(pc)数据见表1.
表1 冲击试验数据
Tab.1 Data of the impact test
试验组E/JD/mmS/cm2pc/MPa132.200.6617.20159.67236.301.0019.65158.26341.401.5218.13132.15447.302.0216.03138.98538.900.9618.88152.35644.401.8221.87124.49740.201.2419.46147.03829.500.4816.00147.86935.200.528.51168.651035.900.6814.48153.671141.100.9010.91165.681246.702.109.81126.251343.301.3419.33146.451444.401.6214.83131.521540.301.2419.34148.941645.401.3820.93136.791765.302.9219.91122.591859.003.0022.95132.351953.601.8018.94127.862055.402.2019.77129.912148.201.3022.00114.482243.901.2824.38149.132343.500.9029.86139.762437.900.5422.13162.60
2 估算模型
2.1 广义回归神经网络
σ对GRNN模型预测性能影响较大，σ取值过大可能使神经网络拟合偏差过大，而取值过小使得神经网络过学习.基础的广义回归神经网络隐含层一般采用同一个σ，样本数据的概率分布对GRNN模型为未知数，所以不能直接从样本中求得最佳光滑因子(σbest)，需使用观测样本对模型进行训练，搜寻求解σbest.搜寻过程一般使用轮询遍历交叉检验法和黄金分割快速搜寻法.
GRNN优秀的泛化能力使其成为一种常用的模型工具，已广泛应用于各领域的预测与估算.本文使用广义回归神经网络技术建立层合板低速冲击能量与冲击后凹痕深度、损伤面积以及低速冲击后剩余压缩强度的映射关系.
2.2 估算模型结构
复合材料层合板低速冲击后剩余压缩强度大小是众多影响因素综合作用的结果.考虑到样本数据信息量的局限性,选取主要影响因素如低速冲击的能量、冲击后凹痕深度与损伤面积作为输入参数，而将冲击后剩余压缩强度作为目标输出.研究表明[18]，对于非线性映射问题，含有4层结构的广义回归神经网络可以实现无限逼近，并且具有柔性的网络结构和较好的容错性，适于解决影响因素间存在高度非线性关系的实际问题.因此，冲击后剩余压缩强度估算模型以广义回归神经网络作为模型框架，其输入层、输出层、隐含层以及求和层等的神经元节点数比较容易确定.根据估算模型设定的输入和输出维数分别为3和1，可以确定输入和输出层神经元节点数为3和1，而隐含层的神经元节点数等于训练样本数，求和层神经元节点数一般为2.
2.3 模型的算法
广义回归神经网络是径向基神经网络的改型，其理论基础是非线性核回归分析.对于假定的随机变量x与y，概率密度函数为f(x,y)，已知x的观察值为X,则y相对于x的回归条件均值为
(1)
当f(x,y)未知时，y相对于x的回归分析实际上是计算具有最大概率值的y.通过对观测样本执行非参数估计，推导出观测样本中自变量和因变量的概率密度函数：
(2)
反归一化处理公式为
若观测样本的第k个样本为xk,隐含层的各个神经元对应一个期望输出yk.加和层的第1个神经元输出S1为隐含层输出yk乘以权值后的和，加和层的第2个神经元的输出S2为隐含层输出的直接求和，输出层的输出为S1/S2.与基于BP神经网络的估算模型相比，广义回归神经网络训练过程通常不需要迭代，且相关参数的确定较为简单，因此神经网络训练过程实际是确定光滑因子的过程.
人工神经网络起源于对神经细胞行为模式在信息处理方面的研究，是智能学科的重要组成部分，为解决复杂问题和智能控制提供了有效的途径.虽然误差反传(Back Propagation,BP)神经网络算法回避了输入相关性问题，但是由于其具有的拓扑结构，估算模型容易陷入局部极小值的困境.另外，隐含层神经元个数依赖于经验公式，因此BP算法对经验与技巧的要求较高.径向基函数(Radical Basis Function,RBF)神经网络解决了局部极小值的问题，但其学习方法采用风险最小化原则，存在一定的经验成分，因此可能出现过学习问题.广义回归神经网络(General Regression Neural Network,GRNN)来源于RBF神经网络的变形，一般由输入层、隐含层、加和层及输出层等4层神经元构成，具有结构简单、收敛速度快、函数逼近能力强等特点.
冲击后压缩强度与冲击能量、凹痕深度以及损伤面积与的关系如图1所示.可以看出，冲击后压缩强度与冲击能量和凹痕深度具有一定的线性关系，而与损伤面积之间不具备线性关系.
(1)轮询遍历交叉检验法：① 假设一个搜索区间，起始搜索点σ0取区间左边界值；② 将观测数据分为N组(N>2)，取出1组作为预测对比数据,其余组作为模型训练数据；③ 用σ构建的广义回归神经模型训练后求方均误差④ 循环选取数据组，使所有数据组都做过预测对比，分别计算均方差；⑤ 取搜索推进步长Δ，使得σ1=σ0+Δ，用σ1重新构建模型，重复步骤3，4，依次循环σi=σ0+iΔ直至达到搜索区间右边界；⑥比较各均方差，最小均方差对应的σ即为σbest.
(2)黄金分割快速搜索法.假设σ1,σ2在σbest的搜索区间(a,b)内，σ1,σ2与a,b两点的距离相同且σ1<σ2，σ1,σ2计算方法为
(3)
式中：约为 1.618.比较Eσ1与Eσ2的大小，若Eσ1<Eσ2，则下个搜索区间设定为(a,σ1)；若Eσ1>Eσ2，则下个搜索区间为(σ2,b).经过多次搜索计算，Eσ达到所需精度时，停止迭代求解，此时Eσ对应的σ即为σbest.
3 模型验证
3.1 数据规范化处理
3.观察组和对照组患者病情复发情况：2组CP患者定期随访1年，观察组复发2例，复发率1.4%；对照组复发15例，复发率10.9%。观察组CP患者病情复发率显著低于对照组，差异有统计学意义(P<0.05)。
大型仪器设备开放共享需要解决的问题 如何增进师生共享意识，提高协作共用？如何协调科研资源，为教学、科研与社会多方服务？如何保证共享实验的规范管理与操作安全？这些都是大型仪器设备开放共享需要解决的问题。
图1 冲击后压缩强度与冲击能量、凹痕深度以及损伤面积与的相关性
Fig.1 Correlation between compressive strength after impact and energy,indentation depth,damage area
(4)
基于GRNN的低速冲击后压缩强度估算模型在使用通常单位的情况下，凹痕深度、冲击能量、损伤面积以及压缩强度数据的数量级分别为100、101、101以及102，量级差较大.为避免该情况对估算模型准确度的影响，将所有观测样本数据归一化转换为[-1,1]之间的数值.对GRNN模型得出的估算数值进行反向处理，将得出的数值反归一化转换至原值域.归一化与反归一化一般使用MATLAB中的PREMNMX、POSTMNMX、TRAMNMX等内置函数进行处理.归一化处理公式为
复合材料层合板低速冲击试验及冲击后压缩强度试验采用的材料体系为CCF300/5428.试件的典型铺层为40层，铺层角顺序为[45/0/-45/90/0]4S，其中45，90以及0表示铺层方向分别为45°，90°以及0°，4为重复次数，S为对称之义.单铺层的名义厚度为 0.2 mm.单层材料的力学性能见表1.E1为纵向弹性模量，E1=148.8 GPa；E2为横向弹性模量，E2=10.02 GPa；G12为剪切模量；υ12为泊松比，υ12=0.31；Xc为纵向压缩强度，Xc=1.44 GPa.
(5)
式中：p为输入变量的维度数；n为观测样本的容量数；σ为光滑因子(高斯函数宽度系数).
(6)
(7)
式中：i=1,2,3,分别为冲击能力、凹痕深度以及损伤面积；k为该维度中元素的坐标；j=1，为输出矩阵的维度.
3.2 参数优化
基于GRNN压缩强度估算模型建模使用的输入、输出参数以及所涉及调优参数较少，只需调整σ且为一维较小的搜索空间，轮询遍历交叉检验法与黄金分割快速搜索法计算量差别不大.因此，下文以均方根误差为目标函数，运用轮询遍历交叉检验法寻找最优光滑因子.
设定高斯函数光滑因子的寻优空间为σ∈[0.01,1]，步长Δ=0.000 5.从观测样本预处理后的数据中用函数随机抽取21组作为GRNN估算模型的训练调参数据，用于搜寻σbest.其余的3组作为验证数据，再将21组训练数据用函数随机分为3分组，每分组含7个数据点，分组依次编号.根据轮询遍历交叉检验法使用MATLAB环境进行编程计算，经过学习训练，得到模型训练过程中目标函数Eσ的最小值仅为 22.330 2 MPa.可见，GRNN估算模型经训练后Eσ较小且训练收敛性很好.可以认为，GRNN估算模型的Eσ已经能满足模型预测的误差要求，最小Eσ对应的σbest为0.415，用此值对基于GRNN的冲击后压缩强度估算模型进行重构，可用于下一步估算模型验证.
3.3 估算结果
GRNN模型重构后，将随机抽取的另外3组观测样本数据输入已完成训练的冲击后压缩强度估算模型进行模型仿真预测，得出估算数值后应用式(6)和(7)计算出模型估算的冲击后压缩强度，见表2.可以看出，模型仿真的3组验证数据的冲击后压缩强度估算值与试验值比较接近，绝对误差区间为[-6.395 7,7.146 8] MPa，相对误差区间为[-4.37%,5.11%].
2） 运行费用估算：氧化药剂费用100元/t，生化药剂费用20元/t，电费30元/t，总计费用150 元 /t。
表2 GRNN模型估算的冲击后压缩强度
Tab.2 Compression strength estimated by GRNN model after impact
E/JD/mmS/cm2pc/MPa试验值估算值绝对误差/MPa相对误差/%43.500.9029.86139.76146.917.155.1155.402.2019.77129.91127.45-2.46-1.8943.301.3419.33146.45140.05-6.40-4.37
由表2可知，使用基于GRNN神经网络的冲击后压缩强度估算模型得到的估算值与试验值能较好地吻合，但个别观测样本的估算数据与试验数据偏差超过5%.如在遭受冲击能量为 43.50 J的低速冲击后，测得凹痕深度为 0.90 mm、损伤面积为 29.86 cm2的情形下，测得压缩强度为 139.76 MPa，而模型估算得到的压缩强度为 146.91 MPa，误差为 7.1 MPa.误差产生的可能原因有：① 复合材料层合板低速冲击试验后，所测损伤面积与实际损伤面积存在一定的误差；② 冲击后产生的表面凹痕有明显回弹效应，测量间隔时间在一定程度上影响凹痕深度精度；③ 基于GRNN估算模型需一定数量的观测数据训练，样本容量较小可能导致精度较低.
4 结语
复合材料层合板遭受低速冲击后，其力学性能的下降程度与损伤严重程度密切相关，而凹痕深度和损伤面积是损伤程度的主要表征参数.本文应用广义回归神经网络，以冲击能量、凹痕深度和损伤面积作为输入参数建立了复合材料层合板受到低速冲击后压缩强度估算模型.使用部分试验样本数据对广义回归神经网络压缩强度估算模型训练后，选取另一部分样本数据进行仿真验证.结果表明：基于广义回归神经网络的压缩强度估算模型计算效率较高且估算结果较为准确，该模型为复合材料层合板低速冲击后压缩强度估算提供了一种有效的新方法.
参考文献：
[1] DANG T D,HALLETT S R.A numerical study on impact and compression after impact behaviour of va-riable angle tow laminates[J].CompositeStructures,2013,96(4):194-206.
[2] SUEMASU H.Analytical approaches to compression after impact (CAI)behavior of carbon fiber-reinforced composite materials[J].AdvancedCompositeMaterials,2014,25(1):1-18.
[3] SUEMASU H,ICHIKI M.Analytical study on low compressive strength of composite laminates with impact damage[J].CompositeStructures,2013,104(5):169-175.
[4] SLATTERY P G,MCCARTHY C T,HIGGINS R M O.Assessment of residual strength of repaired so-lid laminate composite materials through mechanical testing[J].CompositeStructures,2016,147:122-130.
[5] JEFFERSON A J,ARUMUGAM V,SARAVANAKUMAR K,etal.Compression after impact strength of repaired GFRP composite laminates under repeated impact loading[J].CompositeStructures,2015,133:911-920.
[6] KULKARNI M D,GOEL R,NAIK N K.Effect of back pressure on impact and compression-after-impact characteristics of composites[J].CompositeStructures,2011,93(2):944-951.
[7] NASH N H,YOUNG T M,STANLEY W F.An investigation of the damage tolerance of carbon/benzoxazine composites with a thermoplastic tough-ening interlayer[J].CompositeStructures,2016,147:25-32.
[8] HABIB F A.A new method for evaluating the resi-dual compression strength of composites after impact[J].CompositeStructures,2001,53(3):309-316.
[9] BORRELLI R,FRANCHITTI S,CAPRIO F D,etal.A repair criterion for impacted composite structures based on the prediction of the residual compressive strength[J].ProcediaEngineering,2014,88:117-124.
[10] 沈真,杨胜春,陈普会.复合材料抗冲击性能和结构压缩设计许用值[J].航空学报,2007,28(3):561-566.
SHEN Zhen,YANG Shengchun,CHEN Puhui.Behaviors of composite materials to withstand impact and structural compressive design allowableness[J].ActaAeronauticaetAstronauticaSinica,2007,28(3):561-566.
[11] REMACHA M,SNCHEZ-SEZ S,LPEZ-ROMANO B,etal.A new device for determining the compression after impact strength in thin laminates[J].CompositeStructures,2015,127:99-107.
[12] ZHANG A Y,LU H B,ZHANG D X.Effects of voids on residual tensile strength after impact of hygrothermal conditioned CFRP laminates[J].CompositeStructures,2013,95:322-327.
[13] RIVALLANT S,BOUVET C,ABI ABDALLAH E A,etal.Experimental analysis of CFRP laminates subjected to compression after impact:The role of impact-induced cracks in failure[J].CompositeStructures,2014,111:147-157.
[14] TAN W,FALZON B G,CHIU L N S,etal.Predicting low velocity impact damage and compression-after-impact (CAI)behaviour of composite laminates[J].CompositesPartA:AppliedScienceandManufacturing,2015,71:212-226.
[15] MAIO L,MONACO E,RICCI F,etal.Simulation of low velocity impact on composite laminates with progressive failure analysis[J].CompositeStructures,2013,103:75-85.
[16] 贾建东,丁运亮,胡伯仁.复合材料层合板低速冲击后压缩破坏的数值模拟[J].机械科学与技术,2010,29(10):1320-1324.
JIA Jiandong,DING Yunliang,HU Boren.Numerical simulation of the compressive failure of composite laminates under low velocity impact[J].MechanicalScienceandTechnologyforAerospaceEngineering,2010,29(10):1320-1324.
[17] ROMANO F,DI CAPRIO F D,MERCURIO U.Compression after impact analysis of composite panels and equivalent hole method[J].ProcediaEngineering,2016,167:182-189.
[18] SPECHT D F.A general regression neural network[J].IEEETransactionsonNeuralNetworks,1991,2(6):568-576.
AnEstimatingMethodofCompressiveStrengthofCompositeLaminatesAfterLow-VelocityImpact
SHENGMingjian1,2,CHENPuhui1,QIANYibin2
(1.College of Aerospace Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China;2.Commercial Aircraft Corporation of China,Ltd.,Shanghai 200126,China)
Abstract：The compressive strength after low-velocity impact is a key performance index of the composite laminates.By analyzing data of compressive strength experiment of CCF300/5428 composite laminates after low-velocity impact (LVI)under different energy levels,a prediction model based on generalized regression neural network technology of compressive strength after LVI is proposed.This model uses the impact energy,the dent depth and the damage area as input parameters and the Gauss function as the hidden layer excitation function.The model uses part of the experimental data to train for finding the optimal smoothness coefficient.Then,the model is reconstructed with the optimal smoothness coefficient,and is simulated with part of the experimental data.Results show that the model has good generalization ability of experimental data and is feasible to estimate the compressive strength after low-velocity impact.
Keywords：neural network;composite materials;low-velocity impact;compressive strength;estimation
中图分类号：V 258
文献标志码：A
收稿日期：2017-12-26
基金项目：国家自然科学基金(11572152)，江苏高校优势学科建设工程资助项目
作者简介：盛鸣剑(1981-)，男，上海市人，博士生，主要研究方向为复合材料结构设计.
通信作者：陈普会，男，教授，博士生导师，电话(Tel.)：025-84896256;E-mail:phchen@nuaa.edu.cn.
文章编号：1006-2467(2019)10-1182-05
DOI:10.16183/j.cnki.jsjtu.2019.10.006
(本文编辑：陈晓燕)
标签：模型论文; 强度论文; 神经网络论文; 复合材料论文; 损伤论文; 《上海交通大学学报》2019年第10期论文; 国家自然科学基金(11572152)江苏高校优势学科建设工程资助项目论文; 南京航空航天大学航空宇航学院论文; 中国商用飞机有限责任公司论文;