基于麻雀搜索算法考虑病患流的医院诊疗设备维护优化研究

2023-11-16 01:38李岱高刘勤明李佳翔
上海理工大学学报 2023年5期

李岱高 刘勤明 李佳翔

摘要:针对医院诊疗设备出现故障后对病人的影响不可控、产生的维修费用高等问题,提出了一种基于病患流分布不均的预防性维护决策模型。首先,针对医院诊疗设备退化速度受病患流分布影响的特点,采用病患流因子反映其人流量变化规律,并基于设备衰退演化規则构建了设备退化过程模型。其次,通过将量化的诊疗设备故障风险成本以及延迟或提前维护带来的变动成本纳入成本计算体系,提出了病患流分布不均状况下以故障风险值为约束的设备维护策略。最后,采用麻雀搜索算法进行寻优获取关键初始变量,并对模型进行数值仿真分析,结果表明,新维护策略在保障诊疗设备可靠度的同时,也能保持较低的维护成本。

关键词:诊疗设备;病患流分布;故障风险成本;变动成本;麻雀搜索算法

中图分类号:TP 207

文献标志码:A

Maintenance optimization of medical equipment considering patient flow distribution based on the sparrow search algorithm

LI Daigao,LIU Qinming,LI Jiaxiang

(Business School,University of Shanghai for Science and Technology,Shanghai 200093,China)

Abstract:To solve the problems of uncontrollable impact on patients and high maintenance costs after the failure of hospital diagnosis and treatment equipment,a preventive maintenance decision-making model based on uneven distribution of patient flow was proposed.Firstly,in view of the characteristics of the hospital treatment equipment degradation rate affected by the distribution of patient flow,the patient flow factor was used to reflect the changing law of its human flow,and the equipment degradation process model was constructed based on the equipment degradation evolution rule.Secondly,by incorporating the quantified cost of diagnosis and treatment equipment failure risk and the variable cost caused by delayed or early maintenance into the cost calculation system,an equipment maintenance strategy that constrains the failure risk value under the uneven distribution of patient flow was proposed.Finally,through the sparrow search algorithm(SSA)optimization to obtain the key initial variables,and numerical simulation analysis of the model,the results show that the maintenance strategy proposed in this paper can ensure the reliability of the diagnosis and treatment equipment while maintaining alow maintenance cost.

Keywords:medical equipment;distribution of patient flow;failure risk costs;variable costs;sparrow search algorithm

在诊疗活动中,就诊患者人流量受节假日、工作时段等因素影响呈现较大波动,如周末、春节假期尾声等就诊高峰期,经常出现病患集中就诊,导致诊疗设备连续不停长时间运转,而非就诊高峰期内设备运行时间间隔明显增加。病患流时间分布的差异导致诊疗设备运行工况不同;同时,病患流分布不均情况下设备运行时间段的收益不同,进而导致不同时刻停机维护的成本产生较大差异。因此,在设备不同病患流运行工况下,保证诊疗设备可靠度的同时降低维护成本成为亟待解决的问题。

目前,医院各种类型诊疗设备的维护和保养常根据历史经验和维护手册进行[1],易造成维护过剩或维护不足,增加医疗风险。中国《国家医疗器械不良事件监测年度报告(2020年)》显示,2020年共收到医疗器械不良事件报告536055份,相较于2019年增长35.25%。对于医院诊疗设备的预防性维护,不少国内外学者都作了相应的研究。Carmignani[2]提出了基于综合成本的故障影响分析FEMA方法;Abbasgholizadeh Rahimi等[3]用FEMA方法对放射性诊疗设备的故障进行优先排序;Hernández-lópez等[4]采用层次分析法确定医院诊疗设备维护优先级,并给出相应设备权重便于后续评估。然而,他们的研究局限于患者本身的安全性及设备对医院的重要性,忽视了故障维修成本和停机损失成本等客观成本对维护优先级的影响。张宁等[5]建立了基于NSGA-Ⅱ的医疗设备预防性维护与病人调度联合优化模型。但文献[2-5]均未考虑设备运行工况对于设备衰退速率的影响,不能较好地考虑病患流分布不均对诊疗设备维护费用的影响。

研究不同工况下的设备维护,首先需要建立时变工况与设备衰退规律的集成模型。沈南燕等[6]建立了考虑总成本最小化的设备故障率演化模型,并结合支持向量机对模型有效性进行检验;Xia等[7]提出了基于环境变化因子和维修效果的设備预防性维护模型,将外部环境变化对故障率递增的影响纳入模型,但二者并未考虑设备历史工况对设备故障率的影响;Li等[8]提出了考虑故障率刷新因子的设备维护策略,用以求解定期维护和预防性维护相结合的设备故障率递增模型;Vilarinho等[9]将更换部件的成本纳入预防性维护成本中,从而选择更经济的方式对设备进行维护;甘婕等[10]以长期平均维修费用率最小为目标建立了随机规划模型;盖京波等[11]考虑了设备使用时间有限条件下设备预防性维修优化模型;汤乐成等[12]提出了考虑交付时间、产品质量和企业信誉的联合设备维修策略。以上文献均采用了加速失效模型描述设备退化规律,但在医疗设备领域研究较少。

基于上述研究中存在的问题,本文考虑诊疗设备服役过程中病患流分布不均对其退化速度的影响,构建了病患流分布影响下故障率模型;有别于其他研究将故障风险成本与时间变动成本分开研究,本文同时将维护风险、设备重要性等故障风险以及维护时间变动产生的变动成本纳入成本体系,以维护总成本最小为优化目标构建了不同病患流分布下的设备维护调整策略,以满足诊疗设备运行高可靠要求和降低维修费用的需求。最后以某型号的CT机服役数据为例进行仿真分析,并采用麻雀搜索算法对维护策略进行优化,验证了其合理性和有效性。

1问题描述与假设

本文的研究对象主要为大型复杂且不适合长时间连续使用的诊断治疗类医疗设备,例如X射线机、CT机、核磁共振仪等,简称诊疗设备。该类设备使用寿命的衰退趋势会随其连续运行时长的增加而发生变动。例如CT机需要连续曝光产生射线,连续使用多次长时间曝光会导致设备过热,从而影响设备寿命。而就诊病患人流越密集越容易导致诊疗设备的连续运行,就诊排队病患流时域分布不均不仅会造成诊疗设备运行工况的差异,导致设备可靠度的衰退;同时,诊疗设备在病患流高峰期进行维护还将给医院造成较大的经济损失,需要对设备维护计划进行局部修正。如果提前对诊疗设备进行维护,则会造成一定的维护成本浪费,而延迟维护又会造成较大的失效风险。因此,本文基于病患流工况影响建立诊疗设备故障率演化模型,在病患就诊高峰期考虑可靠性及经济性,对诊疗设备制定合理的维护计划。

模型基本假设:

a.诊疗设备以全新状态投入使用,初始役龄为0;

b.在不同病患流量影响下设备失效机理不发生改变,服从加速失效模型;

c.不考虑由于备件不足引发的维修等待,设备发生的各种故障能被及时修复;

d.初始设备故障率服从两参数威布尔分布,形状参数为m,尺度参数为η;

e.每位患者在诊断设备的就诊时间为相同的固定值;

f.病患流区间内病患流因子λ可以表征其真实病患流水平。

2模型建立

2.1病患流分布影响下故障率模型构建

病患流分布密度不同导致各区间设备连续运行时长不同,从而影响医疗设备衰退速率。为更好地表征病患流分布对设备的影响,借鉴文献[13]中对环境应力的等效处理方法,将病患人流量数据进行分段等效处理,由此获得不同病患流区间内等效平均人流量。定义不同病患流区间的人流量因子a i反映其变化规律,即

ai = λi/λ0

(1)

式 中 : λi为 第 i个 病 患 流 区 间 等 效 平 均 人 流 量 ; λ0为该诊疗设备基准人流量, 即在设备设计最佳运 行状态下病患人流量。

诊疗设备在第n个病患流区间的等效役龄tn为

t n=a it I

(2)

式中,ti为第i个病患流区间时长。

诊疗设备在不同病患人流量影响下的等效役龄如图1所示,设备在经历n个病患流区间后的等效役龄t 0n为

在可靠性理论中,可靠度与故障率的关系为[14]

式中:R(t)表示设备可靠度;l(t)表示设备故障率;为设备在当前预防性维护周期内的服役时长。

结合式(2)和式(3),可知在经过n个病患流区间后,诊疗设备的可靠度为

式中,l 0(t)为设备在基准病患流工况下的故障率函数。

执行预防性维护时第n个病患流区间内设备的故障率函数为

虽然诊疗设备运行状态可通过维护来改善,但修复后设备衰退速率较维护前加快,文献[15]使用役龄回退因子和故障率递增因子的概念用以描述这个过程,并验证了其合理性。为刻画维护活动对设备衰退过程的影响,引入役龄回退因子ε和故障率递增因子β描述维护过程对设备故障率的影响,第k次维护后设备的故障率函数为

式中:ε k为第k个预防性维护周期内回退因子,用来描述在役龄层面维护活动对设备故障率的影响;β k为第k个预防性维护周期内故障率递增因子,用来描述不完美维护状况下设备故障率随维护次数增加而变化;εk和β k取值可根据历史维护数据进行估计;T k?1表示第k?1和第k个预防性维护周期时间间隔。

维护策略调整总体流程图如图2所示。

2.2故障风险模型

故障风险的定义为故障发生的概率及故障后果的严重程度,其值由两者量化的结果表示。医院诊疗设备相较于其他普通的生产设备,发生故障造成的后果不仅包含维修成本和停机损失,还包含造成病人安全事故的舆情成本和诉讼成本等。在病患流高峰时,需要对诊疗设备进行提前维护或延迟维护,因此,为辅助决策是否提前或延迟进行维护,需要事先对诊疗设备进行故障风险评估,确定其故障风险系数,使其在保证诊疗活动安全性的情况下尽量兼顾高水平的可靠度。

采用故障风险成本C risk衡量风险后果。

C risk=P dI d

(8)

式中:P d为设备故障发生概率;I d为设备故障发生的后果。

诊疗设备在第n个病患流区间内发生故障的概率[16]为

故障后果由故障维修成本和故障风险因子构成,即

I d=c dφ n

(10)

式中:c d为单位故障维修成本;φ n为故障风险因子,用于衡量设备发生故障后产生的风险。

提前评定医疗设备的故障风险需考虑设备发生故障带来的负面舆情、医院声誉损失、补救成本等多方面因素,参与评估的要素根据设备实际服役过程中关键因素确定[4]。由于涉及模糊环境下考虑多因素的影响问题,因此采用模糊综合评价法提前确定好φ n的值。根据大型诊疗设备的特点分别确定因素集 U = {u1 , u2 , · ·· , ui , ··· , un}和评语集V = {v1 , v2 , · ·· , vi , ··· ,vn},据此构成评价矩阵W为

模糊关系矩阵中各个元素wik′为

通过层次分析法确定因素集中各影响因素的权重[17] ,获得权重集A = [a1 , a2 , · ·· , ai , ··· , an],合成模糊综合评价结果矢量B为

式中: R为模糊关系矩阵; bi为第i个指标的模糊综 合评价结果。

由于医院诊疗设备性质特殊,发生故障的后果相较于其他行业往往更为严重,评判故障风险时需为后果最为严重的情况作准备,故直接取其中最大权重值作为该点故障风险值。

2.3变动成本模型

根据假设 c ,维护活动及时开展, 因此每次维护 的 停 机 时 间 相 同 。图 3 为维 护 时 刻 调 整 示 意图 , 其 中 PM(n?1) , PMn , PM(n+1)分 别 为 不 考 虑 维 护 调 整 时 的 第 n? 1 ,n ,n+ 1次 预 防 性 维 护 时 刻 , PMny为延迟维护时第n次预防性维护时刻, PMnt 为 维护提前时第n次预防性维护时刻。为避免停机损 失过大, 在就诊高峰期时需对诊疗设备进行提前 或延迟维护。提前对诊断设备进行维护会造成其 使用价值的浪费, 而维护活动延迟进行則会产生 一定的设备失效风险, 严重情况下可能导致误诊 或无效诊疗甚至更严重的医疗事故。为保证诊疗 活动的安全性, 延迟维护必须在故障风险限制范 围内进行, 因此, 规定诊疗设备故障风险不能高 于阈值φ0 。

进行设备动态维护时的变动成本ΔC包括设备 提前维护变动成本Ct与延迟维护变动成本Cy 。

提前进行诊疗设备维护可以有效保障诊疗活动的安全性,但会造成其使用价值产生一定的浪费,从而产生一定的经济损失。由于在诊疗活动中对诊疗设备可靠度要求较高,因此将可靠度阈值考虑到维护策略的成本C L建模中十分必要[18]。设备提前维护价值浪费成本的计算公式为

式中: cr为每次维护的可靠度利用价值; Rt为提前 维护时刻设备可靠度; R0为设备维护可靠度阈值。

诊疗设备维护延迟将会导致设备故障风险高于预期,其惩罚成本CY为

式中: cy为延迟维护的基本惩罚成本;Ry为设备延 迟维护时刻可靠度。

维护调整会使诊疗设备在不同周期内的故障 率发生变动, 其故障率变化引起的故障维修变动成本?C d为

式中: Cd为单次故障维修成本; Δt为设备进行维 护调整时其维护时间调整量。

结合式 (15) 、式 (16) 和式 (17),部件提前维护 变动成本Ct与延迟维护变动成本Cy可由以下两式联立计算得出。

假设诊断设备在本次维护周期内进行了J次维护调整,则其由于动态维护调整引起的变动成本为

式中, γ为两种调整成本的决策变量, 当维护提前 时取γ = 1,延迟维护时取γ = 0。

2.4停机成本模型

对于医院诊疗设备来说, 计算停机成本时仅 考虑停机时间是不够的, 其停机成本会随病患流 动态变化。诊疗设备第i次预防性维护的停机损失 成本Csi 为

式中: Tpi 为维护停机时间, 包括故障维修用时和 预防性维护用时; λi为第i个病患流区间内等效平 均病人流量; r为单位病人诊疗平均收益。

若诊疗设备在服役周期内进行了I次预防性维护停机、K次故障停机,则设备维护停机总成本C s为

2.5诊疗设备总维护成本模型

假设诊疗设备在服役周期内进行了I次维护停 机, J次维护调整, K次故障停机, 则其运行过程 中的总成本C包括: 故障风险成本C risk、变动成本 ΔC和停机损失成本Cs ,故考虑病患流分布的诊疗设备维护模型优化目标可表示为

通过求解优化模型获得设备确定合适的最高 故障风险阈值φ0 和诊疗设备最低可靠度阈值R0, 从而得出诊疗设备下一个最优维护时间tk 。

3麻雀搜索算法

3.1算法过程

从上文内容可看出,该优化模型涉及多个非线性函数,理论上难以直接求解,而麻雀搜索算法(sparrow search algorithm,SSA)是受麻雀觅食行为和反捕食行为启发而提出的新型算法,与目标优化领域常用的遗传算法(genetic algorithm,GA)、粒子群算法(particle swarm optimization,PSO)相比,具有参数设置简单、寻优能力强和易于实现的优点[19]。因此,本文采用麻雀搜索算法求解该模型以获取故障风险阈值φ 0和诊疗设备可靠度阈值R 0,算法具体步骤如下:

步骤1 确定种群数量P,初始化麻雀种群

确定种群上界Lub、种群下界Llb和最大迭代次 数ξmax ,设定系统预警值R2和安全值ST;

步骤 2 将总维护成本C(φ, R)作为每只麻雀个 体的适应度函数, 即f = C(φ, R),计算各麻雀个体 的 适 应 度 值 fi并 进 行 排 序 , 确 定 位 置 最 优 个 体 Xbest和位置最差个体Xworst ,更新全局最佳适应度 值fg和全局最差适应度值fw;

步骤 3 通过式 (24) 更新麻雀发现者, 为整个种群寻找食物,并为加入者提供觅食的方向;

式中:x i,j为麻雀加入者个体;α为(0,1)区间随机数;Q为服从正态分布的随机数;L为单位矩阵。

步骤4 在觅食过程中,部分加入者会时刻监视发现者,当发现者发现更好的食物,加入者会与其争夺食物,若成功则获得发现者的食物,若失败则更新加入者位置,通过式(25)更新麻雀发现者和加入者位置;

步骤5 意识到危险的麻雀在种群中的比例为SD,位置由式(26)随机产生;

步骤 6 若当前迭代次数ξ < ξmax ,转到步骤 2,否则输出最优麻雀个体位置信息fg。

3.2基于麻雀搜索算法的维护策略模型求解

考虑病患流分布的诊疗设备维护是一种复杂的多变量优化问题,为将麻雀搜索算法应用于模型求解,需先在麻雀种群矩阵中合理布局变量的位置以便于模型求解。麻雀种群矩阵参数含义如图4所示,矩阵中每一行代表一个单独的麻雀个体。n代表麻雀种群个体数,矩阵k列代表共有k维决策变量。

将麻雀搜索算法应用于优化诊疗设备维护模型的具体步骤如下:

步骤1 读取设置区间内模型所需的诊疗设备所有历史数据与可读信息,并对所得数据进行预处理;

步骤 2 将参数初始化, 输入麻雀种群规模 N、发现者比例 PD、警戒者比例SD、安全值ST、 最大迭代次数ξmax、种群上下限Lub ,Llb。初始化 麻雀种群,生成n只k维麻雀种群矩阵 X。

步骤 3 以诊疗设备总维护成本C(φ, R)最低作 为目标函数, 计算每只麻雀适应度值fi, 排序获得 最优适应度值fg、最差适应度值fw ,及其对应位置 Xbest和Xworst。

步骤 4 根据式 (24) 更新发现者位置, 根据式 (25) 更新加入者位置,根据式 (26) 更新警戒者位置。

步骤 5 计算麻雀个体适应度值, 并更新麻雀 位置,重新确定fg ,fw 以及Xbest和Xworst。

步骤 6 判断迭代次数是否达到最大次数ξmax, 如果是则输出最优适应度值fg作为目标值, 输出最 优麻雀位置Xbest作为最优目标解; 否则转到步骤 4 继续进行运算。

算法求解流程如图5所示。

4算例分析

为验证上述模型的有效性,本文以某医院联影UCT710型号CT机的维护数据为例进行分析。根据假设d,CT机在基准病患流工况下的故障率函数为

4.1算法求解比较

首先,为了验证麻雀搜索算法在考虑病患流分布的医院诊疗设备预防性维护模型中的有效性和优越性,将其与普遍应用在寻优问题上的粒子群算法(PSO)和鲸鱼算法(WOA)进行对比,所有算法的共有参数保持一致,种群规模设置为50,最大迭代次数设置为1000。各算法的参数参照文献[20]设置,如表1所示。表中:W为惯性因子;C 1,C 2为加速度常数;a,b为系数。

比较结果如表2所示,可知SSA算法与WOA算法求解结果相近,且明显优于PSO算法所得结果。各算法迭代如图6所示,可看出PSO算法收敛速度虽比SSA算法和WOA算法快,在迭代到50次左右就达到了最优值,但所得最优值相较于其他两个算法偏差较大。由此可知,PSO算法在此处陷入了局部最优,SSA算法和WOA算法由于较强的全局寻优能力,避免了陷入局部最优。在迭代次数相同且收敛精度相差无几的情况下,SSA算法的收敛速度明显优于WOA算法。

4.2结果分析

通过 SSA 算法寻优得到最佳故障风险阈值Ψ0 = 0.854 5, 最佳设备可靠度阈值R0=0.790 6, 通过数 值仿真求解得到诊疗设备最优维护时间tk等参数如 表 3 所示。

考虑病患流分布情况下诊疗设备故障率演化状况与未考虑病患流分布(即基准病患流)情况下有一定差异。如图7所示,不考虑病患流分布情况下诊疗设备维护时间点延后且维护间隔增加,停机维护成本固然较低,但未考慮病患人流量因素会带来较高的故障风险成本,其区别由是否考虑病患流影响下初始维护方案的不同来反映。将考虑病患流分布记为方案1,基准病患流分布情况记为方案2,对比结果如表4所示。

从表4中可看出,由于不考虑病患流的影响,方案2不产生变动成本,总停机维护成本较低。但忽略病患人流量的影响会埋下更多安全隐患,从而产生更高的故障风险成本,导致总维护成本高,达到40649.8元,高于考虑病患流影响下的总成本。将方案2中维护时间点639~1198代入考虑病患流分布的诊疗设备维护模型中计算故障风险成本和停机维护成本,发现故障风险成本上升,产生惩罚成本15286.5,总成本增加13498.4。因此,方案1相较于方案2能更全面地考虑病患流的影响及其带来的故障风险。考虑病患流的影响在本文框架下有其合理性,建立诊疗设备可靠度和故障率模型时充分考虑病患流造成的影响,有利于更合理地制定设备维护策略。

4.3敏感度分析

为探究考虑病患流分布的诊疗设备维护模型 各参数和模型总维护成本之间的关系, 需要分析 模型对各参数的敏感度。本文采用只改变一个参 数, 其余参数固定不变的方法, 来分析不同参数 取值对总维护成本的影响。需分析参数有故障风 险阈值Ψ0、设备可靠度阈值R0、单位病人平均诊 疗收益 r、可靠度利用价值Cr。参数变化幅度为 ?20%~20%,具体的变化情况如表 5 所示。

由表5结果分析可知:

a. 故障风险成本C risk对故障风险阈值Ψ0和设备 可靠度阈值R0较为敏感, 对单位病人平均诊疗收 益r和可靠度利用价值Cr不敏感。当Ψ0增大时, 对 故障风险值Ψ的要求降低, 制定诊疗设备维护策略 时可接受更高的故障风险, 从而带来更高的故障 风险成本; 同时, 增大的故障风险阈值使维护过 程中的维护变动更少, 从而带来变动成本ΔC下 降。当R0增大时, 对设备可靠度要求更高, 高可 靠度使设备发生故障的风险降低, 进而使故障风 险成本C risk减小; 同时, 高可靠度需要更多次的维 护来保障,因此R0增大导致变动成本ΔC上升。

b. 设备维护停机成本CS对单位病人平均诊疗 收益r最敏感, 对故障风险阈值Ψ0和设备可靠度阈 值R0较为敏感, 对可靠度利用价值Cr不敏感。 r增 大直接影响设备停机维护时产生的病人诊疗费用 损失, 从而增加了停机维护成本。 Ψ0和R0的变动 影响设备停机维护时间点和停机次数, 但停机时 病人诊疗费用损失受该时段病患流密度, 即参数 ai影 响 , 因 此 , CS的 变 动 趋 势 与 Ψ0和 R0无 正 负相关。

c. 设 备 维 护 变 动 成 本 ΔC对 可 靠 度 利 用 价 值 Cr较为敏感, 对单位病人平均诊疗收益 r不敏感。 Cr增大导致设备维护时间提前时产生更多的可靠 度浪费,从而造成维护变动成本增加。

d. 设备维护总成本C对故障风险阈值Ψ0、设备 可靠度阈值R0最敏感, 这两项参数直接影响诊疗 设备预防性维护周期内停机维护与预防性维护次 数; 设备维护总成本C对可靠度利用价值Cr较为敏 感, 因为在模型中参数Cr仅影响维护时间提前时 的维护变动成本ΔC, 而维护变动成本仅为维护总 成本的构成部分, 因此对维护总成本的影响较小。

5结 论

针对诊疗设备这类故障危害大、事故影响坏、故障维修费用和停机损失高昂的特殊设备,将分别考虑安全性、经济性的故障风险方法和考虑设备连续运行状态的病患人流量因子纳入设备预防性维护模型,克服已有模型对于设备故障和设备连续运行对于寿命的影响认识不足的问题。研究结果表明,考量故障风险因素和病患人流量因子使模型更贴合实际,根据模型制定的维护策略更加合理。本文所构建的基于病患流分布不均的预防性维护决策模型较为复杂、求解难度较高,故引入SSA算法工具辅助求解模型,算例分析结果表明,SSA算法能有效提高模型求解速度与精确度。

为诊疗设备制定合理有效的维护策略,对保证诊疗活动安全性和提高其使用效率都具有重大意义。本文提出的基于病患流分布不均的预防性维护决策模型能够帮助医院诊疗设备在运行期间更好地兼顾高可靠度和经济性问题,有效降低设备维护成本,减少设备停机对正常诊疗活动的影响,对诊疗设备维修具有借鉴作用。

未来的研究可考虑多机系统下通过调度进一步降低维护停机成本,降低病人等待时间,使模型实用性更强。另外,鉴于医院诊疗设备是一个复杂的多部件系统,可进一步考虑系统内部件维修之间的关系,使模型更贴合实际。

参考文献:

[1] 韩浩. 医疗设备预防性维护在医疗设备管理中的应用价值分析 [J]. 中国设备工程, 2021(7): 43–44.

[2] CARMIGNANI G. An integrated structural framework to cost-based FMECA: the priority-cost FMECA[J]. Reliability Engineering & System Safety, 2009, 94(4): 861–871.

[3] ABBASGHOLIZADEH RAHIMI S, JAMSHIDI A, AIT- KADI D, et al. Prioritization of failures in radiation therapy delivery[J]. IFAC-PapersOnLine, 2016, 49(12): 1898 – 1903.

[4] HERN?NDEZ-L?PEZ L A, PIMENTEL-AGUILAR A B, ORTIZ-POSADAS M R. An index to prioritize the preventive maintenance of medical equipment[J]. Health and Technology, 2020, 10(2): 399–403.

[5] 張宁, 刘勤明, 叶春明, 等. 基于 NSGA-Ⅱ的医疗设备预 防性维护与病人调度联合优化 [J]. 重庆师范大学学报 (自然科学版), 2021, 38(1): 130–140.

[6] 沈南燕, 武星, 李静, 等. 自动化生产线中关键设备的预 维护策略研究 [J]. 机械工程学报, 2020, 56(21): 231–240.

[7] XIA T B, XI L F, LEE J, et al. Optimal CBPM policy considering maintenance effects and environmental condition[J]. The International Journal of Advanced Manufacturing Technology, 2011, 56(9): 1181–1193.

[8] LI D W, ZHANG Z H, ZHONG Q H, et al. Performance deterioration modeling and optimal preventive maintenance strategy under scheduled servicing subject to mission time[J]. Chinese Journal of Aeronautics, 2014, 27(4): 821–828.

[9] VILARINHO S, LOPES I, OLIVEIRA J A. Preventive maintenance decisions through maintenance optimization models: a case study[J]. Procedia Manufacturing, 2017, 11: 1170–1177.

[10] 甘婕, 曾建潮, 张晓红. 考虑性能可靠度约束的维修决策 模型 [J]. 计算机集成制造系统, 2016, 22(4): 1079–1087.

[11] 盖京波, 孔耀. 有限使用时间内预防性维修策略优化 [J]. 兵工学报, 2015, 36(11): 2164–2172.

[12] 汤乐成, 刘勤明, 叶春明, 等. 考虑批量生产的设备视情 维护策略优化研究 [J]. 上海理工大学学报, 2021, 43(3): 296–302.

[13] 张国龙, 石景岚, 方远, 等. 基于环境应力数据的使用可 靠性预测方法研究 [J]. 电子产品可靠性与环境试验 , 2021, 39(4): 68–71.

[14] 梅嘉健, 刘勤明, 吴擘晟. 基于拆卸序列的多部件设备机 会 维 护 策 略 研 究 [J]. 计 算 机 应 用 研 究 , 2021, 38(12): 3703–3708.

[15] 黄征. 可修系统维修的仿真分析 [D]. 武汉 : 华中农业大 学, 2011.

[16] 陳嘉倩, 刘勤明, 叶春明, 等. 基于故障风险的医院诊疗 设备预防性机会维修策略研究 [J]. 工业工程与管理 , 2021, 26(5): 51–58.

[17] 项伟, 许健, 董芳, 等. 基于模糊数学的甘薯食用品质感 官 评 价 模 型 [J]. 植 物 遗 传 资 源 学 报 , 2021, 22(6): 1624–1634.

[18] 廖雯竹, 潘尔顺, 奚立峰. 基于设备可靠性的动态预防维 护策略 [J]. 上海交通大学学报, 2009, 43(8): 1332–1336.

[19] 刘丽娜, 南新元, 石跃飞. 改进麻雀搜索算法求解作业车 间 调 度 问 题 [J]. 计 算 机 应 用 研 究 , 2021, 38(12): 3634–3639.

[20] 付华 , 刘昊 . 多策略融合的改进麻雀搜索算法及其应 用 [J]. 控制与决策, 2022, 37(1): 87–96.

(编辑:丁红艺)