|
@@ -0,0 +1,699 @@
|
|
|
+import pandas as pd
|
|
|
+import numpy as np
|
|
|
+import datetime
|
|
|
+from LIB.MIDDLE.CellStateEstimation.Common.V1_0_1 import BatParam
|
|
|
+
|
|
|
+class BatDiag:
|
|
|
+ def __init__(self,sn,celltype,df_bms,df_diag_ram_sn,df_bms_ram_sn,df_soh,df_uniform): #参数初始化
|
|
|
+
|
|
|
+ self.sn=sn
|
|
|
+ self.celltype=celltype
|
|
|
+ self.param=BatParam.BatParam(celltype)
|
|
|
+ self.df_bms=df_bms
|
|
|
+ self.df_uniform=df_uniform
|
|
|
+ self.df_soh=df_soh
|
|
|
+ self.df_bms_ram=df_bms_ram_sn.copy()
|
|
|
+ self.df_diag_ram=df_diag_ram_sn.copy()
|
|
|
+ df_bms['时间戳']=pd.to_datetime(df_bms['时间戳'], format='%Y-%m-%d %H:%M:%S')
|
|
|
+
|
|
|
+ if not df_bms_ram_sn.empty:
|
|
|
+ self.df_bms=self.df_bms[self.df_bms['时间戳'] > df_bms_ram_sn.iloc[-1]['time']] #滤除原始数据中的重复数据
|
|
|
+ self.df_bms.reset_index(inplace=True,drop=True) #重置索引
|
|
|
+
|
|
|
+ self.packcrnt=df_bms['总电流[A]']*self.param.PackCrntDec
|
|
|
+ self.packvolt=df_bms['总电压[V]']
|
|
|
+ self.bmstime= df_bms['时间戳']
|
|
|
+ self.bms_soc=df_bms['SOC[%]']
|
|
|
+
|
|
|
+ self.cellvolt_name=['单体电压'+str(x) for x in range(1,self.param.CellVoltNums+1)]
|
|
|
+ # othertemp=['其他温度'+str(x) for x in range(1,self.param.OtherTempNums+1)]
|
|
|
+ self.celltemp_name=['单体温度'+str(x) for x in range(1,self.param.CellTempNums+1)]
|
|
|
+ # self.celltemp_name=celltemp+othertemp
|
|
|
+
|
|
|
+
|
|
|
+ def diag(self):
|
|
|
+ if self.celltype<=50:
|
|
|
+ df_res1,df_res2=self._bat_diag()
|
|
|
+ return df_res1,df_res2
|
|
|
+ else:
|
|
|
+ df_res1,df_res2=self._bat_diag()
|
|
|
+ return df_res1,df_res2
|
|
|
+
|
|
|
+
|
|
|
+ #寻找当前行数据的所有温度值...................................................................................
|
|
|
+ def _celltemp_get(self,num):
|
|
|
+ celltemp = list(self.df_bms.loc[num,self.celltemp_name])
|
|
|
+ return celltemp
|
|
|
+
|
|
|
+ #获取当前行所有电压数据............................................................................................
|
|
|
+ def _cellvolt_get(self,num):
|
|
|
+ cellvolt = list(self.df_bms.loc[num,self.cellvolt_name]/1000)
|
|
|
+ return cellvolt
|
|
|
+
|
|
|
+ #..........................................电池诊断功能..................................................................
|
|
|
+ def _bat_diag(self):
|
|
|
+
|
|
|
+ column_name=['start_time', 'end_time', 'product_id', 'code', 'level', 'info','advice']
|
|
|
+ df_res=pd.DataFrame(columns=column_name)
|
|
|
+
|
|
|
+ end_time='0000-00-00 00:00:00'
|
|
|
+ ah_accum=0 #SOC卡滞初始参数
|
|
|
+ as_chg=0 #过流诊断初始参数
|
|
|
+ as_dis=0 #过流诊断初始参数
|
|
|
+ cellvoltvalid=1
|
|
|
+ voltdsc_time=0
|
|
|
+ voltdsc_count=0
|
|
|
+ voltfail_time=0
|
|
|
+ voltfail_count=0
|
|
|
+ voltloose_time=0
|
|
|
+ voltloose_count=0
|
|
|
+ voltover_time=0
|
|
|
+ voltunder_time=0
|
|
|
+ cov_time=0
|
|
|
+ cuv_time=0
|
|
|
+ cdv_time=0
|
|
|
+ pov_time=0
|
|
|
+ puv_time=0
|
|
|
+ ot_time=0
|
|
|
+ ut_time=0
|
|
|
+ dt_time=0
|
|
|
+
|
|
|
+ for i in range(len(self.df_bms)):
|
|
|
+
|
|
|
+ #电压诊断功能.........................................................................................................................................
|
|
|
+ if i<1:
|
|
|
+ cellvolt2=self._cellvolt_get(i)
|
|
|
+ cellvoltmin2=min(cellvolt2)
|
|
|
+ cellvoltmax2=max(cellvolt2)
|
|
|
+ cellvoltmin_index2=cellvolt2.index(cellvoltmin2)
|
|
|
+ cellvoltmax_index2=cellvolt2.index(cellvoltmax2)
|
|
|
+ if not self.df_bms_ram.empty:
|
|
|
+ time1=self.df_bms_ram.iloc[-1]['time']
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ cellvolt1=self.df_bms_ram.iloc[-1]['cellvolt']
|
|
|
+ cellvoltmin1=min(cellvolt1)
|
|
|
+ cellvoltmax1=max(cellvolt1)
|
|
|
+ cellvoltmin_index1=cellvolt1.index(cellvoltmin1)
|
|
|
+ cellvoltmax_index1=cellvolt1.index(cellvoltmax1)
|
|
|
+ else:
|
|
|
+ cellvolt1=cellvolt2
|
|
|
+ time1=self.bmstime[i]
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ cellvoltmin1=cellvoltmin2
|
|
|
+ cellvoltmax1=cellvoltmax2
|
|
|
+ cellvoltmin_index1=cellvoltmin_index2
|
|
|
+ cellvoltmax_index1=cellvoltmax_index2
|
|
|
+ else:
|
|
|
+ time1=self.bmstime[i-1]
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ cellvolt2=self._cellvolt_get(i)
|
|
|
+ cellvoltmin2=min(cellvolt2)
|
|
|
+ cellvoltmax2=max(cellvolt2)
|
|
|
+ cellvoltmin_index2=cellvolt2.index(cellvoltmin2)
|
|
|
+ cellvoltmax_index2=cellvolt2.index(cellvoltmax2)
|
|
|
+ cellvolt1=self._cellvolt_get(i-1)
|
|
|
+ cellvoltmin1=min(cellvolt1)
|
|
|
+ cellvoltmax1=max(cellvolt1)
|
|
|
+ cellvoltmin_index1=cellvolt1.index(cellvoltmin1)
|
|
|
+ cellvoltmax_index1=cellvolt1.index(cellvoltmax1)
|
|
|
+
|
|
|
+ #电压采样断线..........................................................................................................................................
|
|
|
+ if not 'C308' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if (cellvoltmin2<2 and cellvoltmax2>4.5 and abs(cellvoltmax_index2-cellvoltmin_index2)==1) and (cellvoltmin1<2 and cellvoltmax1>4.5 and abs(cellvoltmax_index1-cellvoltmin_index1)==1): #电压断线故障进入
|
|
|
+ cellvoltvalid=0
|
|
|
+ voltdsc_time=voltdsc_time+(time2-time1).total_seconds()
|
|
|
+ if voltdsc_time>self.param.volt_time: #持续时间
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C308'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电芯{}和{}电压采样断线'.format(cellvoltmin_index2+1, cellvoltmax_index2+1)
|
|
|
+ faultadvice='通知用户更换电池,电池返厂维修'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ elif (cellvoltmin2<2 and cellvoltmax2>4.5 and abs(cellvoltmax_index2-cellvoltmin_index2)==1) and (cellvoltmin1>2.5 and cellvoltmax1<4.3): #连续跳变
|
|
|
+ cellvoltvalid=0
|
|
|
+ voltdsc_count=voltdsc_count+1
|
|
|
+ if voltdsc_count>=3:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C308'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电芯{}和{}电压采样断线'.format(cellvoltmin_index2+1, cellvoltmax_index2+1)
|
|
|
+ faultadvice='通知用户更换电池,电池返厂维修'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ voltdsc_time=0
|
|
|
+ else:
|
|
|
+ cellvoltvalid=0
|
|
|
+ if cellvoltmin2>2 and cellvoltmax2<4.5 and cellvoltmin1>2 and cellvoltmax1<4.5:
|
|
|
+ cellvoltvalid=1
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C308'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #电压采样系统失效.............................................................................................................
|
|
|
+ if not 'C309' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if (cellvoltmin2<1 and cellvoltmin1<1) or (cellvoltmax2>5 and cellvoltmax1>5):
|
|
|
+ cellvoltvalid=0
|
|
|
+ voltfail_time=voltfail_time+(time2-time1).total_seconds()
|
|
|
+ if voltfail_time>self.param.volt_time: #持续时间
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C309'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电芯电压采样系统失效'
|
|
|
+ faultadvice='通知用户更换电池,电池返厂维修'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ elif (cellvoltmin2<1 and cellvoltmin1>2.5) or (cellvoltmax2>5 and cellvoltmax1<4.5): #连续跳变
|
|
|
+ cellvoltvalid=0
|
|
|
+ voltfail_count=voltfail_count+1
|
|
|
+ if voltfail_count>=3:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C309'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电芯电压采样系统失效'
|
|
|
+ faultadvice='通知用户更换电池,电池返厂维修'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ voltfail_time=0
|
|
|
+ else:
|
|
|
+ cellvoltvalid=0
|
|
|
+ if cellvoltmin2>2.5 and cellvoltmax2<4.5 and cellvoltmin1>2.5 and cellvoltmax1<4.5:
|
|
|
+ cellvoltvalid=1
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C309'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #电压采样线松动.................................................................................................................
|
|
|
+ if not 'C208' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellvoltmin2>2 and cellvoltmax2<4.5 and cellvoltmin1>2 and cellvoltmax1<4.5:
|
|
|
+ cellvolt1_std=np.std(cellvolt1)
|
|
|
+ cellvolt1_mean=np.mean(cellvolt1)
|
|
|
+ cellvolt1_3sigma=(np.array(cellvolt1)-cellvolt1_mean)/cellvolt1_std
|
|
|
+ cellvolt2_std=np.std(cellvolt2)
|
|
|
+ cellvolt2_mean=np.mean(cellvolt2)
|
|
|
+ cellvolt2_3sigma=(np.array(cellvolt2)-cellvolt2_mean)/cellvolt2_std
|
|
|
+ if (min(cellvolt1_3sigma)<-3 and max(cellvolt1_3sigma)>3) and (min(cellvolt2_3sigma)<-3 and max(cellvolt2_3sigma)>3): #连续发生
|
|
|
+ cellvoltvalid=0
|
|
|
+ voltloose_time=voltloose_time+(time2-time1).total_seconds()
|
|
|
+ if voltloose_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C208'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电芯电压{}和{}采样松动'.format(cellvoltmin_index2+1, cellvoltmax_index2+1)
|
|
|
+ faultadvice='通知技术人员介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ elif (min(cellvolt1_3sigma)>-3 and max(cellvolt1_3sigma)<3) and (min(cellvolt2_3sigma)<-3 and max(cellvolt2_3sigma)>3): #连续跳变
|
|
|
+ cellvoltvalid=0
|
|
|
+ voltloose_count=voltloose_count+1
|
|
|
+ if voltloose_count>=3:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C208'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电芯电压{}和{}采样松动'.format(cellvoltmin_index2+1, cellvoltmax_index2+1)
|
|
|
+ faultadvice='通知技术人员介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ voltloose_time=0
|
|
|
+ else:
|
|
|
+ voltloose_time=0
|
|
|
+ else:
|
|
|
+ cellvoltvalid=0
|
|
|
+ cellvolt1_std=np.std(cellvolt1)
|
|
|
+ cellvolt1_mean=np.mean(cellvolt1)
|
|
|
+ cellvolt1_3sigma=(np.array(cellvolt1)-cellvolt1_mean)/cellvolt1_std
|
|
|
+ cellvolt2_std=np.std(cellvolt2)
|
|
|
+ cellvolt2_mean=np.mean(cellvolt2)
|
|
|
+ cellvolt2_3sigma=(np.array(cellvolt2)-cellvolt2_mean)/cellvolt2_std
|
|
|
+ if(min(cellvolt1_3sigma)>-3 or max(cellvolt1_3sigma)<3) and (min(cellvolt2_3sigma)>-3 or max(cellvolt2_3sigma)<3):
|
|
|
+ cellvoltvalid=1
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C208'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ if cellvoltvalid==1:
|
|
|
+ #继电器粘连.............................................................................................................................................
|
|
|
+ if not 'C310' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellvoltmax2>self.param.CellOvLv2 and cellvoltmax1>self.param.CellOvLv2 and self.packcrnt[i]<0:
|
|
|
+ voltover_time=voltover_time+(time2-time1).total_seconds()
|
|
|
+ if voltover_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C310'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='充电继电器粘连'
|
|
|
+ faultadvice='通知技术人员介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ voltover_time=0
|
|
|
+ else:
|
|
|
+ if cellvoltmax2<self.param.CellOvLv2-0.05 and cellvoltmax1<self.param.CellOvLv2-0.05:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C310'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ if not 'C311' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellvoltmin2<self.param.CellUvLv2 and cellvoltmin1<self.param.CellUvLv2 and self.packcrnt[i]>0:
|
|
|
+ voltunder_time=voltunder_time+(time2-time1).total_seconds()
|
|
|
+ if voltunder_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C311'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='放电继电器粘连'
|
|
|
+ faultadvice='通知技术人员介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ voltunder_time=0
|
|
|
+ else:
|
|
|
+ if cellvoltmin2>self.param.CellUvLv2+0.05 and cellvoltmin1>self.param.CellOvLv2+0.05:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C311'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #电芯过压.............................................................................................................................................
|
|
|
+ if not 'C501' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellvoltmax2>self.param.CellOvLv2 and cellvoltmax1>self.param.CellOvLv2 and self.packcrnt[i]<0: #二级过压进入
|
|
|
+ cov_time=cov_time+(time2-time1).total_seconds()
|
|
|
+ if cov_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C501'
|
|
|
+ faultlv=5
|
|
|
+ faultinfo='电芯{}过压'.format(cellvolt1.index(cellvoltmax1)+1)
|
|
|
+ faultadvice='联系用户立即停止充电,并通知技术运维人员'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ cov_time=0
|
|
|
+ else: #ram当前故障中有该故障,则判断是否退出该故障
|
|
|
+ if cellvoltmax2<self.param.CellOvLv1-0.05 and cellvoltmax1<self.param.CellOvLv1-0.05: #二级过压故障恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C501'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #欠压诊断.................................................................................................................
|
|
|
+ if not 'C202' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellvoltmin2<self.param.CellUvLv2 and cellvoltmin1<self.param.CellUvLv2: #二级欠压
|
|
|
+ cuv_time=cuv_time+(time2-time1).total_seconds()
|
|
|
+ if cuv_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C202'
|
|
|
+ faultlv=2
|
|
|
+ faultinfo='电芯{}欠压'.format(cellvolt1.index(cellvoltmin1)+1)
|
|
|
+ faultadvice='联系用户询问用车场景,技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ cuv_time=0
|
|
|
+ else:
|
|
|
+ if cellvoltmin2>self.param.CellUvLv1+0.1 and cellvoltmin1>self.param.CellUvLv1+0.1:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C202'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #电芯压差大.....................................................................................................................................................
|
|
|
+ if not 'C104' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if (cellvoltmax2-cellvoltmin2)>self.param.CellVoltDiffLv2 and (cellvoltmax1-cellvoltmin1)>self.param.CellVoltDiffLv2: #二级电芯压差
|
|
|
+ cdv_time=cdv_time+(time2-time1).total_seconds()
|
|
|
+ if cdv_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C104'
|
|
|
+ faultlv=0
|
|
|
+ faultinfo='电芯{}和{}压差大'.format(cellvolt1.index(cellvoltmax1)+1,cellvolt1.index(cellvoltmin1)+1)
|
|
|
+ faultadvice='技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ cdv_time=0
|
|
|
+ else:
|
|
|
+ if (cellvoltmax2-cellvoltmin2)<self.param.CellVoltDiffLv1-0.05 and (cellvoltmax1-cellvoltmin1)<self.param.CellVoltDiffLv1-0.05: #二级欠压恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C104'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ cov_time=0
|
|
|
+ cuv_time=0
|
|
|
+ cdv_time=0
|
|
|
+
|
|
|
+ #电池包诊断.....................................................................................................................................
|
|
|
+ if i<1:
|
|
|
+ if not self.df_bms_ram.empty:
|
|
|
+ packvolt1=self.df_bms_ram.iloc[-1]['packvolt']
|
|
|
+ packvolt2=self.packvolt[i]
|
|
|
+ time1=self.df_bms_ram.iloc[-1]['time']
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ else:
|
|
|
+ packvolt1=self.packvolt[i]
|
|
|
+ packvolt2=self.packvolt[i]
|
|
|
+ time1=self.bmstime[i]
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ else:
|
|
|
+ packvolt1=self.packvolt[i-1]
|
|
|
+ packvolt2=self.packvolt[i]
|
|
|
+ time1=self.bmstime[i-1]
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ if packvolt1<2*self.param.CellVoltNums or packvolt2>4.5*self.param.CellVoltNums: #电池包电压有效性
|
|
|
+ packvoltvalid=0
|
|
|
+ else:
|
|
|
+ packvoltvalid=1
|
|
|
+
|
|
|
+ if packvoltvalid==1:
|
|
|
+ if not 'C502' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if packvolt1>self.param.PackVoltOvLv2 and packvolt2>self.param.PackVoltOvLv2 and self.packcrnt[i]<0: #电池包过压二级进入
|
|
|
+ pov_time=pov_time+(time2-time1).total_seconds()
|
|
|
+ if pov_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C502'
|
|
|
+ faultlv=5
|
|
|
+ faultinfo='电池包过压'
|
|
|
+ faultadvice='联系用户立即停止充电,并通知技术运维人员'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ pov_time=0
|
|
|
+ else:
|
|
|
+ if packvolt1<self.param.PackVoltOvLv1-0.05*self.param.CellVoltNums and packvolt2<self.param.PackVoltOvLv1-0.05*self.param.CellVoltNums: #电池包过压二级恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C502'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ if not 'C203' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if packvolt1<self.param.PackVoltUvLv2 and packvolt2<self.param.PackVoltUvLv2: #电池包二级欠压进入
|
|
|
+ puv_time=puv_time+(time2-time1).total_seconds()
|
|
|
+ if puv_time>self.param.volt_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C203'
|
|
|
+ faultlv=2
|
|
|
+ faultinfo='电池包欠压'
|
|
|
+ faultadvice='联系用户询问用车场景,技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ puv_time=0
|
|
|
+ else:
|
|
|
+ if packvolt1>self.param.PackVoltUvLv1+0.1*self.param.CellVoltNums and packvolt2>self.param.PackVoltUvLv1+0.1*self.param.CellVoltNums: #电池包二级欠压恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C203'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ pov_time=0
|
|
|
+ puv_time=0
|
|
|
+
|
|
|
+ #温度诊断功能.............................................................................................................
|
|
|
+ if i<1:
|
|
|
+ if not self.df_bms_ram.empty:
|
|
|
+ celltemp1=self.df_bms_ram.iloc[-1]['celltemp']
|
|
|
+ celltemp2=self._celltemp_get(i)
|
|
|
+ time1=self.df_bms_ram.iloc[-1]['time']
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ else:
|
|
|
+ celltemp1=self._celltemp_get(i)
|
|
|
+ celltemp2=self._celltemp_get(i)
|
|
|
+ time1=self.bmstime[i]
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ else:
|
|
|
+ celltemp1=self._celltemp_get(i-1)
|
|
|
+ celltemp2=self._celltemp_get(i)
|
|
|
+ time1=self.bmstime[i-1]
|
|
|
+ time2=self.bmstime[i]
|
|
|
+ celltempmin1=min(celltemp1)
|
|
|
+ celltempmin2=min(celltemp2)
|
|
|
+ celltempmax1=max(celltemp1)
|
|
|
+ celltempmax2=max(celltemp2)
|
|
|
+ #温度有效性判断..........................................................................
|
|
|
+ if celltempmax1>self.param.CellTempUpLmt or celltempmin1<self.param.CellTempLwLmt:
|
|
|
+ celltempvalid=0
|
|
|
+ else:
|
|
|
+ celltempvalid=1
|
|
|
+
|
|
|
+ if celltempvalid==1:
|
|
|
+ #过温判断.............................................................................................................
|
|
|
+ if not 'C302' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if celltempmax1>self.param.CellTempHighLv2 and celltempmax2>self.param.CellTempHighLv2: #二级高温进入
|
|
|
+ ot_time=ot_time+(time2-time1).total_seconds()
|
|
|
+ if ot_time>self.param.temp_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C302'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='温度{}高温'.format(celltemp1.index(celltempmax1)+1)
|
|
|
+ faultadvice='技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ ot_time=0
|
|
|
+ else: #ram当前故障中有该故障,则判断是否退出该故障
|
|
|
+ if celltempmax1<self.param.CellTempHighLv1-5 and celltempmax2<self.param.CellTempHighLv1-5: #二级高温恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C302'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #欠温判断.................................................................................................................
|
|
|
+ if not 'C102' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if celltempmin1<self.param.CellTempLowLv2 and celltempmin2<self.param.CellTempLowLv2: #二级低温进入
|
|
|
+ ut_time=ut_time+(time2-time1).total_seconds()
|
|
|
+ if ut_time>self.param.temp_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C102'
|
|
|
+ faultlv=1
|
|
|
+ faultinfo='温度{}低温'.format(celltemp1.index(celltempmin1)+1)
|
|
|
+ faultadvice='技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ ut_time=0
|
|
|
+ else: #ram当前故障中有该故障,则判断是否退出该故障
|
|
|
+ if celltempmax1>self.param.CellTempLowLv1+2 and celltempmax2>self.param.CellTempLowLv1+2: #二级高温恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C102'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #温差判断.............................................................................................................................
|
|
|
+ if not 'C103' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if (celltempmax1-celltempmin1)>self.param.CellTempDiffLv2 and (celltempmax2-celltempmin2)>self.param.CellTempDiffLv2: #二级温差进入
|
|
|
+ dt_time=dt_time+(time2-time1).total_seconds()
|
|
|
+ if dt_time>self.param.temp_time:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C103'
|
|
|
+ faultlv=1
|
|
|
+ faultinfo='温度{}和{}温差大'.format(celltemp1.index(celltempmax1)+1,celltemp1.index(celltempmin1)+1)
|
|
|
+ faultadvice='技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ dt_time=0
|
|
|
+ else: #ram当前故障中有该故障,则判断是否退出该故障
|
|
|
+ if (celltempmax1-celltempmin1)<self.param.CellTempDiffLv1-2 and (celltempmax2-celltempmin2)>self.param.CellTempDiffLv1-2: #二级温差恢复
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C103'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ else:
|
|
|
+ ot_time=0
|
|
|
+ ut_time=0
|
|
|
+ dt_time=0
|
|
|
+
|
|
|
+ #电流过流诊断.......................................................................................................................
|
|
|
+ if i>0.5:
|
|
|
+ step=(self.bmstime[i]-self.bmstime[i-1]).total_seconds()
|
|
|
+ if step<120 and self.packcrnt[i]>self.param.PackDisOc and self.packcrnt[i-1]>self.param.PackDisOc:
|
|
|
+ as_dis=as_dis+(self.packcrnt[i]-self.param.PackDisOc)*step #ah累计
|
|
|
+ elif step<120 and self.packcrnt[i]<self.param.PackChgOc and self.packcrnt[i-1]<self.param.PackChgOc:
|
|
|
+ as_chg=as_chg+(self.param.PackDisOc-self.packcrnt[i])*step #ah累计
|
|
|
+ else:
|
|
|
+ as_dis=0
|
|
|
+ as_chg=0
|
|
|
+
|
|
|
+ if not 'C306' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if as_dis>100:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C306'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电池包放电过流'
|
|
|
+ faultadvice='联系用户询问用车场景,技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if self.packcrnt[i]<self.param.PackDisOc-10:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C306'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ if not 'C305' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if as_chg>100:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C305'
|
|
|
+ faultlv=3
|
|
|
+ faultinfo='电池包充电过流'
|
|
|
+ faultadvice='联系用户询问用车场景,技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if self.packcrnt[i]>self.param.PackChgOc+10:
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C305'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #SOC卡滞、跳变诊断................................................................................................
|
|
|
+ if i<1:
|
|
|
+ bmssoc_st=float(self.bms_soc[i])
|
|
|
+ bmssoc_last=float(self.bms_soc[i])
|
|
|
+ bmssoc_now=float(self.bms_soc[i])
|
|
|
+ else:
|
|
|
+ step=(self.bmstime[i]-self.bmstime[i-1]).total_seconds()
|
|
|
+ if step<120:
|
|
|
+ ah_accum=ah_accum-self.packcrnt[i]*step/3600 #ah累计
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ #SOC卡滞............................................................................................................
|
|
|
+ if abs(ah_accum)>self.param.Capacity*0.1:
|
|
|
+ bmssoc_now=float(self.bms_soc[i])
|
|
|
+ if not 'C106' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if abs(bmssoc_now-bmssoc_st)<self.param.SocClamp: #SOC卡滞故障进入
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C106'
|
|
|
+ faultlv=0
|
|
|
+ faultinfo='电池SOC卡滞'
|
|
|
+ faultadvice='技术介入诊断,检修电池BMS软件'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if abs(bmssoc_now-bmssoc_st)>self.param.SocClamp: #SOC卡滞故障退出
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C106'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ bmssoc_st=bmssoc_now
|
|
|
+ ah_accum=0
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #SOC跳变....................................................................................................................
|
|
|
+ bmssoc_last=float(self.bms_soc[i-1])
|
|
|
+ bmssoc_now=float(self.bms_soc[i])
|
|
|
+ if not 'C107' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if step<70 and abs(bmssoc_now-bmssoc_last)>self.param.SocJump: #SOC跳变进入
|
|
|
+ time=self.bmstime[i]
|
|
|
+ code='C107'
|
|
|
+ faultlv=0
|
|
|
+ faultinfo='电池SOC跳变{}%'.format(bmssoc_now-bmssoc_last)
|
|
|
+ faultadvice='技术介入诊断,检修电池BMS软件'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if abs(bmssoc_now-bmssoc_st)<self.param.SocJump: #SOC跳变故障退出
|
|
|
+ time=self.bmstime[i]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C107'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ #SOC一致性故障报警..........................................................................................................
|
|
|
+ if not self.df_uniform.empty:
|
|
|
+ cellsoc_diff=self.df_uniform.loc[0,'cellsoc_diff']
|
|
|
+ if not 'C201' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellsoc_diff>self.param.SocDiff: #SOC一致性差故障进入
|
|
|
+ time=self.bmstime[0]
|
|
|
+ code='C201'
|
|
|
+ faultlv=0
|
|
|
+ faultinfo='电芯{}和{}SOC差过大:{}'.format(int(self.df_uniform.loc[0,'cellmin_num']),int(self.df_uniform.loc[0,'cellmax_num']),cellsoc_diff)
|
|
|
+ faultadvice='技术介入诊断'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if cellsoc_diff<self.param.SocDiff: #SOC一致性差故障恢复
|
|
|
+ time=self.bmstime[0]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C201'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ cellsoc_diff=3
|
|
|
+
|
|
|
+ #容量过低和一致性故障报警................................................................................................
|
|
|
+ if not self.df_soh.empty:
|
|
|
+ soh=self.df_soh.loc[0,'soh']
|
|
|
+ cellsoh=eval(self.df_soh.loc[0,'cellsoh'])
|
|
|
+ cellsoh=np.array(cellsoh)
|
|
|
+ cellsoh_lowindex=np.argwhere(cellsoh<self.param.SohLow)
|
|
|
+ cellsoh_lowindex=cellsoh_lowindex+1
|
|
|
+ if self.celltype==1 or self.celltype==2 or self.celltype==3 or self.celltype==4:
|
|
|
+ cellsoh_diff=np.max(cellsoh)-np.min(cellsoh)
|
|
|
+ if not 'C204' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if soh<self.param.SohLow: #soh过低故障进入
|
|
|
+ time=self.bmstime[0]
|
|
|
+ code='C204'
|
|
|
+ faultlv=2
|
|
|
+ faultinfo='电池包容量过低:电芯{}'.format(cellsoh_lowindex)
|
|
|
+ faultadvice='检修电池,更换容量过低的电芯或模组'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if soh>self.param.SohLow+2: #soh过低故障恢复
|
|
|
+ time=self.bmstime[0]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C204'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+
|
|
|
+ if not 'C205' in list(self.df_diag_ram['code']): #当前故障中没有该故障,则判断是否发生该故障
|
|
|
+ if cellsoh_diff>self.param.SohDiff:
|
|
|
+ time=self.bmstime[0]
|
|
|
+ code='C205'
|
|
|
+ faultlv=2
|
|
|
+ faultinfo='电池包容量一致性差:电芯{}'.format(cellsoh_lowindex)
|
|
|
+ faultadvice='检修电池,更换容量过低的电芯或模组'
|
|
|
+ self.df_diag_ram.loc[len(self.df_diag_ram)]=[time, end_time, self.sn, code, faultlv, faultinfo, faultadvice]
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ if cellsoh_diff<self.param.SohDiff-2:
|
|
|
+ time=self.bmstime[0]
|
|
|
+ self.df_diag_ram.loc[self.df_diag_ram[self.df_diag_ram['code']=='C205'].index, 'end_time'] = time
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ pass
|
|
|
+ else:
|
|
|
+ cellsoh_diff=5
|
|
|
+
|
|
|
+ #ram更新.......................................................................
|
|
|
+ df_bms_ram=pd.DataFrame(columns=['time', 'sn', 'packvolt', 'cellvolt', 'celltemp'])
|
|
|
+ df_bms_ram.loc[0]=[self.bmstime[0], self.sn, packvolt2, cellvolt2, celltemp2]
|
|
|
+
|
|
|
+ return self.df_diag_ram, df_bms_ram
|
|
|
+
|