DAMAGE LOCALIZATION IN OUTPUT-ONLY SYSTEMS: A FLEXIBILITY BASED APPROACH
|
|
- Felicia Thompson
- 5 years ago
- Views:
Transcription
1 DAAG LOCALIZAION IN OUPU-ONLY SYSS: A LXIBILIY BASD APPROACH Dionisio Bernal an Burcu Gunes Department of Civil an nvironmental nineerin, 47 Snell nineerin Center, Northeastern University, Boston A 0115, U.S.A ABSRAC o assemble flexibility matrices from vibration sinals the input most be measure an that there must be, at least, one co-locate sensor actuator pair. echniques to localize amae that use chanes in flexibility have not, therefore, foun application in the important case of ambient excitation. he central topic of this paper is the iscussion of an approach that allows extension of flexibility base amae localization to cases where only output measurements are available. he funamental iea is that lack of eterministic information on the input can be partially compensate by knowlee of the structure of the mass matrix. In particular, when the inverse of the mass matrix in the coorinates efine by the sensors can be assume iaonal then matrices that iffer from the flexibility by a scalar multiplier can be assemble an these matrices can replace the flexibility with no loss of useful information. NONCLAUR U, D, D DLV nsi WSI K C Ψ Ψ m m n unamae flexibility matrix amae flexibility matrix chane in flexibility matrix amae locatin vector normalize stress inex weihte stress inex mass matrix stiffness matrix ampin matrix isplacement partition of the complex moes complex moes at sensor locations normalization constants number of output sensors number of moes INRODUCION When experimental ata is use to improve a mathematical moel the parameters that are caniates for upatin are typically selecte with uiance from a priori knowlee on what aspects of the moel are most uncertain. In the amae ientification problem, however, uncertainty or sensitivity can not be use as criteria to arrive at a set of free parameters. Since the use of a lare parameter space leas to ill-conitionin an non-uniqueness methos that can provie obective information about amae without the nee to refer to a etaile moel of the structure are of consierable practical importance. Amon the techniques that attempt to localize amae without reference to a moel, those that operate with chanes in moe shapes have receive sinificant attention [1,,3]. A ifficulty often face when usin moe shape chanes to locate amae, however, is the fact that the appropriate pairin of moes from the reference to the potentially amae-state is not always apparent. One way to et aroun this ifficulty is to use all the available moes to assemble flexibility matrices an then focus on the chane of these matrices [4,5]. he fact that the flexibility is ominate by the lower moes (which are typically the ones that can be ientifie experimentally) an that it can be assemble at whatever sensor coorinates are available are convenient features of the approach. An important limitation on the use of flexibility is the fact that these matrices can only be assemble from the ata when the input is measure [6,7]. Since in civil enineerin structures a full characterization of the input is often impractical, the flexibility-base amae localization has not been consiere as a viable option in builins or bries. his paper shows, however, that while it is true that flexibility matrices can
2 not be assemble from vibration ata for output only systems, matrices that iffer from the flexibility by a sinle scalar multiplier can be obtaine when some conitions prevail. Specifically, the paper shows that matrices that are proportional to the flexibility can be compute exclusively from the measure ata if the inverse of the mass matrix can be assume iaonal over the coorinates efine by the available sensors. he obective of this paper is to iscuss the issues associate with the computation of the flexibility proportional matrices an to illustrate how these matrices can be use in amae localization applications. he paper is oranize as follows. he first section presents a summary of a recently evelope technique for interroatin chanes in flexibility about amae localization. he technique, esinate as the Damae Locatin Vector Approach, (DLV) locates the amae by inspectin stress fiels create by vectors that are containe in the null space of the chane in flexibility [8]. he next section illustrates how moal orthoonality can be combine with knowlee of the structure of the inverse of the mass matrix to arrive at matrices that iffer from the flexibility by a sinle scalar multiplier. his evelopment is mae for a viscously ampe system but no limitation on the nature of the ampin (classical or not classical) is introuce. he theoretical part of the paper conclues with a iscussion on how to ensure that the missin scalar in the flexibility matrix, compute for the reference an the amae states, are essentially the same. A numerical example on a 4 DO system illustrates the techniques iscusse. H DLV CHNIQU he Damae Locatin Vector (DLV) approach provies a systematic way for interroatin chanes in flexibility matrices with respect to the localization of the amae. In this section the basic features of the technique are reviewe. A more etaile iscussion of the theoretical backroun, as well as iscussion on robustness an other issues may be foun in [8]. As shall be evient from the results in this section, a missin scalar multiplier in the flexibility matrices is immaterial in the DLV localization, provie the scalar can be mae the same in the reference an the amae states. he basic iea in the DLV approach is that the vectors that span the null-space of the chane in flexibility (from the unamae to the amae states) when treate as static loas on the system, lea to stress fiels that are zero over the amae elements. Depenin on the number an location of the sensors the intersection of the null stress reions ientifie by the DLVs may or may not exclusively contain amae elements. lements that are unamae but which cannot be theoretically iscriminate from the amae ones by chanes in flexibility (for a iven set of sensors) are esinate as inseparable. he steps of the DLV localization can be summarize as follows: 1. Compute the chane in flexibility as; D = U D (1). Obtain a sinular value ecomposition of D, namely; s1 D = U 0 0 s V where s are small sinular values. or ieal conitions s contains zeros an the DLV vectors are simply the columns of V associate with the null space. or the noisy conitions that prevail in practice, however, the values in s are never equal to zero an a cutoff has to be establishe to select the imension of the null space. he vectors in V that can be treate as DLVs for noisy conitions can be selecte as follows (see ref. [8] for the mathematical support) a) Compute the stresses in an unamae moel of the structure usin the columns in V as loas. b) Reuce the internal stresses in every element (for a iven loa vector) to a sinle characterizin stress, σ (strain enery per unit characterizin imension shoul be proportional to σ ). c) Desinate the reciprocal of the maximum value of the characterizin stress as c. Compute the svn inex for every vector in V as; where; q q () s c svn = (3) s c s q c q =max (s c ) for =1: m (4) he vectors for which svn 0.0 can be treate as DLVs. Once the set of DLV vectors has been ientifie the localization proper is carrie out as follows: 3. Compute, for each DLV vector, the normalize stress inex vector as;
3 nsi i σ i = (5) σ i max system has a full set of eienvectors (repeate eienvalues are permitte) the solutions to eq.10 can be oranize as; 4. Compute the vector of weihte stress inices, WSI, as; nlv nsi WSI = i = 1 svn nlv where svn i = max (svn i, ) an nlv is the number of DLV vectors. he potentially amae elements are those havin WSI <1. LXIBILIY ARICS O WIHIN A SCALAR ULIPLIR his section evelops expressions for the inverse of the mass an the stiffness (the flexibility) in terms of the eienvalues of the system an the corresponin eienvectors at the sensor coorinates. It is assume, for enerality, that the ampin is not classical so the moes are irreucibly complex. he expressions obtaine for -1 an K -1 are vali for a particular scalin of the moes which, as shall be shown, can only be enforce to within a missin scalar when the input is not measure. We bein by consierin the homoeneous equation for a time invariant finite-imensional linear system with viscous issipation, namely; i i (6) x& + Cx& + Kx = 0 (7) where, C an K R n x n are the mass, ampin an stiffness matrices. or a state vector efine usin isplacements an velocities the first orer form of eq.7 that preserves symmetry is; where, C =, 0 Assumin a solution eienvalue problem; y& = G y (8) K G = 0 λt 0, x y = (9a,b,c) x & y = φ e eq.8 leas to the φ λ = G φ (10) where Ψ = ΨΛ Ψ Λ, Λ = (11a,b) Ψ Λ Λ Φ Ψ = Φ = ΨΛ 1 1, an Λ = λ [ φ φ L ] λ O an the superscript stans for complex conuate. (1a,b) he symmetry of an G can be exploite to show that; Φ Φ rom where it follows that; G e Φ = (13) e G Φ = (14) = Φ = Φ e e Φ Φ (15) (16) quatin the inverse of the matrices an G in eq.9a an b with the results in eqs.15 an 16 one ets a number of equalities of which we list the ones that are relevant for our purposes: 1 = R ( ΨΛ Λ Ψ 1 ) (17) = R ( Ψ e Λ Ψ ) (18) K 1 = R ( Ψ Ψ ) (19) R ( Ψ Λ Ψ ) = 0 (0) R ( Ψ Ψ ) = 0 (1) rom eqs. 17 an 18 one conclues that; e Λ = e () an thus the previous expressions are reuce to; which, since the matrices are real, yiels real or complex conuate eienvalues λ. Assumin that the = R ( ΨΛ Ψ ) (4)
4 he iaonal matrix 1 K = R ( Ψ Ψ ) (5) R ( Ψ Λ Ψ ) = 0 (6) can be compute from the ata without knowlee of the system matrices in the eterministic input case. In the stochastic case, however, this is not possible an one must introuce some apriori knowlee to procee. Consier as an introuction the case where there is enouh information about the mass matrix to allow the computation of the partition of the inverse over the sensor coorinates. If there are m output sensors then one can use eqs.4 an 6 to set up m(m+1) equations to solve for the complex constants. Assumin all the moes are ientifie one can easily show that the number of sensors require to ientify the normalizin constants has to satisfy m 0. 5( 8n + 1 1) where n is the number of moes (orer is n). or example, if there are 10 moes in the system the normalizin constants can be obtaine from knowlee of the inverse of the mass if m 4. Note that if only a truncate moal space is available then the approach woul lea to some approximation in the constants because one woul be equatin a convere physical quantity (the left sie of eq.4) with a moally truncate approximation. In this paper we purse a less restrictive assumption than that of presumin complete knowlee of -1. In particular, we procee by assumin simply that -1 is iaonal. It is evient, of course, that what we re oin is sayin that we know the off-iaonal part of -1 to be zero. In this case we lose m equations an, because the equations are now all equal to zero we can compute the normalizin factors only up to r unetermine constants, where r is the nullity of the resultin coefficient matrix. It is possible to show that a necessary conition for the nullity to be one (the smallest it can be) is that the number of sensors satisfy m n. So, for n =10, for example, one nees at least 5 sensors. he comment mae previously with respect to equatin a physical quantity (in this case the zeros in the off iaonals of -1 ) to a truncate approximation hols without moification in this case also. he specifics use to arrive at the normalization constants in the case where -1 is assume iaonal are escribe next. Define; H R I m, m, + H i = Ψ Ψ Λ (7) q.4, therefore, can be written in the omain of real numbers as = n = 1 ( H H ) (8) R, R, I, I, where R an I are the real an the imainary components of. akin the upper (lower) trianular portion of -1 (without the main iaonal) an placin it in vector form one can write n 0 = ( Hˆ ˆ R, R, HI, I, = 1 ) (9) where the orer of the entries in the vectors Ĥ is arbitrary as lon as one is consistent in efinin Ĥ R, an Ĥ I,. We now efine the vector β as all the real components of followe by all the imainary components, namely [ R, 1 R, L I, 1 I, L ] β = (30) With the precein notation eq.9 can be written as; where; H β = 0 (31) = [ Ĥ 1,Ĥ,... Ĥ 1, Ĥ... ] (3) H R, R, I, I, It follows then that β (which contains the require normalizin constants orere as per eq.30) is in the null space of the matrix H. One can further restrict the subspace that contains β by takin avantae of the relationship in eq.6. Inee, followin the same approach use to pass from eq.9 to eq.31 one fins that eq.6 ives; S β = 0 (33) where the only ifference is that the assembly of S inclues the iaonal of the matrix in eq.6. Combinin eqs.31 an 33 one ets; H β = Y β = 0 S (34) So β is in the null space of Y. If all the moes are available, the nullity of Y is one an β can be compute to within a sinle scalar. he complex constants (to within a scalar) are then iven by eq.30. an the flexibility proportional matrix is evaluate with eq.5. In practice, of course, one
5 selom obtains all the moes an, as a result, the matrix Y in eq.34 proves to be full rank. An approximate solution can be obtaine by takin β as the sinular vector associate with the smallest sinular value of Y. COPAIBILIY O H SCALAR ULIPLIR he scalar that is missin in the flexibility matrices of the previous section is arbitrary an is not necessarily the same for the unamae an the amae states. In orer to preserve the null space when takin the ifference of the flexibility proportional matrices it is necessary to ensure that the missin constant is consistent. wo proceures have been examine thus far for ensurin compatibility. he first one is base on the iea that the mass matrix has not chane as a result of the amae. If this is true an the moal space is complete then the scalin factor can be auste so that the inverse of the mass in the two states (as iven by eq.8) are the same. Of course, in practice approximations are inevitable an one can not preten that it will be possible to make the two -1 expressions ientical. One can, however, efine a norm an aust the scalin to make this norm the same in the two states. When all the moes are not available the contribution of the available moes to -1 may iffer in the two states an one can not arue that the error in makin the missin scalar compatible erives exclusively from roun off an imprecision. hus far we have looke at two proceures for attainin approximate compatibility between the flexibility proportional matrices of the reference an the amae states. he first approach uses the trace of the -1 as the metric that shoul be equal in the two states an the secon, which oes not use information in -1, is outline next. Assume that the unamae an amae flexibility proportional matrices are U an D while the true moally truncate but, unetermine matrices, are U an D. One can then write; an L U L D = α (37) = η (38) Solvin for L in eq.37 an substitutin the result in eq.38 one ets; U D η = α (39) which shows that the esire η/α ratio is a real eienvalue of the matrix on the left sie. Since there are several real eienvalues, however, the solution oes not immeiately point out which is the correct result. An approach that helps in reucin ambiuity is to first normalize the matrices an in such a way that the η/α ratio is necessarily D U less than 1. his is the case, for example, if the flexibility proportional matrices are normalize to equal trace because the true value of the trace in the amae flexibility is larer (assumin that the effect of amae is more than the inevitable error in the computations). If this scalin is first introuce then one can iscar not only the η/α ratios that are complex but also those that are reater than 1. Unfortunately, there appear to be cases where the solution is still not unique after this is one. While research on the compatibility of the scalar continues, the limite numerical experience aine thus far suests that the simple approach base on the trace of -1 may be sufficiently accurate. NURICAL XAPL he basic steps of amae localization usin the DLV approach for unmeasure input are illustrate usin the simple system shown in fi.1. oal truncation an error in the ientification are not contemplate in this example. an = α (35) U U k 1 = k =1000 k 3 =3000 k 4 =4000 = η (36) D D where α an η are unetermine constants. Assume L is a vector in the effective null space of the moally truncate true chane in flexibility. Reconizin that u L DL one can write (replacin the with = for simplicity); c 1 =40 m 1 m m 3 m 4 c 1 =30 c 1 =0 c 1 =10 k 1 =1500 (after amae) (m 1 =0, m =4, m 3 =8, m 4 =3) iure 1 System consiere
6 Output sensors exist at coorinates 1, an 3 an amae is simulate as 5% reuction in the stiffness of sprin k 1. he first step is to extract arbitrarily scale complex moe shapes Ψ m at the sensor coorinates. his can be one usin any suitable stochastic ientification alorithm [9]. or the system in fi.1 the noise free results in the unamae state are; i i i i i i i i i i i i he secon step is to obtain the normalization constants from eq.34. his requires that one assemble H an S. he computations for the first moe are illustrate in the followin. Assemblin H 1 from eq.7 one ets; i i i H 1 = i i i i i i he upper trianular part of the matrix (without the main iaonal) contains, in this case, 3 numbers. he real an the imainary components of these numbers are use to form the vectors H 1an H 1, namely; Ĥ R,1 = ˆ R, ˆ I, Ĥ I,1 = where the orer selecte is arbitrary. Repeatin these steps for all the moes one obtains the vectors neee to form H (see eq.3). he result is; where the columns of H corresponin to the first moe are hihlihte. he same approach is repeate to form S. One ets; i i i S 1 = i i i i i i Placin the results in vector form ives; Ŝ R,1 = Ŝ I,1 = where we note that the main iaonal is also inclue since it is also zero accorin to eq.6. he matrix S is assemble by repeatin the process for all the moes, the results is; Combinin H an S one obtains the matrix Y. rom inspection of the sinular values, s, it is evient that the nulity is one. he sinular values s an the nullspace β are shown below s = β = rom β one reaily ets the complex constants, from eq.30. he inverse of the mass an stiffness matrices can then be calculate from eqs.4 an 5. One ets; U = α K U = α he exact values for these matrices are; e -1 = K e -1 = x10- from where it can be seen that the ifference is a simple scalar multiplier.
7 ollowin the same steps outline previously the system matrices extracte for the structure in the amae state are; D = α K D = α where the trace of -1 has been normalize to unity in both cases to ensure compatibilty of the missin scalar. Once the flexibility proportional matrices are calculate within a common multiplier, the SVD of the chane in these matrices yiels the sinular values an riht sie sinular vectors shown s= 5.41x10-17 V = x Inspection of the sinular values shows that there are DLVs. Applyin the two vectors ientifie as DLVs as loas at the sensor coorinates one ets the axial force istribution N shown in fi. N N k 1 k k 3 k 4 m 1 m m 3 m 4 DLV1 c 1 c c 3 c DLV iure Axial force istribution 0.50 Both vectors, as expecte, locate the amae correctly at the first sprin. CONCLUSIONS his paper extens the DLV amae localization technique to the important case where the input can not be eterministically characterize. his is one by showin that flexibility proportional matrices can be extracte from the ata without knowlee of the input an that the unetermine missin scalar can be mae consistent for the reference an the amae states. he paper shows that computation of the flexibility proportional matrix is theoretically exact only when the moal basis is complete. xamination of the relationships involve suests, however, that oo approximations of the truncate flexibility shoul result when the truncate space provies a reasonable approximation for -1. he experience erive from numerical experiments has thus far supporte this expectation. RRNCS [1] West, W.., "Illustration of the use of moal assurance criterion to etect structural chanes in an orbiter test specimen." Proc. Air orce Conference on Aircraft Structural Interity, pp.1-6, [] ox, C. H. J., "he location of efects in structures: A comparison of the use of natural frequency an moe shape ata." Proc. of the 10th International oal Analysis Conference, pp.5-58, 199. [3] Salawu, O. S. an Williams, C., "Damae location usin vibration moe shapes." Proc. of the 1 th International oal Analysis Conference, pp , [4] Doeblin, S. C., arrar, C. R., Prime,. B. an Schevitz, D. W., Damae ientification an health monitorin of structural an mechanical systems from chanes in their vibration characteristics: a literature review, Los Alamos, New exico, [5] Paney, A. K. an Biswas,., "Damae etection in structures usin chanes in flexibility." Journal of Soun an Vibration, Vol. 169, No.1, pp.3-17, [6] Bernal, D., "xtractin flexibility matrices from state-space realizations." COS 3 Conference, ari, Spain, Vol.1, pp , 000. [7] Alvin, K.. an Park, K. C., "Secon-orer structural ientification proceure via state-spacebase system ientification", AIAA Journal, Vol.3, No., pp , [8] Bernal, D., "Loa Vectors for Damae Localization", to appear in Journal of nineerin echanics, ASC, January 00. [9] Van Overschee, P., an oor B. L. R., Subspace Ientification for Linear Systems: heory, Implementation, Applications, Kluwer Acaemic Publishers, Boston, 1996.
Applying Axiomatic Design Theory to the Multi-objective Optimization of Disk Brake *
Applyin Aiomatic esin Theory to the Multi-objective Optimization of isk Brake * Zhiqian Wu Xianfu Chen an Junpin Yuan School of Mechanical an Electronical Enineerin East China Jiaoton University Nanchan
More informationUsing Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems
Usin Quasi-Newton Methos to Fin Optimal Solutions to Problematic Kriin Systems Steven Lyster Centre for Computational Geostatistics Department of Civil & Environmental Enineerin University of Alberta Solvin
More information'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 informationDamage 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 informationIncluding the Consumer Function. 1.0 Constant demand as consumer utility function
Incluin the Consumer Function What we i in the previous notes was solve the cost minimization problem. In these notes, we want to (a) see what such a solution means in the context of solvin a social surplus
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationDirect Computation of Generator Internal Dynamic States from Terminal Measurements
Direct Computation of Generator nternal Dynamic States from Terminal Measurements aithianathan enkatasubramanian Rajesh G. Kavasseri School of Electrical En. an Computer Science Dept. of Electrical an
More informationEnergy-preserving affine connections
2 A. D. Lewis Enery-preservin affine connections Anrew D. Lewis 28/07/1997 Abstract A Riemannian affine connection on a Riemannian manifol has the property that is preserves the kinetic enery associate
More informationAn Improvement on Extended Impedance Method towards Efficient Steady-State Analysis of High-Frequency Class-E Resonant Inverters
An Improvement on Extene Impeance Metho towars Efficient Steay-State Analysis of Hih-Frequency Class-E Resonant Inverters Junrui Lian School of Information Science an Technoloy ShanhaiTech University,
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationAnalysis of Halo Implanted MOSFETs
Analysis of alo Implante MOSFETs olin McAnrew an Patrick G Drennan Freescale Semiconuctor, Tempe, AZ, olinmcanrew@freescalecom ABSTAT MOSFETs with heavily ope reions at one or both ens of the channel exhibit
More informationTable 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 informationDamage 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 informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationImplicit 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 informationCode_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 informationarxiv: v1 [math.co] 3 Apr 2019
Reconstructin phyloenetic tree from multipartite quartet system Hiroshi Hirai Yuni Iwamasa April 4, 2019 arxiv:1904.01914v1 [math.co] 3 Apr 2019 Abstract A phyloenetic tree is a raphical representation
More informationDAMAGE 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 informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationStrength 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 informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationAPPROXIMATE 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 informationensembles 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 informationTorque OBJECTIVE INTRODUCTION APPARATUS THEORY
Torque OBJECTIVE To verify the rotational an translational conitions for equilibrium. To etermine the center of ravity of a rii boy (meter stick). To apply the torque concept to the etermination of an
More informationExplicit Formulas for General Integrals of Motion for a Class of Mechanical Systems Subject To Virtual Constraints
Explicit Formulas for General Interals of Motion for a Class of Mechanical Systems Subject To Virtual Constraints Shiriaev, Anton; Perram, John; Robertsson, Aners; Sanber, Aners Publishe in: Decision an
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationTime-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θ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationTHE 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 informationRobustness and Perturbations of Minimal Bases
Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationAPPPHYS 217 Thursday 8 April 2010
APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationQuantum 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 informationAn Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback
Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationOptimal 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 informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationPhysics 251 Results for Matrix Exponentials Spring 2017
Physics 25 Results for Matrix Exponentials Spring 27. Properties of the Matrix Exponential Let A be a real or complex n n matrix. The exponential of A is efine via its Taylor series, e A A n = I + n!,
More informationTEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE
TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,
More informationOptimization 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 informationSTABILITY ESTIMATES FOR SOLUTIONS OF A LINEAR NEUTRAL STOCHASTIC EQUATION
TWMS J. Pure Appl. Math., V.4, N.1, 2013, pp.61-68 STABILITY ESTIMATES FOR SOLUTIONS OF A LINEAR NEUTRAL STOCHASTIC EQUATION IRADA A. DZHALLADOVA 1 Abstract. A linear stochastic functional ifferential
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationRole 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 information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationComputing 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 informationCode_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 informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationA 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 informationAdjoint 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 informationThermal 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 informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationNonlinear 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 informationIntroduction 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 informationPlacement 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 informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationSection 7.1: Integration by Parts
Section 7.1: Integration by Parts 1. Introuction to Integration Techniques Unlike ifferentiation where there are a large number of rules which allow you (in principle) to ifferentiate any function, the
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationA Diagnostic Treatment of Unconstrained Optimization Problems via a Modified Armijo line search technique.
IOSR Journal of Mathematics IOSR-JM e-issn: 2278-5728, p-issn: 2319-765X. Volume 10, Issue 6 Ver. VI Nov - Dec. 2014, PP 47-53 www.iosrjournals.or A Dianostic reatment of Unconstraine Optimization Problems
More informationHarmonic 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 informationUNIFYING PCA AND MULTISCALE APPROACHES TO FAULT DETECTION AND ISOLATION
UNIFYING AND MULISCALE APPROACHES O FAUL DEECION AND ISOLAION Seongkyu Yoon an John F. MacGregor Dept. Chemical Engineering, McMaster University, Hamilton Ontario Canaa L8S 4L7 yoons@mcmaster.ca macgreg@mcmaster.ca
More informationConcentration of Measure Inequalities for Compressive Toeplitz Matrices with Applications to Detection and System Identification
Concentration of Measure Inequalities for Compressive Toeplitz Matrices with Applications to Detection an System Ientification Borhan M Sananaji, Tyrone L Vincent, an Michael B Wakin Abstract In this paper,
More informationLecture 2: Correlated Topic Model
Probabilistic Moels for Unsupervise Learning Spring 203 Lecture 2: Correlate Topic Moel Inference for Correlate Topic Moel Yuan Yuan First of all, let us make some claims about the parameters an variables
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationChapter 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 informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationLQG 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 information3-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 informationParametric Design and Optimization of Cylindrical Gear
Parametric Desin an Optimization of Cylinrical Gear iao Liu, ueyi Li, Collee of Mechanical an Electronic Enineerin, hanon University of cience an echnoloy, Qianwanan Roa, Qinao, hanon, China Abstract Cylinrical
More informationKinematics of Rotations: A Summary
A Kinematics of Rotations: A Summary The purpose of this appenix is to outline proofs of some results in the realm of kinematics of rotations that were invoke in the preceing chapters. Further etails are
More informationHomework 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 informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationFUZZY ARITHMETIC WITHOUT USING THE METHOD OF - CUTS
International Journal o Latest Trens in Computin (E-ISSN: 045-5364) 73 Volume, Issue, December 00 FUZZ ARITHMETIC WITHOUT USING THE METHOD OF - CUTS Supahi Mahanta, Rituparna Chutia, Hemanta K Baruah 3
More informationResilient 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 informationu!i = a T u = 0. Then S satisfies
Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationExperimental 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 informationComparison of Efflux Time between Cylindrical and Conical Tanks Through an Exit Pipe
International Journal of Applie Science an Enineerin 0. 9, : - Comparison of Efflux Time between Cylinrical an Conical Tanks Trou an Exit Pipe C. V. Subbarao * epartment of Cemical Enineerin,MVGR Collee
More informationCharacterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations
Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationa) Identify the kinematical constraint relating motions Y and X. The cable does NOT slip on the pulley. For items (c) & (e-f-g) use
EAMPLE PROBLEM for MEEN 363 SPRING 6 Objectives: a) To erive EOMS of a DOF system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical units (US system) ) To calculate
More information*Agbo, F. I. and # Olowu, O.O. *Department of Production Engineering, University of Benin, Benin City, Nigeria.
REDUCING REDUCIBLE LINEAR ORDINARY DIFFERENTIAL EQUATIONS WITH FUNCTION COEFFICIENTS TO LINEAR ORDINARY DIFFERENTIAL EQUATIONS WITH CONSTANT COEFFICIENTS ABSTRACT *Abo, F. I. an # Olowu, O.O. *Department
More informationLagrangian 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 informationPerfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs
Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite
More information