A Frequency-Domain Approach for Transient Dynamic Analysis using Scaled Boundary Finite Element Method (I): Approach and Validation

Size: px
Start display at page:

Download "A Frequency-Domain Approach for Transient Dynamic Analysis using Scaled Boundary Finite Element Method (I): Approach and Validation"

Transcription

1 COMPUTATIONAL METHODS IN ENGINEERING AND SCIENCE EPMESC X, Aug. 2-23, 26, Sanya, Hanan, Chna 26 Tsnghua Unversty Press & Sprnger A Frequency-Doman Approach for Transent Dynamc Analyss usng Scaled Boundary Fnte Element Method (I): Approach and Valdaton Z. J. Yang, 2 *, A. J. Deeks 2, H. Hao 2 College of Cvl Engneerng and Archtecture, Zhejang Unversty, Hangzhou, 327 Chna 2 School of Cvl and Resource Engneerng, the Unversty of Western Australa, WA 69 Australa Emal: yang@cvl.uwa.edu.au Abstract The sem-analytcal scaled boundary fnte element method (SBFEM) has been successfully appled n elastostatcs. Its applcaton to elastodynamc problems, however, has lagged behnd, prmarly due to the lack of an effectve soluton procedure to the governng equlbrum equaton system. The authors recently developed an easy-to-follow Frobenus soluton procedure to ths equaton system []. Ths study further develops a frequency-doman approach for the general transent dynamc analyss, through combnng the Frobenus soluton procedure and the fast Fourer transform (FFT) technque. The complex frequency-response functons (CFRFs) are frst computed usng the Frobenus soluton procedure. Ths s followed by a FFT of the transent load and an nverse FFT of the CFRFs to obtan the tme hstory of responses. Two transent dynamc problems, subjected to Heavsde step load and trangular blast load, are modelled wth a small number of degrees of freedom usng the new approach. The numercal results agree very well wth the analytcal solutons and those from the FEM. Further applcatons of ths approach to transent dynamc fracture problems are presented n the accompanyng paper [2]. Keywords: transent dynamc analyss, scaled boundary fnte element method, Frobenus soluton procedure, frequency doman, fast Fourer transform INTRODUCTION Transent dynamc analyses are very mportant for understandng responses and falure mechansms of engneerng structures subjected to dynamc loads such as mechancal vbraton, earthquake, mpact, wnd and blast. The fnte element method (FEM) s undoubtedly the domnant method for modellng transent dynamc problems at present, because of ts powerful capablty of smulatng a large varety of problems wth complex structural geometres, complcated materal propertes, and varous loadng and boundary condtons [3]. Although ts success n transent dynamc analyses has been frmly establshed, modellng large-scale lnear and nonlnear structural dynamc systems usng the FEM stll poses great dffcultes for structural analysts [4]. In spte of great efforts, the computatonal cost for the FEM may be stll very hgh for some mportant problems, such as the structural-sol dynamc nteracton nvolvng nfnte domans, solds wth stress sngularty and concentraton (e.g., cracks and corners), and the dynamc problems wth responses domnated by ntermedate and hgh frequency modes, such as those subjected to mpact or blast loads [5]. One alternatve to the FEM s the boundary element method (BEM). The BEM dscretses the boundares only and thus reduces the modelled dmensons by one. The applcaton of the BEM n the transent dynamc analyss s currently very actve, wth many numercal procedures reported [e.g., 6]. BEM/FEM couplng methods have also been developed for lnear elastodynamcs [7]. However, the need of fundamental solutons lmts the applcablty of the BEM consderably. The scaled boundary fnte-element method (SBFEM), developed recently by Wolf and Song [9, ], s a sem-analytcal method combnng the advantages of the FEM and the BEM, and possessng ts own attractons at the same tme. Specfcally, t dscretses the boundares only and thus reduces the modelled spatal dmensons by one as the BEM, but does not need fundamental solutons. The method thus has much wder applcablty than the BEM. In addton, the dsplacement and stress solutons of the SBFEM are approxmate n the crcumferental drecton n the FEM sense, but analytcal n the radal drecton. Ths sem-analytcal nature makes the method an excellent tool for modellng unbounded domans because the radaton condton at nfnty s automatcally satsfed. Ths has been 756

2 demonstrated by a few recent studes of dynamc sol-structure nteracton [, 2]. What s more, the analytcal stress feld explctly represents stress sngulartes at crack tps, whch allows accurate stress ntensty factors to be computed drectly from the defnton. Ths advantage has been used recently by Song to calculate transent dynamc stress ntensty factors [3]. The SBFEM has also a few dsadvantages. For example, ts current formulaton consders only lnear elastc materal behavour. For complex domans that are dvded nto a few subdomans, an egenproblem wth 2n degrees of freedom (DOFs) must be solved for each subdoman of n DOFs, whch may result n hgh computatonal cost. To date, all of these lmted transent dynamc analyses [-3] based on the SBFEM used drect tme-ntegraton schemes, and hence a tme-doman approach. The tme-doman approach requres a mass matrx n the equatons of moton. The mass matrx derved by Song [3] corresponds to the low-frequency expanson of the dynamc stffness matrx [9, ] and thus can only accurately represent the nertal effects at low frequences. Ths requres that the sze of the subdomans or super-elements be small enough to account for the hghest frequency component of nterest [3]. Ths may lead to consderable computatonal cost n solvng a sgnfcant number of egenproblems n each tme step. A frequency-doman approach has certan advantages over the tme-doman approach when modellng elastodynamcs, such as no need for a mass matrx, so that fewer subdomans or coarser meshes may be used, and once a complex frequency-response functon (CFRF) s obtaned, t can be used by the fast Fourer transform (FFT) and the nverse FFT (IFFT) to calculate transent responses for varous forms of dynamc loads. The authors recently developed an easy-to-follow analytcal soluton n the frequency doman usng the Frobenus procedure to the governng equlbrum equatons of the SBFEM[]. Ths study further devses a new frequency-doman approach for general transent dynamc analyss usng the SBFEM. The approach nvolves two steps. The frst step s to employ the Frobenus soluton procedure n the frequency doman to compute the CFRFs. The second step s to use FFT and IFFT to calculate the tme hstores from the CFRFs. The followng sectons descrbe the basc concept of SBFEM and the Frobenus soluton, followed by modellng of a wave propagaton problem and a forced vbraton problem usng the new approach. FUNDAMENTALS OF THE SBFEM A doman of an arbtrary shape s llustrated n Fg. (a). The doman s dvded nto three subdomans. Any scheme of subdvson, wth varous numbers, shapes and szes of subdomans, can be used, as long as a scalng centre for each subdoman can be found from whch the subdoman boundary s fully vsble. Fg. (b) shows the detals of Subdoman. The subdoman s represented by scalng a defnng curve S relatve to a scalng centre. The defnng curve s usually taken to be the doman boundary, or part of the boundary. A normalsed radal coordnate ξ s defned, varyng from zero at the scalng centre and unt value on S. A crcumferental coordnate s s defned around the defnng curve S. A curve smlar to S defned by ξ=.5 s shown n Fg. (b). The coordnates ξ and s form the local coordnate system, whch s used n all the subdomans. (a) Subdomanng of a doman (b) Subdoman Fgure : The concept of the scaled boundary fnte-element method The governng equlbrum equatons of the SBFEM have been derved wthn the context of vrtual work prncple for elastostatcs by Deeks and Wolf [4]. An extenson to elastodynamcs has also been made by Yang et al []. They are not repeated here. In a frequency-doman analyss, the frst step s to compute the CFRFs. The governng equlbrum equatons of the SBFEM under a harmonc exctaton wth constant frequency f (Hz), e.g. p=p e 2πft are [] T p = E u( ξ), ξ + E u ( ξ ) ξ= () 757

3 E ξ u( ξ), + ( E + E E ) ξu( ξ), E u( ξ) + (2 π f ) ξ M u ( ξ) = (2) 2 T ξξ ξ where p are the magntudes of the equvalent nodal loads. E, E, E 2 and M are coeffcent matrces that are dependent on the geometry of the subdoman boundares and the materal propertes. u(ξ) represent the magntudes of the nodal dsplacements and are analytcal functons of the radal coordnate ξ. Eq. () and Eq. (2) are frst appled to the subdomans. The contrbutons of subdomans to the coeffcent matrces and equvalent nodal loads are then assembled n the same way as the contrbutons of fnte elements n FEM. THE FROBENIUS SOLUTION PROCEDURE Eq. (2) s a second-order nonhomogeneous Euler-Cauchy dfferental equaton system wth respect to the radal coordnate ξ. The non-homogenety complcates ts soluton consderably. An easy-to-follow analytcal soluton to Eq. (2) was recently developed by Yang et al [] usng the Frobenus procedure. The soluton wth (k+) seres s n n n 2 k k + ( λ) ( λ) 2 ( λ) k u = cξ φ + ξ c g + ξ + c + g = = = where k+ k+ (3) g = G φ (4a) G = B B B (4b) k+ k+ k 2 k (2 ) ( k + ) ( k + )( T = π f λ + λ ) B E E E E M (4c) d = B ( d ) (4d) k+ k+ k d = cφ (4e) k λ = λ + 2 (4f) k+ where the subscrpt ranges from to n for all the varables and n s the DOFs of the problem. Eq. (3) shows that the soluton to elastodynamcs conssts of (k+) summatons of seres. The frst summaton s the soluton to the correspondng homogeneous equatons of Eq. (2),.e., the governng equatons for elastostatcs. The added k summatons account for dynamc effects. More seres terms are added, more accurately the soluton represents the dynamc effects. λ = λ and φ are the postve egenvalues and correspondng egenvectors of the standard lnear egenproblem [4] formed from the elastostatc governng equatons. c are constants determned by boundary condtons. Eqs. (4a- 4f) descrbe an teratve process n whch the number of added summaton terms k n Eq. (3) s determned by a convergence crteron []: max ( R )< α (5a) 2 n 2 ( k λ ) k (2 π f ) ξ + + R = M ( d ) (5b) = where R can be regarded as the resdual vector and α s the convergence tolerance. Settng ξ= n Eq. (3) leads to k + u = u = Ψ c or = b c= Ψ ub (6) ξ where u b s the nodal dsplacements on the boundary, c={c c 2 c n } T and the matrx k+ k+ k+ j j j Ψ = ϕ+ g ϕ2 + g2 ϕ n + gn (7) j= 2 j= 2 j= 2 Substtutng Eq. (3) nto Eq. () leads to p = E Ψ c+ E Ψ c (8) T 2 where the matrx 758

4 k+ k+ k+ j j j j j j Ψ2 = λϕ+ g λ λ2ϕ2 + g2 λ2 λnϕ n + gn λn (9) j= 2 j= 2 j= 2 Substtutng Eq. (6) nto Eq. (8) gves T p = E ΨΨ + E u = K u () 2 b d b Therefore, the dynamc stffness matrx of the doman wth respect to DOFs on the boundary s K = E ΨΨ + E () T d 2 The nodal dsplacement vector u b can be calculated by Eq. () by applyng boundary condtons on u b and loadng condtons on p. The ntegraton constants c are then obtaned usng Eq. (6). Assumng the (k+)th soluton meets the crteron Eq. (5a), the dsplacement feld s recovered as n n = + + n 2 k + ( λ ) ( λ ) 2 ( λ ) k+ u( ξ, s) N( s) cξ φ ξ c ( g ) L ξ c ( g ) (2) = = = where N(s) s the matrx of shape functons at the crcumferental drecton, whch are constructed as n FEM, typcally usng polynomal functons. The stress feld s then obtaned as n n ( ) n k σ ξ, s = DB s c λ ξ ϕ + c λ ξ ( g ) + + c ( λ ) ξ g k k k k 2 + λ 2 λ 2 + λ + = = = n n ( ) n + DB ϕ + ( g ) + + g λ λ 2 λ + s cξ cξ cξ = = = where D s the elastcty matrx and B (s) and B 2 (s) are relevant matrces []. One can see that the dynamc dsplacement and stress felds from the SBFEM (Eq. (2) and Eq. (3)) are analytcal wth respect to the radal coordnate ξ. Therefore, these state varables at any poston of the doman can be drectly calculated once the nodal dsplacements on the boundary are obtaned by solvng Eq. (). (3) NUMERICAL EXAMPLES, RESULTS AND DISCUSSION Two transent dynamc examples, both subjected to the Heavsde-step load and the trangular blast load, are modelled. For each example, the CFRFs of dsplacements or stresses at desred postons n the doman are frst computed for a wde range of samplng frequences usng the Frobenus soluton procedure presented n the prevous secton. A FFT of the transent load s then carred out to obtan the complex dscrete Fourer coeffcents, followed by an IFFT of the CFRFs to obtan tme hstores of dsplacement or stress at desred postons. The readers are referred to [5] for detals of the dscrete Fourer transform and the fast Fourer transform technques. The functons fft and fft n MATLAB [6] are used to conduct FFT and IFFT respectvely. For each example, the same CFRFs are used to calculate tme hstores for all the forms of transent loadngs that ths example s subjected to. In the frequency-doman SBFEM, the dampng effect s taken nto account by modfyng the elastc modul to ncorporate an nternal dampng coeffcent β, formng the complex Young s modulus E c =E(+2β) and the complex shear modulus G c =G(+2β), where E and G are the elastc Young s modulus and the elastc shear modulus respectvely. A transent dynamc fnte element analyss usng ABAQUS s also carred out for both examples for comparson usng the mplct tme-ntegraton scheme [7]. No physcal dampng (.e., pure elastc materal behavour) s consdered n all the fnte element analyses. The nternatonal standard unt system, namely mass n klograms, length n metres and tme n seconds, s used for both examples and thus no unts are ndcated wth all the physcal enttes. Two-node lnear lne elements and α =. n Eq. (5a) are used for all the analyses usng the SBFEM. A plane stress condton s assumed for both examples.. Example : D rectangular rod The frst example modelled here s the classcal wave propagaton problem wth analytcal solutons, whch was modelled by many, e.g., [6] and [7], to valdate numercal methods. A rectangular rod s subjected to a unform tracton on the rght end wth the left end fxed. The geometry, boundary and loadng condtons are shown n Fg. 2(a). The tme responses at the ponts A, B and C n Fg. 2(a) are examned. Only one subdoman s used to model ths example. Fg. 2(b) shows a typcal mesh wth only 5 nodes. The scalng centre s 759

5 placed at the left-bottom corner so that the two edges connected to t are not dscretsed. β= 6 s used for ths example. (a) Dmensons, materal propertes and loadng condtons Fgure 2: Example, an D rectangular rod (b) A mesh wth 5 nodes Fgure 3: Horzontal dsplacement at the pont A under the Heavsde step load for the st example Fgure 4: Horzontal dsplacement at the pont B under the Heavsde step load for the st example Fg. 3 and Fg. 4 show the tme hstores of the horzontal (axal) dsplacements at the free end of the rod (the pont A) and the mddle of the rod (the pont B) respectvely under the Heavsde step load. The results from both the frequency-doman SBFEM and the tme-doman FEM are n very good agreement wth the analytcal solutons [8]. The axal stress hstores at the mddle and the fxed end of the rod are shown n Fg. 5 and Fg. 6 respectvely. Good agreement between the numercal results and the analytcal soluton can also be observed, except that the FEM results oscllate fercely around the analytcal solutons at moments when the stress jumps. Ths s caused by the sudden applcaton of the step load and cannot be completely corrected by numercal measures. It also happens n other studes [6, 7]. The oscllatons also appear n the results from the new frequency-doman approach, but wth lower magntudes and much less degree. Fg. 7 compares the horzontal dsplacement hstores at the ponts A and B under the trangular blast load. The results from the two numercal methods are vrtually dentcal. Excellent agreement can also be seen n the axal stress hstores at the mddle and the fxed end of the rod, whch are presented n Fg. 8. No oscllaton happens because the trangular blast load functon s contnuous. 76

6 Fgure 5: Horzontal stress at the pont B under the Heavsde step load for the st example Fgure 6: Horzontal stress at the pont C under the Heavsde step load for the st example.75 Horzontal dsplacements FEM (A) Present (A) FEM (B) Present (B) Fgure 7: Horzontal dsplacements at the pont A and B under the blast load for the st example 2.5 Horzontal stress FEM (B) Present (B) FEM (C) Present (C) Fgure 8: Horzontal stress at the pont B and C under the blast load for the st example 76

7 2. Example 2: 2D smply-supported beam subjected to unform bendng force The forced vbraton of a 2D smply-supported beam s modelled as the 2 nd example. The dmensons, boundary and loadng condtons, and materal propertes are llustrated n Fg. 9(a). The beam s subjected to a unformly-dstrbuted bendng tracton at ts upper face. The tme responses at the ponts A and B n Fg. 9(a) are nvestgated. Both the Heavsde step load and the blast load (only the latter s shown n Fg. 9(a)) are modelled. The beam s modelled wth two subdomans consstng of totally 33 nodes usng the SBFEM as shown n Fg. 9(b). The two scalng centres are postoned at the beam lower face so that ths whole face s not dscretsed. A materal dampng coeffcent β= 5 s used. (a) Dmensons, materal propertes and loadng condtons (b) A mesh wth 33 nodes Fg. 9 Example 2: a smply-supported beam subjected to unformly dstrbuted bendng force Fg. and Fg. show the vertcal dsplacement hstores at the pont B and the horzontal stress hstores at the pont A respectvely under the Heavsde step load. Also shown n these fgures are the analytcal solutons [8] takng nto account the effects of rotary nerta and shear deformaton. It s evdent that the results from the tme-doman FEM and the frequency-doman SBFEM almost concde wth each other, and the results from both methods agree very well wth the analytcal solutons, especally for the dsplacement hstores. Because the stresses are dervatves of dsplacements and thus less accurate, the dscrepances between the analytcal stress soluton and the numercal results are understandable. Comparng the results here wth those reported n [9] usng a tme-doman FEM also confrms the effectveness and accuracy of the developed frequency-doman SBFEM. The same concluson can also be drawn from the results under the blast load, whch are shown n Fg. 2 and Fg. 3. Vertcal dsplacement at pont B Analytcal Present FEM Fgure : Vertcal dsplacement at the pont B under the Heavsde step load for the 2 nd example Horzontal stress at pont A 5 45 Analytcal Present FEM Fgure : Horzontal stress at the pont A under the Heavsde step load for the 2 nd example 762

8 Vertcal dsplacement at pont B.5 Present FEM Fgure 2: Vertcal dsplacement at the pont B under the blast load for the 2nd example 4.5 Present FEM Horzontal stress at pont A Fgure 3: Horzontal stress at the pont A under the blast load for the 2nd example CONCLUSIONS Ths study develops a frequency-doman approach usng the SBFEM for modellng general transent dynamc problems. The CFRFs are frst computed by the newly-developed analytcal Frobenus soluton to the governng equatons of SBFEM. The FFT and the IFFT are then carred out to calculate tme hstores of structural responses. Ths new approach does not need a mass matrx as needed by the tme-doman SBFEM, allowng domans to be modelled by fewer subdomans effcently. The new approach has been used to successfully model transent dynamc behavour of two examples wth very small number of DOFs due to the analytcal nature of the Frobenus soluton. ACKNOWLEDGEMENT Ths research s supported by the Australan Research Councl (Dscovery Project No. DP45268) and the Natonal Natural Scence Foundaton of Chna (No ). REFERENCES. Yang ZJ, Deeks AJ, Hao H. A Frobenus soluton to the scaled boundary fnte element equatons n frequency doman. Int. J. Num. Meth. Eng., (under revew). 2. Yang ZJ, Deeks AJ, Hao H. A frequency-doman approach for transent dynamc analyss usng scaled boundary fnte element method (II): applcaton to fracture problems. The th Int. Conf. Comp. Meth. n Engng.&Sc (EPMESC X). Aug. 2-23, 26, Sanya, Hanan, Chna. 3. Zenkewcz OC, Taylor RL. The fnte element method. 5th edton, Butterworth- Henemann, Dokansh MA, Subbaraj K. A survey of drect tme-ntegraton methods n computatonal structural dynamcs I. explct methods. Computers & Structures, 989; 32(6): Zeng LF, Wberg NE. Spatal mesh adaptaton n semdscrete fnte element analyss of lnear elastodynamcs problems. Computatonal Mechancs, 992; 9:

9 6. Carrer JAM, Mansur WJ. Alternatve tme-marchng schemes for elastodynamc analyss wth the doman boundary element method formulaton. Computatonal Mechancs, 24; 34(5): Yu GY, Mansur WJ, Carrer JAM. A more stable scheme for BEM/FEM couplng appled to two-dmensonal elastodynamcs. Computers & Structures, 2; 79(8): Gu YT, Lu GR. A meshless local Petrov Galerkn (MLPG) method for free and forced vbraton analyses for solds. Computatonal Mechancs, 2; Wolf JP, Song CM. Fnte-element modellng of unbounded meda. John Wley and Sons, Chchester, Wolf JP. The scaled boundary fnte element method. Chchester: John Wley and Sons, 23.. Zhang X, Wegner JL, Haddow JB. Three-dmensonal dynamc sol-structure nteracton analyss n the tme doman. Earthquake Engneerng & Structural Dynamcs, 999; 28(2): Yan JY, Zhang CH, Jn F. A couplng procedure of FE and SBFE for sol-structure nteracton n the tme doman. Internatonal Journal for Numercal Methods n Engneerng, 24; 59(): Song CM. A super-element for crack analyss n the tme doman. Internatonal Journal of Numercal Methods n Engneerng, 24; 6(8): Deeks AJ, Wolf JP. A vrtual work dervaton of the scaled boundary fnte-element method for elastostatcs. Computatonal Mechancs, 22; 28(6): Brgham EO. The fast Fourer transform and ts applcatons. Prentce-Hall, MATLAB V7.. The MathWorks, Inc., Hbbtt, Karlsson and Sorensen Inc. ABAQUS/Standard User Manual V6.5, Tmoshenko S, Young DH, Weaver W. Vbraton Problems n Engneerng. 4th edton, John Wley & Sons, Zenkewcz OC, L XK, Nakazawa S. Dynamc transent analyss by a mxed, teratve method. Internatonal Journal for Numercal Methods n Engneerng, 986; 23(7):

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE 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 information

Indeterminate pin-jointed frames (trusses)

Indeterminate 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 information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 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 information

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

Statistical 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 information

Finite Element Modelling of truss/cable structures

Finite 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 information

Numerical Heat and Mass Transfer

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 information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. 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 information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular 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 information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE 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 information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application 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 information

WAVE PROPAGATION, REFLECTION AND TRANSMISSION IN CURVED BEAMS

WAVE PROPAGATION, REFLECTION AND TRANSMISSION IN CURVED BEAMS ICSV4 Carns Australa 9- July, 7 WAVE PROPAGATION, REFECTION AND TRANSMISSION IN CURVED BEAMS Seung-Kyu ee, Bran Mace and Mchael Brennan NVH Team, R&D Centre, Hankook Tre Co., td. -, Jang-Dong, Yuseong-Gu,

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 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 information

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

NUMERICAL 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 information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

More information

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Appendix B. The Finite Difference Scheme

Appendix 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

Inductance Calculation for Conductors of Arbitrary Shape

Inductance 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 information

Lecture Note 3. Eshelby s Inclusion II

Lecture 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 information

Global Sensitivity. Tuesday 20 th February, 2018

Global 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 information

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The 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 information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 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 information

THREE-DIMENSION DYNAMIC SOIL-STRUCTURE INTERACTION ANALYSIS USING THE SUBSTRUCTURE METHOD IN THE TIME DOMAIN

THREE-DIMENSION DYNAMIC SOIL-STRUCTURE INTERACTION ANALYSIS USING THE SUBSTRUCTURE METHOD IN THE TIME DOMAIN The 4 th October 2-7 2008 Bejng Chna THREE-DIENSION DYNAIC SOIL-STRUCTURE INTERACTION ANALYSIS USING THE SUBSTRUCTURE ETHOD IN THE TIE DOAIN X..Yang Y.Chen and B.P.Yang 2 Assocate Professor aculty of Constructon

More information

DYNAMIC BEHAVIOR OF PILE GROUP CONSIDERING SOIL-PILE-CAP INTERACTION

DYNAMIC BEHAVIOR OF PILE GROUP CONSIDERING SOIL-PILE-CAP INTERACTION October 1-17, 8, Bejng, Chna DYNAMIC BEHAVIOR OF PILE GROUP CONSIDERING SOIL-PILE-CAP INTERACTION A. M. Halaban 1 and M. Malek 1 Professor, Faculty of Cvl Engneerng, Isfahan Unversty of Technology, Isfahan,

More information

Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM

Simulation 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 information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-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 information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development

More information

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

OFF-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 information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-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 information

General displacement arch-cantilever element method for stress analysis of arch dam

General 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 information

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,

More information

Lecture 12: Discrete Laplacian

Lecture 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 information

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Visco-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 information

Modeling of Dynamic Systems

Modeling 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 information

The Finite Element Method

The 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 information

Adjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.

Adjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R. th World Congress on Structural and Multdscplnary Optmsaton 7 th - th, June 5, Sydney Australa Adjont Methods of Senstvty Analyss for Lyapunov Equaton Bopng Wang, Kun Yan Department of Mechancal and Aerospace

More information

arxiv:cs.cv/ Jun 2000

arxiv:cs.cv/ Jun 2000 Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC 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 information

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

STATIC 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 information

A new Approach for Solving Linear Ordinary Differential Equations

A 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 information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study 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 information

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

GEO-SLOPE International Ltd, Calgary, Alberta, Canada   Vibrating Beam GEO-SLOPE Internatonal Ltd, Calgary, Alberta, Canada www.geo-slope.com Introducton Vbratng Beam Ths example looks at the dynamc response of a cantlever beam n response to a cyclc force at the free end.

More information

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions IOSR Journal of Mechancal and Cvl Engneerng (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 30-334X, Volume 15, Issue 5 Ver. IV (Sep. - Oct. 018), PP 41-46 www.osrjournals.org Bucklng analyss of sngle-layered FG nanoplates

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION 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 information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor 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 information

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model Identfcaton of Instantaneous Modal Parameters of A Nonlnear Structure Va Ampltude-Dependent ARX Model We Chh Su(NCHC), Chung Shann Huang(NCU), Chng Yu Lu(NCU) Outlne INRODUCION MEHODOLOGY NUMERICAL VERIFICAION

More information

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION VOL. 6, NO., NOVEMBER ISSN 89-668 6- Asan Research Publshng Network (ARPN). All rghts reserved. ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION Adel A.

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer 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 information

A new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems

A new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems Home Search Collectons Journals About Contact us My IOPscence A new ntegrated-rbf-based doman-embeddng scheme for solvng flud-flow problems Ths artcle has been downloaded from IOPscence. Please scroll

More information

In this section is given an overview of the common elasticity models.

In 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 information

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM

NONLINEAR 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 information

Canonical transformations

Canonical 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 information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

Instability Analysis of Damaged Pile Due to Static or Dynamic Overload

Instability Analysis of Damaged Pile Due to Static or Dynamic Overload Geomaterals, 0,, 4-0 http://.do.org/0.436/gm.0.406 Publshed Onlne October 0 (http://www.scrp.org/journal/gm) Instablty Analyss of Damaged Ple Due to Statc or Dynamc Overload P. N. Jk *, J. U. Agber, Cvl

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

One Dimensional Axial Deformations

One 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 information

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites Insttute of Structural Mechancs A comprehensve study: Boundary condtons for representatve volume elements (RVE) of compostes Srhar Kurukur A techncal report on homogenzaton technques A comprehensve study:

More information

A Hybrid Variational Iteration Method for Blasius Equation

A 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 information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 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 information

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

APPROXIMATE 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 information

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

THE 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 information

Lecture 8 Modal Analysis

Lecture 8 Modal Analysis Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal

More information

Computers and Mathematics with Applications

Computers and Mathematics with Applications Computers and Mathematcs wth Applcatons 56 (2008 3204 3220 Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: www.elsever.com/locate/camwa An nnovatve egenvalue

More information

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments. CE7 Structural Analyss II PAAR FRAE EEET y 5 x E, A, I, Each node can translate and rotate n plane. The fnal dsplaced shape has ndependent generalzed dsplacements (.e. translatons and rotatons) noled.

More information

Chapter 3. Estimation of Earthquake Load Effects

Chapter 3. Estimation of Earthquake Load Effects Chapter 3. Estmaton of Earthquake Load Effects 3.1 Introducton Sesmc acton on chmneys forms an addtonal source of natural loads on the chmney. Sesmc acton or the earthquake s a short and strong upheaval

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment 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 information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department

More information

1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave

1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave 014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON 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 information

On Pressure Distributions of Drum Brakes

On Pressure Distributions of Drum Brakes Yuan Mao Huang Professor e-mal: ymhuang@ccmsntuedutw J S Shyr Research Assstant Department of Mechancal Engneerng, Natonal Tawan Unversty, Tape, Tawan, Republc of hna On Pressure Dstrbutons of Drum Brakes

More information

PHYS 705: Classical Mechanics. Canonical Transformation II

PHYS 705: Classical Mechanics. Canonical Transformation II 1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE

GEOSYNTHETICS 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 information

arxiv: v1 [math.co] 12 Sep 2014

arxiv: v1 [math.co] 12 Sep 2014 arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March

More information

Free vibration analysis of a hermetic capsule by pseudospectral method

Free vibration analysis of a hermetic capsule by pseudospectral method Journal of Mechancal Scence and echnology 6 (4) (0) 0~05 wwwsprngerlnkcom/content/78-494x DOI 0007/s06-0-06-y Free vbraton analyss of a hermetc capsule by pseudospectral method Jnhee ee * Department of

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION 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 information

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD SIMUATION OF WAVE POPAGATION IN AN HETEOGENEOUS EASTIC OD ogéro M Saldanha da Gama Unversdade do Estado do o de Janero ua Sào Francsco Xaver 54, sala 5 A 559-9, o de Janero, Brasl e-mal: rsgama@domancombr

More information

High resolution entropy stable scheme for shallow water equations

High 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 information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi 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 information

New Method for Solving Poisson Equation. on Irregular Domains

New 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 information

Numerical Nonlinear Analysis with the Boundary Element Method

Numerical 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 information

Constitutive Modelling of Superplastic AA-5083

Constitutive 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 information

NUMERICAL DIFFERENTIATION

NUMERICAL 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 information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, 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 information

Relaxation Methods for Iterative Solution to Linear Systems of Equations

Relaxation 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 information

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

More information

Construction of Serendipity Shape Functions by Geometrical Probability

Construction of Serendipity Shape Functions by Geometrical Probability J. Basc. Appl. Sc. Res., ()56-56, 0 0, TextRoad Publcaton ISS 00-0 Journal of Basc and Appled Scentfc Research www.textroad.com Constructon of Serendpty Shape Functons by Geometrcal Probablty Kamal Al-Dawoud

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A 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 information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Mathematical Preparations

Mathematical 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 information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

Tensor Smooth Length for SPH Modelling of High Speed Impact

Tensor 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 information

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

A 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 information

Combined Wronskian solutions to the 2D Toda molecule equation

Combined Wronskian solutions to the 2D Toda molecule equation Combned Wronskan solutons to the 2D Toda molecule equaton Wen-Xu Ma Department of Mathematcs and Statstcs, Unversty of South Florda, Tampa, FL 33620-5700, USA Abstract By combnng two peces of b-drectonal

More information

Prediction of Chatter Stability for Milling Process Using an Improved Semi-discretization Method

Prediction of Chatter Stability for Milling Process Using an Improved Semi-discretization Method 5th Internatonal Conference on Advanced Engneerng Materals and Technology (AEMT 05) Predcton of Chatter Stablty for Mllng Process Usng an Improved Sem-dscretzaton Method LI Zhongqun, a, ZHU Fan,b, XIA

More information

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding. Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays

More information