DYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING. University of Bath. Bath BA2 7AY, U.K. University of Cambridge
|
|
- Whitney Stokes
- 6 years ago
- Views:
Transcription
1 DYNAMIC WEIGHT FUNCTIONS FOR A MOVING CRACK II. SHEAR LOADING A.B. Movchan and J.R. Willis 2 School of Mathematical Sciences University of Bath Bath BA2 7AY, U.K. 2 University of Cambridge Department of Applied Mathematics and Theoretical Physics Silver Street, Cambridge CB3 9EW, U.K. ABSTRACT Dynamic weight functions are constructed for general time{dependent shear loading of a plane semi{innite crack propagating with constant speed in an innite isotropic elastic body. The use of Fourier transforms reduces the problem to the analysis of a matrix Wiener{Hopf equation. The solution of the Wiener{Hopf problem is presented. An expression is derived for the rst{order perturbation to the stress intensity factors induced by a small time{dependent deviation from straightness of the crack front. The asymptotic procedure requires consideration of two terms of the asymptotic expansions of the displacement and stress tensor components in a neighbourhood of the crack front.. INTRODUCTION This is a continuation of a study of dynamic weight functions for three-dimensional loading of a plane, semi-innite crack, moving with speed V. Part I Willis and Movchan, 994 presented the basic framework and developed the solution in detail, for the case of Mode I loading. This required the solution of a scalar equation of Weiner-Hopf type. In contrast, even though the crack edge is straight, threedimensional shear loading inevitably produces coupling between Modes II and III. The Fourier transform formulation correspondingly generates, for shear loading, a
2 2 2 matrix Wiener-Hopf problem. This problem has particular structure that permits its solution in closed form; this is obtained. Next, the Fourier transform of the weight function matrix is employed, in a manner similar to that presented in Part I, to obtain a Fourier transformed relation between the rst-order perturbations to the Mode II and Mode III stress intensity factors, induced by a small deviation from straightness of the crack front. Inversion of the transforms is performed employing a technique, alternative to that of Cagniard 939, developed by Willis 973, to obtain the relation in integral form. A by-product of the analysis is a set of explicit relations between second-order terms in the expansions of relative displacement just behind, and traction just ahead of, the crack edge, in the three-dimensional time-dependent case. It is conrmed in an Appendix that, in the limiting case of statics, the formulae derived reproduce results that are already known Gao and Rice, 986. Since the basic formulation has already been presented in Part I, formulae are usually not repeated: reference to I, 3.3, for example, means equation 3.3 of Part I. However, just to set the scene, the reader is reminded that the basic problem is to solve the equations of isotropic elastodynamics, for an innite space which contains a crack occupying the region X := x? V t < 0;? < x 2 < ; x 3 = 0: : The crack faces are loaded, with tractions prescribed as i3 t; X; x 2 ; 0 =?P i t; X; x 2 for X < 0. Attention is focussed, initially, on the relative displacement of the crack faces, [u i ]t; X; x 2 for X < 0 and the traction components on the plane x 3 = 0, i3 t; X; x 2 ; 0 for X > 0. In particular, the weight functions generate the associated stress intensity factors. Later, the crack is taken to occupy the region X < "t; x 2,? < x 2 <, x 3 = 0, and the perturbations to the stress intensity factors are found, to rst order in ". 2. THE MATRIX WIENER{HOPF PROBLEM Part I see Willis and Movchan 994 introduced the notation [U i ] and i3 for the relative displacement and traction components on the plane x 3 = 0, associated with a weight function, and derived, in Appendix A, a relation between their Fourier transforms [ ~ U i ] and ~ i3 with respect to t; x ; and x 2. Mode I uncouples from Modes II and III and was discussed in Part I. For the coupled Modes II and III, the relation 2
3 takes the form given by I, A.6, which is equivalent to ~ 3 ~ 0 ; ; 2 = ? [U R t DR ~ ] 23 [U~ 2 ] where and D denotes the diagonal matrix 0 ; ; 2 ; 2. R = 2 ; 2.2? 2 id D = diag 0 ; ; 2 i 2 0 =b 2 0 ; =b 2? jj 2 =2 0 =b2? jj 2 =2 2 : 2.3 Here, notation is as dened in Part I: in particular, jj = =2 and D is given by I, 4.2. The matrix Wiener{Hopf problem is obtained, as described in Part I, by introducing transforms i3 and [U i ], which are the Fourier transforms of i3 and [U i ] with respect to variables t; X = x? V t and x 2. It is easily veried that, for a function ft; x ; x 2, f; ; 2 = ~ f? V ; ; 2 : 2.4 Thus, 2. relates i3 and [U i ]; i = ; 2; if 0 =? V. Furthermore, since [U i ] = 0 when X < 0 and i3 = 0 when X > 0, it follows that i3 is an analytic function of for each xed with Im > 0 and each real 2 in a half{plane Im < k +, while [U i ] is analytic in a half{plane Im > k?, for some k? ; k + satisfying?k? < 0 < k +. In preparation for the solution of the Wiener{Hopf equation, factorize D? V ; ; 2 as follows: D? V ; ; 2 = D + D? ; 2.5 where D are the diagonal matrices D + = diag [ + V=c 2 + iq c ]T + [ + V=b 2 + iq b ] ; [ =2 + V=b 2 + iq b ] =2 ; D? =? 2 diag + V=c 2? iq c b 2 RV T? ; [ V 2 [ + V=b 2? iq b ] =2 + V=b 2? iq b ] =2 :2.6 Notation is as dened in Part I; D + is analytic in Im > k? and D? is analytic in Im < k +. 3
4 Thus, in the strip k? < Im < k +,? = ? R t D + D? RU + ; 2.7 where? := 3 ; 23 t and U + := [U ]; [U 2 ] t. The objective is to deduce from 2.7 the vector functions?, analytic in Im < k +, and U +, analytic in Im > k?, such that U + is a homogeneous vector function of degree?=2 and, correspondingly,? is homogeneous of degree =2 in ; ; 2 : 3. SOLUTION OF THE WIENER{HOPF EQUATION It is noted that, as 2 0 or, equivalently, as =j 2 j, R I; where I is the identity matrix, and the Wiener{Hopf equation 2.7 uncouples, asymptotically, to give? D + D? U + ; 3. for which a candidate solution with the desired homogeneity is? D? C; U + D? + C; 3.2 for some vector vector function C of and 2, which is homogeneous of degree zero. This provides at least partial motivation for seeking a solution of 2.7 in the form? = D? C + R t? ; U + = D? + C 2 + R t u + ; 3.3 where C and C 2 are vector functions of and 2, homogeneous of degree zero. Substituting 3.3 into 2.7 and multiplying by R from the left give RD? C ? = D + D? RD +? C D + D? u + ; 3.4 or, upon rearrangement, ? C? C F; ; 2 + D +?? = D + u + ; 3.5 4
5 where with F; ; 2 = D?? ED? C? D + ED +? C 2 ; 3.6 : 3.7 E = 0? 0 To construct a general solution, two cases are considered: i C 2 = 0; C 6= 0 and ii C = 0; C 2 6= 0. First, suppose that C 2 = 0: Equation 3.5 takes the form " 2? I + 2 D?? ED? #C + D??? = D + u + : 3.8 The following additive split is obvious: " 2 + ij 2 j +? ij 2 j? sign 2 2i? sign 2 2i # + sign 2 D?? ED? 2i? ij 2 j " D?? ED?? fd?? ED? g;?ij 2 j; 2 fd?? ED? g;?ij 2 j; 2 + ij 2 j The usual Wiener{Hopf argument then gives u + = 2 + ij 2 j D +? + ij 2 j C = D + u +? D??? : 3.9 I + i sign 2 fd?? ED? g;?ij 2 j; 2 # C ; 3.0? =? 2 D?? ij 2 j? i sign 2D?? ED?? ij 2 j and, consequently, U + = Rt D +? 2 + ij 2 j +i sign 2 D?? ED?? fd?? ED? g;?ij 2 j; 2 + ij 2 j C ; 3. I + i sign 2 fd?? ED? g;?ij 2 j; 2 C ; 3.2 5
6 ? = 2 + ij 2 j Rt D? I + i sign 2 fd?? ED? g;?ij 2 j; 2 C : 3.3 It can be checked by direct calculation that the expression 3.3 does represent a \?" function: it does not have a singularity at =?ij 2 j, on account of the particular form of R. It should be mentioned that the rank of the matrix D := I + i sign 2 fd?? ED? g;?ij 2 j; is the same as that of I + i sign 2 E and so is equal to one. Thus, this solution contains, in eect, only one free parameter. For the case when C = 0, equation 3.5 can be rewritten as follows ? " I + 2 D + ED +? #C 2 = D?? which is equivalent to " # 2 + ij 2 j +? sign 2? ij 2 j 2i + sign 2 2i + sign 2 2i " D+?? D + u + ; 3.5 D + ED +? + ij 2 j ED +?? fd + ED +? g; ij 2 j; 2? ij 2 j fd + ED +? g;?ij 2 j; 2? ij 2 j Similarly to 3.2, 3.3 it follows in this case that and? = U + = Rt D +? 2? ij 2 j C 2 = D???? D + u + : 3.6 I? i sign 2 fd + ED +? g; ij 2 j; 2 2? i j 2 j Rt D? I? i sign 2 D + ED +? # ; ij 2 j; 2 C 2 ; 3.7 C 2 : 3.8 The particular form of R ensures that U + has no singularity at = ij 2 j: As in 3.4, the rank of the matrix D 2 := I? i sign 2 D + ED +? ; ij 2 j;
7 is equal to one. The solution of the Wiener{Hopf equation 3.5 is now U + = 2 Rt D +?" + ij 2 j D C +? = 2 Rt D? " + ij 2 j D C + # #? ij 2 j D2 C 2? ij 2 j D2 C 2 ; 3.20 : 3.2 Although C and C 2 are arbitrary vectors independent of, the fact that D and D 2 have rank one implies that this solution is expressible in terms of two arbitrary scalar parameters. Since D can be written in the form D = D?? its eigenvalues are readily evaluated as " I + i sign 2 E = 0; 2 = 2; ; X 2 = # D? ; and the corresponding eigenvectors are D?? X ; D?? X 2 ; where X = i sign2 i sign 2 Similarly, D 2 has eigenvalues = 0; 2 = 2 and corresponding eigenvectors D + X 2 ; D + X : Therefore, in the solution given by 3.20, 3.2 there is no loss in choosing C = c 2 D?? ;?ij 2 j; 2 X 2 ; C 2 = c 2 2 D +; ij 2 j; 2 X ; 3.22 where c ; c 2 are arbitrary scalar functions of ; 2 ; homogeneous of degree zero. Then, 3.20 takes the form U + = R t D +? where a ; b are dened so that 0 B@ i sign 2 + ij 2 ja? a +? ij 2 j + ij 2 jb? i sign 2b +? ij 2 j : CA c c 2 ; 3.23 diagfa? ; b? g := D? ;?ij 2 j; 2 ; diagfa + ; b + g := D + ; ij 2 j; 2 :
8 It is convenient to introduce, starting from 3.23, a \fundamental matrix" solution, with prescribed behaviour as : This is done by expressing c ; c 2 t in the form c a? b? i sign2 b =? +?a + l ; 3.25 c 2 a + a? + b + b??=b? i sign 2 =a? l 2 for some constant l ; l 2. Then U + = R t D +? a + a? + b + b? and, as, 0 B@ b+b? + ij 2 j + a +a?? ij 2 j? 2a?b ?=2 l U + + 0i l 2 2a + b? a + a? + ij 2 j + b +b?? ij 2 j CA l l 2 ; 3.26 : 3.27 The expression for? is similar: it is necessary only to replace D +? by D? on the right side of The matrix multiplying l ; l 2 t in 3.26 is the fundamental matrix that was sought on account of its simple form as : This section is concluded with the observation that, in correspondence with 3.27, the asymptotic form of U + as X 0 is U + t; X; x 2?i l =2 i HXtx 2 : 3.28 l 2 X 4. STRESS{INTENSITY FACTORS It was shown in Part I that [U i ] i3? k3 [u k ] t; X 0 ; x 0 2 = 0; 4. where i3 ; k3 denote stress components for u and U respectively, evaluated at x 3 = 0; as before, [:] denotes the jump across x 3 = 0 of the quantity indicated, and signies convolution with respect to t; X and x 2. Implications of the identity 8
9 4. are now explored in the case that [U ]; [U 2 ] t = U + ; 3 ; 23 t =? and [U 3 ] = 33 = 0. The convolution 4. will rst be evaluated when X 0 > 0. Since [u k ] and k3 are? functions, their convolution has the value zero and 4. reduces to U t + t; X 0 ; x 0 2 = 0; 4.2 where represents 3 ; 23 t, evaluated when x 3 = 0. Now, i3 t; X; x 0 2 ; 0 =?P it; X; x 2 H?X + + i3 t; X; x 2 ; 4.3 where P i t; X; x 2 i = ; 2 are prescribed. Furthermore, as X 0; t; X; x 2 2X?=2 KII K III HX: 4.4 The convolution in 4.2 is now readily evaluated, in the limit as X 0 +0; employing the asymptotic forms 3.28 and 4.4: i =2 i l ; l 2 2 KII t; x 2 K III t; x 2 + fu t + P? gt; 0; x 2 = 0; 4.5 where P? = P ; P 2 t H?X: This result motivates dening a weight function matrix, still called U +, so that it has the Fourier transform 2i U =2 + = R t D +? a + a? + b + b? 0 B@ b+b? + ij 2 j + a +a?? ij 2 j? 2a?b a + b? a + a? + ij 2 j + b +b?? ij 2 j 4.6 This is nothing but the fundamental matrix introduced in the preceding section, multiplied by 2i =2. The corresponding matrix? has Fourier transform? = 0 2i =2 R t D? a + a? + b + b? b+b? + ij 2 j + a +a?? ij 2 j? 2a?b a + b? a + a? + ij 2 j + b +b?? ij 2 j CA : CA : 4.7 9
10 The result 4.5 then gives KII t; x 2 = U t + P? t; 0; x K III t; x 2 2 : 4.8 It is also possible to develop the relation 4.8 by evaluating the convolutions U t + + and t? [u], as X 0; from the corresponding asymptotic forms of the transforms as ; as was done in Part I for the case of Mode-I loading. This will emerge from the more general reasoning presented next. 5. ASYMPTOTIC EXPANSIONS TO HIGHER ORDER As noted in Part I, the identity 4. induces a corresponding identity for its Fourier transform, [U i ] i3? k3 [u k ] = 0; 5. or U t + +? t?[u] =?U t + P? : 5.2 The identity 5.2 will now be expanded as. It is noted rst that, as X 0; + K2X?=2? P + AX =2 HX; 5.3 and [u] V K?2X= =2? B?X 3=2 H?X; 5.4 for some vectors P ; A and B, and some matrix V. Here, K denotes K II ; K III t : Correspondingly, as ; + i 2 =2 K ip? =2 + 0i + 0i? =2 A ; 5.5 2i =2 3=2 + 0i and Next, K [u]? i=2 3i? 3=2 2? 0i 4 i =2 B : 5.6 5=2? 0i U t + P? Ft; x 2? 22= =2 Lt; x 2 X =2 HX + =2 =2 Mt; x 2 X; 5.7 0
11 where Ft; x 2 = U t + P? t; 0; x 2 ; 5.8 Lt; x 2 = Pt; 0; x 2 ; 5.9 and Mt; x 2 = U t + P0? t; 0; x 2 ; 5.0 with P 0? Fourier transforming 5.7 gives, as ; # U t + P? i F" + 0i? L? 2i=2? 0i + 0i # 3=2 +i=2 M" =2 + 0i? : 5. 2? 0i 2 All of these results are direct counterparts of corresponding results in Part I for the case of Mode I loading. Now corresponding expansions of U + and? are required. First, for the diagonal matrices D, as ; A D?? 2? 0i =2? +? 0i A2? ; 5.2 D? + + 0i I?= i B2 + ; 5.3 where A? = diagfb; g; 5.4 A 2? = diag B V c? iq 2 c? 2 V b? iq b + t? ; 2 2 V b? iq b ; with and B 2 + =?diag V c 2 + iq c? 2 V b 2 + iq b + t + ; 2 V b 2 + iq b t =? 2i ; 5.6 B = b 2 RV =V 2 ; 5.7 Z C ln T ; 0 ; 2d 0 : 5.8
12 Notation not explicitly dened is as in Part I. It follows now, from 3.26 and 5.3, that where with and where U + 2i =2 + 0i?=2 q = a +a?? b + b? a + a? + b + b? 0 0? i I + ; + 0i Q 5.9 iq = B 2 + ij 2 jq + 2 2f? E; 5.20 ; 5.2 f = a + a? + b + b?? 0 a + b??a? b + 0?? I? 2 E AV? : 5.22 Inspection of the Wiener{Hopf equation 2. reveals that, correct to order j j?=2 ; " #? I + 2 AV = diag 2V V=b 2 2 ; RV V U + ; 5.23 : 5.24 E Thus, as,??2i =2? 0i =2 [AV ] I? +? i0 A0 V [AV ]? + 2 E? AV E[AV ]? + iq : 5.25 These expressions suce to expand the identity 5.2 to three terms as : ik + 0i? i2i=2 P + 0i 3=2? =2=2 A + Q t K + 0i 2? i=2=2 [AV ]? K? 0i i=2 =2 [AV ]? A 0 V? 2 E? [AV ]? EAV + iq t[av ]? K?? 0i 2 2
13 " # + 3=2=2=2 [AV ]? B = if? 0i 2 + 0i? L? i2i=2? 0i + 0i # 3=2? =2 M" =2 + 0i? : ? 0i 2 It follows, therefore, that K = F; = 2 =2 AV ; A = M? 2 Q t =2 K; B = 2 3 AV A + i 2 =2? 2 E? AV E[AV ]? " t# [AV ]? A 0 V K : 5.27 The rst of relations 5.27 conrms the result 4.8. The second is quoted, for example, by Freund 990. The expressions for A and B were not known previously to the authors. 6. PERTURBATION OF THE CRACK FRONT As in the Mode I analysis of Part I, the identity 5. can also be applied when the stress and displacement elds solve the problem of a crack occupying the region X = x? V t < "t; x 2 ;? < x 2 < ; x 3 = 0; 6. where " is a small positive constant and the function is smooth and bounded. The displacement and stress elds for this problem are written as u + u; + ; where u; solve the problem for the unperturbed crack " = 0: Since 5. is satised for both the perturbed and unperturbed elds, it follows that U t +? t?[u] = 0: 6.2 This identity applies for all and all "; however, information on the rst{order perturbation of the stress{intensity factors will be sought by considering the lowest{ order non{trivial asymptotic expansion of 6.2 as and " 0: 3
14 As in Part I, since corresponds to X 0, it suces to note that, as X 0 and " 0; and [u + u] + K + K [2X? "] =2 + X? "=2 A + A? P + P + " 0 E K 2X =2? P + X =2 A HX? " "? X =2 + K + K? "? X 3=2 B + B +" 0 E?X =2 K??X 3=2 B HX; 6.3 H"? X H?X: 6.4 In comparison with the Mode I problem the expansions 6.3, 6.4 include additional terms which arise from the local rotation of the crack edge. Fourier transforming 6.3, 6.4 and performing formal expansions with respect to " and = gives and i 2=2 + 0i?=2 K + "E 0 K + i" K? =2 i + 0i P + "E 0 P? 2i =2 + 0i 3=2A + i" A + "E 0 A; 6.5 [u]?i=2 K + "E 2? 0i 0 K + K + i" 3=2 K? 3i 4 i =2? 0i 5=2 B + "E 0 B + i" B These formulae are vector analogues of I, 8.6 and I, 8.7. The asymptotic relations 5.9 and 6.5 yield i U t +?"K + + 0i?"Q t K? 2 " K + "E 0 K =2 "A # : 6.6 ; 6.7 4
15 and 5.25 and 6.4 give " i t?[u]?"k + K + "[AV ]? EAV? 0i 0 K + =2 =2 [AV ]? K +" i[av ]? A 0 V? i 2 E? [AV ]? EAV? Q tk? " 3 2 =2 [AV ]? B 2 # : 6.8 The terms that are independent of in 6.7 and 6.8 are the same; hence, they cancel, upon substitution into 6.2. It also follows, from the identities 5.27, that the coecient of + 0i? in 6.7 is the same as the coecient of? 0i? in 6.8. Hence, the identity 6.2 implies to rst order in ", that K = " Q t K? E 0 K + =2 =2 A ; 6.9 where A is given by the third of relations Inverting the transforms yields K = " Q t K? E 0 K + =2 =2 M? Q t K : 6.0 Here, the convolution is over the variables t; x 2 and, since Q; 2 is homogeneous of degree one, it follows from I, 8.2 that " Qt; x 2 2 Q?x2 # + 0i; t + Qx 2 + 0i;?t Ht : 6. 2?x 2 + 0i 2 x 2 + 0i 2 The result 6.0 is the mixed Mode II/III analogue of the Mode I formula I, 8.. The Appendix considers the static limit of the solution. This is obtained by setting V = 0 and allowing to tend to zero. This static limit was addressed directly by Gao and Rice 986; the agreement that is found provides a partial check on the present analysis. 5
16 7. CONCLUDING REMARKS Parts I and II together provide a complete weight function matrix for any loading of the semi{innite crack in uniform motion. The matrix partitions into a scalar component corresponding to Mode I, and a 22 matrix relating to the coupled Modes II and III. Fourier transforms were employed to derive the weight function matrix. Treatment of the Mode I component required the solution of a scalar Wiener{Hopf problem whose solution is routine; the coupled Modes II and III generated a 2 2 matrix Wiener{Hopf problem for which a closed{form solution was derived. Use of a reciprocal theorem permitted the deduction of the perturbation to the stress intensity factors induced by a small time{ and space{dependent deviation from straightness of the crack front. The perturbation is expressed by I, 8. for Mode I, and by 6.0 of the present paper for Modes II and III. Implications of these formulae will be followed up elsewhere. 6
17 REFERENCES Cagniard, L., 939, Reexion et Refraction des Ondes Seismiques Progressives, Gauthier{ Villars, Paris. Freund, L.B., 990, Dynamic Fracture Mechanics, Cambridge Univ. Press. Gao, H. and Rice, J.R., 986, Shear Stress Intensity Factors For a Planar Crack With Slightly Curved Front, ASME Journal of Applied Mechanics, 53, 774{778. Uyand, Y.S., 965, Survey of Articles on the Application of Integral Transforms in the Theory of Elasticity, North Carolina State University, Department of Applied Mathematics, File No. PSR-24/6, Raleigh, N.C. Willis, J.R., 973, Self{Similar Problems in Elastodynamics, Phil. Lond. A274, 435{49. Trans. R. Soc. Willis, J.R. and Movchan, A.B., 994, Dynamic Weight Functions For a Moving Crack. I. Mode I Loading, J. Mech. Phys. Solids. 7
18 APPENDIX The static limit can be obtained from the full dynamic problem by setting V = 0 and allowing to tend to zero. The Wiener{Hopf relation 2. then reduces to where the diagonal matrix is dened as? =? RU ; A. =2Rt = diagf?? ; g; A.2 with denoting Poisson's ratio. Comparison with the full dynamic equation 2.7 shows that D =? =2 ; and, from 2.6, D + + ij 2 j =2 I; D?? 2? ij 2 j =2 : A.3 The static limit of the weight function 4.6 is easily calculated, upon noting that Substituting these values into 4.6 gives U + 22i=2 2? + ij 2 j?=2 Inversion of the transforms then yields U + = p 2x Hx 3=2 x 2 + x 2 2 a + 2ij 2 j =2 ; b + 2ij 2 j =2 ; a??i 2ij 2j =2 2? ; b??i 2 2ij 2j =2 : A.4 0 B@ + 0 B@ 2? + ij 2 j? 2 + ij 2 j 2 2? x 2 =x 2? 4 2?? 2 + ij 2 j 2? + + ij 2 j x 2 =x 4 + x 2 =x 2 2? + x 2 =x 2 x 2 =x? 2? x 2 =x 2 + x 2 =x 2 2? + x 2 =x 2 8 CA : CA : A.5 A.6
19 Allowing for the replacement of x by?x because of the present use of notation that permits the use of convolutions, this agrees with the weight function given by equation 5 of Gao and Rice 986, and attributed by them originally to Uyand 965. Similarly, in the relation 6.9 between K and, the matrix Q reduces to Q =? j 2j 2? i 2 0 : A.7 22? ? 0 Hence, inverting the transform with respect to 2 ; Q = 22? 2? x 2 2 2? the generalized function =x 2 2 part integrals. Substitution of A.8 into 6.0 then gives Kx 2 = " 22? + 2 2? 0 0 x 0 2 ; A.8 being interpreted so as to generate Hadamard nite 2? ?? 0 Z? 0 x 2 Kx Kx 0 2[x 0 2? x 2 ] x 2? x dx 0 2 =2 x 2 Mx 2 ; A.9 This agrees precisely with equation 5 of Gao and Rice 986, except that the term involving M was absent from their work, through their choice of local coordinates. 9
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In
More informationELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY. David O. Olagunju. University of Delaware. Newark, DE 19716
ELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY David O. Olagunju Department of Mathematical Sciences University of Delaware Newark, DE 19716 June 15, 1995 Abstract Consider the axisymmetric,
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationCrack dynamics in elastic media
PHILOSOPHICAL MAGAZINE B, 1998, VOL. 78, NO. 2, 97± 102 Crack dynamics in elastic media By Mokhtar Adda-Bedia and Martine Ben Amar Laboratoire de Physique Statistique de l Ecole Normale Supe  rieure,
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationExplicit evaluation of the transmission factor T 1. Part I: For small dead-time ratios. by Jorg W. MUller
Rapport BIPM-87/5 Explicit evaluation of the transmission factor T (8,E) Part I: For small dead-time ratios by Jorg W. MUller Bureau International des Poids et Mesures, F-930 Sevres Abstract By a detailed
More informationCRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th
ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationFrom Lay, 5.4. If we always treat a matrix as defining a linear transformation, what role does diagonalisation play?
Overview Last week introduced the important Diagonalisation Theorem: An n n matrix A is diagonalisable if and only if there is a basis for R n consisting of eigenvectors of A. This week we ll continue
More informationOn some properties of elementary derivations in dimension six
Journal of Pure and Applied Algebra 56 (200) 69 79 www.elsevier.com/locate/jpaa On some properties of elementary derivations in dimension six Joseph Khoury Department of Mathematics, University of Ottawa,
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
More informationlinearly indepedent eigenvectors as the multiplicity of the root, but in general there may be no more than one. For further discussion, assume matrice
3. Eigenvalues and Eigenvectors, Spectral Representation 3.. Eigenvalues and Eigenvectors A vector ' is eigenvector of a matrix K, if K' is parallel to ' and ' 6, i.e., K' k' k is the eigenvalue. If is
More informationNumerical Linear Algebra Homework Assignment - Week 2
Numerical Linear Algebra Homework Assignment - Week 2 Đoàn Trần Nguyên Tùng Student ID: 1411352 8th October 2016 Exercise 2.1: Show that if a matrix A is both triangular and unitary, then it is diagonal.
More informationContents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces
Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS
ULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS A. S. Cheng Center for QEFP, Northwestern University Evanston, IL 60208-3020 INTRODUCTION A number of researchers have demonstrated
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationWhat is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =
STUDENT S COMPANIONS IN BASIC MATH: THE ELEVENTH Matrix Reloaded by Block Buster Presumably you know the first part of matrix story, including its basic operations (addition and multiplication) and row
More informationBrittle fracture dynamics with arbitrary paths III. The branching instability under general loading
Journal of the Mechanics and Physics of Solids 53 (2005) 227 248 www.elsevier.com/locate/jmps Brittle fracture dynamics with arbitrary paths III. The branching instability under general loading M. Adda-Bedia
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationAsymptotics of generating the symmetric and alternating groups
Asymptotics of generating the symmetric and alternating groups John D. Dixon School of Mathematics and Statistics Carleton University, Ottawa, Ontario K2G 0E2 Canada jdixon@math.carleton.ca October 20,
More informationFracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach
University of Liège Aerospace & Mechanical Engineering Fracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach Ludovic Noels Computational & Multiscale Mechanics of
More informationChapter 1: Systems of Linear Equations
Chapter : Systems of Linear Equations February, 9 Systems of linear equations Linear systems Lecture A linear equation in variables x, x,, x n is an equation of the form a x + a x + + a n x n = b, where
More informationDiagonalization by a unitary similarity transformation
Physics 116A Winter 2011 Diagonalization by a unitary similarity transformation In these notes, we will always assume that the vector space V is a complex n-dimensional space 1 Introduction A semi-simple
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More informationFunctions of Several Variables
Functions of Several Variables The Unconstrained Minimization Problem where In n dimensions the unconstrained problem is stated as f() x variables. minimize f()x x, is a scalar objective function of vector
More informationB. Physical Observables Physical observables are represented by linear, hermitian operators that act on the vectors of the Hilbert space. If A is such
G25.2651: Statistical Mechanics Notes for Lecture 12 I. THE FUNDAMENTAL POSTULATES OF QUANTUM MECHANICS The fundamental postulates of quantum mechanics concern the following questions: 1. How is the physical
More informationProblem Set 9 Due: In class Tuesday, Nov. 27 Late papers will be accepted until 12:00 on Thursday (at the beginning of class).
Math 3, Fall Jerry L. Kazdan Problem Set 9 Due In class Tuesday, Nov. 7 Late papers will be accepted until on Thursday (at the beginning of class).. Suppose that is an eigenvalue of an n n matrix A and
More informationav 1 x 2 + 4y 2 + xy + 4z 2 = 16.
74 85 Eigenanalysis The subject of eigenanalysis seeks to find a coordinate system, in which the solution to an applied problem has a simple expression Therefore, eigenanalysis might be called the method
More information4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial
Linear Algebra (part 4): Eigenvalues, Diagonalization, and the Jordan Form (by Evan Dummit, 27, v ) Contents 4 Eigenvalues, Diagonalization, and the Jordan Canonical Form 4 Eigenvalues, Eigenvectors, and
More informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationTable 1: BEM as a solution method for a BVP dierential formulation FDM BVP integral formulation FEM boundary integral formulation BEM local view is ad
Chapter 1 Introduction to Boundary element Method - 1D Example For reference: Hong-Ki Hong and Jeng-Tzong Chen, Boundary Element Method, Chapter 1 Introduction to Boundary element Method - 1D Example,
More informationMath113: Linear Algebra. Beifang Chen
Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationMAT 2037 LINEAR ALGEBRA I web:
MAT 237 LINEAR ALGEBRA I 2625 Dokuz Eylül University, Faculty of Science, Department of Mathematics web: Instructor: Engin Mermut http://kisideuedutr/enginmermut/ HOMEWORK 2 MATRIX ALGEBRA Textbook: Linear
More informationELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS. 1. Linear Equations and Matrices
ELEMENTARY LINEAR ALGEBRA WITH APPLICATIONS KOLMAN & HILL NOTES BY OTTO MUTZBAUER 11 Systems of Linear Equations 1 Linear Equations and Matrices Numbers in our context are either real numbers or complex
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationFundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved
Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural
More informationSeparation of variables
Separation of variables Idea: Transform a PDE of 2 variables into a pair of ODEs Example : Find the general solution of u x u y = 0 Step. Assume that u(x,y) = G(x)H(y), i.e., u can be written as the product
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationas specic as the Schur-Horn theorem. Indeed, the most general result in this regard is as follows due to Mirsky [4]. Seeking any additional bearing wo
ON CONSTRUCTING MATRICES WITH PRESCRIBED SINGULAR VALUES AND DIAGONAL ELEMENTS MOODY T. CHU Abstract. Similar to the well known Schur-Horn theorem that characterizes the relationship between the diagonal
More informationELASTIC COMPLEX ANALYSIS AND ITS APPLICATIONS IN FRACTURE MECHANICS. N. D. Aparicio. Oxford Brookes University, Gipsy Lane, Oxford OX3 OBP, UK.
EASTIC COMPEX ANAYSIS AND ITS APPICATIONS IN FRACTURE MECHANICS N. D. Aparicio School of Computing & Mathematical Sciences, Oxford Brookes University, Gipsy ane, Oxford OX3 OBP, UK. April 9, 998 Abstract
More informationSolution Set 7, Fall '12
Solution Set 7, 18.06 Fall '12 1. Do Problem 26 from 5.1. (It might take a while but when you see it, it's easy) Solution. Let n 3, and let A be an n n matrix whose i, j entry is i + j. To show that det
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More information(Received 5 May 1980; accepted for print 21 September 1980)
MECHANICS RESEARCH COMMUNICATIONS Vol. 7(6),351-358,1980. Printed in the USA. 0093-6413/80/060351-08502.00/0 Copyright (c) Pergamon Press Ltd ON THE DISPLACEMENT INDUCED IN A RIGID CIRCULAR PUNCH ON AN
More informationPerturbative simulations of crack front waves
Journal of the Mechanics and Physics of Solids 48 () 19±151 www.elsevier.com/locate/jmps Perturbative simulations of crack front waves John W. Morrissey a, James R. Rice b, * a Numerical Algorithms Group,
More informationModels for dynamic fracture based on Griffith s criterion
Models for dynamic fracture based on Griffith s criterion Christopher J. Larsen Abstract There has been much recent progress in extending Griffith s criterion for crack growth into mathematical models
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
More informationYadernaya Fizika, 56 (1993) Soviet Journal of Nuclear Physics, vol. 56, No 1, (1993)
Start of body part published in Russian in: Yadernaya Fizika, 56 (1993) 45-5 translated in English in: Soviet Journal of Nuclear Physics, vol. 56, No 1, (1993) A METHOD FOR CALCULATING THE HEAT KERNEL
More informationTitle: A note on perturbation formulae for the surface-wave speed due to perturbations in material properties
Editorial Manager(tm) for Journal of Elasticity Manuscript Draft Manuscript Number: ELAS248R2 Title: A note on perturbation formulae for the surface-wave speed due to perturbations in material properties
More informationThroughout these notes we assume V, W are finite dimensional inner product spaces over C.
Math 342 - Linear Algebra II Notes Throughout these notes we assume V, W are finite dimensional inner product spaces over C 1 Upper Triangular Representation Proposition: Let T L(V ) There exists an orthonormal
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationLecture Notes in Linear Algebra
Lecture Notes in Linear Algebra Dr. Abdullah Al-Azemi Mathematics Department Kuwait University February 4, 2017 Contents 1 Linear Equations and Matrices 1 1.2 Matrices............................................
More informationNon-Abelian Berry phase and topological spin-currents
Non-Abelian Berry phase and topological spin-currents Clara Mühlherr University of Constance January 0, 017 Reminder Non-degenerate levels Schrödinger equation Berry connection: ^H() j n ()i = E n j n
More information1 Introduction We will consider traveling waves for reaction-diusion equations (R-D) u t = nx i;j=1 (a ij (x)u xi ) xj + f(u) uj t=0 = u 0 (x) (1.1) w
Reaction-Diusion Fronts in Periodically Layered Media George Papanicolaou and Xue Xin Courant Institute of Mathematical Sciences 251 Mercer Street, New York, N.Y. 10012 Abstract We compute the eective
More informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
More informationgr-qc/ Sep 1995
DAMTP R95/45 On the Evolution of Scalar Metric Perturbations in an Inationary Cosmology R. R. Caldwell University of Cambridge, D.A.M.T.P. Silver Street, Cambridge CB3 9EW, U.K. email: R.R.Caldwell@amtp.cam.ac.uk
More informationAbstract. Jacobi curves are far going generalizations of the spaces of \Jacobi
Principal Invariants of Jacobi Curves Andrei Agrachev 1 and Igor Zelenko 2 1 S.I.S.S.A., Via Beirut 2-4, 34013 Trieste, Italy and Steklov Mathematical Institute, ul. Gubkina 8, 117966 Moscow, Russia; email:
More informationIntroduction to Linear Algebra. Tyrone L. Vincent
Introduction to Linear Algebra Tyrone L. Vincent Engineering Division, Colorado School of Mines, Golden, CO E-mail address: tvincent@mines.edu URL: http://egweb.mines.edu/~tvincent Contents Chapter. Revew
More informationSolution to Homework 1
Solution to Homework Sec 2 (a) Yes It is condition (VS 3) (b) No If x, y are both zero vectors Then by condition (VS 3) x = x + y = y (c) No Let e be the zero vector We have e = 2e (d) No It will be false
More informationStatistics 992 Continuous-time Markov Chains Spring 2004
Summary Continuous-time finite-state-space Markov chains are stochastic processes that are widely used to model the process of nucleotide substitution. This chapter aims to present much of the mathematics
More informationMethod for calculating the stress intensity factor for mode-i indentation with eccentric loads
Acta Technica 6 (017), No. 4A, 481488 c 017 Institute of Thermomechanics CAS, v.v.i. Method for calculating the stress intensity factor for mode-i indentation with eccentric loads DUO Yili 1,, XIE Yujun
More informationRATIONAL POINTS ON CURVES. Contents
RATIONAL POINTS ON CURVES BLANCA VIÑA PATIÑO Contents 1. Introduction 1 2. Algebraic Curves 2 3. Genus 0 3 4. Genus 1 7 4.1. Group of E(Q) 7 4.2. Mordell-Weil Theorem 8 5. Genus 2 10 6. Uniformity of Rational
More information2 Eigenvectors and Eigenvalues in abstract spaces.
MA322 Sathaye Notes on Eigenvalues Spring 27 Introduction In these notes, we start with the definition of eigenvectors in abstract vector spaces and follow with the more common definition of eigenvectors
More informationLinear Algebra (part 1) : Matrices and Systems of Linear Equations (by Evan Dummit, 2016, v. 2.02)
Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 206, v 202) Contents 2 Matrices and Systems of Linear Equations 2 Systems of Linear Equations 2 Elimination, Matrix Formulation
More informationMatrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
More informationLocalized transition waves in bistable-bond lattices
Journal of the Mechanics and Physics of Solids 52 (24) 447 479 www.elsevier.com/locate/jmps Localized transition waves in bistable-bond lattices L.I. Slepyan a;, M.V. Ayzenberg-Stepanenko b a Department
More information5 Eigenvalues and Diagonalization
Linear Algebra (part 5): Eigenvalues and Diagonalization (by Evan Dummit, 27, v 5) Contents 5 Eigenvalues and Diagonalization 5 Eigenvalues, Eigenvectors, and The Characteristic Polynomial 5 Eigenvalues
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationDegenerate Perturbation Theory. 1 General framework and strategy
Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying
More informationEXAMPLES OF PROOFS BY INDUCTION
EXAMPLES OF PROOFS BY INDUCTION KEITH CONRAD 1. Introduction In this handout we illustrate proofs by induction from several areas of mathematics: linear algebra, polynomial algebra, and calculus. Becoming
More informationOn the torsion of functionally graded anisotropic linearly elastic bars
IMA Journal of Applied Mathematics (2007) 72, 556 562 doi:10.1093/imamat/hxm027 Advance Access publication on September 25, 2007 edicated with admiration to Robin Knops On the torsion of functionally graded
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationwith angular brackets denoting averages primes the corresponding residuals, then eq. (2) can be separated into two coupled equations for the time evol
This paper was published in Europhys. Lett. 27, 353{357, 1994 Current Helicity the Turbulent Electromotive Force N. Seehafer Max-Planck-Gruppe Nichtlineare Dynamik, Universitat Potsdam, PF 601553, D-14415
More informationMidterm 1 NAME: QUESTION 1 / 10 QUESTION 2 / 10 QUESTION 3 / 10 QUESTION 4 / 10 QUESTION 5 / 10 QUESTION 6 / 10 QUESTION 7 / 10 QUESTION 8 / 10
Midterm 1 NAME: RULES: You will be given the entire period (1PM-3:10PM) to complete the test. You can use one 3x5 notecard for formulas. There are no calculators nor those fancy cellular phones nor groupwork
More informationTopics. Vectors (column matrices): Vector addition and scalar multiplication The matrix of a linear function y Ax The elements of a matrix A : A ij
Topics Vectors (column matrices): Vector addition and scalar multiplication The matrix of a linear function y Ax The elements of a matrix A : A ij or a ij lives in row i and column j Definition of a matrix
More informationChapter 6 Balanced Realization 6. Introduction One popular approach for obtaining a minimal realization is known as Balanced Realization. In this appr
Lectures on ynamic Systems and Control Mohammed ahleh Munther A. ahleh George Verghese epartment of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 6 Balanced
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
More informationLinear Algebra and Matrix Inversion
Jim Lambers MAT 46/56 Spring Semester 29- Lecture 2 Notes These notes correspond to Section 63 in the text Linear Algebra and Matrix Inversion Vector Spaces and Linear Transformations Matrices are much
More informationMatrix operations Linear Algebra with Computer Science Application
Linear Algebra with Computer Science Application February 14, 2018 1 Matrix operations 11 Matrix operations If A is an m n matrix that is, a matrix with m rows and n columns then the scalar entry in the
More informationRegression. Oscar García
Regression Oscar García Regression methods are fundamental in Forest Mensuration For a more concise and general presentation, we shall first review some matrix concepts 1 Matrices An order n m matrix is
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationLinear Algebra, 4th day, Thursday 7/1/04 REU Info:
Linear Algebra, 4th day, Thursday 7/1/04 REU 004. Info http//people.cs.uchicago.edu/laci/reu04. Instructor Laszlo Babai Scribe Nick Gurski 1 Linear maps We shall study the notion of maps between vector
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
More informationAbsolutely indecomposable symmetric matrices
Journal of Pure and Applied Algebra 174 (2002) 83 93 wwwelseviercom/locate/jpaa Absolutely indecomposable symmetric matrices Hans A Keller a; ;1, A Herminia Ochsenius b;1 a Hochschule Technik+Architektur
More information3 (Maths) Linear Algebra
3 (Maths) Linear Algebra References: Simon and Blume, chapters 6 to 11, 16 and 23; Pemberton and Rau, chapters 11 to 13 and 25; Sundaram, sections 1.3 and 1.5. The methods and concepts of linear algebra
More informationSolution: a) Let us consider the matrix form of the system, focusing on the augmented matrix: 0 k k + 1. k 2 1 = 0. k = 1
Exercise. Given the system of equations 8 < : x + y + z x + y + z k x + k y + z k where k R. a) Study the system in terms of the values of k. b) Find the solutions when k, using the reduced row echelon
More information4.3 - Linear Combinations and Independence of Vectors
- Linear Combinations and Independence of Vectors De nitions, Theorems, and Examples De nition 1 A vector v in a vector space V is called a linear combination of the vectors u 1, u,,u k in V if v can be
More informationChapter 6 Balanced Realization 6. Introduction One popular approach for obtaining a minimal realization is known as Balanced Realization. In this appr
Lectures on ynamic Systems and Control Mohammed ahleh Munther A. ahleh George Verghese epartment of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 6 Balanced
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More information1.4 Linear Transformation I
.4. LINEAR TRANSFORMATION I.4 Linear Transformation I MATH 9 FALL 99 PRELIM # 5 9FA9PQ5.tex.4. a) Consider the vector transformation y f(x) from V to V such that if y (y ; y ); x (x ; x ); y (x + x ) p
More informationMoving screw dislocations in piezoelectric bimaterials
phys stat sol (b) 38 No 1 10 16 (003) / DOI 10100/pssb00301805 Moving screw dislocations in piezoelectric bimaterials Xiang-Fa Wu *1 Yuris A Dzenis 1 and Wen-Sheng Zou 1 Department of Engineering Mechanics
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More information