AN APPROACH TO DAMAGE DETECTION IN A BEAM LIKE STRUCTURES FROM EXPERIMENTALLY MEASURED FRF DATA

Size: px
Start display at page:

Download "AN APPROACH TO DAMAGE DETECTION IN A BEAM LIKE STRUCTURES FROM EXPERIMENTALLY MEASURED FRF DATA"

Transcription

1 A APPROACH TO DAMAGE DETECTIO I A BEAM LIKE STRUCTURES FROM EXPERIMETALLY MEASURED FRF DATA Valentina GOLUBOVIĆ-BUGARSKI, Drago BLAGOJEVIĆ 2, 2 University of Bana Luka, Faculty of Mechanical Engineering, BOSIA & HERZEGOVIA ABSTRACT: The metho base on amage etection an relative quantification inicator is investigate in this paper, in orer to etect, locate an quantify the amage of the beam like structure. This metho uses frequency response functions as the characteristics of ynamics response of the mechanical system, from which the amage etection an relative quantification inicators are calculate. The goal of investigation is to etermine effectiveness of the metho base on amage etection an relative quantification inicator to etect, locate an quantify amage of the beam, from the stanpoint of inepenency of the metho to the previously built FEM or analytical moel an for ifferent bounary conitions. KEY WORDS: Damage etection, experimental FRF ata, moal testing.. ITRODUCTIO Structural amage can be efine as changes introuce into a system that aversely affect the current or future performance of that system. If one is focuse on the stuy of amage ientification in structural an mechanical systems, than amage may also be efine as any eviation in the structure s original geometric or material properties that may cause unesirable stresses, isplacements or vibrations on the structure. These weaknesses an eviations may be ue to: cracks, loose bolts, broken wels, corrosion, fatigue, etc. All of them shoul cause a ecrease in the structure s stiffness, an some will also affect its mass an amping properties. Therefore, structural amages shoul always, at a sufficient level of severity, cause a change in a structure s vibration behaviour, escribe by moal properties: natural frequencies, amping loss factor an moe shapes. Since the changes on the ynamic characteristics can be measure an stuie, it is possible to trace what structural changes have cause the ynamic characteristic to change, thus ientifying amage. It is obvious that concept of amage etection is base on comparison between two ifferent states of the system, one of which is assume to represent the initial an often unamage state. The structural amage etection can be ivie into four levels, base on the amount of information provie regaring the amage state. Accoring to Rytter [7] these four levels are: ) ientification of amage existence in a structure; 2) ientification of amage existence in a structure an location of amage; 3) ientification of amage existence in a structure, location of amage an quantification of amage severity; 4) ientification of amage existence in a structure, location of amage an quantification of amage severity an preiction of the remaining service life of the structure. The interest in the ability to monitor a structure an etect amage at the earliest possible stage is pervasive throughout the civil, mechanical an aerospace engineering communities uring last three ecaes. At present, using non-estructive techniques to etect amage status of engineering structures is put into a focus of engineering interests. on-estructive techniques are the means by which structures may be inspecte without isruption or impairment of serviceability. Generally, structural amage etection methos can be classifie as local amage etection methos an global amage etection methos []. Local amage etection techniques refer to non-estructive testing such as acoustic an ultrasonic methos, magnetic fiel methos, raiography, ey-current methos, etc., is mainly use to etect local amage in structure. All these methos require that the vicinity of the amage is known an that the portion of the structure being inspecte is accessible, an can etect amage on or near the surface of the structure. Local amage etection methos utilise only ata obtaine from the amage structure. Baseline ata an theoretical moels of the unamage structure are not use. These are the main avantages of local amage etection, which is very effective for small an regular structures, such as pressure vessels. 296 Tome VIII (year 200), Fascicule 3, (ISS )

2 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) However, for the large an complex structures, for the structures in invisible or close environment, it is very ifficult to etect amage using local amage etection methos, because local amage etection methoology can only be use to inspect some special an accessible components of a structure. In orer to etect amage throughout the whole structure, especially some large, complicate structure, a methoology calle global vibration base structural amage etection has been evelope. Its basic principle is explaine previously, an in a nutshell, when some amage emerges in the structure then structural parameters will change an frequency response functions (FRF) an moal parameters of the structural system will also change. So, the change of the structural FRF s an moal parameters can be taken as the signal of early amage occurrence in the structural system. Such an approach has been introuce for several years in fiels like automotive, aeronautical, mechanical an civil engineering. 2. THE BRIEF SURVEY OF VIBRATIO-BASED DAMAGE DETECTIO APPROACH Structural amage etection using measure vibration ata has been a subect of much practical interests an research efforts an receive consierable attention in the literature. First very etaile survey of the technical literature to etermine the state of the art of the amage-etection fiel using moal analysis proceures was presente by Richarson [6]. The survey focuse on structural integrity monitoring for nuclear power plants, large structures, rotating machinery, an offshore platforms, with by far the largest amount of literature associate with rotating machinery. The author state that while monitoring of overall vibration levels for rotating machinery ha become commonplace, attempts at relating structural amage to measure moal changes was still in its primitive stages. ext survey of the structural health monitoring stuies that have appeare in the technical literature until 996 was presente by Doebling, et al []. That report first categorizes the methos accoring to require measure ata an analysis technique. The analysis categories inclue changes in moal frequencies, changes in measure moe shapes (an their erivatives), an changes in measure flexibility coefficients. Methos that use property (stiffness, mass, amping) matrix upating, etection of nonlinear response, an amage etection via neural networks are also summarize. The applications of the various methos to ifferent types of engineering problems are categorize by type of structure an are summarize. The types of structures inclue beams, trusses, plates, shells, briges, offshore platforms, other large civil structures, aerospace structures, an composite structures. The report escribes the evelopment of the amage-ientification methos an applications an summarizes the current state-of-the-art of the technology. Sohn, et al [9] presents etaile report which is an upate version of the previous literature review report by Doebling, et al. That report contains new technical evelopments publishe between 996. an 200. in the iscipline of structural health monitoring. The authors have organize reviewe articles following the statistical pattern recognition paraigm reporte in [2]. This paraigm can be escribe as a four-part process: () Operational Evaluation, (2) Data Acquisition, Fusion, an Cleansing, (3) Feature Extraction an Information Conensation, an (4) Statistical Moel Development for Feature Discrimination. The reviewe articles are then categorize by type of applications, which inclue beams, truss, plates, briges, aerospace structures, an composite structures. Yan, et al. [] presents general summary an review of state-of-the-art an evelopment of vibration-base structural amage etection methos base on structural ynamics characteristic parameters. They ivie vibration-base structural amage etection methos into traitional-an moern type. The traitional type refers to etection methos for structural amage only utilizing itself mechanical characteristic of structures, such as natural frequencies, moal amping, moal strain energy or moal shapes, etc. This kin of metho generally requires experimental moal analysis or transfer function measure, an the authors fin it not convenient for online etection of structures in service, because these experimental measures often nee multifarious instrument or manual operation. The moern type refers to etection methos for structural amage base on online measure response signal of structures in service. Among the moern type methos for structural amage etection, the representative ones inclue Wavelet analysis, Genetic algorithm an Artificial eural etwork, etc. The cite authors pointe out the main topics for the future research in amage etection base on measure vibration response of a structure: Minimizing the epenence on prior analytical moels an/or prior test ata for the etection an location of amage. Developing of methos that have ability to account for the effects of nonlinear structural response. copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA 297

3 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) Optimization of the position an number of measurement sensors. Construction an extraction of feature inex for small structural amage from structural vibration response. Ability of iscrimination of changes in the moal properties resulting from amage from those resulting from variations in the measurements (resulting from changing environmental an/or test conitions an from the repeatability of the test) Reucing the epenence upon measurable excitation forces, using vibrations inuce by ambient environmental or operating loas for the assessment of structural integrity. Comparison of ifferent amage etection methos by application to a common set of ata. Integration of local nonestructive inspection an global vibration testing. Damage prognosis which estimates the remaining service life of a structure. 3. CLASIFICATIO OF DAMAGE IDETIFICATIO METHODS In orer to etect structural amage from structural ynamic response signal, the first problem is to select amage feature inex to be constructe. The physical variable use to ientify amage may be a global one, but the physical variable use to etermine amage location is better to be local one an must be sensitive to structural local-amage. Determination of structural amage location is equivalent to etermining a region where structural stiffness an loaing capacity ecreases using a measurable quantity. The key factor of vibration base amage etection is to establish the calculation moel an to estimate the vibration parameter to be measure. Doebling [] an Farrar [2] categorize amage ientification methos accoring to: the type of measure ata (feature) use the techniques use to ientify the amage from the measure ata. Common features use in vibration base amage etection stuies are: ) moal frequencies, 2) frequency response functions, 3) moe shapes 4) moe shape curvatures 5) moal strain energy 6) ynamic flexibility, etc. The techniques use to ientify the amage from the measure ata can be classifie as: ) methos base on frequency changes (forwar or inverse problem methos) 2) methos base on moe shape changes 4) methos base on moe shape curvature changes 3) methos base on ynamically measure flexibility: comparison of flexibility changes, stiffness error matrix metho, effects of resiual flexibility, changes in measure stiffness matrix, 4) matrix upate methos, 5) neural-network base methos, 6) time-history an spectral pattern methos, 7) nonlinear methos, 8) statistical pattern recognition methos, etc. He [4] classifies amage etection methos as: Damage etection using experimental ata only. Damage etection without a reaily available spatial moel of the unamage structure is an approach eparte completely from moel upating path. Usually, the ata available are the experimental ata before an after amage occurre. As a result, we are ealing with two sets of moal or FRF ata. The comparison of these two sets shoul yiel the information about the existence an location of amage. The main question is how to relate the ifferences between moal an FRF ata before an after amage to the spatial stiffness changes that resulte in the ifferences. Damage etection using moal ata an analytical ata. Damage etection using moal ata is an approach that was largely aopte from moel upating. Its algorithm aims to etermine amage by using the moal ata from a amage structure an an analytical moel for its counterpart. Although many moel upating methos utilizing moal ata can be aopte irectly for amage etection, metho base on moal force vector is more convenient an efficient. Damage etection using FRF ata. Using measure FRF ata for amage etection has many avantages over the traitional methos using moal analysis ata: ) any numerical errors inherent in moal analysis results cause by inaccurate curve fitting an unavailable resiual terms are avoie; 2) no more efforts is neee to process FRF ata in orer to erive moal ata; 3) the most significant avantages of using measure FRF ata over erive moal analysis ata lies in the fact that FRF ata provie abunant information on the ynamic behavior of a structure. Moal analysis ata lose much of the information that FRF ata have, ue to the necessary numerical process to extract them. 4. OE APPROACH TO DAMAGE DETECTIO USIG FRF DATA Sampaio an Maia [8] present some new evelopment of the Detection an Relative amage Quantification inicator, concerning the etection, localization an the relative severity of amage. 298 copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA

4 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) This metho belongs to the class of methos that use the change in the frequency response functions to etect, locate an relatively quantify the amage. The main avantages of the metho are: ) it is not necessary to perform moal ientification; 2) there is no nee for any analytical or numerical moel of the structure; 3) it uses all measure ata in the form of frequency response functions, without further treatment. Theoretical escription: The equation of motion of a multiple-egree of freeom system with hysteretical amping, which is often use in escribing of complex structure s ynamics [5], is: [ M ]{ x& () t } + i[ D ]{ x ( t )} + [ K ]{ x ( t )} { f ( t )} () If the excitation is harmonic, the realtion between the response an the excitation at each frequency of the analysis is given by: { X} [ α ]{ F} (2) where 2 [ α ( ω )] ([ K + id ] ω [ M ]) (3) is the system receptance matrix, containing all the inforamtion about the ynamic characteristic of the system. Each element α k (ω) of the matrix correspons to an iniviual FRF escribing the relation between the response at a particular coorinate an a single force excitation applie at coorinate k: X α k (ω) F i 0, i... ; i k (4) Fk The column vector, k, of the receptance matrix, { α k }, escribes the shape (in space) exhibite by the structure at each excitation frequency ω, given by the responses normalize by the applie forces. When a structure is amage its stiffness an amping change an, in consequence, so oes the receptance matrix: 2 [ α ( ω )] ([ K + i D ] ω [ M ]) (5) where the superscript stans for amage. It is reasonable to assume that the smaller the egree of correlation between the column vectors,{ α k (ω)} an { α k (ω)}, the larger the amage. To measure the egree of correlation between two vectors, R. Pascual et al. [5,6] propose the Frequency Domain Assurance Criterium (FDAC): where FDAC 2 α k ( ω2 ) α k ( ω ) ( ) k _ ω, ω2 (6) [ α k ( ω2 ) α k ( ω2 )] α k ( ω ) α k ( ω ) [ ] represents the conugate operator an is the total number of co-orinates or measuring points. W. Heylen et al. [7] efine a simplifie form of FDAC, known as RVAC or Response Vector Assurance Criterion (RVAC) (see Figure ), with only one applie force (so that the receptance matrix turns to be ust a vector): ( ω ) α ( ω ) α RVAC _ ( ω ) (7) [ α ( ω ) α ( ω )] α ( ω ) α ( ω ) [ ] The Detection an Relative amage Quantification inicator is formulate as [8]: RVAC ω DRQ (8) ω ω where is the number of frequencies an, so, DRQ will vary between 0 an. 2 copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA 299

5 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) The proceure escribe previously is the strategy for amage etection that is level accoring to Rytter s classification of the structural amage etection. The similar proceure can be further applie to the localization stage, that is level 2 accoring to Rytter. A comparison between shape vectors can be taken as before, but now taking consecutive pairs of coorinates (p an p+) along the structure, instea of oing it for all the coorinates simultaneously (whole structure). The calculation is entirely similar an the local RVAC is efine as Figure. Graphical presentation of the RVAC [8] RVAC_ ( p, ) p + α p p + p ( ω ) α ( ω ) p + [ α ( ω ) α ( ω )] α ( ω ) α ( ω ) ω (9) p 2 [ ] For the sake of comparison with the propose metho, one shall also compute the Generalize Damage Inex, in terms of shape vectors, efine as: β i α i + α i α i ( i i α + i α i α i i i The evelope inicators can also give information about the relative extent of the amage, i.e., it can be use as relative amage quantifiers, which satisfy level 3 of structural amage etection accoring to Rytter. ω) 5. EXPERIMETAL DAMAGE DETECTIO O THE BEAM LIKE STRUCTURE Sampaio an Maia state in [8] that there is no nee for any analytical or numerical moel for structural amage etection base on the Detection an Relative Damage Quantification inicator. They showe numerical performance of the metho an its experimental application on the example of the free-free supporte beam, exciting it by shaker an collecting responses with 2 uniformly istribute accelerometers along the beam. (0) Figure 2. Two bounary conitions: the freely supporte beam an the cantilever beam However, it was intereste to investigate how propose metho operates in the case of impulse excitation an the limite number of accelerometers on isposal. Thus, we ecie to experimentally investigate performance of the propose metho applying it to the freely supporte beam an cantilever beam, Figure 2. Moal testing was performe by means of hammer excitation, using two 300 copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA

6 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) ifferent testing methos: roving hammer metho an roving accelerometers metho. It was necessary to collect FRFs for chosen number of DOFs from the unamage beam an from the amage beam for a few levels of amage. The amage was simulate introucing a cut at a certain location of the beam, an the amage was propagate by eepening the cut for mm at every level of amage. 5.. An equipment use for moal testing Measurement ata were collecte using the Multi-channel Data Acquisition Unit Portable Pulse, Bruel&Kaer type 3560 C, an analyze in the Pulse LabShop 9.0 software, in the frequency range of Hz, Figure 3. An impact hammer Enevco type generates excitation, while the responses were capture by moal accelerometer, B&K type 4507, attache to the structure. Both signals were weighte by some winow functions: the excitation signal by transient winow function an response signal by exponential winow function. Frequency resolution was chosen to be Hz, an the number of averaging was 5 per DOFs. Figure 3. Equipment use for moal testing 5.2. Damage etection of the freely supporte beam The steel beam of imensions mm was use for experimental investigation. The crack of mm with was introuce by wire-cut. The beam was suspene with common strings to simulate free-free conitions. Two test were one with two ifferent moal testing proceure: test : so calle roving hammer test with fixe response measurement position test 2: so calle roving accelerometer test with fixe hammer excitation position. In the roving hammer test the accelerometer was attache to noe 5 to capture the vibration response signals, while the impact hammer generates excitation on the each of 7 noes uniformly arrange along the beam, Figure 4. In the roving accelerometer test the excitation position was fixe at no 8, but the accelerometer was move along the 7 DOF of beam to capture the vibration response signals, Figure 5. To avoi the loaing mass effect ue to local ae transucer mass, 6 ummy transucers were use uring test, Figure 6a an 6b. The ummies were systematically replace by the real transucer for the measurement, until all of the measurement locations have been covere. Thus the systematic error incurre as a result of moving the transucer aroun on the structure uring test was turne into single error by the use of ummy transucer. The beam was cut with 7 epths in the mile of location 4 an 5, Table. As one can see from table the beam was measure in 0 conitions: the unamage (or reference), two unamage but ifferent from reference state, an 7 levels of amage at certain location. Table. Levels of amage of the free-free beam Level of amage Depth of cut [mm] Figure 4. Location of amage an response measurement location for the roving hammer test Figure 5. Location of amage an excitation position for the roving accelerometer test copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA 30

7 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) Figure 6. Dummy transucers use in the roving accelerometer test an the free-free beam equippe by ummy transucers [(m/s²)/] Frequency Response H(ogovor,pobua) - File (Magnitue) After 0 moal tests of the free-free beam were one, all FRFs measure at location 5 were OVERLOAD D:\Valentina\DOKTORAT\G\G bo\frequency Response H(ogovor,pobua) - Mark overlappe, 0k Figure 7. It is obvious from Figure 7, that there is some frequency shift ue to increasing of the amage, that is frequencies move to the left (ecreasing) ue to of ecreasing of the stiffness of the k beam (when level of amage increasing). Moal frequencies for 0 levels of amage an two moal tests are 00 shown in Table 2. a k.6k 2k 2.4k 2.8k 3.2k [Hz] [B/.00 (m/s²)/] Frequency Response H(ogovor 2,pobua) - File7 (Magnitue) Valentina\DOKTORAT\G OVERLOAD va acc\g_5a va acc\frequency Response H(ogovor 2,pobu b k.6k 2k 2.4k 2.8k 3.2k [Hz] Figure 7. Overlappe FRFs measure for 0 levels of amage: a. roving hammer test, b. roving accelerometer test Test Test 2 Table 2. Values of moal frequencies for unamage an amage beam Moal frequencies (Hz) Level of amage "" f f f f f f G - inikator DRQ Test :DRQ G inikator DRQ Test 2:DRQ DRQ DRQ Level ivo of amage oštećena Level ivo of oštećena amage Figure 8. Detection an quantification of the amage using the DRQ inicator for two moal tests After calculation the Response Vector Assurance Criterion (eq.7.) an the Detection an Relative amage Quantification inicator (eq.8.), results are graphically interprete as follow on Figure 8. It is obvious that DRQ inicator shows ecreasing tren as increasing the level of amage, for both moal tests. However, roving accelerometer test shows better tenency of ecreasing the DRQ inicator 302 copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA

8 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) value as the level of amage increasing. Generally, it can be state that the DRQ inicator is able to etect an relatively quantify amage, which satisfy the level an the level 3 of amage etection accoring to Rytter. For the purpose to locate amage, much better results showe the General Damage Inex - GDI (eq. 0) than the local RVAC inicator (eq.9). To calculate GDI, some correction has been one to the inex efine with the equation as follows: α ( ) ( ) ( ) ( ) f + α f α ( ) ( f ) ( ) β p, f () α ( ) ( ) + ( ) ( ) ( ) f ( f ) α However, GDI efine from equation was still not enough sensitive to a low level of amage, so that the new inex, name cumulative GDI was propose [3]. The Cumulative GDI was calculate by successive aing the values of GDI for the each level of amage. Some measurement inaccuracies occurre on the certain locations uring testing coul be average, but it is suppose that GDI shoul increase continuously on the location on amage. Figure 9 shows cumulative GDI inicating the location of amage between measurement location 4 an 5, for both moal tests. a.05 G Inikator GDI - Cumulative suma sukcesivnih GDI nivoa oštećena 5 4 b Cumulative G inikator GDI_4 4 5 GDI GDI Location Lokacia (noe (bro čvora) number) Location Lokacia (noe (bro čvora) number) Figure 9. Location of amage element using the GDI inicator: a) roving hammer test, b) roving accelerometer test 5.3. Damage etection of the cantilever beam The steel beam of imensions mm was use for experimental investigation. The beam was clampe at one en, forming bracket of 300 mm length, Figure 0. An impact hammer generates excitation on the each of 4 noes uniformly arrange along the beam. An accelerometer was attache to noe to capture the vibration response signals. The amage was simulate by cut of 0.5 mm with introuce the mile of location 4 an 5, with 5 level of epths, Table 2. As one can see from Table 3 the cantilever beam was measure in 7 conitions: the unamage (or reference), one unamage but ifferent from reference state, an 5 levels of amage at certain location Figure 0. Location of amage an response measurement location of the cantilever beam Table 3. Levels of amage of cantilever beam Level of amage Depth of cut [mm] After 7 moal tests of the cantilever beam were one, all FRFs measure at location were overlappe, Figure. It is obvious from figure, that there is some frequency shift ue to increasing of the amage, that is frequencies move to the left (ecreasing) ue to of ecreasing of the stiffness of the beam (when level of amage increasing). Moal frequencies for 7 levels of amage are shown in Table 4. copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA 303

9 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) Table 4. Values of moal frequencie for unamage an amage beam Moal frequencies (Hz) Level of amage "" Moe: f Moe 2: f Moe 3: f Moe 4: f [(m/s²)/] Frequency Response H(ogovor,pobua) -50% (Magnitue) OVERLOAD D:\Valentina\DOKTORAT\U k3\k3 50%\Frequency Response H(ogovor,pobua) - Mar k m k.6k 2k 2.4k 2.8k 3.2k [Hz] Figure. Overlappe FRFs measure for 7 levels of amage DRQ GDI K3 inikator DRQ DRQ Damage ivo level oštećena Figure 2. Detection an quantification of amage using the DRQ inicator K3 suma sukcesivnih GDI nivoa oštećena Location (noe Lokacia number) Figure 3. Location of amage element using GDI inicator After calculation the Detection an Relative amage Quantification inicator, results are graphically interprete as follow on Figure 2. It is obvious that DRQ inicator shows ecreasing tren with increasing the level of amage. Therefore, the DRQ inicator is able to etect an relatively quantify amage. As previously mentione for the freely supporte beam, to locate the amage much better results showe the General Damage Inex (GDI) than the local RVAC inicator. The Cumulative GDI was calculate by successive aing the values of GDI for the each level of amage. Figure 3 shows cumulative GDI inicating the location of amage between measurement location 4 an 5 of cantilever beam. 6. COCLUSIOS This paper presents one approach in amage etection using frequency response functions ata. It is point out that using measure FRF ata for amage etection has many avantages over the traitional methos using moal analysis ata, especially that FRF ata provie abunant information on the ynamic behavior of a structure. By reviewing accomplishment in amage etection methos base on vibration responses, one of the observe tasks for future evelopment an avancing of amage etection methos is minimizing the epenence on prior analytical moels an/or prior test ata for the etection an location of amage. So, the metho base on Detection an Relative amage Quantification inicator is investigate an approve. Experimental investigation was conucte on the free-free beam an the cantilever beam, using hammer excitation an roving hammer metho of moal testing. The results of experiments show that DRQ inicator is able to etect an relatively quantify the amage that is to recognize the pattern of amage variation. It was foun that Generalize Damage Inex gave better location of amage than local RVAC inicator. As certain measurement inaccuracies happene on the certain locations uring testing it coul be average by successive aing the values of GDI for the each level of amage. However, it is suppose that GDI shoul increase continuously on the location on amage. So, some improvement of the GDI inicator is propose that is the Cumulative GDI. Propose amage etection metho showe 304 copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA

10 AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG. Tome VIII (Year 200). Fascicule 3 (ISS ) goo performance even for the hammer excitation an one response transucer available, which is important consiering the practical implementation of the metho in the frugally equippe laboratories. REFERECES [.] Doebling S.W., Farrar C.R., Prime M.B., Shevitz DW., Damage Ientification an Health Monitoring of Structural an Mechanical Systems from Changes in their Vibration Characteristics: A Literature Review, Los Alamos ational Laboratory report LA-3070-MS.,996. [2.] Farrar C.R., Doebling S.W., Vibration-Base Damage Detection, Los Alamos ational Laboratory Report, 999. [3.] Golubović-Bugarski V., Moels of correlation between structural amages an ynamics response of mechanical structure, Ph.D.Dissertation, University of Bana Luka, B&H, 200. [4.] He J., Damage Detection an Evaluation I, Moal Analysis an Testing, ATO Science Series, 999. [5.] Maia,.M.M, Silva, J.M.M. et al., Theoretical an Experimental moal Analysis, Research stuies Press Lt., 997. [6.] Richarson M.H., Detection of Damage in Structures from Changes in their Dynamic (Moal) Properties- A survey, UREG/CR-43, U.S. uclear Regulatory Commission, Washington, D.C., 980. [7.] Rytter A., Vibration base inspection of civil engineering structures, Ph.D.Dissertation, Department of Builing Technology an Structural Engineering, Aalborg University, Denmark, 993. [8.] Sampaio R.P.C., Maia.M.M., Strategies for an efficient inicator of structural amage, Mechanical system an Signal processing, Vol.22, [9.] Sohn H., Farrar C.R., A Review of Structural Health Monitoring Literature: , Los Alamos ational Laboratory report LA-3976-MS, [0.] Farrar C.R., Doebling S.W., Damage Detection an Evaluation II, Moal Analysis an Testing, ATO Science Series, 999. [.] Yan Y.J., Cheng L., Wu Y.Y., Yam L.H., Development in vibration-base structural amage etection technique, Mechanical system an signal processing 2, , AALS OF FACULTY EGIEERIG HUEDOARA ITERATIOAL JOURAL OF EGIEERIG copyright University Politehnica Timisoara, Faculty of Engineering Huneoara, 5, Revolutiei, 3328, Huneoara, ROMAIA copyright FACULTY of EGIEERIG - HUEDOARA, ROMAIA 305

Damage detection of shear building structure based on FRF response variation

Damage detection of shear building structure based on FRF response variation , pp.18-5 http://x.oi.org/10.1457/astl.013.3.05 Damage etection of shear builing structure base on FRF response variation Hee-Chang Eun 1,*, Su-Yong Park 1, Rae-Jung im 1 1 Dept. of Architectural Engineering,

More information

Damage identification based on incomplete modal data and constrained nonlinear multivariable function

Damage identification based on incomplete modal data and constrained nonlinear multivariable function Journal of Physics: Conference Series PAPER OPEN ACCESS Damage ientification base on incomplete moal ata an constraine nonlinear multivariable function To cite this article: S S Kourehli 215 J. Phys.:

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

Optimum design of tuned mass damper systems for seismic structures

Optimum design of tuned mass damper systems for seismic structures Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty

More information

Placement and tuning of resonance dampers on footbridges

Placement and tuning of resonance dampers on footbridges Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version

More information

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints EUROSTEEL 2, August 3 September 2, 2, Buapest, Hungary A NUMERICAL EVALUATION OF CHS T JOINTS UNDER AXIAL LOADS Raphael S. a Silva a, Luciano R. O. e Lima b, Pero C. G. a S. Vellasco b, José G. S. a Silva

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1

DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

More information

VIBRATION CONTROL AND FULL-SCALE MEASUREMENT OF A STEEL TV TOWER WITH A DAMPER DEVICE OF PTTMD

VIBRATION CONTROL AND FULL-SCALE MEASUREMENT OF A STEEL TV TOWER WITH A DAMPER DEVICE OF PTTMD 13 th Worl Conference on Earthquake Engineering Vancouver, B.C., Canaa August 1-6, 24 Paper No. 1439 VIBRATION CONTROL AND FULL-SCALE MEASUREMENT OF A STEEL TV TOWER WITH A DAMPER DEVICE OF PTTMD Renle

More information

Code_Aster. Detection of the singularities and computation of a card of size of elements

Code_Aster. Detection of the singularities and computation of a card of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size

More information

Monitoring of Mechanical Properties of Gypsum Using the Impulse Excitation Method

Monitoring of Mechanical Properties of Gypsum Using the Impulse Excitation Method Monitoring of Mechanical Properties of Gypsum Using the Impulse Excitation Metho RICHARD ŤOUPEK, TOMÁŠ PLACHÝ, MICHAL POLÁK, PAVEL TESÁREK Department of Mechanics, Faculty of Civil Engineering Czech Technical

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

A new identification method of the supply hole discharge coefficient of gas bearings

A new identification method of the supply hole discharge coefficient of gas bearings Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Active control based on finite element analysis. Li Hongying

Active control based on finite element analysis. Li Hongying 5th International Conference on ucation, Management, Information an Meicine (MIM 5) Active control base on finite element analysis Li Hongying Automotive ngineering institute, Jiangxi University of Technology,

More information

Level Construction of Decision Trees in a Partition-based Framework for Classification

Level Construction of Decision Trees in a Partition-based Framework for Classification Level Construction of Decision Trees in a Partition-base Framework for Classification Y.Y. Yao, Y. Zhao an J.T. Yao Department of Computer Science, University of Regina Regina, Saskatchewan, Canaa S4S

More information

EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL

EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL 386 EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL Lien-Kwei CHIEN 1, Yan-Nam OH 2 An Chih-Hsin CHANG 3 SUMMARY In this stuy, the fille material in Yun-Lin

More information

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor

More information

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,

More information

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić Sensors an Actuators A: Physical 168 (211) 149-154 oi: 1.116/j.sna.211..7 211 Elsevier B.V. CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS Jože Kutin *, Ivan Bajsić Laboratory of

More information

Design and Application of Fault Current Limiter in Iran Power System Utility

Design and Application of Fault Current Limiter in Iran Power System Utility Australian Journal of Basic an Applie Sciences, 7(): 76-8, 13 ISSN 1991-8178 Design an Application of Fault Current Limiter in Iran Power System Utility M. Najafi, M. Hoseynpoor Department of Electrical

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL 6th European Signal Processing Conference EUSIPCO 28, Lausanne, Switzerlan, August 25-29, 28, copyright by EURASIP A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL Leonaro Tomazeli

More information

Code_Aster. Detection of the singularities and calculation of a map of size of elements

Code_Aster. Detection of the singularities and calculation of a map of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of

More information

Modelling dowel action of discrete reinforcing bars in cracked concrete structures

Modelling dowel action of discrete reinforcing bars in cracked concrete structures Title Moelling owel action of iscrete reinforcing bars in cracke concrete structures Author(s) Kwan, AKH; Ng, PL; Lam, JYK Citation The 2n International Symposium on Computational Mechanics an the 12th

More information

An Approach for Design of Multi-element USBL Systems

An Approach for Design of Multi-element USBL Systems An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Final. Report. Finite Element Model Updating and Damage Detection for Bridges Using Vibration Measurements. December 2013

Final. Report. Finite Element Model Updating and Damage Detection for Bridges Using Vibration Measurements. December 2013 University Transportation Research Center - Region 2 Final Report Finite Element Moel Upating an Damage Detection for Briges Using Vibration Measurements Performing Organization: Columbia University December

More information

Essential Considerations for Buckling Analysis

Essential Considerations for Buckling Analysis Worlwie Aerospace Conference an Technology Showcase, Toulouse, France, Sept. 24-26, 2001 Essential Consierations for Buckling Analysis 2001-120 Sang H. Lee MSC.Software Corporation, 2 MacArthur Place,

More information

Non-Conservative Stability Analysis of Hauger Types of Columns with Different Boundary Conditions

Non-Conservative Stability Analysis of Hauger Types of Columns with Different Boundary Conditions Proceeings of the Worl Congress on Engineering 08 Vol II WCE 08, July 4-6, 08, onon, U.K. on-conservative Stability nalysis of Hauger Types of Columns with Different Bounary Conitions S.. Fazelzaeh,. Tashakorian,

More information

Optimal LQR Control of Structures using Linear Modal Model

Optimal LQR Control of Structures using Linear Modal Model Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

Ductility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2

Ductility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2 Ductility an Failure Moes of Single Reinforce Concrete Columns Hiromichi Yoshikawa 1 an Toshiaki Miyagi 2 Key Wors: seismic capacity esign, reinforce concrete column, failure moes, eformational uctility,

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating droplets tracking by holographic high speed video in turbulent flow Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es

More information

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope

Comparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel

More information

Crack-tip stress evaluation of multi-scale Griffith crack subjected to

Crack-tip stress evaluation of multi-scale Griffith crack subjected to Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Real-time arrival prediction models for light rail train systems EDOUARD NAYE

Real-time arrival prediction models for light rail train systems EDOUARD NAYE DEGREE PROJECT IN TRANSPORT AND LOCATION ANALYSIS STOCKHOLM, SWEDEN 14 Real-time arrival preiction moels for light rail train systems EDOUARD NAYE KTH ROYAL INSTITUTE OF TECHNOLOGY SCHOOL OF ARCHITECTURE

More information

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study

Resilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:

More information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive Control of a Laboratory Time Delay Process Experiment Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results

More information

Experiment 2, Physics 2BL

Experiment 2, Physics 2BL Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

International Conference on Advances in Energy, Environment and Chemical Engineering (AEECE-2015)

International Conference on Advances in Energy, Environment and Chemical Engineering (AEECE-2015) International Conference on Avances in Energy, Environment an Chemical Engineering (AEECE-2015) Stuy on Damage Characteristic of Unergroun Cavern Blasting Excavation Base on Dynamic Damage Constitutive

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

INFLUENCE OF SURFACE ROUGHNESS THROUGH A SERIES OF FLOW FACTORS ON THE PERFORMANCE OF A LONGITUDINALLY ROUGH FINITE SLIDER BEARING

INFLUENCE OF SURFACE ROUGHNESS THROUGH A SERIES OF FLOW FACTORS ON THE PERFORMANCE OF A LONGITUDINALLY ROUGH FINITE SLIDER BEARING ANNALS of Faculty Engineering Huneoara International Journal of Engineering Tome XIV [2016] Fascicule 2 [May] ISSN: 1584-2665 [print; online] ISSN: 1584-2673 [CD-Rom; online] a free-accessmultiisciplinarypublication

More information

Lecture 3: Development of the Truss Equations.

Lecture 3: Development of the Truss Equations. 3.1 Derivation of the Stiffness Matrix for a Bar in Local Coorinates. In 3.1 we will perform Steps 1-4 of Logan s FEM. Derive the truss element equations. 1. Set the element type. 2. Select a isplacement

More information

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d

SYNTHESIS ON THE ASSESSMENT OF CHIPS CONTRACTION COEFFICIENT C d SYNTHESIS ON THE ASSESSMENT OF HIPS ONTRATION OEFFIIENT Marius MILEA Technical University Gheorghe Asachi of Iași, Machine Manufacturing an Inustrial Management Faculty, B-ul D Mangeron, nr 59, A, 700050,

More information

Time-Optimal Motion Control of Piezoelectric Actuator: STM Application

Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Yongai Xu, Peter H. Mecl Abstract This paper exaes the problem of time-optimal motion control in the context of Scanning Tunneling

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

Image Based Monitoring for Combustion Systems

Image Based Monitoring for Combustion Systems Image Base onitoring for Combustion Systems J. Chen, ember, IAENG an.-y. Hsu Abstract A novel metho of on-line flame etection in vieo is propose. It aims to early etect the current state of the combustion

More information

Radar Sensor Management for Detection and Tracking

Radar Sensor Management for Detection and Tracking Raar Sensor Management for Detection an Tracking Krüger White, Jason Williams, Peter Hoffensetz Defence Science an Technology Organisation PO Box 00, Einburgh, SA Australia Email: Kruger.White@sto.efence.gov.au,

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Impact Experimental Analysis and Computer Simulation Yucheng Liu Department of Mechanical Engineering, University of Louisville

Impact Experimental Analysis and Computer Simulation Yucheng Liu Department of Mechanical Engineering, University of Louisville Impact Experimental Analysis an Computer Simulation Yucheng Liu Department of Mechanical Engineering, University of Louisville Abstract In this paper, an automotive bumper system (a bumper connecte to

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

Leaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes

Leaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes Leaving Ranomness to Nature: -Dimensional Prouct Coes through the lens of Generalize-LDPC coes Tavor Baharav, Kannan Ramchanran Dept. of Electrical Engineering an Computer Sciences, U.C. Berkeley {tavorb,

More information

Improved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1

Improved Geoid Model for Syria Using the Local Gravimetric and GPS Measurements 1 Improve Geoi Moel for Syria Using the Local Gravimetric an GPS Measurements 1 Ria Al-Masri 2 Abstract "The objective of this paper is to iscuss recent research one by the author to evelop a new improve

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

Optical wire-grid polarizers at oblique angles of incidence

Optical wire-grid polarizers at oblique angles of incidence JOURNAL OF APPLIED PHYSICS VOLUME 93, NUMBER 8 15 APRIL 003 Optical wire-gri polarizers at oblique angles of incience X. J. Yu an H. S. Kwok a) Center for Display Research, Department of Electrical an

More information

arxiv:quant-ph/ v1 29 Jun 2001

arxiv:quant-ph/ v1 29 Jun 2001 Atomic wave packet basis for quantum information Ashok Muthukrishnan an C. R. Strou, Jr. The Institute of Optics, University of Rochester, Rochester, New York 14627 (March 15, 2001) arxiv:quant-ph/0106165

More information

Finite element analysis of electromagnetic bulging of sheet metals

Finite element analysis of electromagnetic bulging of sheet metals International Journal of Scientific & Engineering Research Volume 3, Issue 2, Febraury-212 1 Finite element analysis of electromagnetic bulging of sheet metals Ali M. Abelhafeez, M. M. Nemat-Alla, M. G.

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper

More information

Modelling machine tool dynamics using a distributed parameter tool holder joint interface

Modelling machine tool dynamics using a distributed parameter tool holder joint interface International Journal of Machine Tools & Manufacture 47 (27) 96 928 www.elsevier.com/locate/ijmactool Moelling machine tool ynamics using a istribute parameter tool holer joint interface Keivan Ahmai,

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Sparse Reconstruction of Systems of Ordinary Differential Equations

Sparse Reconstruction of Systems of Ordinary Differential Equations Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA

More information

Part I: Web Structure Mining Chapter 1: Information Retrieval and Web Search

Part I: Web Structure Mining Chapter 1: Information Retrieval and Web Search Part I: Web Structure Mining Chapter : Information Retrieval an Web Search The Web Challenges Crawling the Web Inexing an Keywor Search Evaluating Search Quality Similarity Search The Web Challenges Tim

More information

Adjoint Transient Sensitivity Analysis in Circuit Simulation

Adjoint Transient Sensitivity Analysis in Circuit Simulation Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear

Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Sujecte to In-Plane Shear Husam Al Qalan 1)*, Hasan Katkhua 1) an Hazim Dwairi 1) 1) Assistant Professor, Civil Engineering

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information