BOUNDARY ELEMENT ANALYSIS OF THIN STRUCTURAL PROBLEMS IN 2D ELASTOSTATICS
|
|
- Hester Wood
- 5 years ago
- Views:
Transcription
1 Journal of Marne Scence and Technology, Vol. 19, No. 4, pp (011) 409 BOUNDARY ELEMENT ANALYSIS OF THIN STRUCTURAL PROBLEMS IN D ELASTOSTATICS Yao-Mng Zhang,, Yan Gu, and Jeng-Tzong Chen Key words: boundary element method, nearly sngular ntegrals, thn structural problems, elastcty. ABSTRACT Thn structures, such as thn flms and coatngs, have been wdely desgned and utlzed n many ndustres recently. However, the wdespread expermental research n thn structural problems underles a general lack of modelng efforts whch can accurately and effcently predct ther performance. In ths paper, the boundary element method (BEM) based on elastcty theory s developed for two-dmensonal (D) thn structures wth the thckness to length rato n the mcro (10-6 ) or nano (10-9 ) scales. The BEM-based approach proposed n ths paper s constructed usng a combnaton of the regularzed ndrect boundary ntegral equatons (BIEs) and a general nonlnear transformaton whch can elmnate the nearly sngular propertes of the ntegral kernels. For the test problems studed, very promsng results are obtaned when the thckness to length rato s n the orders between 10-6 and 10-9, whch s suffcent for modelng most thn structures n the mcro- or nano-scales. I. INTRODUCTION Structures made from plate and shell assembles, also known as thn-walled structures, are mportant n many areas of engneerng. The effcent soluton of such structures has wde applcaton, from the desgn of more conventonal cvl engneerng systems to more sophstcated structures employed, for nstance, by the arcraft ndustry. However, the wdespread expermental research n thn structures underles a general lack of modelng efforts whch can accurately and Paper submtted 0/05/10; revsed 05/04/10; accepted 05/08/10. Author for correspondence: Yao-Mng Zhang (e-mal: zymfc@163.com). Insttute of Appled Mathematcs, Shandong Unversty of Technology, Zbo, Shandong Provnce, Chna. Industral Equpment State Key Laboratng of Structures Analyss, Dalan Unversty of Technology, Dalan 11604, Chna. Department of Harbor and Rver Engneerng, Natonal Tawan Ocean Unversty, Keelung 04, Tawan. effcently predct ther performance. For computatonal models of thn structures or thn shapes n structures, two numercal methods can be employed, namely, the fnte element method (FEM) and the boundary element method (BEM). In the last three decades, the FEM based on plate and shell theores has been a successful tool for the analyss of 3-D thn structures such as plates, shells and layered composte structures to study ther deformaton and stress n macro-scale. However, the FEM element count ncreases dramatcally for thn structures due to aspect rato lmtatons, and t was found out that, as shown n Ref. [14], the number of the fnte elements were so large that the task quckly exceed the capacty of the computer used as the thckness decreases. The BEM s a powerful and effcent computatonal method f ntegrals are evaluated accurately, and the BEM based on the elastcty theory s n general more accurate n stress analyss of structures. Ths accuracy wll be mantaned n the analyss of thn structures as well, f the BEM s mplemented correctly to deal wth the dffcultes assocated wth thn structures. Studes show that the conventonal boundary element method (CBEM) usng the standard Gaussan quadrature fals to yeld relable results for thn-walled structures. The major reason for ths falure s that the kernels ntegraton presents varous orders of near sngulartes, owng to the mesh on one sde of the thn-body beng too close to the mesh on the opposte sde. Nearly sngular ntegrals are not sngular n the sense of mathematcs. However, from the pont of vew of numercal ntegratons, these ntegrals can not be calculated accurately by usng the conventonal numercal quadrature snce the ntegrand oscllates very fercely wthn the ntegraton nterval. Other than the nearly sngular ntegral, many drect and ndrect algorthms for sngular ntegral have been developed and used successfully [1-4, 6, 7, 10, 11, 19, 3, 4]. Therefore, the key pont n achevng the requred accuracy and effcency of the BEM s not the sngular ntegral but the nearly sngular ntegral. Although that dffculty can be overcome by usng very fne meshes, the process requres too much preprocessng and CPU tme. In the past decades, tremendous effort s devoted to derve convenent ntegral forms or sophstcated computatonal technques for calculatng the nearly sngular ntegrals. These
2 410 Journal of Marne Scence and Technology, Vol. 19, No. 4 (011) proposed methods can be dvded on the whole nto two categores: ndrect algorthms and drect algorthms. The ndrect algorthms, whch beneft from the regularzaton deas and technques for the sngular ntegrals, are manly to calculate ndrectly or avod calculatng the nearly sngular ntegrals by establshng new regularzed BIE [5, 13, 17, 0,, 7]. The drect algorthms are calculatng the nearly sngular ntegrals drectly, such as nterval subdvson method [9], specal Gaussan quadrature method [15], exact ntegraton method [18, 8], and varous nonlnear transformaton methods [8, 16, 5, 6]. In a recent study, the above methods have been revewed n detal by Zhang and Sun [6]. Wth the development of the numercal technques for calculaton nearly sngular ntegrals, consderable progress has been made n the applcaton of the BEM to the analyss of thn walled structures n the past few years. Sladek et al. have obtaned amount of the orgnal results n ths feld [1]. Nonsngular ntegral equatons for thn-walled structures are proposed based on a subtracton technque and mathematcal regularzaton. Lu et al. [1, 14] have undertaken a lot of researches n ths feld. The nearly sngular surface ntegrals are transformed nto a sum of weakly sngular ntegrals, and a nonlnear coordnate transformaton s developed for nearly weakly sngular ntegrals. The theory s also appled to nterfacal stress analyss for mult-coatng systems, thermal stress analyss of mult-layer thn flms and coatngs, and also thn pezoelectrc solds. However, n Lu s work only some boundary unknowns can be computed, for example, ths approach only gves the results of radal stresses on the boundary nodes. The tangental stresses on the boundary and the physcal quanttes at nteror ponts need further nvestgaton. Zhou et al. proposed sem-analytcal or analytcal ntegral algorthms to solve -D or 3-D thn body problems [30, 31]. Both boundary and nteror unknowns are computed n ther works whch the geometry boundary s depcted by usng lnear shape functons when nearly sngular ntegrals need to be calculated. Although great progresses have been acheved of each of the above works, t should be pont out that most of these earler methods are nether general used nor provdng accurate results when the thckness of the thn structure s smaller than 1.0E-6. The objectve of ths paper s to develop a general BEM-based smulaton for predctng the physcal quanttes n thn structures when the thckness to length rato s n the orders of 1.0E-6 to 1.0E-9. It s well-known that the geometres of many problems of practcal nterest are created from crcular or ellptc arcs. Arc boundary elements can represent crcular and ellptc boundares exactly, and consequently, errors caused by representng such geometres usng polynomal shape functons can be removed usng exact geometrcal representatons. Therefore, the exact geometrcal representaton s expected to gve more accurate results than lower-order boundary element analyss when nearly sngular ntegrals need to be calculated. Another more mportant reason for usng exact geometrcal representaton can be found n the followng condtons. If the boundary geometry s depcted by usng lower order element, the lnear element of the outer surface wll attach or even pass through the nner boundary f the thckness of the consdered structure s very small. Consequently, the actual geometry of consdered doman can not be descrbed lvely, and thus lowerorder geometry approxmaton wll fal to yeld relable results for such problems. In order to avod ths phenomenon, very fne meshes mush be used n ths stuaton, but ths s bound to requre more preprocessng and CPU tme. Most mportantly, a great number of meshes wll produce a lot of artfcal corners whch wll lead to the uncontnuty of the tangent dervatve of the boundary unknowns. Ths s fatal to many engneerng problems. Obvously, the utlzaton of exact geometrcal representaton can be a good choce to avod ths problem. In vew of these reasons, ths paper wll gve an effcent strategy to treat thn body problems wth crcular and ellptc boundares. Ths paper s an extenson of our prevous work [5] where a new nonlnear transformaton method was proposed and appled to treat the thn body effect occurrng n D potental problems. Heren, we derve a general BEM-based approach to treat thn body problems n D elastostatcs. The BEMbased approach proposed n ths paper s constructed usng a combnaton of the regularzed ndrect BIEs and a general nonlnear transformaton whch can remove or damp out the nearly sngular propertes of the ntegral kernels. It s shown that ths combned BIE formulaton can provde stable results for the analyss of thn body problems. For the test problems studed, very promsng results are obtaned when the thckness to length rato of the coatngs s n the orders between 1.0E-1 and 1.0E-9, whch s suffcent for modelng most thn coated cuttng tools n the mcro- or nano- scales. In concluson, the seemng dffcult task of determnng the physcal quanttes n thn structures can be dealt wth effectvely and effcently. II. NON-SINGULAR BOUNDARY INTEGRAL EQUATIONS (BIES) It s well known that the doman varables can be computed by ntegral equatons only after all the boundary quanttes have been obtaned, and the accuracy of boundary quanttes drectly affects the valdty of the nteror quanttes. However, when calculatng the boundary quanttes, we have to deal wth the sngular boundary ntegrals, and usng the regularzed boundary ntegral equatons (BIEs). In ths paper, we always assume that Ω s a bounded doman n R ; Ω c s ts open complement, and denotes the boundary; t(x) and n(x) (or t and n) are the unt tangent and outward normal vectors of to the doman Ω at the pont x, respectvely. For D elastc problems, the non-sngular BIEs wth ndrect varables are gven n Ref. [9]. Regardless of to the rgd body dsplacement and the body forces, the non-sngular BIEs on Ω c can be expressed as
3 Y.-M. Zhang et al.: Boundary Element Analyss of Thn Structural Problems n D Elastostatcs 411 u y = x u y x d y (1) ( ) φk( ) k(, ), u = u d ( y) [ φk( x) φk( y)] k( y, x) φk( y ) u (, ) [( ) ( )] k yx t x t y d+ [ ( ) ( )] n x n y t uk ( yx, ) k0 d+ ny ( ) [ nk( ) nk( )] G x y n ln r ln r d+ nk( y) [ t( ) t( )] d x x y t ln r + nk( y) [ n( x) n( y)] d, y n For the doman Ω, the non-sngular BIEs are gven as ( ) φk( ) k(, ), ( () u y = x u x y d y (3) 1 nk( ) n( ) u ( ) φ k( ) ( ) [ k ] [ k( ) G δ y y y = y n y + (1 v) φ x φ x y k( y)] uk( y, x) d φk( y) [ t( ) t ( )] uk ( yx, ) uk ( yx, ) d+ [ ( ) ( )] d nx ny t n ( k0 ln r + n( y) [ nk( ) nk( )] d+ nk( ) G x y y x ln r [ t( x) t( y)] d+ nk( y) t ln r [ n( x) n( y)] d, y n For the nternal pont y, the ntegral equatons can be wrtten as ( ) ( ) (, ), ˆ φk k { (4) u y = x u y x d y Ω (5) ( ) ( ) (, ), ˆ φk k u y = x u y x d y Ω (6) In Eqs. (1)~(6),, k = 1, ; k 0 = 1/4π(1 v); G s the shear modulus; φ k (x) s the densty functon to be determned; u k (y, x) denotes the Kelvn fundamental soluton. In Eqs. (5) and (6) ˆΩ = Ω or Ω c. Note that, the Gaussan quadrature s drectly used to calculate the ntegrals n dscretzed equatons n the CBEM. However, f the doman of a consdered problem s thn, some boundares wll be very close to each other. Thus, the dstance r between some boundary nodes and boundary ntegral elements probably approaches zero. Ths causes the ntegrals n dscretzed Eqs. (1)~(4) nearly sngular, and the results of the Gaussan quadrature become nvald. Therefore, the densty functons cannot be calculated accurately, needless to say, to calculate the physcal quanttes at nteror ponts. Moreover, almost all the nteror ponts of thn bodes are very close to the ntegral elements. Thus, there also exst nearly sngular ntegrals n Eqs. (5) and (6). These nearly sngular ntegrals can be expressed as I1 = ψ ( x)lnrd e 1 I = ψ ( ) d x e r α where α > 0, ψ(x) s a well-behaved functon ncludng the Jacoban, the shape functons, and a fnte sum of polynomals dvded by r n. Under such a crcumstance, ether a very fne mesh wth massve ntegraton ponts or a specal ntegraton technque needs to be adopted. In the last two decades, numerous research works have been publshed on ths subject n the BEM lterature. Most of the works have been focused on the numercal approaches, such as subdvsons of the element of ntegraton, adaptve ntegraton schemes, exact ntegraton methods and so on. However, most of these earler methods are ether neffcent or can not provde accurate results when the thckness of the thn structure s smaller than 1.0E-6. In ths work, a very effcent transformaton method s proposed to avod the ntegraton dffculty for thn coatngs wth the thckness to length rato n mcro- or nano- scales. III. THE APPROXIMATION OF GEOMETRY BOUNDARIES The quntessence of the BEM s to dscretze the boundary nto a fnte number of segments, not necessarly equal, whch are called boundary elements. Two approxmatons are made over each of these elements. One s about the geometry of the boundary, whle the other has to do wth the varaton of the unknown boundary quantty over the element. 1. Lnear Element Approxmaton In ths secton, the geometry segment s modeled by a contnuous lnear element. Assumng x 1 = ( x 1 1 1, x), x = ( x 1, x) are the two extreme ponts of the lnear element j, then the element j can be expressed as 1 k 1 k k (7) x ( ξ) = N ( ξ) x + N ( ξ) x, ξ [ 1,1], k = 1, (8) where N 1 (ξ) = (1 ξ)/, N (ξ) = (1 + ξ)/. Lettng s = x 1 1 x, w = y ( x x )/, one has
4 41 Journal of Marne Scence and Technology, Vol. 19, No. 4 (011) r y x sξ + w r, = = = (9) r r r = = = ξ + ξ + = ξ η + r x y rr A B E L [( ) d ] (10) where A = s s /4, B = s w, E = w w, ηη = B/A, L = A, d = 4 AE B / A. Wth the ad of the Eq. (10), the nearly sngular ntegrals n Eq. (7) can be rewrtten as η { η } η 1 { η } 1 1 I1 = + ( )ln[( ) ] ln ( ) 1 g ξ ξ η + d dξ + L g ξ dξ (11) g( ξ ) I = + dξ α α L [( ξ η) + d ] where g( ) s a regular functon that conssts of shape functon and Jacoban.. Arc Element Approxmaton It s well-known that the geometres of many problems of practcal nterest are created from crcular or ellptc arcs. Arc boundary elements can represent crcular and ellptc boundares exactly, and consequently, errors caused by geometres usng polynomal shape functons whch can be removed by usng exact geometrcal representatons. Therefore, the exact geometrcal representaton s expected to gve more accurate results than lower-order or even hgh-order boundary element analyss when nearly sngular ntegrals need to be calculated. An exact geometrcal representaton, termed arc element, for crcular and ellptc boundares has been proposed by Zhang [9] n 004. Consder a crcular are element of radus of curvature R wth ts centre of curvature located at ( x 1, x ). Suppose (R, θ 1 ), (R, θ ) are the coordnates of the two extreme ponts of the arc element j, respectvely. Then the exact boundary elements can be expressed as: x1( ξ ) = x 1+ Rcosθ x( ξ ) = x + Rsnθ (1) where θ = 1 ξ θ θ, ( ξ 1). ξ The Jacoban of transformaton from arc element to ntrnsc coordnate ξ reduces to: J( ) R θ ξ θ 1 = (13) For the nteror pont y = (R 0 cos θ 0, R 0 sn θ 0 ), we can suppose θ 1 < θ 0 < θ and ( x 1, x ) = (0, 0). Then 1 η 1+ η θ0 = θ1+ θ ( 1< η < 1) Usng the procedure descrbed above, the dstance between the source pont and the feld pont can be expressed as: { γ } = = 4 0 sn + r x y RR d (14) R R0 where γ = β(ξ η), β = (θ θ 1 )/4, d =. 4 RR0 Substtutng r nto Eq. (7), then the ntegrals I 1 and I can be dvded nto two parts at pont η as follows η { η } η 1 α { η } 1 1 I1 = + ( )ln(sn ) ln ( ) 1 g ξ γ + d dξ + L g ξ dξ 1 g( ξ ) I = + dξ α L (sn γ + d ) (15) where L = RR0 and g( ) s a regular functon that conssts of shape functon, the Jacoban and a fnte sum of polynomals dvded by r n. IV. THE TRANSFORMATION FOR NEARLY SINGULAR INTEGRAL In Eqs. (11) and (15), f d s very small, the above ntegrals would present varous orders of near sngularty. The key to achevng hgh accuracy s to fnd a method to calculate these ntegrals accurately for a small value of d. The ntegrals I 1 and I n Eqs. (11) and (15) can be reduced to the followng ntegrals by smple deducton = + A gx ( ) T dx 0 = ( x + d ) α A T1 g( x)ln( x d ) dx 0 Introducng the followng nonlnear transformaton x d e t (16) k(1 + t) = ( ), [ 1,1] (17) where k = ln 1+ A d. Substtutng (17) nto (16), then the ntegrals I 1 and I can be rewrtten as follows: + dk g t e + e dt 1 k(1 + t) T1 = dk g()ln t d e dt 1 k(1 + t) k(1 + t) ()ln[( 1) 1] k(1 + t) 1 1 α gte () T = d k dt k(1 + t) α [( e 1) + 1] (18) We can observe that (e k(1+t) 1) Thus, the nte-
5 Y.-M. Zhang et al.: Boundary Element Analyss of Thn Structural Problems n D Elastostatcs 413 p A L = m Fg. 1. A thn plate under constant pressure p (Shear modulus μ = , Posson s rato v = 0.). h y x Normal stresses Present CBEM Exact grand s fully regular even f the value of d s very small. By followng the procedures descrbed above, the near sngularty n the BIEs has been fully regularzed. The fnal ntegral formulatons are obtaned as shown n Eq. (18), and be computed straghtforwardly by usng the standard Gaussan quadrature. V. NUMERICAL EXAMPLES To verfy the method developed above, two smple test problems are studed n whch BEM solutons are compared wth the exact solutons. 1. Test Problem 1: A Thn Plate Frst, a thn plate under external pressure p shown n Fg. 1 s studed. We assume the length of the plate n z drecton s large so that ths problem can be smplfed as a plane stran problem. The length L of the plate n x drecton s constant n ths study, whle the thckness h changes from L to 10-9 L. Note that the thckness s changng from macro-scale to mcroscale, and eventually to nano-scale. On the boundary y = 0, dsplacement components n both x and y drectons are constraned. The fnte element analyss of ths problem was attempted n Ref. [14], but t was soon found out that the number of the D fnte element was so large that the task quckly exceed the capacty of the computer used. In the BEM model, the boundary of the plate s dscretzed wth only 0 lnear boundary elements, 16 elements wth length L = m, and 4 other elements wth thckness h. When the thckness h ranges from 1.0E-1 to 1.0E-10, the results of the normal stresses σ xx and σ yy at nteror pont A(0.8, h/ are shown n Fg. and Fg. 3, respectvely. It s obvous that the results calculated by usng the CBEM deterorate quckly as the thckness less than 1.0E-. In contrast, the results calculated by usng the proposed method are very consstent wth the exact solutons even for the thckness as small as 1.0E-10. Ths proves that the developed transformaton work n the BEM procedure s effectve. For h = 1.0E 9, the stress results at the nteror ponts on the lne x = 1.8 are lsted n Table 1 and Table. These results calculated by the present BEM are all n good agreement wth the exact solutons wth the largest relatve error less than.0e-4 (%). However, the CBEM are nvald to calculate both the radal and tangental stresses at these ponts E-0 1E-04 1E-06 1E-08 1E-10 Thckness h Fg.. Normal stresses σ xx at the nteror pont A of Fg. 1. Normal stresses Present CBEM Exact E-0 1E-04 1E-06 1E-08 1E-10 Thckness h Fg. 3. Normal stresses σ yy at the nteror pont A of Fg. 1. The example problem studed here s a very smple one n thn structures. The purpose of usng ths example s to verfy the correctness of the present transformaton method n dealng wth nearly sngular ntegrals.. Test Problem : Thn Coatng on a Shaft In ths secton, the method developed n ths paper wll be used to solve a problem of a shaft wth a thn coatng [14], as shown n Fg. 4. The shaft and coatng have outer rad r a and r b respectvely, wth ther centre of curvature located at the pont o(0, 0). In ths example, the coated system s loaded by a unform pressure p = 1, and the shaft s consdered to be rgd when compared to the coatng, so the boundary condtons are u x = u y = 0 for all nodes at the shaft/coatng nterface. In ths example, δ = (r b r a )/r a s defned as the thckness to length rato. As r a s held constant at 1, the rato reduces as r b decreases. There are totally 0 dscontnuous arc elements dvded along the shaft and coatng surfaces, regardless of the thck-
6 414 Journal of Marne Scence and Technology, Vol. 19, No. 4 (011) Table 1. Normal stresses σ xx at the nteror ponts on the lne x = 1.8, where h = 1.0E 9. y Exact CBEM Relatve error (%) Present Relatve error (%) 1.0E E E E E-03.0E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E-03 Table. Normal stresses σ yy at the nteror ponts on the lne x = 1.8, where h = 1.0E 9. y Exact CBEM Relatve error (%) Present Relatve error (%) 1.0E E E E E-04.0E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E-04 r a r b p B o Fg. 4. A shaft wth a thn unform coatng. It s obvous that the results calculated by usng the proposed method, shown n Table 3, are very consstent wth the exact solutons even for the thckness to length rato as small as 1.0E-10. We can see from Table 4 that the results of tangental stresses calculated by usng the proposed method are very consstent wth the exact solutons, wth the largest relatve error less than 0.%, even when the thckness-to-length rato as small as 1.0E 9. For δ = 1.0E 9, convergence curves of the computed tangental stresses at nteror ponts B by usng the presented method are shown n Fg. 5, whch we can observe that the convergence speeds are stll fast even when the thcknessto-length rato s as small as 1.0E 9. In Fg. 5, only the errors of the present method are gven snce the errors of the CBEM are relatvely too large. ness of the structure. The elastc shear modulus s G = Pa, Posson s rato s v = 0.. In 1998, Luo et al. [14] have handled ths coatng system. However, n ther work only boundary unknown radal stresses σ r are computed. The physcal quanttes at nteror ponts need further nvestgaton. For dfferent thckness-to-length ratos, the results of the radal and tangental stresses at nteror pont B((r a + r b )/, 0) are lsted n Table 3. For δ = 1.0E 9, the results of tangental stresses on the lne y = 0 are lsted n Table 4. Both the CBEM and the proposed method are employed for the purpose of comparson. VI. CONCLUSION In ths paper, a BEM-based approach s presented and appled to deal wth -D elastc problems of thn bodes. Usng the transformaton technque demonstrated n ths paper, the seemngly dffcult task of evaluatng the nearly sngular ntegrals n the BIE for -D thn structures can be dealt wth effectvely and effcently. Two numercal examples of elastc thn-walled structures wth thckness-to-length ratos rangng between 1.0E 1 and 1.0E 9, whch s suffcent for modelng most thn structures n ndustral applcatons, are presented to test the proposed method. In concluson, the
7 Y.-M. Zhang et al.: Boundary Element Analyss of Thn Structural Problems n D Elastostatcs 415 Table 3. Radal and Tangental stresses at the nteror pont B of Fg. 4. δ Radal stresses σ r Tangental stresses σ θ Exact Present Exact Present 1.0E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+00 Table 4. Tangental stress σ θ at the nteror pont on the lne y = 0, where δ = 1.0E 9. x Exact CBEM Relatve error (%) Present Relatve error (%) E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+00 Relatve error (%) tangental stresses Number of elements Fg. 5. Convergence curves of σ θ at the nteror pont B of Fg. 4. thn-body problem has been overcome successfully by usng the proposed strategy, whch ndcates that BEM s especally accurate and effcent for numercal analyss of thn boy problems. The developed method for analyzng -D thn structures can be extended to model layered coatngs, thn flms or other layered structures, and some work along ths lne for thn structures wll be dscussed later. ACKNOWLEDGMENTS The research s supported by the Natonal Natural Scence Foundaton of Chna ( ), the Openng Fund of the State Key Laboratory of Structural Analyss for Industral Equpment (GZ1017) and the Natonal Natural Scence Foundaton of Shangdong Provnce of Chna (ZR010AZ003). REFERENCES 1. Brebba, C. A., Tells, J. C. F., and Wrobel, L. C., Boundary Element Technques, Sprnger-Verlag, Berln, Germany (1984).. Chen, J. T., Chen, P. Y., and Cha, T. C., Surface moton of multple alluval valleys for ncdent plane SH-waves by usng a sem-analytcal approach, Sol Dynamcs and Earthquake Engneerng, Vol. 8, No. 3, pp (008). 3. Chen, J. T. and Hong, H. K., Revew of dual boundary element methods wth emphass on hypersngular ntegrals and dvergent seres, Appled Mechancs Revews, Vol. 5, No. 1, pp (1999). 4. Chen, J. T. and Shen, W. C., Degenerate scale for multply connected laplace problems, Mechancs Research Communcatons, Vol. 34, No. 1, pp (007b). 5. Granados, J. J. and Gallego, G., Regularzaton of nearly hypersngular ntegrals n the boundary element method, Engneerng Analyss wth Boundary Elements, Vol. 5, No. 3, pp (001). 6. Hong, H.-K. and Chen, J. T., Dervatons of ntegral equatons of elastcty, Journal of Engneerng Mechancs, ASCE, Vol. 114, No. 6, pp (1988). 7. Hsao, S. S., Ln, M. C., and Hu, N. C., Drbem analyss of combned wave refracton and dffracton n the presence of current, Journal of Marne Scence and Technology, Vol. 10, No. 1, pp (00). 8. Huang, Q. and Cruse, T. A., Some notes on sngular ntegral technques n boundary element analyss, Internatonal Journal for Numercal
8 416 Journal of Marne Scence and Technology, Vol. 19, No. 4 (011) Methods n Engneerng, Vol. 36, No. 15, pp (1993). 9. Jun, L., Beer, G., and Meek, J. L., Effcent evaluaton of ntegrals of order usng Gauss quadrature, Engneerng Analyss, Vol., No. 3, pp (1985). 10. Lu, C. S., A modfed collocaton Trefftz method for the nverse Cauchy problem of laplace equaton, Engneerng Analyss wth Boundary Elements, Vol. 3, No. 9, pp (006). 11. Lu, C. S., Chang, C. W., and Chang, C. Y., A regularzed ntegral equaton method for the nverse geometry heat conducton problem, Internatonal Journal of Heat and Mass Transfer, Vol. 51, No. 1, pp (008). 1. Lu, Y. J., Analyss of shell-lke structures by the boundary element method based on 3-D elastcty: formulaton and verfcaton, Internatonal Journal for Numercal Methods n Engneerng, Vol. 41, No. 3, pp (1998). 13. Lu, Y. J., On the smple soluton and non-sngular nature of the BIE/BEM-a revew and some new results, Engneerng Analyss wth Boundary Elements, Vol. 4, No. 10, pp (000). 14. Luo, J. F., Lu, Y. J., and Berger, E. J., Analyss of two-dmensonal thn structures (from mcro- to nano-scales) usng the boundary element method, Computatonal Mechancs, Vol., pp (1998). 15. Lutz, E. L., Exact Gaussan quadrature methods for near-sngular ntegrals n the boundary element method, Engneerng Analyss wth Boundary Elements, Vol. 9, No. 3, pp (199). 16. Ma, H. and Kamya, N., Dstance transformaton for the numercal evaluaton of near sngular boundary ntegrals wth varous kernels n boundary element method, Engneerng Analyss wth Boundary Elements, Vol. 5, No. 4, pp (00a). 17. Mukerjee, S., Chat, M. K., and Sh, X. L., Evaluaton of nearly sngular ntegrals n boundary element contour and node methods for threedmensonal lnear elastcty, Internatonal Journal of Solds and Structures, Vol. 37, No. 51, pp (000). 18. Nu, Z. R., Wang, X. X., and Zhou, H. L., A general algorthm for calculatng the quanttes at nteror ponts close to the boundary by the BEM, Chnese Journal of Theoretcal and Appled Mechancs, Vol. 33, No., pp (001). (n Chnese) 19. Sladek, V. and Sladek, J., Sngular Integrals n Boundary Element Methods, Computatonal Mechancs Publcatons, Southampton, UK and Boston (1998). 0. Sladek, V., Sladek, J., and Tanaka, M., Regularzaton of hypersngular and nearly sngular ntegrals n the potental theory and elastcty, Internatonal Journal for Numercal Methods n Engneerng, Vol. 36, No. 10, pp (1993). 1. Sladek, V., Sladek, J., and Tanaka, M., Nonsngular BEM formulatons for thn-walled structures and elastostatc crack problems, Acta Mechanca, Vol. 99, pp (1993).. Sun, H. C., Nonsngular Boundary Element Method, Dalan Unversty of Technology Press (1999). (n Chnese) 3. Tanaka, M., Sladek, V., and Sladek, J., Regularzaton technques appled to BEM, Appled Mechancs Revews, Vol. 47, No. 10, pp (1994). 4. Wang, W. H., Lou, J. H., Robert, S., and Ben, D., Machne vbraton nduced underwater acoustc radaton, Journal of Marne Scence and Technology, Vol. 8, No. 1, pp (000). 5. Zhang, Y. M. and Gu, Y., An effectve method n BEM for potental problems of thn bodes, Journal of Marne Scence and Technology, Vol. 18, No. 1, pp (010). 6. Zhang, Y. M. and Sun, C. L., A general algorthm for the numercal evaluaton of nearly sngular boundary ntegrals n the equvalent nonsngular BIEs wth ndrect unknowns, Journal of the Chnese Insttute of Engneers, Vol. 31, No. 3, pp (008). 7. Zhang, Y. M. and Sun, H. C., Theoretc analyss on vrtual boundary element, Chnese Journal of Computatonal Mechancs, Vol. 17, No. 1, pp (000). (n Chnese) 8. Zhang, Y. M. and Sun, H. C., Analytcal treatment of boundary ntegrals n drect boundary element analyss of plane potental elastcty problems, Appled Mathematc and Mechancs, Vol. 6, pp (001). (n Chnese) 9. Zhang, Y. M., Wen, W. D., Wang, L. M., and Zhao, X. Q., A knd of new nonsngular boundary ntegral equatons for elastc plane problems, Acta Mech, Vol. 36, No. 3, pp (004). (n Chnese) 30. Zhou, H. L., Nu, Z. R., Cheng, C. Z., and Guan, Z. W., Analytcal ntegral algorthm n the BEM for orthotropc potental problems of thn body, Engneerng Analyss wth Boundary Elements, Vol. 31, pp (007). 31. Zhou, H. L., Nu, Z. R., Cheng, C. Z., and Guan, Z. W., Analytcal ntegral algorthm appled to boundary layer effect and thn body effect n BEM for ansotropc potental problems, Computers and Structures, Vol. 86, pp (008).
Numerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationSimulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM
Smulaton of 2D Elastc Bodes wth Randomly Dstrbuted Crcular Inclusons Usng the BEM Zhenhan Yao, Fanzhong Kong 2, Xaopng Zheng Department of Engneerng Mechancs 2 State Key Lab of Automotve Safety and Energy
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationThe Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationEVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationNumerical Nonlinear Analysis with the Boundary Element Method
Blucher Mechancal Engneerng Proceedngs May 2014, vol. 1, num. 1 www.proceedngs.blucher.com.br/evento/10wccm Numercal Nonlnear Analyss wth the Boundary Element Method E. Pneda 1, I. Vllaseñor 1 and J. Zapata
More informationAPPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS
6th ICPT, Sapporo, Japan, July 008 APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS James MAINA Prncpal Researcher, Transport and Infrastructure Engneerng, CSIR Bult Envronment
More informationCHAPTER 9 CONCLUSIONS
78 CHAPTER 9 CONCLUSIONS uctlty and structural ntegrty are essentally requred for structures subjected to suddenly appled dynamc loads such as shock loads. Renforced Concrete (RC), the most wdely used
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationANALYSIS OF CONTACT PROBLEM USING IMPROVED FAST MULTIPOLE BEM WITH VARIABLE ELEMENTS LENGTH THEORY
Journal of Marne Scence and Technology, Vol., No., pp. -7 () DOI:.69/JMST--7- NLYSIS OF CONTCT PROBLEM USING IMPROVED FST MULTIPOLE BEM WITH VRIBLE ELEMENTS LENGTH THEORY Ha-Lan Gu, Qang L, Qng-Xue Huang,
More informationAssessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion
Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationConstitutive Modelling of Superplastic AA-5083
TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationInstituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245
Insttuto Tecnológco de Aeronáutca FIITE ELEMETS I Class notes AE-5 Insttuto Tecnológco de Aeronáutca 5. Isoparametrc Elements AE-5 Insttuto Tecnológco de Aeronáutca ISOPARAMETRIC ELEMETS Introducton What
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationGeneral displacement arch-cantilever element method for stress analysis of arch dam
Water Scence and Engneerng, 009, (): 3-4 do:0.388/j.ssn.674-370.009.0.004 http://kkb.hhu.edu.cn e-mal: wse@hhu.edu.cn General dsplacement arch-cantlever element method for stress analyss of arch dam Hao
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationSpeeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem
H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationStudy on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI
2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationAdaptive sliding mode reliable excitation control design for power systems
Acta Technca 6, No. 3B/17, 593 6 c 17 Insttute of Thermomechancs CAS, v.v.. Adaptve sldng mode relable exctaton control desgn for power systems Xuetng Lu 1, 3, Yanchao Yan Abstract. In ths paper, the problem
More informationNUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST
Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationA boundary element method with analytical integration for deformation of inhomogeneous elastic materials
Journal of Physcs: Conference Seres PAPER OPEN ACCESS A boundary element method wth analytcal ntegraton for deformaton of nhomogeneous elastc materals To cte ths artcle: Moh. Ivan Azs et al 2018 J. Phys.:
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationStudy on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component
Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationOFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES
ICAMS 204 5 th Internatonal Conference on Advanced Materals and Systems OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES VLAD LUPĂŞTEANU, NICOLAE ŢĂRANU, RALUCA HOHAN, PAUL CIOBANU Gh. Asach Techncal Unversty
More informationOne Dimensional Axial Deformations
One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationA Computational Viewpoint on Classical Density Functional Theory
A Computatonal Vewpont on Classcal Densty Functonal Theory Matthew Knepley and Drk Gllespe Computaton Insttute Unversty of Chcago Department of Molecular Bology and Physology Rush Unversty Medcal Center
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationFast Simulation of Pyroshock Responses of a Conical Structure Using Rotation-Superposition Method
Appled Mathematcs & Informaton Scences An Internatonal Journal 211 NSP 5 (2) (211), 187S-193S Fast Smulaton of Pyroshock Responses of a Concal Structure Usng Rotaton-Superposton Method Yongjan Mao 1, Yulong
More informationGEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE
GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE Prof. J. N. Mandal Department of cvl engneerng, IIT Bombay, Powa, Mumba 400076, Inda. Tel.022-25767328 emal: cejnm@cvl.tb.ac.n Module - 9 LECTURE - 48
More informationNONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM
Advanced Steel Constructon Vol. 5, No., pp. 59-7 (9) 59 NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM M. Abdel-Jaber, A.A. Al-Qasa,* and M.S. Abdel-Jaber Department of Cvl Engneerng, Faculty
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationIntegrals and Invariants of
Lecture 16 Integrals and Invarants of Euler Lagrange Equatons NPTEL Course Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng, Indan Insttute of Scence, Banagalore
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationTHE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES
THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor
More informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More information