Apparatus and method for damage detection

Claim 3

The method of claim 1, wherein said indicator function indicates a change in energy stored in each said location.

Claim 4

The method of claim 1, wherein said indicator function comprises a function of a pre-damage sensitivity matrix, F 0 , and of a post-damage sensitivity matrix, F * , wherein said pre-damage sensitivity matrix relates changes in said pre-damage modal parameters to a fractional change in stiffness of one or more elements of the structure, and wherein said post-damage sensitivity matrix relates changes in said post-damage modal parameters to each said fractional change.

Claim 5

The method of claim 3, wherein said indicator function comprises a function of a difference between a first element of said pre-damage sensitivity matrix and a corresponding first element of said post-damage sensitivity matrix, wherein said first and corresponding first elements correspond to a particular mode shape and location in the structure.

Claim 6

The method of claim 4, wherein said indicator function comprises said difference raised to a power of N, where N is a real number.

Claim 7

The method of claim 4, wherein said indicator function comprises said difference multiplied by N, where N is a real number.

Claim 8

The method of claim 3, wherein said indicator function comprises a function of a ratio between a first element, F ik , of said pre-damage sensitivity matrix and a corresponding first element, F * ik , of said post-damage sensitivity matrix, wherein said first and corresponding first elements correspond to a particular mode shape i and location k in the structure.

Claim 9

The method of claim 7 wherein said indicator function further comprises a function of a ratio between (a) a pre-damage differential displacement (Φ u -Φ iC ) k of said location k with respect to said particular mode shape i, where said location k is connected to a first node A and to a second node C, and where Φ u and Φ iC are elements of a pre-damage modal quantities mode shape vectors matrix Φ, for said particular mode shape i and (b) a post-damage differential displacement (Φ iA * -Φ iC * ) k of said location k with respect to said particular mode shape i, where Φ u * and Φ iC * are elements of a post-damage modal quantities mode shape vectors matrix Φ i * for said particular mode shape i.

Claim 10

The method of claim 8, wherein said indicator further comprises a function of a ratio between a pre-damage modal stiffness, K i , for said particular mode shape i and a post-damage modal stiffness, K * i , for said particular mode shape i.

Claim 11

The function of claim 9 wherein said indicator function is a function of ##EQU26##

Claim 12

The method of claim 8 wherein said indicator function is a function of ##EQU27##

Claim 13

The method of claim 3 wherein said indicator function is a function of an angle between a first mathematical surface and a second mathematical surface, wherein said first and second surfaces correspond respectively to said pre-damage and said post-damage sensitivity matrices, and wherein said angle is determined at a first point on said first surface corresponding to a particular mode shape and location in the structure and at a corresponding first point on said second surface corresponding to said particular mode shape and said location.

Claim 14

The method of claim 12, wherein said angle is a function of gradient responses of said first and second surfaces at said first and corresponding first points, respectively.

Claim 15

The method of claim 1, wherein said indicator function is a function of a difference between (a) a first ratio between a first fractional change in modal parameters corresponding to a first mode shape and a second fractional change in modal parameters corresponding to a second mode shape, and (b) a second ratio between a first element of a pre-damage sensitivity matrix corresponding to said first mode shape and a second element of said pre-damage sensitivity matrix corresponding to said second mode shape, wherein said pre-damage sensitivity matrix relates changes in said pre-damage modal parameters to a fractional change in stiffness of one or more elements of the structure.

Claim 16

The method of claim 1, wherein said indicator function is a function of a difference between (a) a first ratio between a first fractional change in modal parameters corresponding to a first mode shape and a second fractional change in modal parameters corresponding to a second mode shape, and (b) a second ratio between a first element of a post-damage sensitivity matrix corresponding to said first mode shape and a second element of said post-damage sensitivity matrix corresponding to said second mode shape, wherein said post-damage sensitivity matrix relates changes in said post-damage modal parameters to a fractional change in stiffness of one or more elements of the structure.

Claim 17

The method of claim 15, wherein said first and second fractional changes in modal parameters comprise fractional changes in damped eigenvalue frequencies.

Claim 18

The method of claim 1, wherein said first data comprise data indicative of theoretical modal parameters of the structure.

Claim 19

The method of claim 17, and further comprising the step of adapting said theoretical modal parameters to correspond with said second data.

Claim 20

The method of claim 1, and further comprising the step of determining a magnitude of the damage at each said location based upon each said identified location and upon said first and second data.

Claim 21

The method of claim 1, wherein said second data determining step comprises the step of measuring a post-damage vibration frequency response of the structure.

Claim 22

The method of claim 1, wherein said pre-damage and post-damage modal parameters respectively comprise at least one pre-damage resonant frequency and at least one post-damage resonant frequency.

Claim 23

The method of claim 1, wherein said pre-damage and post-damage mode shapes each comprise at least one maximum deflection of at least one location of the structure.

Claim 24

The method of claim 22, wherein said pre-damage and post-damage modal parameters respectively comprise at least one pre-damage damping value and at least one post-damage damping value.

Claim 25

The method of claim 1, wherein said indicator function is a function of a pre-damage sensitivity and of a post-damage sensitivity.

Claim 26

The method of claim 24, wherein said indicator function is a function of a difference between said pre-damage sensitivity and said post-damage sensitivity.

Claim 27

The method of claim 24, wherein said indicator function is a function of a gradient response of said pre-damage sensitivity and of said post-damage sensitivity.

Claim 28

The method of claim 1 wherein said identifying step comprises the step of determining one or more Z matrices of fractional changes in modal parameters for one or more vibrational modes, based upon said first and second data.

Claim 29

The method of claim 27 wherein one of said Z matrices comprises a matrix of fractional changes in damped eigenvalue frequencies Z d .

Claim 30

The method of claim 27 wherein one of said Z matrices comprises a matrix of fractional changes in modal masses Z m .

Claim 31

The method of claim 27 wherein one of said Z matrices comprises a matrix of fractional changes in modal damping ratios Z damp .

Claim 32

The method of claim 27 wherein one of said Z matrices comprises a matrix of fractional changes in the norm of mode shapes Z modes .

Claim 33

The method of claim 27 and further comprising the steps of: generating a stiffness matrix of the structure based on one or more material and geometric properties of the structure; using said stiffness matrix and said first data, computing a pre-damage sensitivity matrix F o which relates said Z matrices to a fractional change in stiffness of one or more elements of the structure; and using said stiffness matrix and said second data, computing a sensitivity change matrix ΔF which accounts for changes in said pre-damage sensitivity matrix F o resulting from damage to the structure.

Claim 34

The method of claim 32 wherein said indicator function relates an assumed location of damage and an assumed severity of damage to a damage detection equation.

Claim 35

The method of claim 33 wherein said indicator function is given by ##EQU28## wherein a represents said assumed location, b represents said assumed severity, α represents a matrix of fractional changes in stiffness of one or more elements in the structure, k represents one of said elements, n is a real number, and wherein said damage detection equation is given by ΣZ matrices=(Fo+ΔF)α(k).

Claim 40

The method of claim 36, wherein said location identifying step comprises the step of applying said first and second data to an indicator function.

Claim 43

The method of claim 38 wherein the number of elements B exceeds the number of vibrational modes N, the method further comprising the steps of: identifying at least one location of the damage using an indicator function; eliminating elements of the damage detection equation based upon each said location to form a reduced set damage detection equation; and solving for an element of α corresponding to each said location using the reduced set damage detection equation.

Claim 44

The method of claim 38, wherein said element of α is solved for by dividing an element of said quantity F 0 +ΔF into an element of at least one said matrix of modal parameters, wherein said element of said matrix F 0 +ΔF corresponds to a particular mode and to a particular location, wherein said element of each said matrix of modal parameters corresponds to said particular mode, and wherein said element of a corresponds said particular location.

Claim 45

The method of claim 38, wherein the number of elements B exceeds the number of vibrational modes N, the method further comprising the steps of: inverting the matrix F 0 +ΔF using one or more pseudo inverse techniques; and solving for α using the damage detection equation.

Claim 46

The method of claim 38, wherein the number of vibrational modes N exceeds the number of elements B, the method further comprising the steps of: inverting the matrix F 0 +ΔF using one or more exact inverses; and solving for α using the damage detection equation.

Claim 47

The method of claim 38, wherein the number of vibrational modes N exceeds the number of elements B, the method further comprising the steps of: inverting the matrix F 0 +ΔF using at least one over-determined least squares approach; and solving for α using the damage detection equation.

Claim 48

The method of claim 43, wherein said least squares approach is performed according to the expression: F.sup.1 =F.sup.T (FF.sup.T).sup.-1.

Claim 49

The method of claim 38, wherein at least one element of the matrix F 0 +ΔF is eliminated if, after inverting the matrix F 0 +ΔF, a corresponding element of α has a value outside a predetermined range.

Claim 50

The method of claim 45, and further comprising the step of inverting the matrix F 0 +ΔF after eliminating each said element of the quantity matrix F 0 +ΔF.

Claim 51

The method of claim 38, wherein each of the B elements is modelled as a spring with element stiffness matrix given by ##EQU29## wherein i is an element number, E is Young's modulus of elasticity, I is a second moment of area, and L is a length of element number i.

Claim 52

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer operable to measure a first signal set indicative of pre-damage modal parameters of the structure in response to a first excitation force, said transducer further operable to determine a second signal set indicative of post-damage modal parameters of the structure in response to a second excitation force; a converter operable to convert said first signal set to first data indicative of pre-damage modal parameters of the structure, said pre-damage modal parameters including at least one measured mode shape, said converter further operable to convert said second signal set to second data indicative of post-damage modal parameters of the structure, said post-damage modal parameters including at least one measured mode shape; and means coupled to said converter operable to identify at least one location of the damage based upon said first and second data applied to an indicator function.

Claim 53

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer operable to measure a first signal set indicative of pre-damage modal parameters of the structure in response to a first excitation force, said transducer further operable to determine a second signal set indicative of post-damage modal parameters of the structure in response to a second excitation force; a converter operable to convert said first signal set to first data indicative of pre-damage modal parameters of the structure, said pre-damage modal parameters including at least one measured mode shape, said converter further operable to convert said second signal set to second data indicative of post-damage modal parameters of the structure, said post-damage modal parameters including at least one measured mode shape; and means coupled to said converter operable to identify at least one location of the damage based upon said first and second data applied to an indicator function.

Claim 54

The apparatus of claim 28, wherein said indicator function indicates a change in energy stored in each said location.

Claim 55

The apparatus of claim 48, wherein said indicator function comprises a function of a pre-damage sensitivity matrix, F 0 , and of a post-damage sensitivity matrix, F * , wherein said pre-damage sensitivity matrix relates changes in said pre-damage modal parameters to a fractional change in stiffness of one or more elements of the structure, and wherein said post-damage sensitivity matrix relates changes in said post-damage modal parameters to each said fractional change.

Claim 56

The apparatus of claim 50, wherein said indicator function comprises a function of a difference between a first element of said pre-damage sensitivity matrix and a corresponding first element of said post-damage sensitivity matrix, wherein said first and corresponding first elements correspond to a particular mode shape and location in the structure.

Claim 57

The apparatus of claim 51, wherein said indicator function comprises said difference raised to a power of N, where N is a real number.

Claim 58

The apparatus of claim 51, wherein said indicator function comprises said difference multiplied by N, where N is a real number.

Claim 59

The apparatus of claim 50, wherein said indicator function comprises a function of a ratio between a first element, F ik , of said pre-damage sensitivity matrix and a corresponding first element, F * ik , of said post-damage sensitivity matrix, wherein said first and corresponding first elements correspond to a particular mode shape i and location k in the structure.

Claim 60

The apparatus of claim 54, wherein said indicator function further comprises a function of a ratio between (a) a pre-damage differential displacement (Φ u -Φ iC ) k of said location k with respect to said particular mode shape i, where said location k is connected to a first node A and to a second node C, and where Φ u and Φ iC are elements of a pre-damage modal quantities mode shape vectors matrix Φ i for said particular mode shape i and (b) a post-damage differential displacement (Φ u * -Φ iC * ) k of said location k with respect to said particular mode shape i, where Φ u * and Φ iC * are elements of a post-damage modal quantities mode shape vectors matrix Φ i * for said particular mode shape i.

Claim 61

The apparatus of claim 55, wherein said indicator further comprises a function of a ratio between a pre-damage modal stiffness, K i , for said particular mode shape i and a post-damage modal stiffness, K * i , for said particular mode shape i.

Claim 62

The apparatus of claim 56, wherein said indicator function is a function of ##EQU30##

Claim 63

The apparatus of claim 55 wherein said indicator function is a function of ##EQU31##

Claim 64

The apparatus of claim 50 wherein said indicator function is a function of an angle between a first mathematical surface and a second mathematical surface, wherein said first and second surfaces correspond respectively to said pre-damage and said post-damage sensitivity matrices, and wherein said angle is determined at a first point on said first surface corresponding to a particular mode shape and location in the structure and at a corresponding first point on said second surface corresponding to said particular mode shape and said location.

Claim 65

The apparatus of claim 59, wherein said angle is a function of gradient responses of said first and second surfaces at said first and corresponding first points, respectively.

Claim 66

The apparatus of claim 48, wherein said indicator function is a function of a difference between (a) a first ratio between a first fractional change in modal parameters corresponding to a first mode shape and a second fractional change in modal parameters corresponding to a second mode shape, and (b) a second ratio between a first element of a pre-damage sensitivity matrix corresponding to said first mode shape and a second element of said pre-damage sensitivity matrix corresponding to said second mode shape, wherein said pre-damage sensitivity matrix relates changes in said pre-damage modal parameters to a fractional change in stiffness of one or more elements of the structure.

Claim 67

The apparatus of claim 48, wherein said indicator function is a function of a difference between (a) a first ratio between a first fractional change in modal parameters corresponding to a first mode shape and a second fractional change in modal parameters corresponding to a second mode shape, and (b) a second ratio between a first element of a post-damage sensitivity matrix corresponding to said first mode shape and a second element of said post-damage sensitivity matrix corresponding to said second mode shape, wherein said post-damage sensitivity matrix relates changes in said post-damage modal parameters to a fractional change in stiffness of one or more elements of the structure.

Claim 68

The apparatus of claim 62, wherein said first and second fractional changes in modal parameters comprise fractional changes in damped eigenvalue frequencies.

Claim 69

The apparatus of claim 48, wherein said identifying means comprises means for determining one or more Z matrices of fractional changes in modal parameters for one or more vibrational modes, based upon said first and second data.

Claim 70

The apparatus of claim 64, wherein one of said Z matrices comprises a matrix of fractional changes in damped eigenvalue frequencies Z d .

Claim 71

The apparatus of claim 64, wherein one of said Z matrices comprises a matrix of fractional changes in modal masses Z m .

Claim 72

The apparatus fo claim 64, wherein one of said Z matrices comprises a matrix of fractional changes in modal damping ratios Z damp .

Claim 73

The apparatus of claim 64, wherein one of said Z matrices comprises a matrix of fractional changes in the form of mode shapes Z modes .

Claim 74

The apparatus of claim 64, and further comprising: means for generating a stiffness matrix of the structure based on one or more material and geometric properties of the structure; means coupled to said stiffness matrix generating means and to said first data determining means for using said stiffness matrix and said first data to compute a pre-damage sensitivity matrix F o which relates said Z matrices to a fractional change in stiffness of one or more elements of the structure; and means coupled to said stiffness matrix generating means and to said second data determining means for using said stiffness matrix and said second data to compute a sensitivity change matrix ΔF which accounts for changes in said pre-damage sensitivity matrix F o resulting from damage to the structure.

Claim 75

The apparatus of claim 69, wherein said indicator function relates an assumed location of damage and an assumed severity of damage to a damage detection equation.

Claim 76

The apparatus of claim 70, wherein said indicator function is given by ##EQU32## wherein a represents said assumed location, b represents said assumed severity, α represents a matrix of fractional changes in stiffness of one or more elements in the structure, k represents one of said elements, n is a real number, and wherein said damage detection equation is given by ΣZ matrices=(Fo+ΔF)α(k).

Claim 77

The apparatus of claim 48, wherein said first data comprise data indicative of theoretical modal parameters of the structure.

Claim 78

The apparatus of claim 72, and further comprising means for adapting said theoretical modal parameters to correspond with said second data.

Claim 79

The apparatus of claim 48, and further comprising means coupled to said identifying means for determining a magnitude of the damage at each said location based upon each said identified location and upon said first and second data.

Claim 80

The apparatus of claim 48, wherein said transducer is operable to measure a post-damage vibration frequency response of the structure.

Claim 81

The apparatus of claim 48, wherein said pre-damage and post-damage modal parameters respectively comprise at least one pre-damage resonant frequency and at least one post-damage resonant frequency.

Claim 82

The apparatus of claim 34, wherein said pre-damage and post-damage mode shapes each comprise at least one maximum deflection of at least one location of the structure.

Claim 83

The apparatus of claim 77, wherein said pre-damage and post-damage modal parameters respectively comprise at least one pre-damage damping value and at least one post-damage damping value.

Claim 84

The apparatus fo claim 48, wherein said indicator function is a function of a pre-damage sensitivity and of a post-damage sensitivity.

Claim 85

The apparatus of claim 79, wherein said indicator function is a function of a difference between said pre-damage sensitivity and said post-damage sensitivity.

Claim 86

The apparatus of claim 79, wherein said indicator function is a function of a gradient response of said pre-damage sensitivity and of said post-damage sensitivity.

Claim 87

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer coupled to the structure operable to measure a pre-damage vibration frequency response of the structure in response to a first excitation force, said transducer further operable to measure a post-damage vibration frequency response of the structure in response to a second excitation force; a converter operable to convert said pre-damage vibration frequency response to a first data indicative of pre-damage modal parameters of the structure, said pre-damage modal parameters including at least one measured mode shape, said converter further operable to convert said post-damage vibration frequency response to second data indicative of post-damage modal parameters of the structure, said post-damage modal parameters including at least one measured mode shape; and means coupled to said converter for identifying at least one location of the damage based upon said first and second data.

Claim 88

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer coupled to the structure operable to measure a pre-damage vibration frequency response of the structure in response to a first excitation force, said transducer further operable to measure a post-damage vibration frequency response of the structure in response to a second excitation force; a converter operable to convert said pre-damage vibration frequency response to a first data indicative of pre-damage modal parameters of the structure, said pre-damage modal parameters including at least one measured mode shape, said converter further operable to convert said post-damage vibration frequency response to second data indicative of post-damage modal parameters of the structure, said post-damage modal parameters including at least one measured mode shape; and means coupled to said converter for identifying at least one location of the damage based upon said first and second data.

Claim 89

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer operable to measure a first signal set in response to a first excitation force, said transducer further operable to measure a second signal set in response to a second excitation force; a converter operable to convert said first signal set to first data indicative of pre-damage modal parameters of the structure, said pre-damage modal parameters including at least one measured mode shape, said converter further operable to convert said second signal; set to second data indicative of post-damage modal parameters of the structure, said post-damage modal parameters including at least one measured mode shape; means coupled to said converter for identifying at least one location of the damage based upon said first and second data; and means coupled to said location identifying means for determining a magnitude of the damage at each said location based upon each said identified location and upon said first and second data.

Claim 90

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer operable to measure a first signal set in response to a first excitation force, said transducer further operable to measure a second signal set in response to a second excitation force; a converter operable to convert said first signal set to first data indicative of pre-damage modal parameters of the structure, said pre-damage modal parameters including at least one measured mode shape, said converter further operable to convert said second signal; set to second data indicative of post-damage modal parameters of the structure, said post-damage modal parameters including at least one measured mode shape; means coupled to said converter for identifying at least one location of the damage based upon said first and second data; and means coupled to said location identifying means for determining a magnitude of the damage at each said location based upon each said identified location and upon said first and second data.

Claim 91

The apparatus of claim 83, wherein said location identifying means comprises means for applying said first and second data to an indicator function.

Claim 92

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer coupled to the structure operable to acquire a pre-damage vibration frequency response for the structure; in response to a first excitation force, said transducer further operable to determine a post-damage vibration frequency response for the structure as damaged; a converter operable to convert said pre-damage vibration frequency response to first data indicative of pre-damage modal parameters of the structure, said converter further operable to convert said post-damage vibration frequency response to second data indicative of post-damage modal parameters of the structure; means for generating a stiffness matrix of the structure based on one or more material and geometric properties of the structure; means coupled to said converter for extracting from the pre-damage and post-damage frequency response data a plurality of N×1 matrices of modal parameters for one or more vibrational modes, the parameters including fractional changes in damped eigenvalue frequencies Z d , fractional changes in modal masses Z m , fractional changes in modal damping ratios Z damp , and fractional changes in the norm of mode shapes Z modes , N being selected as a number of said vibrational modes; means coupled to said stiffness matrix generating means and to said converter for computing an N×B pre-damage sensitivity matrix F 0 using the stiffness matrix and the pre-damage frequency response data, wherein said sensitivity matrix relates changes in said modal parameters to a fractional change in stiffness of one or more elements of the structure, B being selected as a number of said elements; means coupled to said stiffness matrix generating means and to said converter for computing an N×B sensitivity change matrix ΔF using the stiffness matrix and the post-damage frequency response data, wherein said sensitivity change matrix accounts for changes in the sensitivity matrix F 0 resulting from damage to the structure; and means coupled to said extracting means, to said sensitivity matrix computing means and to said sensitivity change matrix computing means for using a damage detection equation: Z.sub.d +Z.sub.m +Z.sub.damp +Z.sub.modes =(F.sub.0 +ΔF)α to solve for α, a B×1 matrix containing the fractional change in stiffness of each said element.

Claim 93

Apparatus for detecting damage in a structure, comprising: an exciter operable to provide an excitation force to the structure; a transducer coupled to the structure operable to acquire a pre-damage vibration frequency response for the structure; in response to a first excitation force, said transducer further operable to determine a post-damage vibration frequency response for the structure as damaged; a converter operable to convert said pre-damage vibration frequency response to first data indicative of pre-damage modal parameters of the structure, said converter further operable to convert said post-damage vibration frequency response to second data indicative of post-damage modal parameters of the structure; means for generating a stiffness matrix of the structure based on one or more material and geometric properties of the structure; means coupled to said converter for extracting from the pre-damage and post-damage frequency response data a plurality of N×1 matrices of modal parameters for one or more vibrational modes, the parameters including fractional changes in damped eigenvalue frequencies Z d , fractional changes in modal masses Z m , fractional changes in modal damping ratios Z damp , and fractional changes in the norm of mode shapes Z modes , N being selected as a number of said vibrational modes; means coupled to said stiffness matrix generating means and to said converter for computing an N×B pre-damage sensitivity matrix F 0 using the stiffness matrix and the pre-damage frequency response data, wherein said sensitivity matrix relates changes in said modal parameters to a fractional change in stiffness of one or more elements of the structure, B being selected as a number of said elements; means coupled to said stiffness matrix generating means and to said converter for computing an N×B sensitivity change matrix ΔF using the stiffness matrix and the post-damage frequency response data, wherein said sensitivity change matrix accounts for changes in the sensitivity matrix F 0 resulting from damage to the structure; and means coupled to said extracting means, to said sensitivity matrix computing means and to said sensitivity change matrix computing means for using a damage detection equation: Z.sub.d +Z.sub.m +Z.sub.damp +Z.sub.modes =(F.sub.0 +ΔF)α to solve for α, a B×1 matrix containing the fractional change in stiffness of each said element.

Claim 94

The apparatus fo claim 85, wherein the number of elements B exceeds the number of vibrational modes N, the apparatus further comprising: means coupled to said converter for identifying at least one location of the damage using an indicator function; means coupled to said identifying means and to said damage detection equation means for eliminating elements of the damage detection equation based upon each said location to form a reduced set damage detection equation, wherein said damage detection equation means comprises means for solving for an element of α corresponding to each said location using the reduced set damage detection equation.

Claim 95

The apparatus of claim 85, wherein said damage detection equation means comprises means for solving for said element of α by dividing an element of said quantity F 0 +ΔF into an element of at least one said matrix of modal parameters, wherein said element of said matrix F 0 +ΔF corresponds to a particular mode and to a particular location, wherein said element of each said matrix of modal parameters corresponds to said particular mode, and wherein said element of α corresponds said particular location.

Claim 96

The apparatus of claim 85, wherein the number of elements B exceeds the number of vibrational modes N, said damage detection equation means further comprising: means for inverting the matrix F 0 +ΔF using one or more pseudo inverse techniques; and means coupled to said inverting means for solving for α using the damage detection equation.

Claim 97

The apparatus of claim 85, wherein the number of vibrational modes N exceeds the number of elements B, said damage detection equation means further comprising: means for inverting the matrix F 0 +ΔF using one or more exact inverses; and means coupled to said inverting means for solving for α using the damage detection equation.

Claim 98

The apparatus of claim 85, wherein the number of vibrational modes N exceeds the number of elements B, said damage detection equation means further comprising: means for inverting the matrix F 0 +ΔF using at least one over-determined least squares approach; and means coupled to said inverting means for solving for α using the damage detection equation.

Claim 99

The apparatus of claim 90, wherein said inverting means comprises means for performing said least squares approach according to the expression: F.sup.-1 =F.sup.T (FF.sup.T).sup.-1.

Claim 100

The apparatus of claim 85, wherein said damage detection equation means comprises means for eliminating at least one element of the matrix F 0 +ΔF if, after inverting the matrix F 0 +ΔF, a corresponding element of α has a value outside a predetermined range.

Claim 101

The apparatus of claim 92, wherein said damage detection equation means comprises means coupled to said eliminating means for inverting the matrix F 0 +ΔF after each said element of the quantity matrix F 0 +ΔF is eliminated.

Claim 102

The apparatus of claim 85 and further comprising means for modelling each of the B elements as a spring with element stiffness matrix given by ##EQU33## wherein i is an element number, E is Young's modulus of elasticity, I is a second moment of area, and L is a length of element number i.