Summary of the invention
The present invention proposes a kind of structure health monitoring failure sensor data reproduction method of anthropomorphic memory, this method comprises: the multiple configuration state information that sensor acquisition is arrived, be saved in the long-term memory storage area chronologically and form memory mapping, set up the clue relation of configuration state information and memory mapping; According to the clue relation subsequent treatment is discerned and done to the new construction status information: perhaps, according to the clue relation new construction status information is screened to realize the data yojan, perhaps, concern the memory mapping that finds new construction status information correspondence according to clue, the disappearance in the new construction status information is partly compensated, reappears.
Set up the clue relation of configuration state information and memory mapping according to following formula,
In the formula, M is the quantity of memory mapping; N is the quantity of memory clue;
P
j, P
NBeing respectively j kind and the pairing memory clue of N kind configuration state information, is the current data value of j kind and N kind configuration state information also, and j is a numerical sequence, and j=1,2 ... N;
S
kBe k memory mapping of configuration state information; S
MBeing M memory mapping of configuration state information, also is the memory mapping of last position; S
kAnd S
MIn all comprised all original memory clue corresponding historical data values, i.e. the P in the identical moment
1, P
2P
NCorresponding configuration state information is kept in the same memory mapping, and k is a numerical sequence, and k=1, and 2 ... M; { S
1, S
2, S
3S
MThe data acquisition formed is the memory mapping storehouse, initiate memory mapping is designated as S
M+1
W
jIt is the weight index;
R (P
j, S
k) be P
jAnd S
kStrength of joint;
S
k[P
j] be that the memory clue is P
jConfiguration state information be kept at S
kIn the historical data value;
S (P
j, S
k[P
j]) be memory clue P
jWith S
k[P
j] the middle P that preserves
jThe relationship degree formula of historical data value.
S (the P that the present invention adopts
j, S
k[P
j]) expression formula is as follows:
In the formula, P
jBeing the pairing memory clue of j kind configuration state information, also is the current data value of j kind configuration state information;
S
k[P
j] be that the memory clue is P
jConfiguration state information be kept at S
kIn the historical data value;
Method to the new construction status information that collects is discerned comprises: definition new construction status information is the current data value (as strain, pressure, amount of deflection, vibration, temperature etc.) of multiple configuration state information; Current data value P with various configuration state information
j, the clue relational expression F (P of substitution configuration state information and memory mapping
1, P
2P
N):
If F (P
1, P
2P
N) 〉=θ then is identified as " oldState " with this new construction status information, and promptly this new state information is Already in the memory mapping storehouse;
If F (P
1, P
2P
N)<θ then is identified as " new state " with this new construction status information, this new construction status information need be added in the memory mapping storehouse, forms S
M+1
Wherein, F (P
1, P
2P
N) be the clue relational expression of configuration state information and memory mapping information;
θ is set thresholding.
The method that the new construction status information is screened to realize the data yojan is: if a certain new construction status information is identified as " oldState ", then with this new construction status information deletion.
The method that disappearance in the new construction status information is partly compensated, reappears is:
If the current data value of the part-structure status information in a certain new construction status information can't obtain because of sensor failure, and the current data value of another part configuration state information is identified as " oldState ", then according to the pairing memory clue of discernible configuration state information, from the memory mapping storehouse, seek the memory mapping of maximal phase like coupling, according to the memory mapping of the maximal phase that finds the data value that can't obtain because of sensor failure is compensated, reappears like coupling, simultaneously, the pairing P of memory mapping that the maximal phase that finds is seemingly mated
jAnd S
kThe strength of joint fractional increments that adds up makes the time of the memory mapping that search repeatedly visited shorter and shorter.
Compensation, the method for reappearing can be summarized as following steps: be P if the current data of disappearance is worth pairing memory clue 1)
T, then ignore formula
In S (P
T, S
k[P
T]) its result of back calculating; T is a numerical sequence, and 1≤T≤N;
If F (P
1, P
2P
N) satisfy F (P
1, P
2P
N) 〉=θ promptly, removes P in the new construction status information
TThe outer pairing current data value of all the other memory clues has been kept in the data base, judges that then the missing data in this new construction status information has reproducibility;
2) have under the prerequisite of reproducibility at missing data, utilize memory to recall the similarity formula and search maximal phase in the data base like the memory mapping of coupling, and with in the memory mapping that finds with memory clue P
TCorresponding S
k[P
T] give memory clue P to current data value disappearance
TThereby, realize compensation, reproduction to missing data;
It is as follows that the similarity formula is recalled in memory:
Wherein, S
kBe k memory mapping of configuration state information;
P
jIt is the pairing memory clue of j kind configuration state information; J is a numerical sequence, and j=1,2 ... N, N is the quantity of memory clue;
W
jIt is the weight index;
S
k[P
j] be that the memory clue is P
jConfiguration state information be kept at S
kIn the historical data value; K is a numerical sequence, and k=1,2 ... M, M are the quantity of memory mapping;
Searching the method that maximal phase seemingly mates memory mapping is:
If it is P that the current data of disappearance is worth pairing memory clue
T, T is a numerical sequence, and 1≤T≤N, then ignores formula
In S (P
T, S
k[P
T]), calculate all J (P according to formula
1, P
2P
N| S
q), q=1,2 ... M gets minimum J (P
1, P
2P
N| S
q) value (minimum J (P
1, P
2P
N| S
q) the corresponding memory mapping of value is maximal phase and seemingly mates memory mapping); The corresponding ordinal number of the memory mapping that definition finds is K, and the memory mapping that then finds is S
K[P
T], with S
K[P
T] give memory clue P to current data value disappearance
T
Wherein, S (P
T, S
k[P
T]) be memory clue P
TWith S
k[P
T] the middle P that preserves
TThe relationship degree formula of historical data value;
S
K[P
T] be that the memory clue is P
TConfiguration state information be kept at S
KIn the historical data value; 1≤K≤M
S
qBe q memory mapping of configuration state information, S
qIn comprised memory clue P
TPairing historical data value; Q=1,2 ... M;
To P
jAnd S
kThe add up method of a fractional increments of strength of joint is: by memory clue P
1, P
2P
NSearch a certain memory mapping S
kSuccessfully be identified as the memory mapping of maximal phase, then to this memory mapping and corresponding memory clue P like coupling
1, P
2P
NThe pairing strength of joint fractional increments v that adds up, that is:
R(P
j,S
k)=R(P
j,S
k)+v,k=1,2…M
Wherein, P
jBeing the pairing memory clue of j kind configuration state information, is the current data value of j kind configuration state information also, and j is a numerical sequence, and j=1,2 ... N;
S
kBe k memory mapping of configuration state information; K=1,2 ... M, M are memory mapping quantity.
Useful technique effect of the present invention: 1) the configuration state monitoring information is designated " oldState " or " new state " by configuration state memory recognition function, can abandon in a large number the information that is designated " oldState " according to system's needs, thus the data yojan of implementation structure health monitoring; 2) when the sensor failure of long-term work, recall function by the configuration state memory and from data base, recover the Monitoring Data of failure sensor, thereby saved later stage system maintenance cost.
Embodiment
The inventor has proposed a kind of structure health monitoring failure sensor data reproduction method of anthropomorphic memory, the method is the simulating human memory pattern, at the sensor failure problem that occurs in the later maintenance of SHM system, the data rule that model is lived according to memory when certain dimension or several dimension data disappearance is carried out information playback, provide correct assessment result, thereby the monitoring system erroneous judgement of avoiding causing because of the sensor faults itself is broken, and reduces SHM system later maintenance cost simultaneously.This method core is: with the multiple configuration state information that sensor acquisition arrives, be saved in the long-term memory storage area chronologically and form memory mapping, set up the clue relation of configuration state information and memory mapping; According to the clue relation subsequent treatment is discerned and done to the new construction status information: perhaps, according to the clue relation new construction status information is screened to realize the data yojan, perhaps, concern the memory mapping that finds new construction status information correspondence according to clue, the disappearance in the new construction status information is partly compensated, reappears.
Adopt the system of this method can be divided into two working stages: the early stage learning phase, the operation phase in later stage;
Referring to Fig. 3, system is in early stage during learning phase: at a certain works, with sensor acquisition to different types of configuration state information be labeled as P in turn
1, P
2P
N, P
NThe memory clue that just is called N kind configuration state information, P
NNumerical value be the current data that sensor acquisition arrives; At T
1Constantly, with the P that collects
1, P
2P
NData are recorded in memory mapping S
1In, with T
2The P that constantly collects
1, P
2P
NData are recorded in memory mapping S
2In, mode like this is with the P in the follow-up identical moment
1, P
2P
NData are recorded in the memory mapping S of serial number
kIn, after memory mapping quantity forms on a fairly large scale, can think that memory mapping storehouse (data base) basically forms; It is worthy of note, be not that a dirt is constant after data base forms, and also may need to add new memory mapping in the operational process.
Operation phase in later stage: structure is after long-time the use, and the sensor that has is owing to expire fatigue lifetime, or other reason damage, causes the operative sensor can't image data, at this moment, and the P of the sign structural health conditions that collects
1, P
2P
NData lack with regard to some, and under most of sensors can also the situation of operate as normal, the data that we can arrive according to the sensor acquisition of operate as normal, search in the data base the similar memory mapping of part normal data therewith, historical data according to corresponding missing data in the memory mapping is reappeared missing data, for the evaluation of structural health conditions provides foundation.Data reproduction is the aspect that the present invention uses, and the another application of the invention aspect is, according to the P that newly collects
1, P
2P
NData are searched and whether are preserved the similar memory mapping that similarity meets the demands in the data base, if having then delete the P that newly collects
1, P
2P
NData realize the yojan of system data.
During set forth summary of the present invention the front, the Several Key Problems that relates to was: how to set up configuration state information and concern with the clue of memory mapping, how to search memory mapping similar to normal data in (identification) data base, how to judge the similarity of memory mapping and new data; Further describe with regard to above 3 problems below.
(1) memory identification rule:
GillundG and ShiffrinRM have proposed a fairly simple human brain memory familiarity formula:
This formula is by introducing strength of joint R (P
j, S
k) index has been set up between memory clue and memory mapping has the mapping relations that intensity quantizes, for analyst's brain memory function provides the mathematical model basis.Use for reference the characteristics of this technological thought and integrated structure health monitoring, the present invention introduces relationship degree computing function in preceding formula, proposes following memory identification familiarity formula:
In the formula, M is the quantity of memory mapping; N is the quantity of memory clue;
P
j, P
NBeing respectively j kind and the pairing memory clue of N kind configuration state information, is the current data value of j kind and N kind configuration state information also, and j is a numerical sequence, and j=1,2 ... N;
S
kBe k memory mapping of configuration state information; S
MBeing M memory mapping of configuration state information, also is the memory mapping of last position; S
kAnd S
MIn all comprised all original memory clue corresponding historical data values, i.e. the P in the identical moment
1, P
2P
NCorresponding configuration state information is kept in the same memory mapping, and k is a numerical sequence, and k=1, and 2 ... M; { S
1, S
2, S
3S
MThe data acquisition formed is the memory mapping storehouse, initiate memory mapping is designated as S
M+1
W
jBe reflection memory clue P
jTo memory mapping S
kThe weight index of importance;
R (P
j, S
k) be P
jAnd S
kStrength of joint;
S
k[P
j] be that the memory clue is P
jConfiguration state information be kept at S
kIn the historical data value;
S (P
j, S
k[P
j]) be memory clue P
jWith S
k[P
j] the middle P that preserves
jThe relationship degree formula of historical data value, can be according to the demand structure of SHM application system; S (the P that the present invention adopts
j, S
k[P
j]) expression formula is as follows:
In the formula, P
jBeing the pairing memory clue of j kind configuration state information, also is the current data value of j kind configuration state information;
S
k[P
j] be that the memory clue is P
jConfiguration state information be kept at S
kIn the historical data value;
To set thresholding θ, memory identification rule is as follows:
1) if F is (P
1, P
2P
N) 〉=θ, data processing model is one " oldState " to the configuration state of this memory clue representative of system prompt, Already in the memory mapping storehouse;
2) if F is (P
1, P
2P
N)<θ, data processing model is one " new state " to the configuration state of this memory clue representative of system prompt, needs to add to advance in the memory mapping storehouse, produces memory mapping S
M+1, S
M+1Be at learning phase and finish the new memory mapping (correlating mutually) that adds in the operational process of back with preamble.
Above-mentioned memory identifying is effectively managed the sensor acquisition data, and mass of redundancy data is identified as " oldState " and rejects, and what enter the long-term memory district is the reduced set of memory mapping information.
(2) rule is recalled in memory:
Configuration state is identified as the memory clue P of " oldState "
1, P
2P
N, if seek its similar record in the memory mapping storehouse, this process is promptly remembered and is recalled.Consider to remember to recall to come down to a process of seeking maximal phase like occurrence, the present invention proposes following memory and recalls the similarity formula:
Wherein, S
kBe k memory mapping of configuration state information;
P
jIt is the pairing memory clue of j kind configuration state information; J is a numerical sequence, and j=1,2 ... N, N is the quantity of memory clue;
W
jIt is the weight index;
S
k[P
j] be that the memory clue is P
jConfiguration state information be kept at S
kIn the historical data value; K is a numerical sequence, and k=1,2 ... M, M are the quantity of memory mapping;
Searching the method that maximal phase seemingly mates memory mapping is:
If it is P that the current data of disappearance is worth pairing memory clue
T, T is a numerical sequence, and 1≤T≤N, then ignores formula
In S (P
T, S
k[P
T]), calculate all J (P according to formula
1, P
2P
N| S
q), q=1,2 ... M gets minimum J (P
1, P
2P
N| S
q) value; The corresponding ordinal number of the memory mapping that definition finds is K, and the memory mapping that then finds is S
K[P
T], with S
K[P
T] give memory clue P to current data value disappearance
T
Wherein, S (P
T, S
k[P
T]) be memory clue P
TWith S
k[P
T] the middle P that preserves
TThe relationship degree formula of historical data value;
S
K[P
T] be that the memory clue is P
TConfiguration state information be kept at S
KIn the historical data value; 1≤K≤M
S
qBe q memory mapping of configuration state information, S
qIn comprised memory clue P
TPairing historical data value; Q=1,2 ... M;
The human brain memory theory thinks that repetitive learning can enhance memory, sets each other off therewith, and the present invention is by promoting P
jAnd S
kStrength of joint R (P
j, S
k) value simulates this memory characteristic: recall by once successful configuration state memory, find the memory mapping S of maximal phase like coupling
kWith memory clue P
1, P
2P
NBe complementary, then be provided with:
R(P
j,S
k)=R(P
j,S
k)+v,k=1,2…M
Wherein, P
jBeing the pairing memory clue of j kind configuration state information, is the current data value of j kind configuration state information also, and j is a numerical sequence, and j=1,2 ... N;
S
kBe k memory mapping of configuration state information; K=1,2 ... M, M are memory mapping quantity.
Formula R (P
j, S
k)=R (P
j, S
kThe implication of)+v is: in observation process when long, if the data of certain memory mapping are recalled (recalling) continually, strength of joint between this memory mapping and the memory clue will be enhanced, this means that this memory mapping recalls and have more probability preferentially to be wandered back in the process in memory next time, also promptly reduce search time.According to probability theory, S
kThe chance probability of preferentially being recalled can calculate by following formula:
In the formula, formula
Reflection S
kAnd P
1, P
2P
NBetween the strength of joint comprehensive effect.
Value big more, memory mapping data S
kPossiblely more wandered back to by priority access.
(3) failure sensor data reproduction:
Memory identification rule provides the criterion of failure sensor data reproducibilitys, possesses the reproducibility of data through the memory identification decision for the configuration state of " oldState " that is:, can further recall rule by memory and realize data reproduction.With above-mentioned two rules as shown in Figure 4 in conjunction with the data reproduction process that is applied to failure sensor.Its method is as follows:
If the current data value of the part-structure status information in a certain new construction status information can't obtain because of sensor failure, and the current data value of another part configuration state information is identified as " oldState ", then according to the pairing memory clue of discernible configuration state information, from the memory mapping storehouse, seek the memory mapping of maximal phase, the data value that can't obtain because of sensor failure is compensated, reappears according to the memory mapping of the maximal phase that finds like coupling like coupling.
Specific practice is:
1) be P if the current data of disappearance is worth pairing memory clue
T, then ignore formula
In S (P
T, S
k[P
T]) its result of back calculating; T is a numerical sequence, and 1≤T≤N;
If F (P
1, P
2P
N) satisfy F (P
1, P
2P
N) 〉=θ promptly, removes P in the new construction status information
TThe outer pairing current data value of all the other memory clues has been kept in the data base, judges that then the missing data in this new construction status information has reproducibility;
2) have under the prerequisite of reproducibility at missing data, utilize memory to recall the similarity formula and search maximal phase in the data base like the memory mapping of coupling, and with in the memory mapping that finds with memory clue P
TCorresponding S
k[P
T] give memory clue P to current data value disappearance
TThereby, realize compensation, reproduction to missing data;
Further set forth method of the present invention below in conjunction with embodiment.
Embodiment:
A free beam strain monitoring experiment porch, experimental system is as shown in Figure 5; The experiment porch main body by
1200mm strides the little bent steel plate and the moment of torsion force application apparatus in footpath and forms, equidistant six groups of strain measurement full-bridge circuits that distribute on little bent steel plate based on the metal forming resistance strain gage, and numbered sequence is 6,4,2,1,3,5 from left to right.Turn to pressure and act on little bent steel plate center by the torque axis that will manually apply, cause the steel plate structure strain, see shown in Figure 6.The strain of each monitoring location of foil gauge full-bridge circuit sensing also uploads to host computer through data collector and is for data processing.Data collector adopts NI CRIO of company 9237 and embedded controller CRIO-9014, is responsible for gathering and handles obtaining strain value.Upload to industrial computer through Ethernet interface and carry out the data analysis processing.Each strain bridge is typical Huygens's full-bridge circuit, and circuit theory and installation diagram are seen shown in Fig. 7,8.
Corresponding with memory models, the function class of NI CRIO-9237 is like immediate memory, and the structural strain state is converted to electric signal by the collection of full-bridge strainometer, and CRIO-9237 gathers this signal and passes to controller CRIO-9014 as the immediate memory memory block.After this, controller CRIO-9014 is strain data as the short-term memory district in the memory models with the electric signal conversion processing and the data combination of diverse location is the memory clue.Industrial computer is forever preserved useful information as the long-term memory district in the memory models and is carried out works of treatment such as configuration state data reproduction.
The experiment of configuration state data reproduction is divided into two parts: configuration state memory identification and configuration state memory are recalled.Degenerative process by driving torque handle model configuration.Because free beam is to have flexible steel, so during the revolution of moment of torsion handle, free beam can return to its original state.This specific character makes whole experiment have repeatability.We are at first to structure memory input priori, and then carry out configuration state memory identification and process is recalled in memory.
(1) configuration state memory identification experiment
At first, obtain basic data memory by a learning process.The moment of torsion handle is by Rotate 180 degree successively clockwise, 7 times like this.Six groups of strain datas of a sequence are stored the into long-term memory district of industrial computer, are shown in Table 1.
Table 1
This experiment is with relationship degree formula S (P
j, S
k[P
j]) be defined as follows:
If P
j=S
k[P
j], S (P then is set
j, S
k[P
j])=10
-6In addition, all R (P
j, S
k) initial value all be made as 0.1, k=1,2 ... M; W
j(j=1,2 ... N) all be made as 1.
Simulated data input memory mapping S with the memory clue
T1, S
T2... S
T8And substitution formula
To determine familiarity thresholding θ, obtain following memory mapping sequence:
S
t1=[1 50?16 130?50?30]*(-1.0000e-006)
S
t2=[27?30?147?26 23?30]*(-1.0000e-006)
S
t3=[27?41?147?30 23?25]*(-1.0000e-006)
S
t4=[27?41?166?40 25?38]*(-1.0000e-006)
S
t5=[27?41?166?150?26?40]*(-1.0000e-006)
S
t6=[27?41?166?150?57?40]*(-1.0000e-006)
S
t7=[27?41?166?150?57?18]*(-1.0000e-006)
S
t8=[26?42?164?151?58?20]*(-1.0000e-006)
With the memory mapping S in the table 1
4In the memory clue and the memory clue in the above-mentioned memory mapping compare: S
T1And S
4Different fully, S
T2A clue element and S are arranged
4The same, S
T3Two clue elements and S are arranged
4The same, S
T4, S
T5, S
T6, S
T7By that analogy.S
T8Each clue element all and S
4Approaching.
The result of calculation of familiarity is seen Fig. 9, and as shown in the figure, the memory clue has been divided into different grades with familiarity between memory mapping, and this shows certain particular system, can distinguish the familiarity grade by defining a thresholding θ.
This experiment is provided with thresholding θ=20.After this with the moment of torsion handle from 0 the degree rotate to 180 * n (n=1,2 ...) and the degree.To rotation each time, the data reproduction model can both be that " oldState " or " new state " judge to the memory mapping data according to thresholding indicating structure strain regime.
For example, when the moment of torsion handle is transferred to 180 * 3 when spending, the clue data in the memory mapping are [23 32 148 129 4312] * (1.0000e-006), familiarity result of calculation F (P
1, P
2P
N)=3935>θ, this shows that this strain regime is one " oldState " once took place; When the moment of torsion handle is transferred to 180 * 9 when spending, the clue data in the memory mapping are [47 103 268,251 118 38] * (1.0000e-006), F (P
1, P
2P
NIt is one " new state " that)=0.0059<θ, data reproduction model identify this strain regime thus, should store in the long-term memory district.New memory mapping data storage is seen shown in Figure 10.
(2) experiment is recalled in the configuration state memory
The memory clue that is identified as " oldState " can be recalled its matched data of searching in memory mapping storage data by the configuration state memory.For example, be [23 32 148 129 43 12] * (1.0000e-006), application of formula to memory clue data
Similarity result of calculation J (P
1, P
2P
N| S
k)=[7.9763 1.6262 0.1880 0.8710 1.3751 1.7296 1.9918 2.3833].Because min (J (P
1, P
2P
N| S
k))=0.1880, so the memory mapping S that wanders back to
kSequence number be 3.This memory clue data and memory mapping data S accurately
3=[24 30 147 130 44 13] * is (1.0000e-006) close.This shows that configuration state memory recalls model and can be used as to disturb and suppress to filter, even the Monitoring Data of obtaining has certain noise, recalls process by memory and can eliminate noise effect.
(3) configuration state failure sensor data reproduction experiment
But utilize the memory identification of configuration state and the data recovery that function implementation structure monitors failure sensor is recalled in memory.For example, strainometer 4 and 1 test plug forward the moment of torsion handle to about 180 * 9 positions in pulling out Fig. 6, and * (1.0000e-006) for [46 X, 268 X 119 38] for the memory clue data of obtaining.If P
2=P
4=0, by formula
Calculate S (P
2, S
k[P
2])=1, S (P
4, S
k[P
4])=1, F (P
1, P
2P
N)=56>θ, clue is imperfect although this shows memory, and it is one " oldState " for memory mapping, has reproducibility.By formula
Can calculate J (P
1, P
2P
N| S
k)=[14.8458 8.0691 5.3024 3.8335 3.1072 2.6952 2.4910 2.3577].Because min (J (P
1, P
2P
N| S
k))=2.3577, therefore remember the memory mapping S that wanders back to
kSequence number be 8.Can from the memory mapping storehouse, find in view of the above and lack big data value P
2=S
8[P
2]=103, P
4=S
8[P
4]=251.The data of acquired original and data recovered are seen shown in Figure 11.
Above-mentioned experiment supposition R (P
j, S
k) (k=1,2 ... value M) is a constant basis.If consider the enhance memory strength of joint, formula R (P can be set
j, S
k)=R (P
j, S
kFractional increments v=0.001 among the)+v is as memory mapping data S
3Recalled and visited when reaching 10 times memory mapping data S
k(k=1,2 ... 8) when being visited once more by priority access to Probability Distribution see Figure 12.From Figure 12 as seen, by recalled the visit 10 times after, S
3Have higher probability preferentially to be wandered back to, this phenomenon and human brain memory characteristic are approximate.