BEM Formulation plane orthotropic bodies a modification for exterior regions and its proof
|
|
- Kerry Flowers
- 5 years ago
- Views:
Transcription
1 Ŕ periodica polytechnica Civil Engineering 5/ 7 3 web: ci c Periodica Polytechnica 7 BEM Formulation plane orthotropic bodies a modification for exterior regions its proof RESEARCH ARTICLE Received 7--7 Abstract Assuming linear displacements constant strains stresses at infinity our aim is to reformulate the equations of the direct boundary element method for plane problems of elasticity. We shall consider a body made of orthotropic material. The equations we have reformulated make possible to attack plane problems on exterior regions without replacing the region in question by a bounded one. Keywords Boundary element method Somigliana relations direct method orthotropic body plane problems exterior regions Acknowledgement The support provided by the Hungarian National Research Foundation within the framework of the project OTKA T46834 is gratefully acknowledged. Introduction As is well known a large literature studies plane problems for orthotropic bodies including for instance as well as the books the references therein. However as explained by Schiavone the stard formulation for exterior regions has the disadvantage that it is impossible to prescribe a constant stress state at infinity. The reason is that an assumption about the far field pattern of the displacements is needed in order to establish an appropriate Betti s formula to prove uniqueness existence for the exterior Dirchlet Neuman problems. Unfortunately this assumption excludes those problems from the theory for which the displacements are linear while the strains stresses are constant at infinity. To make progress on plane problems with such displacements we note that if the direct formulation reproduces this displacement field then the resulting strain stress conditions must also be constant at infinity. Consequently plane problems for the exterior regions can be attacked without replacing the region by a bounded one. The work in presents such direct formulations by assuming constant strains stresses at infinity for an isotropic body. For exterior regions reformulates the classical approach to plane problems. For the same class of problems but in a dual formulation sets up the equations of the direct method in terms of stress functions of order one. The present paper is an attempt to clarify how the formulation changes if we apply the ideas presented in paper to orthotropic bodies. György Szeidl Department of Mechanics University of Miskolc H-355 Miskolc- Egyetemváros Hungary Gyorgy.SZEIDL@uni-miskolc.hu Judit Dudra Department of Mechanics University of Miskolc H-355 Miskolc- Egyetemváros Hungary dudrajudit@fre .hu Preliminaries Throughout this paper x x are rectangular Cartesian coordinates referred to an origin O. Greek subscripts are assumed to have the range summation over repeated subscripts is implied. The doubly connected exterior region under consideration is denoted by A e is bounded by the contour. We stipulate that the contour admits a nonsingular parametrization in terms of its arc length s. The outer normal is denoted by n π. In accordance with the notations introduced δ κλ is the Kronecker BEM Formulation plane orthotropic bodies 7 5
2 symbol α sts for the derivatives taken with respect to x α ɛ 3κλ is the permutation symbol. Assuming plane problems let u κ e κλ t κλ be the displacement field the in plane components of strain stress respectively. For orthotropic bodies the material constants are denoted by s s = s s. For homogenous orthotropic material the plane problem of classical elasticity is governed by the kinematic equations Hook s law e ρλ = ρu λ + λ u ρ t = c e + c e t = c e + c e t = t = c 66 e c = s d c = c = s d c = s d the equilibrium equations c 66 = d = s s s ; 3 t ρλ λ + b ρ = 4 which should be complemented with appropriate boundary conditions not detailed here since they play no role in the present investigations. The basic equation for u λ takes the form D ρλ u λ + b ρ = the differential operator D ρλ has the form Dρλ = c +c 66 c + c 66 c + c 66 c +c 66 5a. 5b Let Qξ ξ M x x be two points in the plane the source point the field point. We shall assume temporarily that the point Q is fixed. The distance between Q M is R the position vector of M relative to Q is r κ. The small circle as a subscript for instance M or Q indicates that the corresponding points i.e. Q or M are taken on the contour. It is obvious that r α M Q = x α M ξ α Q = x α ξ α. 6 Let us introduce the following notations λ + λ = s + /s 7 λ λ = s /s 8 A α = s λ α s 9 ρ α = λ αr + r D = π λ λ s. For our later considerations we note that Eqs. 7 8 imply λ = s + s ± s + s s s. The well-known singular fundamental solutions for the basic Eq. 5a 7 9 are given by the formulas U M Q = D λ A ln ρ λ A ln ρ λ λ r r U M Q = D A A arctan λ λ r + r 3 U M Q = U M Q A U M Q = D ln ρ A ln ρ λ λ λ A λ A T M Q = D r n + r n ρ λ A λ A T M Q = D r n T M Q = D ρ λ A ρ λ λ A ρ ρ ρ ρ λ A ρ r n } λ λ A r n ρ } λ A λ A r n λ A λ A T M Q = D r n + r n 4 ρ u λ M = U λκ M Qe κ Q t λ M = T λκ M Qe κ Q are the displacement vector stress vector on a line element with a normal n λ = n λ M to it at the point M due to the force e κ = e κ Q at Q. 3 Basic formulas for exterior regions Fig. depicts a triple connected region A e bounded by the contours L ε the circle with radius center at O. Here L ε is the contour of the neighborhood A ε of Q with radius R ε while is sufficiently large so that the region bounded by covers both L L ε. If R ε then clearly A e A e. Let u κ M g κ M be sufficiently smooth continuously differentiable at least twice but otherwise arbitrary displacement fields on A e. The stresses obtained from these displacement fields are denoted by t λκ uρ M t λκ gρ M respectively. ρ ρ Per. Pol. Civil Eng.
3 O x x Lo ds ds Q e n L A R A e we obtain A e u λ M µ D M λσ U σ κ M Q µ D M λσ u σ M U λκ M Q da M e κ Q = = u λ M T λκ M Q t λ M U λκ M Q ds M e κ Q+ u λ M T λκ M Q t λ M U λκ M Q ds M e κ Q L ε u λ M T λκ M Q t λ M U λκ M Q ds M e κ Q 8 since t λκ uρ M n κ M = t λ M is the stress on the contour obviously ds t λκ gρ M n κ M = T λκ M Qe κ Q. Fig.. Equation A e u λ M µ D M λσ g σ M g λ M µ D M λσ u σ M da M = uλ M t λκ gρ M n κ M g λ M t λκ uρ M n κ M ds M + L ε uλ M t λκ gρ M n κ M g λ M t λκ uρ M n κ M ds M + uλ M t λκ gρ M n κ M g λ M t λκ uρ M n κ M ds M 5 in which M over a letter denotes derivatives taken with respect to the point coordinates M n κ M is the outward normal is the primal Somigliana identity applied to the triple connected region A e. Let g λ Q = U λκ M Qe κ Q which is a non singular elastic state of the plane in A e. We regard u λm as a different elastic state in the region A e. Further we assume that u λ M has the far field pattern asymptotic behaviour ũ κ M = c κ + ε 3ρκ x ρ ω + e κβ x β 6 as x β or equivalently M tends to infinity. Here c κ is a translation ω is a rotation in finite c κ + ε 3ρκ x ρ ω is the corresponding rigid body motion e κβ is a constant strain tensor at infinity e κβ x β is the corresponding displacement field. The stresses induced by the strains e κβ can be obtained by the Hooke law: In the sequel we shall assume that there are no body forces. This assumption has no effect on the result we will obtain. It is clear that one can omit e κ Q. Regarding the equation obtained by omitting e κ Q our goal is to compute its limit as R ε. As is well known the left h side vanishes under the conditions detailed abovesee hence + lim = R ε L ε u κ Q + u λ MT o λκ M Q t λ M U λκ M Q ds M. Consequently u κ Q = lim t λ M U λκ M Q u λ M T λκ M Q ds M + L R t λ M U λκ M Q u λ M T λκ M Q ds M. 9 In order to establish the first Somigliana formula for the exterior region we need to find the limit of the first integral on the right h side. 4 Somigliana formulas modified for exterior regions In this section our main objective is to prove that I κ = lim t λ M U λκ M Q u λ M T λκ M Q ds M = c κ + ε 3ρκ ξ ρ ω + e κβ ξ β = ũ κ Q. t = c e + c e t = c e + c e t = t = c 66 e 7 Using x β o M = n β o M Substituting the above quantities into the Somigliana identity BEM Formulation plane orthotropic bodies 7 5 3
4 we substitute t λρ n ρ M for t λ M c λ + ε 3ρλ x ρ ω + e κβ x β for u λ M. This implies I κ = I κ + I κ + 3 I κ + 4 I κ = lim c λ T λκ M Q ds M lim ε 3ρλω n ρ MT o λκ M Q ds M + L R lim t λρ n ρ M U λκ M Q ds M lim e λβ n β MT o λκ M Q ds M. REMARK.: In accordance with Eq. 6 the stresses strains are taken as constant quantities which are therefore independent of the arc coordinate s. Since T λκ M Q ds M = δ κλ ε 3ρλ r ρ T λκ M Q ds M = which is the moment about the origin of the stresses due to a unit force applied at Q one can write I κ + I κ = c κ lim ε 3ρλω ξ ρ +r ρ T λκ M Q ds M = L R c κ lim ε 3ρλ ωξ ρ T λκ M Q ds M + L R ε 3ρλ ω r ρ T λκ M Q ds M = c κ + ε 3ρκ ωξ ρ which clearly shows that reproduces the rigid body motion. Determining the limits 3 I κ 4 I κ requires long formal transformations. For this reason we confine ourselves to basic argument the results of the most important steps. First we have to exp U λκ T λκ into series in terms of to the power etc. These transformations are based on the relations: n α M o = n α 3a r α = x α M o ξ α Q = x α ξ α = n α ξ α 3b ρ α = λ α r + r λ α n + n λ α n ξ + n ξ λ α n + + λ α ξ + ξ } n λ α n + 3c n ln ρ α ln + λ ln α n + n λ α n ξ + n ξ λ α n + n ρα λ α n + λ α n ξ + n ξ + n λ α n + n arctanx + εa = arctan x + + λ α ξ + ξ λ α n + 3d n λ α ξ ξ λ α n + n εa + x 3e a is constant ε is a very small quantity. Consequently λ λ r r λ λ n n arctan λ λ r + arctan r λ λ n + + n λ λ λ λ n n n 3 λ ξ λ n 3 n n ξ n 4 λ λ + n 4 + λ + λ n. 3f n Using relations 3a... 3f for a sufficiently large we have U M Q λ D A ln + λ ln n + n λ n ξ + n ξ λ n + n λ A ln + λ ln n + n λ n ξ + n ξ λ n + 4 n U M Q A D ln + λ λ ln n + n λ n ξ + n ξ λ n + n A ln + λ λ ln α n + n λ α n ξ + n ξ λ α n + 5 n λ U M Q D A A arctan λ n n λ λ λ λ n + + n λ λ n n n 3 λ ξ λ n 3 n n ξ n 4 λ λ + n 4 + λ + λ n 6 n T M Q = D λ A λ n + + λ n ξ + n ξ n λ n + n D λ A λ n + + λ n ξ + n ξ n λ n + n ξ n + ξ n ξ n + ξ n 7 T M Q = D λ A λ n + n λ n n + λ n ξ + n ξ λ n + n λ ξ n ξ n D λ A λ n + n λ n n + λ n ξ + n ξ λ n + n λ ξ n ξ n 8 T M Q = D A λ λ n + n n λ n + λ n ξ + n ξ λ n + n D A λ λ n + n n λ n + λ n ξ + n ξ λ n + n ξ n ξ λ n ξ n ξ λ n 9 4 Per. Pol. Civil Eng.
5 T M Q = D λ A λ n + + λ n ξ + n ξ n λ n + n ξ n ξ n D λ A λ n + + λ n ξ + n ξ n λ n + n ξ n ξ n. 3 Knowledge of the series of U λκ T λκ allows us to calculate the integrals 3 I κ 4 I κ. The following observations are used in the calculations: The outward unit normal on is given by equation ϕ is the polar angle. The arc element on admits the form n a = sin ϕ cos ϕ 3 ds o M = dϕ. 3 3 As the coefficients of is are always an integral integrals of zero value. The integrals we have used are all listed in the Appendix. 4 The structure of the coefficients of to the power zero is similar; but these terms involve ξ α the trigonometric integrals that constitute these coefficients are not necessarily equal to zero. After performing the integrations for 3 I κ we have 3 I κ = lim t λρ n ρ M U λκ M Q ds M = = t ξ I κ + ξ I κ + t ξ I κ + ξ I κ + + t ξ 33 I κ + ξ 33 I κ + t } 3 I = π D A λ + A λ λ + λ } 3 I = π D A λ + A λ λ + λ ξ I κ + ξ I κ 33 3 I = 3 I = 33 λ λ I = π D A A λ + λ + 33 I = 34 I = 33 I 3 I = 33 I 3 I = 34 I 33 I = π D A λ + A } λ + } 34 A I = π D λ λ + A λ λ + In the same way we obtain for 4 I κ that 34 I = 3 I = 3 I = 33 I = 34 I =. 4 I κ = lim e λβ n β MT o λκ M Q ds M = = e ξ I κ + ξ I κ + e ξ I κ + ξ I κ + + e ξ 43 I κ + ξ 43 I κ + e 34a 34b 34c 34d 34e 34f 34g 34h ξ I κ + ξ I κ 35 4 I = π D A λ λ + A 4 λ I = π D A λ A 43 λ I = π D A λ λ 44 λ I = π D A λ λ 4 λ I = π D A λ 4 λ I = π D A 43 I = π D A } λ λ + } λ λ } A λ λ λ λ A } λ λ λ λ A λ } λ λ λ } λ λ A λ } λ λ λ + A λ + 44 λ I = π D A λ A } λ λ 4 I = 36a 4 I = 36b 43 I = 36c 44 I = 36d 4 I = 4 I = 36e 36f 43 I = 36g 44 I =. 36h Substituting the integrals with zero value from Eqs. 34a...34h 36a...36h into Eqs then utilizing the Hooke law 7 3 we obtain 3 I + 4 I = s e s e s s s e ξ s e s e s s s e ξ 3 I + 3 I + 33 I ξ I e ξ I + ξ I 4 ξ I I + ξ I + 43 I 3 I + 4 I = s e s e s s s e ξ + s e s e s s s e e ξ I + ξ I ξ 3 I + 3 I + 33 I + ξ 34 I e 37a ξ 4 I + ξ 4 I 4 ξ I I + ξ I + 43 I e ξ I + ξ I. 37b Using Eqs. 7 to one can prove that the following equations hold the formal calculations are presented in Appendix B.: s s s s s s s s 3 I + s 34 I 4 I = s s s 3 I + 33 I 4 I 43 I = 3 I s s s s 34 I 44 I = 38a BEM Formulation plane orthotropic bodies 7 5 5
6 s s s s s s s s 3 I + s 34 I 4 I = s s s 3 I + 33 I 4 I 43 I = 3 I s s s s 34 I 44 I =. 38b If we substitute Eqs. 38ab into equtions 37ab we find that 3 I κ + 4 I κ = e κβ ξ β. 39 Neglecting the rigid body motion i.e. setting I κ + I κ to zero utilizing 39 we obtain I κ = I κ = I κ + I κ + 3 I κ + 4 I κ = e κβ ξ β. Then the first modified Somigliana formula immediately follow from Eqs. 9 : u κ Q = e κβ ξ β Q + t λ M U λκ M Q u λ M T λκ M Q ds M Q A e 4 If Q = Q is on nothing changes concerning the limit of the integral taken on. Consequently C κρ u ρ Q = e κβ ξ β Q + t λ M U λκ M Q u λ M T λκ M Q ds M Q = Q 4 C κρ = δ κρ / if the contour is smooth at Q. This integral equation is that of the direct method or the second Somigliana formula for exterior regions. If Q is inside the contour this region is referred to as A i then it is easy to show that = e κβ ξ β Q+ t λ M U λκ M Q u λ M T λκ M Q ds M Q = Q A i 4 which is the third Somigliana formula for exterior regions. Using the formulae set up for the strains in Appendix C the Hooke law 7 we can calculate the stresses: t αβ Q = t αβ + t λ M ˆD λαβ M Q ds M L o u λ M S λαβ M Q ds M Q A e 43 ˆD λ = c D λ + c D λ ˆD λ = c D λ + c D λ ˆD λ = c 66 D λ = ˆD λ ˆ S λ = c S λ + c S λ ˆ S λ = c S λ + c S λ Sˆ λ = c 66 S λ = Sˆ λ Behaviour at infinity Our goal in this section is to compute the limit of representation 4 as Q. This will lead to a characterization of the asymptotic behaviour of u κ Q. If this behaviour is the same what we have assumed i.e. if the limit coincides with 6 provided that in the latter the rigid body motion is neglected then we confirm that the results of the previous section are correct. It is clear from representation 4 that it is sufficient to show that the following relations hold lim Qξ ξ t λ M U λκ M Qds M = lim Qξ ξ u λ M T λκ M Qds M =. 45 In order to find the limit of the above integrals we have to set up asymptotic relations for the fundamental solutions U λκ M Q T λκ M Q if Q. Using the notations introduced in Fig. as well as Eqs. 6 we have r α M o Q = x α M o ξ α Q = x α ξ α = 46a ˆR ˆn α x α ˆR ˆn α ˆn α = 46b ˆR ρ α = λ α r + r ˆR λ α ˆn + ˆn 46c ln ρ α ln ˆR + λ ln α ˆn + ˆn. 46d Substituting equations 46bbc into 3 then performing Fig.. x n O x ds M o r QM o r OQ = R some manipulations we obtain the following asymptotic relations for the fundamental solution of order one: U M Q D λ A λ A ln ˆR 47a U M Q = U M Q λ λ ˆn ˆn D A A arctan λ λ ˆn + 47b ˆn U M Q D Q n A A ln ˆR. λ λ 47c It is obvious that asymptotically U λκ M Q U λκ Q i.e. the kernel in integral 45 is independent of M. 6 Per. Pol. Civil Eng.
7 Consequently lim t λ M U λκ M Qds M = Qξ ξ = lim U λκq t λ M ds M Qξ ξ L } o } resultant =. 48 By repeating the line of thought leading to the asymptotic relations 47abc for the fundamental solutions of order two we obtain T M Q Ḓ λ A λ A R λ ˆn + ˆn λ ˆn + = T Q 49a ˆn ˆR T M Q 4 Iintegral Equations of the First Kind in Plane Elasticity Quarterly Ḓ λ A R λ λ ˆn + λ A ˆn λ λ ˆn + ˆn ˆn = T of Applied Mathematics LIII 995 no Q 5 Gradstein IS Ryzhik IM Table of Integrals Series Products Academic New York 98. ˆn ˆR 49b 6 Table of Integrals Series Products Moscow Nauka Pub. in Ḓ λ λ A λ λ A R λ ˆn + ˆn λ ˆn + ˆn T M Q ˆn ˆn = T Q ˆR 49c T M Q Ḓ λ A λ A R λ ˆn + ˆn λ ˆn + = T Q 49d ˆn ˆR T λκ = T λκ Q is defined by the above equations. Since for large ˆR the kernel T λκ M Q is independent of M tends to zero it follows that lim u λ M T λκ M Qds M = Qξ ξ T λκ Q = lim u λ M ds M =. 5 Qξ ξ ˆR This verifies that the asymptotic behaviour of the displacement representation 4 is as expected. Under this condition the strain energy density is bounded although the strain energy is not there is no need to replace the exterior region by a finite one if a constant stress condition is prescribed at infinity. This can be an advantage if one considers an infinite plane with holes or cracks subjected to constant stresses at infinity an attempt is made to determine the stresses in finite. Existing codes can be modified easily to perform computations. References Banarjee PK Butterfield R Boundary Element Methods in Engineering Science Mir Moscow 984. Banarjee PK The Bondary Element Methods in Engineering McGraw Hill New York Consta C The Boundary Iintegral Equation Method in Plane Elasticity Proceedings of the American Mathematical Society no. Russian Huang L Sun X Liu Y Cen Z Parameter identification for twodimensional orthotropic material bodies by the boundary element method Engineering Analysis with Boundary Elements 8 4 no Vable M Sikarskie DL Stress analysis in plane orthotropic material by the boundary element method Int. J. Solids Structures no.. 9 Rizzo RJ Shippy DJ A method for stress determination in plane anisotropic elastic bodies J. Composite Materials Schiavone P Chong-Quing Ru On the exterior mixed problem in plane elasticity Mathematics Mechanics of Solids Szeidl Gy Boundary integral equations for plane problems remark to the formulation for exterior regions Publications of the University of Miskolc Series D Natural Sciences Mathematics no Boundary integral equations for plane problems in terms of stress functions of order one Journal of Computational Applied Mechanics no Concluding remarks For the sake of completeness note that Consta 3 gives an asymptotic expansion for the displacements at infinity which ensures the validity of the Betti formula for exterior regions isotropic bodies. Under this condition the total strain energy stored in the region is bounded. In addition uniqueness existence proofs are easy to give. We have modified the Somigliana formulas for exterior regions of orthotropic bodies by assuming that the strains are constants accordingly the displacements are linear at infinity. BEM Formulation plane orthotropic bodies 7 5 7
8 A Trigonometric integrals When determining the integrals 3 I κ + 4 I κ we have used the following trigonometric integrals: n ln ds = ln n ln n ln π λ α n + n ds = λ α n + n ds = n λ α n + ds = n n n λ α n + ds = n n λ α n + ds = n π π π n3 n π π π cos ψ dψ = π π λ n + n λ n + ds = n λ n + n λ n + ds = n n n3 λ n + n λ n + ds = n n4 n ln ds = ln π sin ψ dψ = 5 cos ψ ln λ α cos ψ + sin ψ dψ = 5 sin ψ ln λ α cos ψ + sin ψ dψ = 53 cos ψ λ α cos ψ + sin ψ dψ = π λ α + λ α 54 cos ψ sin ψ λ α cos ψ + sin dψ = 55 ψ sin ψ λ α cos ψ + sin ψ dψ = π π π λa + 56 cos 3 ψ sin ψ λ cos ψ + sin ψ λ cos ψ + sin dψ = 57 ψ cos ψ sin 3 ψ λ cos ψ + sin ψ λ cos ψ + sin dψ = 58 ψ cos 4 ψ λ cos ψ + sin ψ λ cos ψ + sin ψ dψ = π λ + λ + λ λ λ + λ λ + λ + 59 λ n + n λ n + ds = n n n cos ψ sin ψ λ cos ψ + sin ψ λ cos ψ + sin ψ dψ = π λ + λ λ + λ + 6 n4 n 4 λ λ + n 4 + λ + λ n n ds = sin 4 ψ λ cos ψ + sin ψ λ cos ψ + sin ψ dψ = π λ + λ + λ λ λ + λ λ + λ +. 6 The integrals detailed above can be checked either by using Maple 9.5 or the Table of Integrals 5 6 by Gradstein Ryzhik. B Proof of equation 38a B.. First consider equation 38a. Substituting the integrals 3 I 34 I 4 I see equations 34a 34d 36a we have π s s s D s A λ + A λ λ + λ } λ + s D A A λ λ + λ + π D A λ λ + + π D A λ λ + =. 6 8 Per. Pol. Civil Eng.
9 After multiplying throughout by the common denominator s s s λ λ some further manipulations we obtain s A λ λ + A λ λ + + s A A λ λ A λ λ + } λ λ If we divide throughout by λ λ substitute 9 then we find + s s s + λ λ s s s s A λ λ + λ λ =. 63 s λ s λ λ + s λ s s λ + s s s + s λ s λ λ + s λ s s λ + s s s λ λ s s s s λ + λ + s s λ s s λ s λ λ =. 64 Here we can again divide throughout by λ λ. Hence λ 3 λ 3 s 3 + λ s s + λ s s λ λ s3 λ λ s 3 λ λ 3 s 3 λ 3 λ s 3 + λ s s + λ s s s s s + λ λ s s + λ λ s s λ λ s s = 65 Consequently λ s s + λ s s = s s λ + λ = s s s + = s + s s s λ λ s3 λ λ s 3 = λ λ s 3 λ + λ = s s 3 s + = s + s s. s s Here λ 3 λ 3 s 3 λ λ 3 s 3 λ 3 λ s 3 + λ s s + λ s s s s s + λ λ s s + λ λ s s =. λ λ 3 s 3 λ 3 λ s 3 = s3 Therefore λ λ 3 + λ 3 λ = s 3 λ λ λ + λ = s 3 s λ + λ = s s s λ + λ λ s s + λ s s = s s λ + λ. λ 3 λ 3 s 3 s s s + λ λ s s + λ λ s s = λ 3 λ 3 s 3 + s s λ λ s λ 3 λ 3 + λ λ s s = s s s λ λ s λ 3 λ 3 s 3 + λ λ s s = s s s λ λ s s s s λ λ s =. B.. Consider equation 38a. If we substitute the integrals 3 I 33 I 4 I 43 I see equations 34b 34c 36b 36c perform some manipulations we obtain 4π D A λ λ + A λ λ + A λ λ A λ + λ + + A λ λ + A After substituting D multiplying throughout by the common denominator we have s λ λ λ + λ + A λ λ + A λ λ + A A λ λ A λ + λ A λ + λ λ λ + =. 66 s λ λ λ + λ + =. 67 BEM Formulation plane orthotropic bodies 7 5 9
10 If we make use of equation 9 perform some transformations we find s λ s λ + s λ s λ λ s λ s λ + s λ s λ + λ + s λ s s λ s λ λ s λ λ λ + λ + =. If we factor the above expression we have λ λ s + s λ λ + λ + λ s λ + λ s =. The first factor is obviously different from zero. Hence we should investigate the other two. We shall start with the last one: s λ + λ s = s + s s + s = s + s + =. Since the above expression is zero equation 38a is truly satisfied. B.3. The third equation to consider is 38a 3. If we substitute the integrals 3 I 34 I 4 I - see equations 34a 34d 36d we obtain π D s s s s A λ λ + A λ λ λ + s A A λ + λ + λ + + π D λ A λ + λ A λ + =. 68 Let us multiply throughout by the common denominator s s s λ + λ + substitute D. We have s A λ λ + A λ λ + + s A A λ λ + If we make use of equations 9 we find s s s A λ λ + A λ λ + =. 69 s s λ s λ λ + s s λ s λ λ + + s s λ s s λ s λ λ + s s s s λ s λ λ + s s s s λ s λ λ + =. 7 After factoring the above expression we obtain the product: λ λ s + s λ + λ + λ + λ + λ λ λ λ s s =. The first factor is apparently not equal to zero. Therefore we have to consider the other two. The last one clearly vanishes: s + λ λ s = s + s s = s + s =. s Consequently equation 38a 3 is also satisfied. C Proof of equation 38b C.. Third consider equation 38b. Upon substitution of the integrals 3 I 34 I 4 I see equations 34e 34h 36e we have π D s A A λ λ A s s s λ + λ + + s λ λ + A λ λ + A π D λ λ + A λ λ + =. 7 Per. Pol. Civil Eng.
11 After multiplying throughout by the common denominator s s s λ λ λ + λ + making use of equations 9 we obtain λ λ s 3 λ + λ s s + s 3 + λ 3 λ 3 s 3 λ s s λ s s s s s λ λ s s + λ λ 3 s s λ λ 3 s s } } + λ 3 λ s s λ 3 λ s s } } +λ λ s s λ s s s λ s s s + +λ λ 3 s s + λ 3 λ s s + λ λ s s λ λ s s } } Let us divide throughout by λ λ cancel the terms braced. If in addition we take into account that we shall have s s + s 3 = s3 + s3 = =. λ s 3 + λ s 3 + λ 3 λ 3 s 3 λ s s in which λ s s s s s λ λ s s + λ λ s s λ λ 3 s s + + λ λ 3 s s λ 3 λ s s + λ 3 λ s s + λ λ s s λ s s s λ s s s + λ λ 3 s s + λ 3 λ s s λ λ s s s = λ 3 λ 3 s 3 s s s λ λ s s + λ λ s s λ s s s λ s s s + λ λ 3 s s + λ 3 λ s s = λ λ 3 s s + λ 3 λ s s λ s s s λ s s s = Therefore = s s λ λ λ + λ s s s λ + λ = s s λ + λ s s s λ + λ = s λ 3 λ 3 s 3 + λ λ s s s s s λ λ s s = s = s s s λ + λ s s s λ + λ =. = s λ λ s λ λ + s s s s + s λ λ = s s λ λ + s s s s + s λ λ = s s = s s λ λ + s s s s s s s s λ λ =. C.. Consider equation 38b. If we substitute the integrals 3 I 33 I 4 I 43 I see equations 34b 34c 36b 36c perform some manipulations we obtain π D A λ + A λ + A A λ λ λ + λ + π D A λ + A λ + π D A λ λ + A λ =. 7 λ + BEM Formulation plane orthotropic bodies 7 5
12 Let us multiply throughout by the common denominator λ + λ + make some manipulations. In this way we have A λ + A λ + If simplify?? then we have A A λ λ A λ + A λ + A λ λ + A λ λ + λ λ s λ + λ + =. 73 s λ s λ + s λ s λ + s λ s s λ s λ λ + s λ s λ + λ + + s λ s λ + λ λ λ s λ + λ + =. 74 If we factoring this expression we shall get λ λ s + λ s + λ s + λ λ s s + λ s + λ s =. The first factor is evidently different from zero. Hence we have to investigate the other two. We shall start with the last one s + λ s + λ s = s +s λ +λ = s +s s + s = s +s + =. C.3. Third consider equation 38b 3.Upon substitution of the integrals 3 I 34 I 4 I see equations 34e 34h 36e we have π D s A A λ λ s s s λ + λ + + A λ λ + A λ λ + A π D λ + A =. 75 λ + Then multiplying throughout by the common denominator s s s λ λ λ + λ + making some manipulations we obtain s A A λ λ λ λ s A λ λ + A λ λ + s s s λ λ A λ + A λ + λ λ s s s s λ λ λ + λ + =. 76 If we substitute equation 9 we have s s λ s s λ + s s λ λ s s λ λ + s s λ s s λ + + s s λ λ s s λ λ + s s λ λ 3 s s λ 3 λ s λ λ 3 + s λ 3 λ =. 77 After factoring this equations we shall find λ λ s + λ s + λ s + λ λ s s λ λ s =. One can check easily that the third factor is zero: Consequently 38b 3 is fulfilled. s λ λ s = s s s = s s =. s Per. Pol. Civil Eng.
13 D Formulae for strains Making use of equations from the kinematic equations we obtain e αβ = t λ M D λαβ M Q ds M u λ M S λαβ M Q ds M Q A e 78 D λ = U λ S λ = T λ If we introduce the notations D λ = U λ + U λ = D λ D λ = U λ S λ = T λ + T λ = D λ S λ = T λ. 79 b = λ ρ k = λ 3 A ρ b = λ ρ k = λ 3 A ρ c = λ λ r + r d = λ λ r A A f = f λ 3 = r + r λ 3 r + r for the derivatives in equation 79 we have U = D b λ A r b λ A r U = D A A + d r λ λ dr cd = U U = D b A r b A r U = D b A r b A r U = D A A dr + d r r c λ λ r = U A U = D r ρ A r λ ρ λ 8 8 T = D T = D T = D T = D b λ A r ρ b λ A r ρ b λ A r ρ b λ A r ρ b λ A r ρ r n + r n D b A b A n b λ A r ρ b λ A r ρ b λ A r ρ r n D b A b A n D f λ r f λ r r n r n D b λ A b λ A n D r n + r n D b A b A n b A r T = D b A r r n + r n D b A b A n T = D T = D T = D ρ b A r ρ b λ A r ρ b A r ρ ρ b A r ρ b λ A r r n D f r f r r n + D f f n b λ A r r n D ρ b A r ρ b A r ρ r n + r n D b A b A n. ρ b λ A r ρ r n b A r ρ r n + D b A b A n 8a 8b BEM Formulation plane orthotropic bodies 7 5 3
2 Basic Equations in Generalized Plane Strain
Boundary integral equations for plane orthotropic bodies and exterior regions G. Szeidl and J. Dudra University of Miskolc, Department of Mechanics 3515 Miskolc-Egyetemváros, Hungary Abstract Assuming
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
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 informationE&M. 1 Capacitors. January 2009
E&M January 2009 1 Capacitors Consider a spherical capacitor which has the space between its plates filled with a dielectric of permittivity ɛ. The inner sphere has radius r 1 and the outer sphere has
More informationEFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES
EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC
More informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationComputable representations of influence functions for multiply connected Kirchoff plates
Computable representations of influence functions for multiply connected Kirchoff plates Y.A. Melnikov Department of Mathematical Sciences Middle Tennessee State University, USA Abstract A modification
More informationThe Solution of Torsion Problem for the Bars with Irregular Cross Sections by Boundary Element Method
International Journal of Mechanics and Applications 04, 4(): 9-49 DOI: 0.59/j.mechanics.04040.0 The Solution of Torsion Problem for the Bars with Irregular Cross Sections by Boundary Element Method Hakan
More informationTHE ANALYTICAL SOLUTIONS OF THE BOUNDARY-VALUE PROBLEMS BY THE METHOD OF FACTORIZATION
Materials Physics and Mechanics (5) 47-5 Received: March 7 5 THE ANALYTICAL SOLUTIONS OF THE BOUNDARY-VALUE PROBLEMS BY THE METHOD OF FACTORIZATION VA Babeshko * OV Evdokimova OM Babeshko Kuban State University
More informationMathematical Preliminaries
Mathematical Preliminaries Introductory Course on Multiphysics Modelling TOMAZ G. ZIELIŃKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Vectors, tensors, and index notation. Generalization of the concept
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationRocking behaviour of a rigid foundation with an arbitrary embedment
Rocking behaviour of a rigid foundation with an arbitrary embedment H. Moghaddasi Islamic Azad University, Qazvin, Iran. M. Kalantari University of Tehran, Tehran, Iran N. Chouw The University of Auckland,
More informationDegenerate scale problem for plane elasticity in a multiply connected region with outer elliptic boundary
Arch Appl Mech 00) 80: 055 067 DOI 0.007/s0049-009-0357-3 ORIGINAL Y. Z. Chen X. Y. Lin Z. X. Wang Degenerate scale problem for plane elasticity in a multiply connected region with outer elliptic boundary
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 informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationMOHR S CIRCLES FOR NON-SYMMETRIC STRESSES AND COUPLE STRESSES
CRNOGORSKA AKADEMIJA NAUKA I UMJETNOSTI GLASNIK ODJELJENJA PRIRODNIH NAUKA, 16, 2005. QERNOGORSKAYA AKADEMIYA NAUK I ISSKUSTV GLASNIK OTDELENIYA ESTESTVENNYH NAUK, 16, 2005. THE MONTENEGRIN ACADEMY OF
More informationClassical Oscilators in General Relativity
Classical Oscilators in General Relativity arxiv:gr-qc/9709020v2 22 Oct 2000 Ion I. Cotăescu and Dumitru N. Vulcanov The West University of Timişoara, V. Pârvan Ave. 4, RO-1900 Timişoara, Romania Abstract
More informationWeak solutions to anti-plane boundary value problems in a linear theory of elasticity with microstructure
Arch. Mech., 59, 6, pp. 519 539, Warszawa 2007 Weak solutions to anti-plane boundary value problems in a linear theory of elasticity with microstructure E. SHMOYLOVA 1, S. POTAPENKO 2, A. DORFMANN 1 1
More informationACCELERATION WAVES IN THE VON KARMAN PLATE THEORY 1
ACCELERATION WAVES IN THE VON KARMAN PLATE THEORY V. VASSILEV and P. DJONDJOROV Institute of Mechanics Bulgarian Academy of Sciences Acad. G. Bontchev St., Block 4, 3 Sofia, BULGARIA arxiv:math-ph/000606v2
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationBeam Models. Wenbin Yu Utah State University, Logan, Utah April 13, 2012
Beam Models Wenbin Yu Utah State University, Logan, Utah 843-4130 April 13, 01 1 Introduction If a structure has one of its dimensions much larger than the other two, such as slender wings, rotor blades,
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationChapter 9 Uniform Circular Motion
9.1 Introduction Chapter 9 Uniform Circular Motion Special cases often dominate our study of physics, and circular motion is certainly no exception. We see circular motion in many instances in the world;
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationOn some exponential decay estimates for porous elastic cylinders
Arch. Mech., 56, 3, pp. 33 46, Warszawa 004 On some exponential decay estimates for porous elastic cylinders S. CHIRIȚĂ Faculty of Mathematics, University of Iaşi, 6600 Iaşi, Romania. In this paper we
More informationAPPENDIX 2.1 LINE AND SURFACE INTEGRALS
2 APPENDIX 2. LINE AND URFACE INTEGRAL Consider a path connecting points (a) and (b) as shown in Fig. A.2.. Assume that a vector field A(r) exists in the space in which the path is situated. Then the line
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
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 information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationStress of a spatially uniform dislocation density field
Stress of a spatially uniform dislocation density field Amit Acharya November 19, 2018 Abstract It can be shown that the stress produced by a spatially uniform dislocation density field in a body comprising
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 08,000.7 M Open access books available International authors and editors Downloads Our authors
More informationJournal of Solid Mechanics and Materials Engineering
and Materials Engineering Analysis of In-Plane Problems for an Isotropic Elastic Medium with Many Circular Holes or Rigid Inclusions Mutsumi MIYAGAWA, Jyo SHIMURA, Takuo SUZUKI and Takanobu TAMIYA Dept.of
More information3 Chapter. Gauss s Law
3 Chapter Gauss s Law 3.1 Electric Flux... 3-2 3.2 Gauss s Law (see also Gauss s Law Simulation in Section 3.10)... 3-4 Example 3.1: Infinitely Long Rod of Uniform Charge Density... 3-9 Example 3.2: Infinite
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationPROPAGATION OF A MODE-I CRACK UNDER THE IRWIN AND KHRISTIANOVICH BARENBLATT CRITERIA
Materials Science, Vol. 39, No. 3, 3 PROPAGATION OF A MODE-I CRACK UNDER THE IRWIN AND KHRISTIANOVICH BARENBLATT CRITERIA I. I. Argatov, M. Bach, and V. A. Kovtunenko UDC 539.3 We study the problem of
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationAnalytical solutions for some defect problems in 1D hexagonal and 2D octagonal quasicrystals
PRAMANA c Indian Academy of Sciences Vol. 70 No. 5 journal of May 008 physics pp. 911 933 Analytical solutions for some defect problems in 1D hexagonal and D octagonal quasicrystals X WANG 1 and E PAN
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationAppendix A: Data completeness
Appendix A: Data completeness The fraction of relps in a given constant-resolution sphere that are lost in the blind region lie in two polar caps intersecting the rotation axis (z). That is, all points
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More information1.4 LECTURE 4. Tensors and Vector Identities
16 CHAPTER 1. VECTOR ALGEBRA 1.3.2 Triple Product The triple product of three vectors A, B and C is defined by In tensor notation it is A ( B C ) = [ A, B, C ] = A ( B C ) i, j,k=1 ε i jk A i B j C k =
More informationResearch Article Divergent Integrals in Elastostatics: General Considerations
International Scholarly Research Network ISRN Applied Mathematics olume 2011, Article ID 726402, 25 pages doi:10.5402/2011/726402 Research Article Divergent Integrals in Elastostatics: General Considerations..
More informationDr. Parveen Lata Department of Basic and Applied Sciences, Punjabi University, Patiala, Punjab, India.
International Journal of Theoretical and Applied Mechanics. ISSN 973-685 Volume 12, Number 3 (217) pp. 435-443 Research India Publications http://www.ripublication.com Linearly Distributed Time Harmonic
More informationCorrections to the First Printing: Chapter 6. This list includes corrections and clarifications through November 6, 1999.
June 2,2 1 Corrections to the First Printing: Chapter 6 This list includes corrections and clarifications through November 6, 1999. We should have mentioned that our treatment of differential forms, especially
More informationThe Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density
Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationBohr & Wheeler Fission Theory Calculation 4 March 2009
Bohr & Wheeler Fission Theory Calculation 4 March 9 () Introduction The goal here is to reproduce the calculation of the limiting Z /A against spontaneous fission Z A lim a S. (.) a C as first done by
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationElectric fields in matter
Electric fields in matter November 2, 25 Suppose we apply a constant electric field to a block of material. Then the charges that make up the matter are no longer in equilibrium: the electrons tend to
More informationHooke s law and its consequences 1
AOE 354 Hooke s law and its consequences Historically, the notion of elasticity was first announced in 676 by Robert Hooke (635 73) in the form of an anagram, ceiinosssttuv. He explained it in 678 as Ut
More informationarxiv:gr-qc/ v1 15 Apr 1997
Indeterministic Quantum Gravity and Cosmology VII. Dynamical Passage through Singularities: Black Hole and Naked Singularity, Big Crunch and Big Bang Vladimir S. MASHKEVICH 1 arxiv:gr-qc/9704038v1 15 Apr
More informationThe Basic Properties of Surface Waves
The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationEffect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities
Effect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities Hany El Naggar, Ph.D., P. Eng. and M. Hesham El Naggar, Ph.D., P. Eng. Department of Civil Engineering
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
More information2. Geometrical Preliminaries
2. Geometrical Preliminaries The description of the geometry is essential for the definition of a shell structure. Our objective in this chapter is to survey the main geometrical concepts, to introduce
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker
THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092
More informationThe double layer potential
The double layer potential In this project, our goal is to explain how the Dirichlet problem for a linear elliptic partial differential equation can be converted into an integral equation by representing
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationIntroduction to the J-integral
Introduction to the J-integral Instructor: Ramsharan Rangarajan February 24, 2016 The purpose of this lecture is to briefly introduce the J-integral, which is widely used in fracture mechanics. To the
More informationLinear elastic analysis of thin laminated beams with uniform and symmetric cross-section
Applied and Computational Mechanics 2 (2008) 397 408 Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section M. Zajíček a, a Faculty of Applied Sciences, UWB in Pilsen,
More informationEFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL
AMS Subject Classification Index: 74Q20 EFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL Irini DJERAN-MAIGRE and Sergey V.KUZNETSOV INSA de Lyon, URGC 34 Avenue des Arts, 69621
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationAnalytical formulation of Modified Upper Bound theorem
CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are
More informationLECTURE 10: REVIEW OF POWER SERIES. 1. Motivation
LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the
More informationTHE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng. 46, 1883}1895 (1999) THE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY ZHOU SHENJIE*, CAO ZHIYUAN AND SUN
More informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
More informationSIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas
More information1 The Lagrange Equations of Motion
1 The Lagrange Equations of Motion 1.1 Introduction A knowledge of the rudiments of dynamics is essential to understanding structural dynamics. Thus this chapter reviews the basic theorems of dynamics
More informationChapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING )
Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) 5.1 DEFINITION A construction member is subjected to centric (axial) tension or compression if in any cross section the single distinct stress
More informationCE-570 Advanced Structural Mechanics - Arun Prakash
Ch1-Intro Page 1 CE-570 Advanced Structural Mechanics - Arun Prakash The BIG Picture What is Mechanics? Mechanics is study of how things work: how anything works, how the world works! People ask: "Do you
More informationChapter 1. Condition numbers and local errors in the boundary element method
Chapter 1 Condition numbers and local errors in the boundary element method W. Dijkstra, G. Kakuba and R. M. M. Mattheij Eindhoven University of Technology, Department of Mathematics and Computing Science,
More informationJun-ichi Takada
Wave Theory II Numerical Simulation of Waves 2) Representation of Wave Function by Using Green Function I) Green s Theorem and Free Space Green Function Jun-ichi Takada takada@ide.titech.ac.jp)) This lecture
More informationCHAPTER -6- BENDING Part -1-
Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationIntroduction to fracture mechanics
Introduction to fracture mechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 6-9 November, 2017 Institute of Structural Engineering, ETH Zu rich November 9, 2017 Institute
More informationEXAMINATION OF AN OPTIMIZED REPLACEABLE CUTTING TOOTH OF EXCAVATOR
Geosciences and Engineering, Vol. 1, No. (01), pp. 337 34. EXAMINATION OF AN OPTIMIZED REPLACEABLE CUTTING TOOTH OF EXCAVATOR ZOLTÁN VIRÁG 1 SÁNDOR SZIRBIK 1 Department of Geotechnical Equipment, University
More informationCoordinate Systems and Canonical Forms
Appendix D Coordinate Systems and Canonical Forms D.1. Local Coordinates Let O be an open set in R n. We say that an n-tuple of smooth realvalued functions defined in O, (φ 1,...,φ n ), forms a local coordinate
More informationGeneralization of the Dick Model
Generalization of the Dick Model arxiv:hep-ph/0092v 7 Sep 200 M. Ślusarczyk and A. Wereszczyński Institute of Physics, Jagellonian University, Reymonta, Kraków, Poland July 5, 208 Abstract We discuss a
More information1. (3) Write Gauss Law in differential form. Explain the physical meaning.
Electrodynamics I Midterm Exam - Part A - Closed Book KSU 204/0/23 Name Electro Dynamic Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to tell about
More informationLecture notes Models of Mechanics
Lecture notes Models of Mechanics Anders Klarbring Division of Mechanics, Linköping University, Sweden Lecture 7: Small deformation theories Klarbring (Mechanics, LiU) Lecture notes Linköping 2012 1 /
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationSSPH basis functions for meshless methods, and comparison of solutions with strong and weak formulations
Comput Mech 8 4:57 545 DOI.7/s466-7-9- ORIGINAL PAPER SSPH basis functions for meshless methods, and comparison of solutions with strong and weak formulations R. C. Batra G. M. Zhang Received: 7 February
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More information, is the time to fatigue failure. According to [1]. The equation of plastic flow theory accept in the form (1)
Fatigue Failure of an Annular Plate Under the Action of Pulsating Moment Pressure along the internal contour LatifKh. T., Nagiyeva N. M. National Academy of Sciences of Azerbaijan, Institute of Mathematics
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More information