美国核管会(NRC)2018年9月17日宣布批准为印第安角2号和3号机组换发运行许可证，允许这两台机组延寿20年，即可分别运行至2024年4月30日和2025年4月30日。
基于对照经验公式的构件本构离线模型更新混合试验方法的主要思想是：钢筋混凝土框架中柱构件因构件参数如轴压比、配筋率、截面尺寸等异同，导致构件的恢复力模型参数各不相同。在整体结构中，所有柱构件都可以利用混凝土规范或经验公式计算获得一组非精确恢复力模型参数，将通过这种方式计算获得的结果称为对照组。同时，对于混合试验中的物理子结构，通过参数识别的方式可以获取柱构件的真实恢复力模型参数。将物理子结构的参数识别结果同该构件的对照组计算结果进行正向对比，以正向对比误差作为修正基准，反向修正待更新数值子结构中柱构件对照组恢复力模型参数，从而完成数值子结构中构件参数不同柱构件的修正更新，以此来达到提高混合试验的精度的目的。基于对照经验公式的离线模型更新混合试验方法流程如图1所示。
图1基于对照经验公式的离线模型更新混合试验流程图Fig.1Flowchartofhybridsimulationforoff-linemodelupdatebasedoncomparativeempiricalformula
该方法与传统模型更新混合试验不同的地方是增加了对照经验公式修正模块，用以修正构件参数不同的数值子结构构件的恢复力模型参数。对照经验公式修正模块主要包含3个子模块：对照组建立模块、正向对比分析模块、反向修正更新模块。在本文中对照组的含义是作为建立对比和修正基准的一个恢复力模型参数经验公式计算结果集合。对照组的建立需要选择合适的恢复力模型参数经验公式，这将对混合试验的精度造成影响。正向对比分析子模块是由物理子结构试验的恢复力模型参数识别结果及该构件对照组恢复力模型参数组成。对于混合试验，物理子结构试验恢复力模型参数仅能通过参数识别获取一组样本，故无法进行线性或非线性的回归分析，因此本文选用正比例线性修正方法，如式(1)所示，计算得到参数修正系数。根据正向对比分析结果，对待更新数值子结构构件的对照组恢复力模型参数进行线性修正更新。通过上述修正模块的修正计算即可完成对构件参数不同的数值子结构构件更新。
式中：xi为通过经验公式计算的构件各个恢复力模型参数初始值；yi为物理子结构参数识别获得的恢复力模型参数真实值；ki为各恢复力模型参数的修正系数。
此外，以杀人偿命观为例，大众观念没有一个具体的衡量标准。如果杀人后需要一命抵一命，那么无论是故意杀人、故意伤害致人死亡、过失致人死亡、意外事件等都需要行为人承担责任，明显地违反了主客观相统一的原则。若是一贯地践行杀人偿命观，在具体认定犯罪时也不免出现主观归罪的倾向。
在线模型更新的混合试验采用拟动力试验的方式，试验流程较为复杂，且受本文修正更新算法及经验公式选择的不确定性影响，可能会使混合试验无法准确的反映结构的抗震性能。而离线模型更新混合试验对数据进行批量后处理，使经验公式及修正算法具有一定选择空间，且试验流程简单，便于操作。因此本文采用离线模型更新验证对照经验公式修正的可行性。
本文修正方法需对柱构件恢复力模型参数按照经验公式进行计算并建立对照组，因此选用IMK模型作为混凝土框架中待更新数值子结构柱构件的恢复力模型，IMK模型由IBARRA等[1213]提出，该模型将恢复力骨架曲线简化为三线性，即弹性段、强化段及软化段，它能够模拟混凝土梁柱的强度、刚度退化等现象，如图2所示。其骨架曲线形状由5个参数EIy，θcap,pl，θpc，My和Mc/My确定。HASELTON等[14]通过大量钢筋混凝土柱的试验数据建立了IMK恢复力模型中骨架曲线参数的经验公式，当钢筋混凝土构件为弯曲破坏时，计算公式如式(2)～(10)所示。
图2IMK模型骨架曲线Fig.2SkeletoncurveofIbarra-Medina-Krawinklermodel
式中：EIy为割线模量；θcap,pl为塑性转角；θpc为峰值点后软化段塑性转角；My为屈服弯矩值；Mc/My为峰值弯矩与屈服弯矩的比值；φy为屈服曲率；ky，A和B为构件基本信息计算出的参数；P/Agf′c，v为轴压比；b和d分别为柱横截面的宽和高；δ′=d′/d，d′为受压钢筋中心与受压区边缘之间的距离；n=Es/Ec，Es和Ec分别为钢筋和混凝土的弹性模量；asl为纵向钢筋滑移系数(考虑取1，不考虑取0)；f′c为混凝土圆柱体抗压强度；fy为受拉钢筋屈服强度，MPa；ρsh为柱塑性铰区横向钢筋面积比；ρ为纵向钢筋配筋率；ρ′为受压钢筋配筋率；ρv为抗剪钢筋的配筋率。
为验证基于对照经验公式的构件本构离线模型更新混合试验方法的可行性，本节以中柱配筋率与边柱配筋率不同的1榀2层2跨弹性钢筋混凝土框架为例，通过OpenSees有限元软件进行离线模型更新混合试验，钢筋混凝土框架结构模型如图3所示。
图3弹性RC框架结构计算模型Fig.3CalculationmodeloftheelasticRCframestructure
(a)框架结构原型；(b)柱截面尺寸及配筋图
框架建模采用二维杆系模型，各节点包含2个平动和1个转动共3个自由度。框架层高3.0m，跨度为6.0m，梁柱截面尺寸分别为400mm×400mm，200mm×500mm。框架梁柱构件采用IMK模型进行建模，第1刚度段参数设置为合理范围使构件始终处于弹性阶段。因需采用对照经验公式的更新修正方法，故对照组参数取为梁柱单元的弹性模量，对照经验公式取为框架柱换算弹性模量，计算公式如式(11)所示：
式中：As和Ac分别为钢筋和混凝土的截面面积。
表1钢筋和混凝土弹性模量Table1ElasticmodulusofconcreteandrebarMPa
混合试验取底层边柱为物理试验子结构，进行拟静力试验，对试验结果进行参数识别，并将公式(11)计算所得的弹性参数值作为对照组，对数值子结构中中柱刚度进行线性修正，柱弹性模量预测值如表2所示。
表2柱弹性模量预测值Table2PredictedcolumnelasticmodulusMPa
弹性钢筋混凝土框架地震时程分析采用以下4种工况：1)初始值；2)真实值；3)传统构件本构离线模型更新混合试验，简称为传统构件更新；4)基于对照经验公式的构件本构离线模型更新混合试验，称之为对照构件更新。框架时程分析结果如图4所示。
图4弹性框架结构时程响应Fig.4Timehistoryresponsesofelasticframestructure
(a)底部剪力时程局部放大图；(b)顶点位移时程局部放大图
表3不同工况下RMSE值Table3RMSEunderdifferentworkingconditions
式中：Dref，D分别为真实值和各工况仿真结果的向量，T为向量转置。
由各个工况底部剪力及顶点位移RMSE值可以明显的看出，基于对照经验公式的构件本构模型更新混合试验方法均方根误差小于传统模型更新混合试验方法，相较而言本文试验方法可以从构件本构层面显著提升混合试验精度。
不同于传统构件本构模型更新混合试验仅能对构件恢复力模型参数相同的构件进行更新，这种新的模型更新方法可以对不同恢复力模型参数的相似构件进行更新。对于复杂结构的模型更新混合试验，本文方法可以扩大更新的构件数量及范围，进而提升混合试验的精度。
为检验基于对照经验公式的构件本构模型更新混合试验方法的有效性及离线模型更新数据重复利用的可靠性，本节通过OpenSees有限元软件进行1榀5层3跨钢筋混凝土弹塑性框架(框架A)离线模型更新混合试验数值仿真，并重复利用物理子结构试验数据完成1榀3层2跨钢筋混凝土框架(框架B)离线模型更新数值仿真。
表4材料本构模型参数Table4Materialparametersofconstitutivemodel
注：fc，εo为混凝土的抗压强度及对应的应变；fu，εu为混凝土的残余强度及对应的应变；fy为混凝土的屈服强度；E为钢筋的弹性模量；b钢筋的应变硬化率；R0为控制从弹性段向塑性段过度的参数。
取框架A底层边柱作为物理子结构，物理子结构模型建立采用4.1节相同的构件及截面参数。对子结构构件进行足尺拟静力数值仿真试验，构件底端约束，顶端有限转动，柱顶施加恒定轴力P=247.6kN。以构件顶端位移控制加载，所加水平位移依次为5，10，15，20，25，30，40，45，55mm。
物理子结构拟静力试验数值仿真结果如图5(a)所示，其中真实值及初始值分别对应表4中材料的不同数值仿真结果。提取构件拟静力试验结果真实值的骨架曲线，利用Matlab编程的差分进化算法[15]对骨架曲线进行参数识别。拟静力试验构件未进入第3刚度段及倒塌阶段，故IMK模型参数识别仅识别弹性段和强化段。根据参数识别结果，对物理子结构IMK恢复力模型参数进行计算。利用OpenSees完成物理子结构的IMK模型建模，并施加以相同加载方式，IMK模型拟静力计算结果如图5(b)所示，可见使用参数识别结果的IMK恢复力模型可以较好地反应构件的滞回性能。
图5物理子结构拟静力结果及参数识别Fig.5Pseudo-staticresultsandparametersidentificationofphysicalsubstructure
(a)拟静力试验数值仿真结果；(b)IMK模型数值仿真结果
因框架结构中柱构件的轴压比各不相同，即恢复力模型参数不同，所以采用对照经验公式的构件本构离线模型更新混合试验对钢筋混凝土框架进行地震动作用下的时程分析数值仿真。
T·S·艾略特早期诗歌中的反英雄群像及认知价值…………………………………………………江群（2.106）
图6框架A结构时程响应Fig.6FrameAstructuretimehistoryresponse
(a)底部剪力时程曲线；(b)顶点位移时程曲线
表5框架A中IMK恢复力模型参数预测Table5IMKhystereticmodelparameterpredictionofframeA
对比几种工况下的RMSE值，如表6所示。可以明显看出，对照经验公式修正可以提升模型更新混合试验的精度。但由于受到修正算法及恢复力模型经验公式选择的限制，修正精度还有待提高。
表6框架A各工况RMSE值Table6FrameARMSEforeachworkingcondition
根据4.3节物理子结构对照组正向对比修正子模块，对框架B进行底层柱离线模型更新混合试验，以此来完成离线模型更新数据的重复利用。时程分析参数及工况分类与框架A相同。对各构件恢复力模型参数的修正如表7所示。
对挖沟机进行了原理样机设计，并分析了需要考虑的基本性能指标,提出了轮腿式挖沟机的构型方案,最后通过Solidworks进行基本构型参数设计,用ADAMS仿真软件对挖沟机进行运动特性仿真分析,得到了挖沟机各部件的运动特性参数曲线,包括车体、车轮的位移、速度、加速度和角速度、角加速度和动能等特性曲线.仿真试验证明:该构型完全满足设计要求,并显示出轮腿式挖沟机在复杂地形地貌下具有良好的稳定性和适应性,为挖沟机的整车结构设计提供了理论依据.
表7框架B中IMK恢复力模型参数预测Table7IMKhystereticmodelparameterpredictionofframeB
对比几种工况下的RMSE值，如表8所示。可以看出对于构件及截面参数相同的不同结构，离线模型更新的试验数据可以通过本文提出的对照经验公式方法进行重复利用，并提高混合试验的精度。
表8框架B各工况RMSE值Table8FrameBRMSEforeachworkingcondition
1)提出了基于对照经验公式的构件本构离线模型更新混合试验方法，通过弹性及弹塑性钢筋混凝土框架离线模型更新混合试验算例，证明该方法可以提高模型更新混合试验的精度。
痊愈为NIHSS评分减少≥90%，脑电图正常；显效为90%＞NIHSS评分减少≥50%，脑电图显著改善；有效为49%＞NIHSS分恢复≥15%，脑电图有所好转；无效为以上指标均未满足。
2)修正更新方法的精度受恢复力模型选择及修正算法等因素影响。使用更精确的恢复力模型和更可靠的修正算法可以进一步提升该方法的修正精度。
这个人喝下热水之后，觉得自己可以站起来了，甚至还可以像想象中一个快死的人那样走路了。他每走一两分钟，就不得不停下来休息一会。他的步子软弱无力，很不稳，就像跟在他后面的那只狼一样又软又不稳；这天晚上，等到黑夜笼罩了光辉的大海的时候，他知道他和大海之间的距离只缩短了不到四哩。
3)离线模型更新混合试验流程精简，本文所提出的对照经验公式修正方法亦可应用于离线模型更新数据的重复利用。
参考文献：
[1]郭玉荣,曾东,肖岩.桥梁结构远程协同拟动力试验平台开发[J].湖南大学学报(自然科学版),2009,36(9):16.GUOYurong,ZENGDong,XIAOYan.Developmentofaremotelycollaborativepseudo-dynamictestingplatformforbridgestructures[J].JournalofHunanUniversity(NaturalSciences),2009,36(9):16.
[2]王涛,吴斌.基于约束UKF模型更新的混合试验方法[J].地震工程与工程振动,2013,33(5):100109.WANGTao,WUBin.HybridtestingmethodbasedonmodelupdatingwithconstrainedunscentedKalmanfilter[J].EarthquakeEngineeringandEngineeringDynamics,2013,33(5):100109.
[3]郭玉荣,刘钟真.基于OpenSees的子结构拟动力试验方法研究[J].湖南大学学报(自然科学版),2015,42(9):17.GUOYurong,LIUZhongzhen.Researchonsubstructurepseudo-dynamictestmethodbasedonOpenSees[J].JournalofHunanUniversity(NaturalSciences),2015,42(9):17.
[4]YANGY,TSAIK,ELNASHAIA,etal.Preliminarystudyononlineupdatinghybridsimulation[R].Taipei,China,2009.
[5]张健.自适应子结构拟动力试验方法[D].哈尔滨:哈尔滨工业大学,2010.ZHANGJian.Adaptivesubstructurepseudo-dynamictestingmethod[D].Harbin:HarbinInstituteofTechnology,2010.
[6]陈凡,郭玉荣.基于AUKF的框架结构离线模型更新混合试验方法[J].地震工程与工程振动,2020,40(1):162170.CHENFan,GUOYurong.AdaptiveunscentedKalmanfilterbasedoff-linemodelupdatinghybridsimulationofframestructures[J].EarthquakeEngineeringandEngineeringDynamics,2020,40(1):162170.
[7]HASHEMIM,MASROORA,MOSQUEDAG.Implementationofonlinemodelupdatinginhybridsimulation[J].EarthquakeEngineering&StructuralDynamics,2014,43(3):395412.
[8]WUBin,WANGTao.ModelupdatingwithconstrainedunscentedKalmanfilterforhybridtesting[J].SmartStructuresandSystems,2014,14(6):11051129.
[9]WUBin,CHENYongsheng,XUGuoshan,etal.Hybridsimulationofsteelframestructureswithsectionalmodelupdating[J].EarthquakeEngineering&StructuralDynamics,2016,45(8):12511269.
[10]ELANWARHH,ELNASHAIAS.Frameworkforonlinemodelupdatinginearthquakehybridsimulations[J].JournalofEarthquakeEngineering,2016,20(1):80100.
[11]梅竹,吴斌,杨格,等.钢筋混凝土结构的模型更新混合试验方法[J].工程科学与技术,2018,50(6):6570.MEIZhu,WUBin,YANGGe,etal.Modelupdatinghybridsimulationofreinforcedconcretestructures[J].AdvancedEngineeringSciences,2018,50(6):6570.
[12]IBARRALF,KRAWINKLERH.Globalcollapseofframestructuresunderseismicexcitations[R].PEER,2005.
[13]IBARRALF,MEDINARA,KRAWINKLERH.Hystereticmodelsthatincorporatestrengthandstiffnessdeterioration[J].EarthquakeEngineering&StructuralDynamics,2005,34(12):14891511.
[14]HASELTONCB,LIELAB,LANGEST.Beam-columnelementmodelcalibratedforpredictingflexuralresponseleadingtoglobalcollapseofRCframebuildings[R].PEER,2007.
[15]余波,李长晋,吴然立.钢筋混凝土柱的非对称恢复力模型与参数识别[J].工程力学,2017,34(2):153161.YUBo,LIChangjin,WURanli.Asymmetricrestoringforcemodleandparameteridentificationofreinforcedconcretecolumn[J].EngineeringMechanics,2017,34(2):153161.