结构准静态解耦载荷识别方法与结构随机载荷谱
载荷谱是进行结构可靠性设计和可靠性试验验证的基础。建立载荷谱的前提是获取结构在实际运用条件下的完整载荷历程。工程结构运用载荷历程的获取，一般是通过测量载荷在结构上产生的响应进行反求，称之为载荷识别问题[1]。这方面的研究工作较早始于航空领域 [2]，后来逐步扩展至航天、车辆、船舶、海洋平台等研究领域[3-9]。
载荷识别问题属于结构动力学中的反问题，这类问题解决起来通常要比结构动力学中的正问题更为复杂[1,9]。载荷识别方法广义上分为频域法和时域法[10]。经过四十余年的发展，已建立起直接求逆方法、正则化方法、小波分析法、遗传算法、逆虚拟激励法和神经网络算法等多种方法[1,3,11-14]。这些方法总体上还属于探索性的，成功应用于获取工程结构运用载荷历程的案例还有限[1,9-10,15]。
载荷识别问题解决起来更为困难的主要原因在于必须对载荷响应传递矩阵进行求逆运算。识别误差过大常常导致载荷识别失败，其根源在于载荷响应传递矩阵往往具有病态特性，求逆过程中会导致严重的误差放大[16-18]。
本文以转向架构架这一轨道车辆结构中可靠性问题最为突出的典型结构为对象，研究构架在实际运用条件下的载荷识别问题。结合构架实际运用情况下的动态特性分析，针对载荷响应传递矩阵的病态特性控制和矩阵求逆过程中的误差控制问题，探索针对性的解决方案。在此基础上，进一步研究确立构架基本载荷系。结合构架运用载荷历程的获取和构架基本载荷系的确立，可为建立构架随机载荷谱创造充分条件。
例10：校园歌舞表演片段，可以在后期剪辑时适当添加多种多样的转场效果，还可以利用“画中画”滤镜，设置移动路径，从无到有、大小渐变、旋转翻转、倒影阴影、定格动画等，产生较强的艺术效果，给观众以新颖、活泼、明快、变幻无穷的感觉，使表演气氛更加热烈浓厚。
1 载荷识别方法的基本原理
对于任一结构的载荷识别问题，载荷与响应之间的传递关系均可以表示为
ε(ω)=H(ω)F(ω)
(1)
式中:ε(ω)为响应向量;F(ω)为载荷向量;H(ω) 为载荷响应传递矩阵;ω为向量和矩阵元素的载荷频率。
在实际测试响应时，总会产生测量误差，由上式求解后，产生的载荷误差可以通过下式估计[16-19]
(2)
式中:‖F‖为载荷向量范数;‖δF‖为载荷误差向量范数;‖ε‖为响应向量范数;‖δε‖为响应误差向量范数;cond(H)为传递矩阵的条件数，可以表示为
近年来，网络技术发展也非常的迅速，而在智能电网中的继电保护装置，其实本质上也是一个计算机装置，作用就是对系统进行全面的监控工作，对整个网络的运行进行指导，让智能电网的系统运行能够高效和安全。同时，继电保护装置在对电力信息的数据和故障信息进行获取的时候，会及时反馈信息到系统的网络控制中心。另外，随着我国智能电网的发展越来越迅速，在自动化程度方面，电力系统也做得越来越好，这样就能够保障了电力系统的安全和高效。
(3)
其中，λmax和λmin分别为传递矩阵的最大特征值和最小特征值。显然传递矩阵的条件数不小于1。
屠呦呦从中药青蒿中提取的倍半萜内酯药物，说明了是通过干扰细胞表膜--线粒体，使疟原虫氨基酸饥饿，从而抑制蛋白合成，大量胞浆快速死亡而治疗疟疾的机理，被国际公认的获得了诺贝尔医学奖[2] ；陈竺发现白血病细胞分化的RIG基因，揭示了三氧化二砷（中药砒霜）通过诱导癌细胞程序性死亡，靶向性作用于早幼粒细胞白血病的关键致病基因编码蛋白质，使中药砒霜治疗白血病的原理得到了科学的证实。这种应用现代科学研究方法取得成果的范例，开启了中草药治病原理研究与临床应用，步入现代医学科学殿堂的大门。
综合式(2)和式(3)可知，通过式(1)求解载荷时，产生的载荷相对误差一般会得到放大。实际上，大量的载荷识别案例表明，当传递矩阵条件数很大时，载荷识别会产生无法接受的误差。载荷识别过程中，通常情况下传递矩阵的条件数难于时时把握和控制，这使得载荷识别误差的控制也难以有效，成为导致载荷识别失败的主要原因[1,9-10,15]。如何在载荷识别过程中，可靠控制传递矩阵的条件数不至于过大是保证载荷识别取得成功的关键因素。
对于转向架构架而言，由于置于一系悬挂之上，构架动载荷与动应变响应主要集中在10 Hz以下频段[20-22]，而构架弹性模态的第一阶固有频率通常在30 Hz左右[23-25]；在载荷识别的范畴中，构架动载荷识别基本上属于低频范围。目前大部分的动载荷识别方法基本上是针对中、高频动载荷的，却很难准确识别出低频范围的动载荷[9]。鉴于构架动载荷的低频特征和动载荷识别方法的发展现状，本研究尝试采用准静态标定方法来确定构架传递矩阵，并以构架动应变响应作为测试响应，此时式(1)可以表示为
动物实验表明，CTLA-4基因敲除的小鼠会迅速发展为淋巴细胞浸润性心肌炎，进一步导致严重的多器官衰竭[44-45]。小鼠心肌炎模型中CTLA-4(-/-)细胞毒性T淋巴细胞比CTLA-4(+/+)更易使小鼠罹患心肌炎[46]。因此，CTLA-4单抗引起心脏毒性的机制可能是其降低了心脏中T细胞特异性自身抗原活性的阈值。心脏毒性已经与其他免疫疗法如自身抗MAGE-A3 T淋巴细胞与心肌细胞的免疫反应导致细胞死亡相关联。
ε=KP
(4)
式中：ε为构架动应变响应向量；P为构架动载荷向量；K为构架载荷应变传递矩阵，在准静态标定方式下，K为常数矩阵。
为实现对传递矩阵K的条件数的可靠把控，有效抑制载荷识别误差，本研究采用解耦方法标定构架传递矩阵，通过最大限度地消除应变响应中载荷之间的相互影响，使传递矩阵K实现准对角化。当传递矩阵中非对角元素远小于对角元素时，对角线上的元素可以作为传递矩阵K的特征值的可靠估计，进而通过式(3)能够可靠估计传递矩阵K的条件数。因此，当采用解耦方法标定构架传递矩阵时，可以方便地确定传递矩阵的条件数，通过标定方案的适当调整可以实现对于条件数的有效把控，确保传递矩阵具备良好特性，进而有效控制载荷识别误差。此外，传递矩阵求逆过程是载荷识别产生误差的主要环节，当传递矩阵高度实现对角化后，矩阵求逆的大量运算可以省略，这也十分有利于减小载荷识别误差。
综上所述，消除载荷系之间的响应耦合，实现载荷系解耦的原则应包括两点：①载荷响应传递矩阵实现准对角化，即非对角元素远远小于对角元素；②载荷响应传递矩阵的条件数应尽量小。
2 构架载荷传递矩阵的解耦标定
以提速客车209P 转向架构架为研究对象。该构架承受的载荷有18个，分别为轮轨激扰产生的各轴端垂向载荷、横向载荷和纵向载荷(12个)，牵引产生的纵向载荷(2个)，制动时产生的制动吊座垂向载荷(4个)。各载荷作用位置及作用方式见图1。
图1 构架载荷示意图
2.1 构架基本载荷系
在实际运用过程中，转向架构架承受的载荷与构架的基本运动方式和运用方式相适应，以若干组相对独立的方式成组作用于构架，每一载荷组合均以自平衡方式呈现，这些载荷组合称之为构架基本载荷系。构架的基本运动方式主要包括浮沉、侧滚、扭转、横移、菱形等，基本运用方式主要包括牵引、制动等。每种基本运动和运用方式所对应的构架基本载荷系由UIC 515-4—2003[26]和BS EN 13749—2011[27]等国际规范详细描述，分别为浮沉载荷系Fv、侧滚载荷系Fr、扭转载荷系Fw、横向载荷系Ft、菱形载荷系Fo、纵向牵引载荷系Fl和制动载荷系Fb，见图2。
图2 构架各基本载荷系示意
在北京交通大学结构强度实验室多年开展构架载荷测试研究[28-29]的基础上，依据构架的运动特征，构架上还有可能存在现行国际规范中尚不包含的蛇行载荷系Fa，构架蛇行载荷系示意见图3。图3中，py1～py2为蛇行载荷，pRx1～pRx4为蛇行载荷反作用力。本研究将通过测试蛇行载荷系和对比分析由该载荷系产生的损伤情况，来判断是否有必要将该载荷系增加到构架基本载荷系中。
图3 构架蛇行载荷系示意图
2.2 构架基本载荷系的解耦标定
构架载荷以具有自平衡特征的基本载荷系方式成组出现。相应地构架的载荷解耦传递矩阵标定也是以基本载荷系为对象开展的。
在具体标定过程中，将构架基本载荷系分成2组：一组为浮沉载荷系、侧滚载荷系、扭转载荷系、横向载荷系、菱形载荷系和制动载荷系，以构架应变响应作为测试响应，通过构架标定载荷应变传递系数；另一组为纵向牵引载荷系和蛇行载荷系，以牵引拉杆应变响应作为测试响应，通过牵引拉杆标定载荷应变传递系数。
构架基本载荷系的解耦标定是一项十分细腻的工作。首先要对各基本载荷系在构架上产生的应变分布特性进行有限元分析，确定与每一基本载荷系对应的强响应区和弱响应区，并分析其应变分布具有的对称性特征。在确定某一基本载荷系的解耦测量方案时，测点区域选取方面一般优先考虑该基本载荷系对应的强响应区，最好同时又属于其他基本载荷系对应的弱响应区；测点应变组合选取方面则主要利用各基本载荷系产生的应变分布的对称性特征，必要时可进一步利用同一测点相互垂直方向上的应变通常符号相反的特点。构架各基本载荷系之间具有相对独立性，各自产生的应变分布情况和对称性特征差异显著，实践表明按上述方式确定构架基本载荷系的解耦测量方案是可行的。要获得构架基本载荷系之间比较理想的解耦效果，还必须对选定测点区域的测点应变组合进行优化，这需要对测点应变组合进行大量、细致的比选。
在解耦标定过程中，为了有效控制构架基本载荷系传递矩阵可能出现的病态状况，还需要对各基本载荷系对于各自测点应变组合的传递系数进行比较，实时估计传递矩阵的条件数，必要时调整解耦标定方案。
下面介绍2个典型的构架基本载荷系的解耦标定情况，分别是构架浮沉载荷系和菱形载荷系，其中浮沉载荷系的解耦实现相对容易，而菱形载荷系的解耦实现则相对困难。构架浮沉载荷系和菱形载荷系的解耦标定见表1、表2。
表1 构架浮沉载荷系的解耦标定
浮沉载荷系FvFtFbFwFrFo应变传递系数/(10-6·kN-1)12.40.040.060.150.640.27响应比/%1000.40.51.25.22.2
表2 构架菱形载荷系的解耦标定
菱形载荷系FoFvFtFbFwFr应变传递系数/(10-6·kN-1)32.40.630.490.570.400.25响应比/%1001.91.51.71.20.7
由表1和表2可见：对于选定的构架浮沉载荷系优化测点应变组合，与构架浮沉载荷系自身影响相比，其他基本载荷系的影响不超过6%；对于选定的构架菱形载荷系优化测点应变组合，与构架菱形载荷系自身影响相比，其他基本载荷系的影响不超过2%。总体来看，解耦标定后的构架基本载荷系传递矩阵的非对角元素与对角元素相比，在量值上小两个数量级，比较好地实现了构架基本载荷系传递矩阵解耦。
解耦标定后，构架基本载荷系相对于各自测点应变组合的传递系数，居于构架基本载荷系传递矩阵的对角线上，见表3。由于该传递矩阵已实现高度解耦，这些对角线元素可以作为构架基本载荷系传递矩阵特征值的可靠估计。对于采用构架应变作为测试响应的基本载荷系，最大特征值对应菱形载荷系，最小特征值对应制动载荷系，作为相应传递矩阵条件数的良好估值，两者之比仅为3.6，显示出基于构架测试的基本载荷系传递矩阵状态优良。对于采用牵引拉杆应变作为测试响应的基本载荷系，两者对应的特征值相同，条件数估值为1，显示出基于拉杆测试的基本载荷系传递矩阵处于最佳状态。
表3 构架基本载荷系与各自测点应变组合的传递系数 10-6·kN-1
基本载荷系基本载荷系-应变组合传递系数浮沉载荷系12.4侧滚载荷系11.3扭转载荷系11.1横向载荷系18.1制动载荷系9.4菱形载荷系33.4纵向牵引载荷系2.4蛇行载荷系2.4
综合上述，本研究比较好地实现了构架基本载荷系传递矩阵解耦，获得的传递矩阵状态优良，能够有效控制传递矩阵求逆产生的误差，为开展构架载荷识别工作奠定了坚实的基础。
3 试验验证
开展构架运用载荷识别的目标是建立构架载荷谱，进而服务于构架定量化可靠性设计、提升和试验验证。要实现这一目标，识别出的构架运用载荷必须能够正确再现构架的运用损伤，这是检验构架载荷识别是否成功的根本标准。构架的运用损伤直接取决于构架在运用过程中所承受的动应力历程，可由对应运用里程的等效动应力来度量，具体计算方法详见文献[30-31]。因此，识别出的构架运用载荷能否正确再现相同运用区段对应的动应力，特别是整个运用过程对应的等效动应力，就成为检验构架载荷识别是否成功的直接标准。要进行上述对比，需要开展两项工作：①在实验室内标定构架基本载荷系与构架关键部位应力(或应变)之间的传递关系;②同时测试构架载荷和这些关键部位的动应力。
近几十年来，新的教学方法不断涌现，其中PBL教学法的发展引人注目，有大量研究对PBL教学法进行了探讨[1～5]。很多研究发现，PBL教学法在培养学生自主学习、解决问题和创新能力等方面取得了良好效果。但也有研究发现，PBL教学法存在一些问题[4，5]，比如许多课程直接套用经典PBL教学法困难重重，或在应用过程中因受学生、教师或教学条件等方面的限制导致教学效果欠佳等。PBL教学法在临床和基础医学的很多课程中得到了广泛应用，但在大学生心理健康课程中的研究报道很少，PBL教学法是否适用于大学生心理健康课程，经典PBL教学法是否需要调整以便取得更好效果等问题有待深入探讨。
3.1 构架载荷应变传递关系
在构架上选取典型测点，①覆盖构架疲劳控制部位，②充分体现与构架所有基本载荷系的关联性。本研究共选取了56个典型测点，其中D1～D24测点为构架侧梁区域，D25～D44测点为构架横侧梁连接区域, D45～D48测点为构架牵引座区, D49～D56测点为构架制动座区，具体见图4。
从表5的数据可知，2011-2015年台湾对大陆农产品贸易额增加了约6.7～10亿美元，而其生产效果则增加了约19.61～26.96亿美元，约为贸易增加值的2～3倍。由此可见，台湾对大陆贸易附加价值GDP效果逐年增加，最高的一年增加1亿美元。尽管2015年的生产效果下降，但附加价值GDP效果仍是增加的，也就是说两岸农产品贸易有助于促进台湾提供高附加价值的产品。
图4 构架应变典型测点示意
采用类似于构架运用状态的约束方式，按照与构架各基本载荷系对应的加载方式，分别对构架分级施加载荷，在确认构架约束状况和载荷接触区域稳定后，测试构架典型测点的应变，计算出构架基本载荷系与典型测点应变之间的传递系数，具体见表4。
表4 构架基本载荷系与典型测点应变之间的传递系数 10-6·kN-1
测点FoFwFvFrFtFbFlFa测点FoFwFvFrFtFbFlFaD12.1-2.8-1.7-1.6-3.22.80.60.5D29-9.013.63.82.71.62.3-0.4-0.7D2-1.21.9-1.61.53.21.1-0.8-0.8D307.9-13.33.8-2.7-1.62.23.20.5D3-1.22.0-1.7-1.6-3.3-1.20.8-0.5D318.9-13.33.82.71.8-2.3-3.80.6D42.1-2.8-1.61.63.3-3.0-0.60.5D32-8.213.23.6-2.6-1.8-2.00.2-0.7D50.10.3-0.50.42.31.8-5.44.7D33-0.7-4.6-3.7-7.83.10.60.7-0.4D60.10.3-0.5-0.4-2.5-2.16.04.7D340.75.5-4.97.8-3.10.80.50.6D7-3.95.34.14.31.5-1.41.5-3.2D350.77.3-5.2-7.83.1-0.8-0.50.6D83.3-5.94.7-4.8-1.5-1.5-0.33.9D36-0.7-3.9-3.27.8-3.1-0.51.0-0.3D93.1-6.45.05.21.61.60.14.1D376.1-10.69.88.03.010.63.03.3D10-3.36.95.4-5.6-1.61.81.3-3.2D38-6.210.19.3-7.7-3.010.35.6-2.7D110.43.94.24.1-1.80.31.52.0D39-6.710.09.37.62.5-10.2-4.9-3.4D12-0.4-3.54.3-4.01.80.9-0.7-2.0D406.1-10.79.9-8.1-2.5-10.8-2.43.0D13-0.4-3.84.64.3-2.9-0.90.9-2.0D410.75.1-4.7-4.13.4-6.1-1.0-2.7D140.44.85.2-5.02.9-0.3-2.02.0D42-0.7-5.1-4.74.2-3.4-6.22.2-3.2D15-0.8-2.96.35.40.45.01.4-0.4D43-0.7-5.0-4.7-4.23.76.0-2.1-3.2D160.93.67.5-6.5-0.46.11.5-1.0D440.75.1-4.74.1-3.76.11.4-2.7D170.53.57.36.30.3-5.9-1.4-1.0D45-0.9-4.53.2-8.14.00.72.9-2.1D18-0.9-3.57.9-6.9-0.3-6.1-1.7-0.4D46-0.9-4.53.28.1-4.2-0.7-2.7-2.1D19-3.02.15.24.8-1.81.11.0-0.6D47-0.10.2-0.10.20.1-0.17.9-7.8D203.3-2.45.7-5.21.81.1-1.94.5D48-0.10.2-0.1-0.2-0.10.1-9.5-9.2D213.1-1.94.84.5-1.7-1.01.64.3D49-9.74.84.71.11.610.31.1-1.0D22-3.22.15.2-4.81.7-1.1-1.0-0.9D508.8-3.45.6-1.1-1.611.84.03.4D231.22.0-4.93.8-0.7-9.4-3.62.4D518.8-3.45.61.01.9-11.9-4.94.1D241.22.3-5.6-4.40.710.74.13.8D52-8.32.84.7-0.9-1.9-9.8-0.6-0.8D25-8.35.64.14.66.2-2.70.3-7.1D53-10.315.64.12.41.4-2.6-0.4-0.9D266.6-5.13.9-4.3-6.2-2.20.46.9D5411.2-15.74.2-2.5-1.4-2.82.80.6D275.2-5.13.94.3-6.32.2-0.46.8D559.7-15.74.22.51.52.5-3.00.6D28-6.95.64.1-4.66.32.70.3-7.1D56-11.815.64.2-2.5-1.52.70.2-1.0
3.2 线路试验简介
提速客车209P转向架构架于2015-10-13—2015-10-19，在大连至赤峰区段实际运用条件下进行测试。经运用部门确认，该线路路况较一般线路恶劣，具备典型性，适合于构架运用可靠性考核。测试共包含3个往返，累计运行里程5 826 km，累计运行时长近96 h。
构架测滚载荷系和菱形载荷系的测试情况见图5。构架典型测点的动应力测试情况见图6。限于数据量的庞大，为较清楚地显示工况特征，图中仅展示了时长600 s典型工况的测试结果。
（3）必须相信用户之间不会共谋。如果一个群体或个人控制了PoW系统中51%的算力，或者PoS系统中51%的投票权，整个区块链的防篡改根本就不成立了。
图5 构架基本载荷系典型工况时域图
图6 构架典型测点实测动应力时域图
3.3 构架载荷识别分析评价
前已说明，对于构架载荷识别的效果，应当以能否正确再现构架等效动应力来评价。为此，首先需要采用线路实测获得的构架各基本载荷系的时间历程，借助构架各基本载荷系与构架典型测点应变传递系数，通过叠加各基本载荷系产生的动应力响应，来预测构架典型测点动应力的时间历程。
构架代表性测点实测动应力和载荷预测动应力在典型工况下的时间历程见图7。由图7可见，预测动应力与实测动应力相比，大小相近，变化趋势基本一致。上述测点与一个完整测试往返相对应的实测动应力和预测动应力的频谱图见图8。从图8中可以看出：实测动应力频域响应主要集中在10 Hz以下；采用实测构架载荷预测的动应力与实测动应力的频谱变化趋势基本相同，总体吻合良好。
客观、科学衡量县域经济发展水平是一项具有系统性、动态性和相对性的复杂工作.构建经济综合评价体系需要依据科学性、全面性和可获取性原则，参考前人相关研究成果[15-16]并结合成都平原城市群经济发展的实际情况，从经济规模、经济效益、经济结构和政府调控能力4个方面，选取了17个具体评价指标，从而构建了县域经济综合评价指标体系(见表1).
图7 构架典型测点实测动应力与预测动应力的时域对比情况
图8 构架测点预测动应力与实测动应力频谱图
检验构架载荷识别成效的综合性标准，也是更具根本性的标准，在于实测构架载荷预测的等效动应力与实测等效动应力是否一致。本文中，各类等效动应力均对应于设计寿命里程(全寿命运用里程)，实测载荷预测的动应力和实测动应力均扩展至全寿命运用里程进行等效。对于209P转向架构架，全寿命运用里程取1 200万km。构架关键测点实测载荷预测等效动应力与实测等效动应力之间的比较情况见图9。由图9可见，对于构架上所有的关键测点，预测等效动应力与实测等效动应力的比值均在0.8～1.3之间，误差不超过30%。
再然后是夏小春约温简逛街，温简买了一条五百块钱的裙子，那是温简最贵的一条裙子，其实也是心疼的，但看着夏小春买了好些东西，她也有些浮躁了。夏小春说：“温简，这裙子好像是为你量身定做的，不如我买来送你。”温简急急地说：“不用，我自己买。”是有些故做的镇静和骄傲。
图9 构架关键测点载荷预测等效动应力与实测等效动应力之比
为更进一步评价构架实测载荷的效果，图10还给出了构架可靠性设计现行国际规范对应的用于可靠性设计评价的应力幅(等价于本文中的等效动应力)与构架实测等效动应力的比较情况。由图10可见，构架可靠性设计现行国际规范给出的应力幅与实测等效动应力的比值在1.3～4.0 之间，误差可高达400%。按本研究实测构架载荷预测的等效动应力，比国际现行规范给出的等效动应力精度提高一个数量级。
图10 构架现行国际规范等效动应力与实测等效动应力之比
根据文献[30]，运用寿命与等效动应力之间存在如下关系
(5)
式中：Dp为基于实测载荷预测的运用寿命；Dm为基于实测动应力预测的运用寿命；σp为实测载荷预测的等效动应力；σm为实测等效动应力；m为S-N曲线斜率，对于209P转向架焊接构架，m取值为3.5。
根据式(5)可知，本研究基于实测载荷预测的运用寿命与基于实测动应力预测的运用寿命的比值在0.46～2.5倍之间，预测精度基本满足定量预测构架运用寿命需要，可用于既有构架可靠性提升；而基于现行国际规范载荷预测的运用寿命与基于实测动应力预测的运用寿命的比值在2.5～128倍之间，表明构架可靠性设计现行国际规范还处于定性层面，不适合定量预测构架运用寿命。
于是除了沉潜于作画（仅1958年就创作了200多幅作品，可谓是一个人的大跃进）潘天寿再次有了辞职的念头，却未能如愿，他不知道，在毛泽东那个讲话后，即使有人真对他有意见，也不敢轻易动他，所以越请辞越升职。
上述比较显示：对于载荷响应主要集中于低频区域的转向架构架，采用本研究建立的结构准静态解耦载荷识别方法，获得的构架载荷识别结果可以比较好地预测构架运用等效动应力，也即构架运用损伤情况，已能够满足既有构架可靠性提升需要。从适应构架可靠性设计要求这一根本目标来看，本研究开展的构架载荷识别工作是相当成功的。
4 构架随机载荷谱
建立构架随机载荷谱需要两方面条件：①确立构架完备的基本载荷系；②获取构架实际运用条件下各基本载荷系具有代表性的载荷时间历程。
在构架基本载荷系的确立方面，在构架可靠性设计现行国际规范的基础上，补充了蛇行载荷系。构架各测点与蛇行载荷系对应的等效动应力同全部基本载荷系对应等效动应力的对比见图11。由图11可见，在构架牵引拉杆座区域，蛇行载荷系对等效动应力有显著影响；对于其他区域则影响不大；总的来看，增加蛇行载荷系是有意义的。
图11 蛇行载荷系对构架关键测点预测等效动应力的响应比
由图9可知，基于本研究确立的8个构架基本载荷系及其实测运用载荷时间历程，可以比较高的精度预测构架运用等效动应力；一方面表明这8个基本载荷系对于拖车构架具有完备性；另一方面也表明结合这些基本载荷系的运用载荷时间历程，可以直接建立构架随机载荷谱。
随机载荷谱是最接近于实际运用条件的载荷谱，应用时全部基本载荷系同步施加于构架之上，每个基本载荷系对应着各自的载荷时间历程。
所谓自主阅读，就是充分调动学生的主观能动性，让学生带着一定的阅读目的，运用一定的阅读方法，让学生通过探索、实践、质疑、创造等方法独立提出问题，分析问题和解决问题，从而实现阅读目标。
十一届全国人大常委会第二十六次会议听取国务院关于农田水利建设工作情况的报告，并针对报告向水利部、国家发改委、财政部、农业部等部门展开专题询问。面对常委会委员们的提问，国务院相关部门负责人坦诚回答，使得诸多疑问在询问中化解，共识在互动中形成，未来农田水利建设的思路得到进一步厘清。
当应用于构架可靠性设计或者构架可靠性试验验证时，随机载荷谱还不够方便，对试验设备的要求也往往过高，实际应用时一般还需要在此基础上编制出方便应用的载荷谱。
5 讨论
在应用本研究建立的结构准静态解耦载荷识别方法时，载荷传递矩阵的解耦标定过程至关重要，需要把控好每一个关键细节。由图8可知，在5 Hz以下的低频区域，构架关键测点动应力的载荷预测结果与实测结果还存在较为明显的偏差。具体来看，测点14的预测结果大于实测值，而测点31的预测结果则小于实测值。不仅两者的偏离趋势不同，而且在低频区域偏差更为显著，也与准静态识别方法的基本原理相抵触。经过细致地分析整个测试过程，认为问题的主要原因在于：在标定过程中，对于构架的约束处理还不够到位，特别是在侧滚载荷系和扭转载荷系的标定过程中，垂向约束刚度显著大于构架运用时的垂向在位刚度，这对于构架横侧梁连接区域的载荷响应传递关系可以产生较为明显的影响。
对于城市轨道车辆和高速动车组，构架动应力测试情况显示，构架的1～2个弹性模态可能被有效激发[21,29,32-33]。对于弹性模态能够被有效激发的构架，载荷识别时应在本研究建立的结构准静态解耦载荷识别方法的基础上，进一步考虑弹性模态的动态放大效应。就构架而言，激发出的弹性模态一般呈明显分离状态，可以在某一模态频率附近频段，按单自由度方式识别模态参数，进而计算并消除弹性模态的动态放大效应。
6 结论
本研究在分析结构载荷识别误差产生机理的基础上，针对低频响应结构，建立了结构准静态解耦载荷识别方法，通过确保载荷响应传递矩阵的良性状态和省略传递矩阵主要的求逆运算过程来控制载荷识别误差。采用提速客车转向架构架进行了测试验证，结果表明本研究建立的载荷识别方法可以有效控制长周期服役条件下的载荷识别误差，成功解决了现有结构载荷识别方法不适用于低频响应结构载荷识别的问题。
在构架可靠性设计现行国际规范的基础上，补充确立了蛇行载荷系，完善了拖车构架基本载荷系，进一步结合构架各基本载荷系的实测运用载荷时间历程，可建立构架随机载荷谱。依据实测载荷预测的构架等效动应力与实测等效动应力相比，误差不超过30%，运用寿命预测误差不超过2.5倍；比构架可靠性设计现行国际规范预测的等效动应力精度提高一个数量级，运用寿命预测精度提高2个数量级。运用上述体系建立的载荷谱，适用于对既有构架进行可靠性升级，并具备将现有国际规范由定性层次提升至定量层次的潜力。
参考文献：
[1] 张方, 秦远田. 工程结构动载荷识别方法[M]．北京：国防工业出版社,2011.
[2] BARTLETT F D, FLANNELLY W G. Model Verification of Force Determination for Measuring Vibratory Loads[J]. Journal of the American Helicopter Society, 1979, 24(2):10-18.
[3] 贾有, 杨智春. 一种飞机垂尾抖振载荷识别的新方法[J]. 航空学报, 2013,34(10):2333-2340.
JIA You, YANG Zhichun. A New Approach to Identify Buffet Loads for Aircraft Vertical Tail[J]. Acta Aeronautica et Astronautica Sinica, 2013,34(10):2333-2340.
[4] 张玉良, 杨飞, 岳洪浩, 等. 基于频域法的星箭连接分离装置的冲击载荷识别[J]. 振动与冲击, 2018,37(17):79-85.
ZHANG Yuliang, YANG Fei, YUE Honghao, et al. Impact Load Identification of Connection-separation Device between Satellite and Rocket with Frequency Domain Method Based on EEMD[J]. Journal of Vibration and Shock, 2018,37(17):79-85.
[5] 邓江华, 孙健颖, 李奥飞. 基于主成分分析的运行工况车轮轮心载荷识别方法研究[J]. 噪声与振动控制, 2019,39(4):145-149.
DENG Jianghua, SUN Jianying, LI Aofei. Research on Operational Six-Axis Wheel Force Identification Method Based on Principal Component Analysis[J]. Noise and Vibration Control, 2019, 39(4):145-149.
[6] 王明猛,朱涛,王小瑞,等.一种逆结构滤波法的轨道车辆轮轨力识别[J].振动工程学报,2019,32(4):602-608.
WANG Mingmeng, ZHU Tao, WANG Xiaorui, et al. An Inverse Structural Filter Method for Wheel-rail Contact Forces Identification of Railway Vehicles[J].Journal of Vibration Engineering, 2019,32(4):602-608.
[7] 杨志荣, 戴乐阳, 于洪亮, 等. 船舶中间轴承的载荷识别技术研究[J]. 噪声与振动控制, 2018,38(S1):462-464.
YANG Zhirong, DAI Leyang, YU Hongliang, et al. Research on Marine Intermediate Bearing Load Identification Technology[J]. Noise and Vibration Control, 2018, 38(S1):462-464.
[8] 张韶光, 范勇, 马汝建, 等. 海洋平台振动载荷识别研究进展[J]. 济南大学学报(自然科学版), 2004,18(4):332-336.
ZHANG Shaoguang, FAN Yong, MA Rujian, et al. Progress of Vibration Loading Recognition of Offshore Platforms[J]. Journal of Jinan University(Science and Technology), 2004,18(4):332-336.
[9] 杨智春，贾有．动载荷的识别方法 [J]．力学进展，2015，45(1)：29-54.
YANG Zhichun，JIA You．The Identification of Dynamic Loads [J]．Advances in Mechanics，2015，45(1)：29-54.
[10] 毛玉明, 林剑锋, 刘靖华, 等. 动载荷反演分析技术研究综述[J]. 动力学与控制学报, 2014, 12(2):97-104.
MAO Yuming, LIN Jianfeng, LIU Jinghua, et al. Recent Advances of Dynamic Force Estimation Technique[J]. Journal of Dynamics and Control, 2014, 12(2):97-104.
[11] 薛晓峰, 陈雪峰, 耿佳, 等. 基于Hermitian小波有限元的叶轮叶片载荷识别[J]. 振动与冲击, 2018,37(10):71-79.
XUE Xiaofeng, CHEN Xuefeng, GENG Jia, et al. Load Identification of Impeller Blades Based on the Hermitian Wavelet Finite Element Method[J]. Journal of Vibration and Shock, 2018,37(10):71-79.
[12] 常晓通, 闫云聚. 一种基于遗传算法的时域内稳态动载荷识别方法[J].应用力学学报, 2012,29(3):274-277.
CHANG Xiaotong, YAN Yunju. A Genetic Algorithm Based Method for Steady-state Dynamic Force Identification in Time Domain[J]. Chinese Journal of Applied Mechanics, 2012,29(3):274-277.
[13] 郑敏，任芳，杨兆建，等．自适应BP神经网络的转子系统载荷识别的研究 [J]．机械设计与制造，2016(6)：85-88，92.
ZHENG Min，REN Fang，YANG Zhaojian，et al．Study on the Load Identification of the Rotor System with Self-adaptive BP Neural Networks [J]．Machinery Design & Manufacture，2016(6)：85-88，92.
[14] 田燕, 邢士勇, 张志斌, 等. 基于径向基网络的变速箱载荷识别研究[J].振动、测试与诊断, 2004,24(2):100-102.
TIAN Yan, XING Shiyong, ZHANG Zhibin, et al. Load Identification of Gearbox Based on Radial Basis Function Neural Networks[J]. Journal of Vibration Measurement & Diagnosis, 2004,24(2):100-102.
[15] 周盼, 张权, 率志君, 等. 动载荷识别时域方法的研究现状与发展趋势[J]. 噪声与振动控制, 2014,34(1):6-11.
ZHOU Pan, ZHANG Quan, SHUAI Zhijun, et al. Review of Research and Development Status of Dynamic Load Identification in Time Domain[J]. Noise And Vibration Control, 2014,34(1):6-11.
[16] 秦领,刘西拉.关于结构刚度方程病态问题的疑难分析[J].力学季刊,2005,26(3):370-376.
QIN Ling,LIU Xila.Puzzle Analysis of Ill-conditioned Structure Stiffness Equation[J].Chinese Quarterly of Mechanics, 2005,26(3):370-376.
[17] 郭志松．基于矩阵条件数的载荷识别方法研究 [D]．哈尔滨：哈尔滨工程大学，2015.
[18] 杨承民.矩阵求逆条件数在误差估计中的最优性[J].华东师范大学学报(自然科学版),1986(3):35-38.
YANG Chengmin.Optimum of K(A)=‖A‖‖A-1‖ in Error Estimate [J].Journal of East China Normal University(Natural Science), 1986(3):35-38.
[19] GOLUB G H, LOAN C F V. Matrix Computations [M]. 2nd ed.Baltimore: Johns Hopkins University Press, 1989.
[20] 张立. 北京地铁2号线车辆转向架构架载荷测试与研究[D]. 北京:北京交通大学, 2017.
[21] 王斌杰,孙守光,王曦,等.地铁转向架构架运用载荷与疲劳损伤特征研究[J].铁道学报,2019,41(6):53-60.
WANG Binjie,SUN Shouguang,WANG Xi,et al.Research on Characteristics of Operation Loads and Fatigue Damage of Metro Train Bogie Frame[J]. Journal of the China Railway Society,2019,41(6):53-60.
[22] 王萌. 焊接转向架构架线路载荷的特征与应用研究[D]. 北京:北京交通大学, 2016.
[23] 杨萌,李强,彭永明,等.30 t轴重电力机车转向架构架的静强度与疲劳强度分析及模态分析[J].铁道机车车辆,2014,34(3):18-22.
YANG Meng,LI Qiang,PENG Yongming,et al. Strength and Modal Analysis of 30 t Electric Locomotive Bogie Frame[J]. Railway Locomotive & Car,2014,34(3):18-22.
[24] 吴丹,郭富强,商跃进,等.CW-200K转向架构架疲劳强度及模态分析研究[J].兰州交通大学学报,2011,30(3):51-54.
WU Dan,GUO Fuqiang,SHANG Yuejin,et al.Fatigue Strength and Modal Analysis of CW-200K Bogie Frame[J]. Journal of Lanzhou Jiaotong University,2011,30(3):51-54.
[25] 李凯. 出口伊朗地铁车辆转向架构架疲劳可靠性研究[D]. 北京:北京交通大学, 2013.
[26] International Union of Railways. Passenger Rolling Stock-Trailer Bogies Running Gear Bogie Frame Structure Strength Tests: UIC 515-4—2003[S]. Paris, France: International Union of Railways, 2003.
[27] European Commiftee for Standardication.Railway Applications—Wheelsets and Bogies—Method of Specifying the Structural Requirements of Bogie Frames: BS EN 13749—2011[S], London:British Standards Institution, 2011.
[28] 朱宁. 高速列车转向架构架结构损伤一致性载荷谱理论研究 [D]. 北京：北京交通大学, 2016.
[29] 邹骅. 城际动车组转向架构架载荷谱研究[D]. 北京：北京交通大学, 2016.
[30] DIETER R．焊接结构疲劳强度[M]．郑朝云，张式程,译．北京：机械工业出版社，1994.
[31] 王文静,王燕,孙守光,等. 高速列车转向架载荷谱长期跟踪试验研究[J]. 西南交通大学学报, 2015, 50(1): 84-89.
WANG Wenjing, WANG Yan, SUN Shouguang, et al. Long-term Load Spectrum Test of High-speed Train Bogie [J]. Journal of Southwest Jiaotong University, 2015, 50(1): 84-89.
[32] 张丽,任尊松,孙守光,等.构架弹性振动对疲劳寿命影响研究[J].铁道机车车辆,2015,35(2):115-119.
ZHANG Li,REN Zunsong,SUN Shouguang,et al. Research on the Influence of Railway Bogie Elastic Vibration to Fatigue Life[J]. Railway Locomotive & Car,2015,35(2):115-119.
[33] 金新灿,孙守光,陈光雄.列车通过道岔时转向架结构振动特性的研究[J].力学学报,2006,38(4):530-539.
JIN Xincan,SUN Shouguang,CHEN Guangxiong. Study on Characteristics of Railway Bogie’s Structural Vibrations while Train Passing Through a Switch[J]. Chinese Journal of Theoretical and Applied Mechanics, 2006,38(4):530-539.