OPTIMUM DESIGN OF BRIDGE SYSTEM SUBJECTED TO DEVASTATING EARTHQUAKE CONSIDERING PERFORMANCE AT ULTIMATE STATE

Size: px
Start display at page:

Download "OPTIMUM DESIGN OF BRIDGE SYSTEM SUBJECTED TO DEVASTATING EARTHQUAKE CONSIDERING PERFORMANCE AT ULTIMATE STATE"

Transcription

1 OTIMUM DESIGN OF BRIDGE SYSTEM SUBJECTED TO DEVASTATING EARTHQUAKE CONSIDERING ERFORMANCE AT ULTIMATE STATE Kzuiro Tniw, Kent Tnk Associte rofessor,, Deprtment of Civil nd Environmentl Eng., Fukui University of Tecnology, Gkuen 3-6-, Fukui , Jpn E-mil: Telepone: Cief Exminer, Codi Co. Ltd., ootemci -8-4, Nk-ku, Hirosim , Jpn E-mil: Telepone: Abstrct: In tis study, rtionl nd efficient optiml seismic design metod for bridge system subjected to devstting ertqukes considering performnce t ultimte stte is proposed. Te bridge system consists of superstructure, rubber berings, RC piers nd cst-in-plce concrete pile foundtion. In te proposed optimum design metod, te optimum solutions for te eigts of rubber berings, cross-sectionl dimensions nd mount of steel reinforcements for RC piers nd te detil of concrete pile foundtion re determined for severl llowble ductile fctors of RC piers considering te constrints on te reltive orizontl displcements of rubber berings to te bot bridge nd trnsverse directions, te ductile fctor of RC piers, nd te constrint on te cst-in-plce concrete pile foundtion. From te prcticl design te eigts of rubber berings cn tke continuous vlues, but te oter vribles must be selected from discrete vrible sets. Terefore, te construction cost minimiztion problem cn be expressed s mixed discrete-continuous problem. Tis problem is trnsformed into convex pproximtion problem wit te estimtion formule by using te experimentl design, nd te dynmic beviors nd tose sensitivities re clculted nlyticlly by using te estimtion formule witout nlyzing te structures. Te optimum design problem is solved by clssicl brnc nd bound metod wit dul lgoritm. In te numericl design exmples, it is empsized tt te optimum solutions cn be obtined efficiently by using te experimentl design. It is lso demonstrted tt te reductions of te eigts of rubber berings nd cross-sectionl dimensions of RC piers cn be observed by incresing te llowble ductility fctor. Keywords: Bridge system, Optimiztion, Seismic design, Design of experiments, erformnce t ultimte stte. Introduction After te Hyogoken Nnbu Ertquke in 995, te seismic design code for igwy bridges, JSHB [], s been revised in order to ensure sufficient ultimte dynmic cpcities in te bridge systems for lrge displcements cused by devstting ertqukes. Recently, te performnce-bsed design metod s been introduced for te seismic design t ultimte stte in te JSHB. According to te JSHB, te bridge members re not llowed to yield for te frequent ertqukes (Level), nd te bridge members must ve te sufficient ultimte dynmic cpcities to be ble to repir tose rpidly fter te excittions due to devstting ertqukes (Level). Tis tsk ccompnies wit tremendous complexity in te process of design of te bridge system. In generl, Level design is criticl for determintions of member sizes for lrge scle bridge systems. Terefore, te estblisment of rtionl nd efficient optiml seismic design metod, wic cn determine te optimum member sizes considering performnce t ultimte stte in te Level design process, s been wited expectntly in te prcticl design. From tis point of view, one of utors proposed n optiml seismic design metod using te design of experiments nd suboptimiztion tecnique [,3]. In tis reserc works, utors mde effort to introduce severl reltions between construction cost nd design vribles to mke te optimiztion problem simple. Te design vribles for bridge members re delt wit s continuous vribles nd te optimum solutions considering te displcement constrints for bridge direction re determined. In tis study, n rtionl nd efficient optiml performnce-bsed seismic design metod for bridge system subjected to devstting ertqukes is proposed. In te design of bridge system, te Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00

2 dimensions of superstructure re ssumed to be given, nd te eigts of rubber berings, cross-sectionl dimensions nd mount of steel reinforcements for RC piers, nd numbers of piles nd te dimeters of piles in te cst-in-plce concrete pile foundtion re tken into ccount s design vribles. Te dynmic nonliner beviors of te bridge system re nlyzed precisely by using te generl purpose nonliner nlysis softwre (TDA-III) wit te ccelertion specified in te JSHB. Te reltive orizontl displcements between superstructure nd piers to te bot bridge nd trnsverse directions re delt wit s design constrints for te rubber berings. Te ductile fctor, wic is given by te rtio of worng curvture to te yield curvture, is delt wit s te design constrints for te RC piers so s to ensure te performnce specified t te ultimte stte. Furtermore, te constrint on te cst-in-plce concrete pile foundtion is lso delt wit to ensure te sufficient ultimte dynmic cpcity in te RC pile foundtion. However, te constrint on te RC pile foundtion is not treted in te optimiztion process to simplify te optimiztion lgoritm. After determintion of optimum solution te constrint on te RC pile foundtion is exmined, nd te RC pile foundtion is replced wit te lrger one so s to stisfy te constrint. From te prcticl design te eigts of rubber berings cn tke continuous vlues, but te oter vribles must be selected from discrete vrible sets. Terefore, te construction cost minimiztion problem cn be expressed s mixed discrete-continuous problem, nd it is solved by clssicl brnc nd bound metod [4] wit dul lgoritm nd convex pproximtion [5,6] in tis study. Te sensitivities of te design constrints need in te optimiztion process nd we encounter te difficulty to obtin tose in utilizing te generl purpose nonliner nlysis softwre. To overcome tis problem te design of experiments is pplied successfully in order to clculte te dynmic beviors nd tose sensitivities in te optimiztion process. In te design of experiments, te estimtion formule for dynmic beviors re introduced in te expression of qudrtic functions of te design vribles. Te dynmic beviors nd tose sensitivities re clculted nlyticlly by using te estimtion formule witout nlyzing te structures. After te determintion of optimum solution te design constrints wit te estimtion formule re exmined by re-nlyzing te bridge system using TDA-III. In cse tt te design constrints violte te llowble limits, te estimtion formule for dynmic beviors re improved nd te minimum cost design problem is re-solved. Tis optimiztion process is iterted until te reltive errors between te estimted design constrins nd te exct ones stisfy te llowble limits. Te proposed optiml design metod is pplied to five-spn continuous steel girder bridge system, nd te optiml solutions t vrious llowble ductility fctors of RC pier re compred. In te numericl results, it is demonstrted tt te reductions of te eigts of rubber berings nd cross-sectionl dimensions cn be observed by incresing te llowble ductility fctor. It is lso empsized tt te optimum solutions cn be obtined efficiently t few itertions of improvements of te estimtion formule for dynmic beviors. Te ccurcy of te estimtion formule is excellent H H H H B/ oint of ppliction to te force of inerti of superstructure 00 B 6000-B/ H A 3 4 A Fig. Five-spn continuous steel girder bridge system () Front view (b) Side view Fig. Front nd side views of piers nd RC pile foundtion Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00

3 M H My H H-H 5 Mc Kd H 5 φ y φ C 0 6 φ C 3 φ y 9 φ B Fig.3 Cross section of pier witin 7 percent of reltive errors between te exct beviors nd estimted ones.. OTIMUM DESIGN FORMULATION AND OTIMIZATION ALGORITHM. Design Model In tis study, te five-spn continuous steel girder bridge system sown in Fig. is considered in wic te superstructure is supported by six rubber berings, RC piers nd te cst-in-plce concrete pile foundtion. Te front nd side views of pier nd RC pile foundtion re described in Fig.. Te lengts of piles re 5m nd five types of soil conditions in strtum re considered to clculte spring constnts. Te reinforcements in te cross section of piers re rrnged in two lyers for te bridge direction nd one lyer for te trnsverse direction, nd te intervl of ec reinforcement re fixed t 5mm s sown in Fig.3. Following n enlrgement of cross section te numbers of reinforcements increse so s to keep te intervls of reinforcements. Te stiffness of RC pier is tken into ccount s te triliner rigidity reduction type model (Tked model) sown in Fig.4. Te nonliner beviors of te bridge system for te bot bridge nd trnsverse directions subjected to devstting ertqukes re nlyzed precisely by using TDA-III in wic te Type II stndrd strong ccelertion wve motion model t te Type II soil ground specified in te JSHB is pplied. In te time-istory response nlysis te spring constnts of rubber berings, pile foundtions nd superstructure re elstic, nd bot te superstructure nd butment re ssumed s rigid body. Te piers re divided into 50 segments in order to clculte te nonliner dynmic beviors ccurtely.. Optimum Design Formultion In te design of te bridge system, te dimension of superstructure is ssumed to be given nd widts of rectngulr rubber berings re ssumed to be 70cm nd 80cm t butment nd piers, respectively. Te design vribles for rubber berings re te eigts of tose t butment nd piers, B nd B. For te cst-in-plce concrete pile foundtions te numbers of piles nd dimeters of piles re intensively summrized s te properties of orizontl nd rottion spring constnts. In tis study te orizontl spring constnts of RC pile foundtion, K, wic cn be commonly used for te time-istory response nlysis to te bot bridge nd trnsverse directions, re considered s te design vribles. Te widts to te bridge nd trnsverse directions nd te mount of steel reinforcements in cross section, H, B nd A s, re tken into ccount s te design vribles for RC piers. Te bridge system sown in Fig. is symmetricl to te centerline nd te totl number of design vribles is six of B, B, K, H, A, B S. B-B B B Fig.4 Triliner rigidity ysteresis model for RC pier (Tked model) Engineers ve to design te bridge system wic ve sufficient ultimte dynmic cpcities for lrge displcements cused by devstting ertqukes. Terefore, te reltive orizontl displcements between superstructure nd piers to te bot bridge nd trnsverse directions re delt wit s te design Mc φ, Μ: Crcng curvture nd moment -My c φ, Μ: Initil yieldcurvture nd moment y c y M C + M y φmx kd = φ + φ φ c y y 0.4 Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00 3

4 Tble roperties of tree types of RC piles Dimeter φ Number of piles Widt of footing B Widt of footing H Heigt of footing Construction cost(0 3 yen) K(kN/m) K θ (knm/rd) (bridge direction) K θ (knm/rd) (trnsverse direction) weigt(kn).0m 9 7.0m 7.0m.5m 3, m 9 8.4m 8.4m.5m 6, m 7.0m 9.5m.5m 7, constrints, g, g, gt, gt, for te sfety of te rubber berings. Furtermore, te ductile fctors re lso delt wit s te design constrints for te RC piers, g μ, so s to ensure te performnce specified t te ultimte stte. In te design of RC pile foundtion, for te cse tt te orizontl ultimte dynmic bering cpcity for te RC pier is not enoug lrge for te orizontl force clculted using te specified design seismic coefficient, RC pile foundtion is not llowed to yield wen te equivlent lods corresponding to te orizontl ultimte dynmic bering cpcity for te RC pier re pplied to te RC pile foundtion. For te cse tt te orizontl ultimte dynmic bering cpcity for te RC pier is sufficient, RC pile foundtion is llowed to yield up to te ductile fctor 4.0. Tis constrint is quite complex to del wit in te optimiztion process. Furtermore, we need to consider tt te design vrible for RC pile foundtion is dependent on te design vribles for RC piers. To simplify te optimum design problem, terefore, it is ssumed tt te design vrible for RC pile foundtion is independent, nd te constrint on te RC pile foundtion is not delt wit in te optimiztion process. After te determintion of optimum solution te constrint on te RC pile foundtion is exmined. Te totl construction cost minimiztion problem, wic is expressed s te summtion of bering construction cost, COST B ( B,B ),pier construction cost, COST F ( K ), nd pier construction cost, COST ( H, B ), cn be formulted s find B, B, K, H, B wic minimize COST ( B, B, K, H, B ) = COST B ( B, B ) + COSTF ( K ) + COST ( H, B ) () subject to g δ δ 0 (), g δ δ 0 (3), g δ δ 0 (4) = = t = t gt = δ t δ 0 (5), g μ = μ μ 0 (6), were δ nd δ re te llowble reltive orizontl displcements of berings t butment nd piers, wic re given s te products of te eigts of berings B, B multiplied by.5. μ is te ductile fctor of pier, wic is given by te rtio of worng curvture to te yield curvture for te bridge direction. In te optimum design problem B nd B cn tke continuous vlues, but te oters must be selected from list of discrete vlues. In tis study, K, H nd B re selected from te following discrete sets in wic tree types of pile foundtions summrized in Tble re considered to clculte K. { 657( kn / m), 76477, 9500} K H { 000( mm), 00, 00, 300, 400, 500, 600, 700, 800, 900,3000} A S { 98.6( mm ), 86.5, 387., 506.7, 64.4, 794., 956.6, 40} B { 3000( mm), 3500, 4000, 4500, 5000, 5500, 6000, 6500} Terefore, te construction cost minimiztion problem cn be expressed s mixed discrete-continuous problem. Severl types of optimiztion tecniques ve been developed, nd Hung nd Aror [4] investigted te efficiency nd relibility of tose for discrete nd mixed discrete-continuous problems. Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00 4

5 Replcement of pile foundtion Fig.5 Mcro-flow of te proposed optimum design metod In tis study te optimiztion problem is solved by te clssicl brnc nd bound metod wit dul lgoritm nd convex pproximtion [5,6] for te reson tt te pproc is efficient nd relible for mixed discrete-continuous problem witout ny prmeters..3 Optimiztion Algoritm In tis optimiztion process, in generl, number of nonliner seismic response nlyses nd sensitivity nlyses re necessry to determine te optiml solutions. To void tese complexity nd difficulties nd mke te optimum design process tremendously efficient, te design of experiments [7] is pplied to introduce te estimtion formule for te dynmic beviors. Te dynmic beviors nd tose sensitivities re clculted by using te estimtion formule witout nlyzing te structure. In te design of experiments, ccording to te ortogonl rry tble L ( ) [7] given in Tble, te tree levels for ll design vribles re ssumed nd te twenty seven runs of nonliner seismic response nlyses re crried out in usge of TDA-III for te bot bridge nd trnsverse directions, respectively. Te first six fctors mong tirty fctors in Tble re ssigned to te design vribles B, B, K, H, B, respectively. Assuming tt te intended vrible for te kt fctor is x k nd te men vlue of tree levels ( xˆ, i =, L,3) for te kt fctor is x k, te generl form of estimtion formul is introduced in te expression of qudrtic functions of te design vribles given s eqs.(7)-(0). y = b m m 0 + bkzk + bk Mk Mk3zk + Mkzk ) k= k= (, (7) i i i were M ( ˆ ˆ ˆ = zk + zk + + zkn) ( k =, L, m), (8) n zˆ ˆ = x xk ( i =, L, n) ( k =, L, m) (9), zk = xk xk, (0) m nd n re respectively te number of fctors, i.e. te number of design vribles (= 6), nd te number of levels for ec fctor (= 3). Te estimtion vlues of b 0, b k nd b k in eq.(7) re given s n b ˆ 0 = T i rs, n b ˆ kl = WT ( l =,) rs () were = n S k W, () i= Initil tree levels Clcultion of eigenvlue Nonliner time-istory response nlysis Estimtion formule of dynmic beviors by te experimentl design Formultion of optimum design problem Optimiztion by te brnc nd bound metod NO NO Convergence? YES Clcultion of eigenvlue Nonliner time-istory response nlysis NO Accurte? YES Stisfy te pile foundtion constrints? YES stop k i= Improvement of tree levels Tble Ortogonl rry tble L ( 3 3 ) No No. No. No No No No No No No No No No No No No No No No No No No No No No No No No Fctor:, Fctor: B, Fctor3: K, Fctor4: H, Fctor5: A, Fctor6: B B r is te number of runs wit te level xˆ (= 9). T is te summtion of results by te design of experiments wit te level of xˆ. W is te vlue of function of coefficient f k ( z k ) in eq.(7) wit respect to b k nd b k were zk = zˆ, nmely, W = zˆ nd W = M M zˆ + M z. Fctor ˆ k k3 k i= 7 S Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00 5

6 After te determintion of optimum solutions te design constrints wit te estimtion formule re exmined by re-nlyzing te bridge system. In cse tt te design constrints violte te llowble limit, te tree levels for ll design vribles nd estimtion formule for dynmic beviors re improved nd te minimum cost design problem is re-solved. Tis optimiztion process is iterted until te reltive errors between te estimted design constrins nd te exct ones stisfy te llowble limit. After te determintion of optimum solution te constrint on te RC pile foundtion is exmined. In te cse tt te constrint is violted te RC pile foundtion is replced wit te lrger one nd te bridge system is re-optimized. Te mcro-flow of te proposed optimiztion lgoritm is depicted in Fig.5. Levels Itertion Itertion Tble 3 Improvements of tree levels in te optimiztion process B (cm) (spring constnt(kn/m)) B (cm) (spring constnt(kn/m)) spring constnt of pile K (kn/m) K θ (knm/rd) K θ (knm/rd) H (mm) A s(mm ) B (mm) 6.0(533) 4.0(857) (7500).0(6667) (047) 0.0(3000) (3065) 8.0(40000) (4500) 9.0(35556) (047) 0.0(3000) Allowble ductile fctors B (spring constnt) B B Tble 4 Optimum solutions for μ =.0, 3.0 nd 4.0 μ (spring constnt) K (φ, n).0 4.8cm (773kN/m).8cm (8366kN/m) (φ=.m, n=9) 700mm cm (8035kN/m) 8.0cm (40000kN/m) (φ=.m, n=9) 800mm 8.37cm (957kN/m) 8.0cm (40000kN/m) 657 (φ=.0m, n=9) 700mm As 40mm 387.mm 387.mm H 4000mm 4500mm 3000mm δ / δ δ / δ δ t / δ δ t / δ μ / μ yield of pile foundtion D.exp. *.000 D.exp. *.000 D.exp. * 0.99 Anl ** 0.97 Anl **.0 Anl **.03 D.exp. * 0.9 D.exp. * D.exp. * 0.43 Anl ** Anl ** 0.60 Anl ** D.exp. * D.exp. * D.exp. *.000 Anl ** Anl **.0 Anl ** D.exp. * D.exp. * D.exp. * Anl ** Anl **.038 Anl ** 0.83 D.exp. *.000 D.exp. * D.exp. * 0.97 Anl **.0 Anl **.000 Anl **.008 Bridge dir: μfr=.5 Trnsverse dir:μfr=.6 Bridge dir: not yield Trnsverse dir:μfr=.6 not yield Totl cost (0 3 yen) D.exp.* : Fesibility of design constrints wit te estimtion formule by te design of experiments Anl** : Fesibility of design constrints using exct beviors by nlysis Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00 6

7 3. DESIGN EXAMLES Te proposed optiml design metod is pplied to te five-spn continuous steel girder bridge system sown in Fig. nd te optiml solutions for severl llowble ductile fctors μ re compred. Te unit cost of rubber is s 45yen/cm 3. Te construction costs of pile re ssumed s 6500yen/m 3 for te dimeter.0m nd 73800yen/m 3 for te dimeter.m. Te construction costs of footing nd form for pile foundtion re ssumed s 33500yen/m 3 nd 8000yen/m. Te construction costs of concrete, form nd reinforcement for piers re ssumed s 8500yen/m 3, 8000yen/m nd 0000yen/tf, respectively. Following te flow-crt in Fig.5 te optimiztion processes for μ =. 0, 3.0 nd 4.0 re initited wit te levels of itertion sown in Tble 3. In te optimiztion process, te lower nd upper limits for discrete design vribles re set t te djcent discrete vlues of te minimum nd mximum vlues of te tree levels to limit improvements of design vribles. Te optimum solution for μ =. 0 cn be obtined quite efficiently witout ny improvements of te tree levels for ll design vribles. Te optimum solutions for μ = 3. 0 nd 4.0 determined by te lower limits set s te move limits. After ten, te tree levels re improved to te vlues of itertion in Tble 3 referring to te optimum solutions wit te previous tree levels. Te optimum solutions for μ = 3. 0 nd 4.0 cn be obtined efficiently t tis stge witout dditionl improvements of te tree levels. Te optimum solutions for μ =.0, 3.0 nd 4.0 re summrized in Tble 4. Te orizontl spring constnts of pile foundtions for ll cses re determined by te lower limit wic indictes te lowest cost. Ten, te RC pile foundtions for μ =. 0 nd 3.0 re replced wit te lrger one so s to stisfy te constrint on te RC pile foundtion. In cse of μ =. 0 te lrgest dimensions of cross section, H nd B, nd reinforcement in te piers A S re required in order to stisfy te llowble ductile fctor. By incresing te eigts of rubber berings B nd B, nmely reducing te vlues of spring constnt, te period of bridge system is mde longer nd te effect from superstructure is minimized. As te result te totl construction cost is minimized. In cse of μ = 3. 0 B, B nd A S re remrkbly reduced compred wit tose in cse of μ =. 0, nd B nd A S re determined by te lower limits. Te totl cost is reduced to 8 percent of tt in cse of μ =. 0. In cse of μ = 4. 0 B, K, AS nd H re determined by te lower limits. Te totl cost is reduced to 7 percent of tt in cse of μ =. 0. As clerly seen from te vlues of fesibility of design constrints using exct beviors by nlysis in Tble4, bot te constrints on reltive orizontl displcements t butment to te bridge direction g nd ductile fctors g μ re ctive for ll cses simultneously. Te displcements t butment nd piers to te trnsverse direction g t re lso ctive for μ = Te displcements t piers to te bridge direction g re inctive for ll cses. Te ccurcy of te estimtion formule is excellent witin 7 percent of reltive errors. Te exct constrints re enoug fesible witin 3.8 percent of violtion for ll cses. 4. CONCLUSIONS Te following conclusions cn be drwn from tis study: ) Te proposed optiml design metod cn determine te eigts of rubber berings, cross-sectionl dimensions nd mount of steel reinforcements for RC piers, nd numbers nd dimeters of piles rigorously nd efficiently. Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00 7

8 ) By pplying te design of experiments, te estimtion formule for te ductile fctor in piers nd te mximum orizontl displcements to te bridge nd trnsverse directions cn be introduced ccurtely wit smll number of nonliner seismic response nlyses. Te ccurcy of te estimtion formule is excellent witin 7 percent of reltive errors between te exct beviors nd estimted ones. 3) A few itertions of improvements for tree levels re required to obtin te optimum solutions in te proposed design metod. 4) In te cse tt te llowble ductile fctor is set t smll vlue, te eigts of rubber berings increse in order to mke te period of bridge system longer, nd te effect from superstructure is minimized. As incresing te vlue of llowble ductile fctor te eigts of rubber berings re reduced nd te dimension of cross section nd reinforcement in te piers re lso reduced. 5) In te proposed design process, te constrint on te RC pile foundtion is not delt wit in te optimiztion process nd, ten, te RC pile foundtion is replced wit te lrger one so s to stisfy te constrint on te RC pile foundtion. Tis design process cn simplify te optimiztion lgoritm gretly. 6) Te constrints on reltive orizontl displcements t butment to te bridge direction g nd ductile fctors g μ re ctive t te optimum solutions simultneously. Acknowledgements rt of tis work s been supported by FUT Reserc romotion Fund. References. Jpn Rod Assocition, Specifiction for Higwy bridges, rt V seismic design, Tokyo, Mruzen (in Jpnese), 00.. Okubo,S., Tnk,K., Wtnbe,T., Yositke,R., Optuml design metod for seismic-isoltion bridge systems subject to uge ertqukes, Journls of JSCE., No.703/I-59, pp.67-8 (in Jpnese), Okubo,S., Tnk,K., Kdot,K., Minimum cost design metod for seismic-isoltion bridge systems considering initil construction cost nd repir cost cused by uge ertquke, Journls of JSCE., No.70/I-60, pp.9-08 (in Jpnese), Hung,M., Aror,J.S., Optiml design wit discrete vribles: some numericl exmples, Int. J. Num. Met. Eng., 40, pp.65-88, Fleury,C. nd Bribnt,V., Structurl optimiztion: A new dul metod using mixed vribles, Int. J. Num. Met. Eng., 3, pp , Tniw,K. nd Okubo,S.: Optiml syntesis metod for trnsmission tower truss structures subjected to sttic nd seismic lods, Int. J. Struct. Multidisc. Optim. Vol.6, pp , Tguci,G., System of Experimentl Design, New York, UNIUB/Krss Interntionl ublictions, Volume &, 987. Interntionl Conference on Sustinble Built Environment (ICSBE-00) Kndy, 3-4 December 00 8

A Discussion on Formulas of Seismic Hydrodynamic Pressure

A Discussion on Formulas of Seismic Hydrodynamic Pressure Interntionl Forum on Energy Environment Science nd Mterils (IFEESM 2017) A Discussion on Formuls of Seismic Hydrodynmic Pressure Liu Himing1 To Xixin2 1 Cin Mercnts Congqing Communiction Reserc & Design

More information

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy.

Key words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy. O SPECTRA ACCURACY OF QUADRATURE FORMUAE BASED O PIECEWISE POYOMIA ITERPOATIO A KURGAOV AD S TSYKOV Abstrct It is well-known tt te trpezoidl rule, wile being only second-order ccurte in generl, improves

More information

Chapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS

Chapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS S Cpter Numericl Integrtion lso clled qudrture Te gol of numericl integrtion is to pproximte numericlly. f(x)dx Tis is useful for difficult integrls like sin(x) ; sin(x ); x + x 4 Or worse still for multiple-dimensionl

More information

New data structures to reduce data size and search time

New data structures to reduce data size and search time New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute

More information

Design Data 1M. Highway Live Loads on Concrete Pipe

Design Data 1M. Highway Live Loads on Concrete Pipe Design Dt 1M Highwy Live Lods on Concrete Pipe Foreword Thick, high-strength pvements designed for hevy truck trffic substntilly reduce the pressure trnsmitted through wheel to the subgrde nd consequently,

More information

Shear and torsion interaction of hollow core slabs

Shear and torsion interaction of hollow core slabs Competitive nd Sustinble Growth Contrct Nº G6RD-CT--6 Sher nd torsion interction of hollow core slbs HOLCOTORS Technicl Report, Rev. Anlyses of hollow core floors December The content of the present publiction

More information

Experimental Study, Stiffness of Semi-Rigid Beam-to-Column Connections Using Bolts and Angles

Experimental Study, Stiffness of Semi-Rigid Beam-to-Column Connections Using Bolts and Angles nd rd Interntionl Conference on Electricl, Electronics nd Civil Engineering (ICEECE'1) Jnury 4-5, 1 Bli (Indonesi) Experimentl Study, Stiffness of Semi-Rigid -to- s Using Bolts nd Angles Khled M. Amtered

More information

IMPULSE-BASED SIMULATION OF INEXTENSIBLE CLOTH

IMPULSE-BASED SIMULATION OF INEXTENSIBLE CLOTH IMPULSE-BASED SIMULATION OF INEXTENSIBLE CLOTH Jn Bender nd Dniel Byer Institut für Betries- und Dilogssysteme Universität Krlsrue Am Fsnengrten 5 76128 Krlsrue Germny ABSTRACT In tis pper n impulse-sed

More information

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University

Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur Module Anlysis of Stticlly Indeterminte Structures by the Mtrix Force Method Version CE IIT, Khrgpur esson 8 The Force Method of Anlysis: Bems Version CE IIT, Khrgpur Instructionl Objectives After reding

More information

Chapter 18 Two-Port Circuits

Chapter 18 Two-Port Circuits Cpter 8 Two-Port Circuits 8. Te Terminl Equtions 8. Te Two-Port Prmeters 8.3 Anlysis of te Terminted Two-Port Circuit 8.4 nterconnected Two-Port Circuits Motivtion Tévenin nd Norton equivlent circuits

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Topic 6b Finite Difference Approximations

Topic 6b Finite Difference Approximations /8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) 747 6958 E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions?

More information

Multi-objective optimization of dielectric layer photonic crystal filter

Multi-objective optimization of dielectric layer photonic crystal filter Optic Applict, Vol. XLVII, No. 1, 017 DOI: 10.577/o170103 Multi-objective optimiztion of dielectric lyer photonic crystl filter HONGWEI YANG *, CUIYING HUANG, SHANSHAN MENG College of Applied Sciences,

More information

The Active Universe. 1 Active Motion

The Active Universe. 1 Active Motion The Active Universe Alexnder Glück, Helmuth Hüffel, Sš Ilijić, Gerld Kelnhofer Fculty of Physics, University of Vienn helmuth.hueffel@univie.c.t Deprtment of Physics, FER, University of Zgreb ss.ilijic@fer.hr

More information

Quantum transport (Read Kittel, 8th ed., pp )

Quantum transport (Read Kittel, 8th ed., pp ) Quntum trnsport (Red Kittel, 8t ed., pp. 533-554) Wen we ve structure in wic mny collisions tke plce s crriers trnsport cross it, te quntum mecnicl pse of te electron wvefunctions is essentilly rndomized,

More information

12. DYNAMIC ANALYSIS. Force Equilibrium is Fundamental in the Dynamic Analysis of Structures 12.1 INTRODUCTION

12. DYNAMIC ANALYSIS. Force Equilibrium is Fundamental in the Dynamic Analysis of Structures 12.1 INTRODUCTION 12. DYNAMIC ANALYSIS Force Equilibrium is Fundmentl in the Dynmic Anlysis of Structures 12.1 INTRODUCTION { XE "Newton's Second Lw" }All rel physicl structures behve dynmiclly when subjected to lods or

More information

Math Week 5 concepts and homework, due Friday February 10

Math Week 5 concepts and homework, due Friday February 10 Mt 2280-00 Week 5 concepts nd omework, due Fridy Februry 0 Recll tt ll problems re good for seeing if you cn work wit te underlying concepts; tt te underlined problems re to be nded in; nd tt te Fridy

More information

Comparison of the Design of Flexural Reinforced Concrete Elements According to Albanian Normative

Comparison of the Design of Flexural Reinforced Concrete Elements According to Albanian Normative ISBN 978-93-84422-22-6 Proceedings of 2015 Interntionl Conference on Innovtions in Civil nd Structurl Engineering (ICICSE'15) Istnbul (Turkey), June 3-4, 2015 pp. 155-163 Comprison of the Design of Flexurl

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

CHAPTER 4a. ROOTS OF EQUATIONS

CHAPTER 4a. ROOTS OF EQUATIONS CHAPTER 4. ROOTS OF EQUATIONS A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering by Dr. Ibrhim A. Asskk Spring 00 ENCE 03 - Computtion Methods in Civil Engineering II Deprtment

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

A P P E N D I X POWERS OF TEN AND SCIENTIFIC NOTATION A P P E N D I X SIGNIFICANT FIGURES

A P P E N D I X POWERS OF TEN AND SCIENTIFIC NOTATION A P P E N D I X SIGNIFICANT FIGURES A POWERS OF TEN AND SCIENTIFIC NOTATION In science, very lrge nd very smll deciml numbers re conveniently expressed in terms of powers of ten, some of wic re listed below: 0 3 0 0 0 000 0 3 0 0 0 0.00

More information

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical Analysis: Trapezoidal and Simpson s Rule nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =

More information

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge

More information

Stiffness Reduction Factor for Flat Slab Structures under Lateral Loads

Stiffness Reduction Factor for Flat Slab Structures under Lateral Loads TECHNICAL NOTES Stiffness Reduction Fctor for Flt Slb Structures under Lterl Lods Sng-Whn Hn, Ph.D., P.E. 1 ; Young-Mi Prk 2 ; nd Seong-Hoon Kee 3 Abstrct: Effective bem width model EBWM hs been widely

More information

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

INTRODUCTION. The three general approaches to the solution of kinetics problems are: INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The

More information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

R. I. Badran Solid State Physics

R. I. Badran Solid State Physics I Bdrn Solid Stte Physics Crystl vibrtions nd the clssicl theory: The ssmption will be mde to consider tht the men eqilibrim position of ech ion is t Brvis lttice site The ions oscillte bot this men position

More information

Fuzzy transform to approximate solution of boundary value problems via optimal coefficients

Fuzzy transform to approximate solution of boundary value problems via optimal coefficients 217 Interntionl Conference on Hig Performnce Computing & Simultion Fuzzy trnsform to pproximte solution of boundry vlue problems vi optiml coefficients Zr Alijni Mtemtics nd Sttistics Dept University of

More information

Kiel Probes. General Information

Kiel Probes. General Information Kiel s Generl Informtion Aerodynmic Properties Kiel probes re used to mesure totl pressure in fluid strem were te direction of flow is unknown or vries wit operting conditions. Teir correction fctor is

More information

Job No. Sheet 1 of 8 Rev B. Made by IR Date Aug Checked by FH/NB Date Oct Revised by MEB Date April 2006

Job No. Sheet 1 of 8 Rev B. Made by IR Date Aug Checked by FH/NB Date Oct Revised by MEB Date April 2006 Job o. Sheet 1 of 8 Rev B 10, Route de Limours -78471 St Rémy Lès Chevreuse Cedex rnce Tel : 33 (0)1 30 85 5 00 x : 33 (0)1 30 5 75 38 CLCULTO SHEET Stinless Steel Vloristion Project Design Exmple 5 Welded

More information

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of

More information

Finite Element Determination of Critical Zones in Composite Structures

Finite Element Determination of Critical Zones in Composite Structures Finite Element Determintion of Criticl Zones in Composite Structures Alexey I. Borovkov Dmitriy V. Klimshin Denis V. Shevchenko Computtionl Mechnics Lb., St. Petersburg Stte Polytechnicl University, Russi

More information

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data A07 Uncertinties in Locl Anisotropy Estimtion from Multi-offset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.

More information

2. General formula for Runge-Kutta methods

2. General formula for Runge-Kutta methods NOLTA, IEICE Pper Equivlent circuits for implicit unge-kutt metods in circuit simultors for nonliner circuits Ysuiko Toym 1), Ttusy Kuwzki 1, nd Jun Sirtki 1 Mkiko Okumur 1 1 Kngw Institute of Tecnology,

More information

Tests for the Ratio of Two Poisson Rates

Tests for the Ratio of Two Poisson Rates Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson

More information

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations AMATH 731: Applied Functionl Anlysis Fll 2009 1 Introduction Some bsics of integrl equtions An integrl eqution is n eqution in which the unknown function u(t) ppers under n integrl sign, e.g., K(t, s)u(s)

More information

Effect of soil profile modulus distribution on pile head lateral stiffness

Effect of soil profile modulus distribution on pile head lateral stiffness Proc. 18 th NZGS Geotechnicl Symposium on Soil-Structure Interction. d. CY Chin, Aucklnd Michel Pender University of Aucklnd, New Zelnd eywords: pile hed stiffness, effect of pile shft size, soil modulus

More information

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

PKK Energy chains. PKK 320 / 200 x 2340 / 200; 5Pz, 1PT55

PKK Energy chains. PKK 320 / 200 x 2340 / 200; 5Pz, 1PT55 Energy cins Te order for plstic energy cin sould contin te following dt: Type / dius x Lengt / Widt "Arrngement"; Seprtor rrngement 1 3 Type selection is in ccordnce wit dimeter nd quntity of te lines

More information

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below. Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we

More information

LABYRINTH SEALS DYNAMIC COEFFICIENTS AND CRITICAL SPEED RESEARCH BASED ON CHILDS MODEL

LABYRINTH SEALS DYNAMIC COEFFICIENTS AND CRITICAL SPEED RESEARCH BASED ON CHILDS MODEL Abstrct LABYINTH SALS DYNAMIC COFFICINTS AND CITICAL SPD SACH BASD ON CHILDS MODL M Wensheng, Hung Hi, Feng Guoqun AVIC Shenyng Aeroengine eserch Institute, Shenyng, Chin Keywords:lbyrinth sels, dynmics

More information

Flutter frequency based on bending - torsion coupling theory

Flutter frequency based on bending - torsion coupling theory Flutter frequency bsed on bending - torsion coupling theory *ZHENG Xin ¹, LIU Yu-Bin¹, CHEN Pu, SHEN Feng 3,ZHANG Sheng-Jun 3,nd FU Xing-Rong ¹ 1 College of Wter Conservncy nd Civil Engineering, Chin Agriculturl

More information

Deteriorating Inventory Model for Waiting. Time Partial Backlogging

Deteriorating Inventory Model for Waiting. Time Partial Backlogging Applied Mthemticl Sciences, Vol. 3, 2009, no. 9, 42-428 Deteriorting Inventory Model for Witing Time Prtil Bcklogging Nit H. Shh nd 2 Kunl T. Shukl Deprtment of Mthemtics, Gujrt university, Ahmedbd. 2

More information

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM ROMAI J., v.9, no.2(2013), 173 179 THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM Alicj Piseck-Belkhyt, Ann Korczk Institute of Computtionl Mechnics nd Engineering,

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

Effects of peripheral drilling moment on delamination using special drill bits

Effects of peripheral drilling moment on delamination using special drill bits journl of mterils processing technology 01 (008 471 476 journl homepge: www.elsevier.com/locte/jmtprotec Effects of peripherl illing moment on delmintion using specil ill bits C.C. Tso,, H. Hocheng b Deprtment

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

Study on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method

Study on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method Avilble online t www.sciencedirect.com Physics Procedi 4 () 9 97 Interntionl Conference on Applied Physics nd Industril Engineering Study on the Clcultion of Mgnetic Force Bsed on the Equivlent Mgnetic

More information

Explain shortly the meaning of the following eight words in relation to shells structures.

Explain shortly the meaning of the following eight words in relation to shells structures. Delft University of Technology Fculty of Civil Engineering nd Geosciences Structurl Mechnics Section Write your nme nd study number t the top right-hnd of your work. Exm CIE4143 Shell Anlysis Tuesdy 15

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model: 1 2 MIXED MODELS (Sections 17.7 17.8) Exmple: Suppose tht in the fiber breking strength exmple, the four mchines used were the only ones of interest, but the interest ws over wide rnge of opertors, nd

More information

Solutions to Supplementary Problems

Solutions to Supplementary Problems Solutions to Supplementry Problems Chpter 8 Solution 8.1 Step 1: Clculte the line of ction ( x ) of the totl weight ( W ).67 m W = 5 kn W 1 = 16 kn 3.5 m m W 3 = 144 kn Q 4m Figure 8.10 Tking moments bout

More information

99/105 Comparison of OrcaFlex with standard theoretical results

99/105 Comparison of OrcaFlex with standard theoretical results 99/105 Comprison of OrcFlex ith stndrd theoreticl results 1. Introduction A number of stndrd theoreticl results from literture cn be modelled in OrcFlex. Such cses re, by virtue of being theoreticlly solvble,

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

ME311 Machine Design

ME311 Machine Design ME11 Mchine Design Lecture 10: Springs (Chpter 17) W Dornfeld 9Nov018 Firfield University School of Engineering A Free Body Digrm of coil spring (cutting through nywhere on the coil) shows tht there must

More information

IDENTIFICATION AND MODIFICATION OF FRAME STRUCTURE

IDENTIFICATION AND MODIFICATION OF FRAME STRUCTURE The 14 th World Conference on Erthque Engineering IDENTIFICATION AND MODIFICATION OF FRAME STRUCTURE P. Roso 1 1 Ass. Professor, Center of Mechnics nd Structurl Dynmics, Vienn University of Technology,

More information

Evaluation of Allowable Hold Loading of, Hold No. 1 with Cargo Hold No. 1 Flooded, for Existing Bulk Carriers

Evaluation of Allowable Hold Loading of, Hold No. 1 with Cargo Hold No. 1 Flooded, for Existing Bulk Carriers (997) (Rev. 997) (Rev.2 ept. 2000) (Rev.3 July 2004) Evlution of Allowble Hold Loding of Crgo, Hold No. with Crgo Hold No. Flooded, for Existing Bulk Crriers. - Appliction nd definitions These requirements

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

Advanced Computational Analysis

Advanced Computational Analysis Advnced Computtionl Anlysis REPORT REPORT NO: S2149-2 Revision A Title: Closed-Form Anlysis Of Forces And Moments In Bungee Trmpoline Structure Client: Mr Jmes Okey Author: Dr M Lcey BSc PhD CEng F I Mech

More information

ANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS

ANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS ANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS A. N. Anoshkin *, V. Yu. Zuiko, A.V.Glezmn Perm Ntionl Reserch Polytechnic University, 29, Komsomolski Ave., Perm, 614990, Russi

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE Determintion RevAdvMterSci of mechnicl 0(009) -7 properties of nnostructures with complex crystl lttice using DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES

DOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m

More information

Rolling Contact Bearings (pg 599)

Rolling Contact Bearings (pg 599) Bering V9.xmcd [Pg / 6] Title [234] The Units used s stndrd: m, kg, N, P, sec, wtts N, kg, m, P, sec/min, wtts/kw Rolling Contct Berings (pg 599) This note is only guideline for using the text book. Detiled

More information

New Fuzzy Numerical Methods for Solving Cauchy Problems

New Fuzzy Numerical Methods for Solving Cauchy Problems Article New Fuzzy Numericl Metods for Solving Cucy Problems Hussein ALKssbe, *, Irin Perfiliev, Mummd Zini Amd nd Zinor Ridzun Yy Institute of Engineering Mtemtics, Universiti Mlysi Perlis, Kmpus Tetp

More information

Tangent Line and Tangent Plane Approximations of Definite Integral

Tangent Line and Tangent Plane Approximations of Definite Integral Rose-Hulmn Undergrdute Mthemtics Journl Volume 16 Issue 2 Article 8 Tngent Line nd Tngent Plne Approximtions of Definite Integrl Meghn Peer Sginw Vlley Stte University Follow this nd dditionl works t:

More information

APPROXIMATE INTEGRATION

APPROXIMATE INTEGRATION APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be

More information

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer 5 Proceedings of the Interntionl Conference on Informtion nd Automtion, December 58, 5, Colombo, Sri Ln. Appliction of Exct Discretiztion for Logistic Differentil Equtions to the Design of Discrete-ime

More information

13.4 Work done by Constant Forces

13.4 Work done by Constant Forces 13.4 Work done by Constnt Forces We will begin our discussion of the concept of work by nlyzing the motion of n object in one dimension cted on by constnt forces. Let s consider the following exmple: push

More information

An inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation

An inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation Americn Journl of Engineering Reserch (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-02, Issue-10, pp-276-281 www.jer.org Reserch Pper Open Access An inverse stedy stte therml stresses in thin clmped

More information

System Identification with Noisy Data

System Identification with Noisy Data System Identifiction with Noisy Dt U K Dewngn & S V Bri 2 Dept of Civil Engineering, Ntionl Institute of Technology, Ripur, Indi 2 Dept of Civil Engineering, IIT Khrgpur, Khrgpur -72 302, Indi E-mil :

More information

Damage of Houses and Residential Areas by Niigata Prefecture Earthquakes (Part2)

Damage of Houses and Residential Areas by Niigata Prefecture Earthquakes (Part2) Proceedings of the Eighteenth () Interntionl Offshore nd Polr Engineering Conference Vncouver, BC, Cnd, July -11, Copyright by The Interntionl Society of Offshore nd Polr Engineers (ISOPE) ISBN 97-1-53-7-

More information

Stepsize Control for Cubic Spline Interpolation

Stepsize Control for Cubic Spline Interpolation Nonme mnuscript No. (will be inserted by te editor) Stepsize Control for Cubic Spline Interpoltion Dvod Kojste Slkuye Received: dte / Accepted: dte Abstrct Numericl lgoritms wic include stepsize re ffected

More information

Entropy ISSN

Entropy ISSN Entropy 006, 8[], 50-6 50 Entropy ISSN 099-4300 www.mdpi.org/entropy/ ENTROPY GENERATION IN PRESSURE GRADIENT ASSISTED COUETTE FLOW WITH DIFFERENT THERMAL BOUNDARY CONDITIONS Abdul Aziz Deprtment of Mechnicl

More information

Credibility Hypothesis Testing of Fuzzy Triangular Distributions

Credibility Hypothesis Testing of Fuzzy Triangular Distributions 666663 Journl of Uncertin Systems Vol.9, No., pp.6-74, 5 Online t: www.jus.org.uk Credibility Hypothesis Testing of Fuzzy Tringulr Distributions S. Smpth, B. Rmy Received April 3; Revised 4 April 4 Abstrct

More information

ENERGY-BASED METHOD FOR GAS TURBINE ENGINE DISK BURST SPEED CALCULATION

ENERGY-BASED METHOD FOR GAS TURBINE ENGINE DISK BURST SPEED CALCULATION 28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES ENERGY-BASED METHOD FOR GAS TURBINE ENGINE DISK BURST SPEED CALCULATION Anton N. Servetnik Centrl Institute of Avition Motors, Moscow, Russi servetnik@cim.ru

More information

Reliable Optimal Production Control with Cobb-Douglas Model

Reliable Optimal Production Control with Cobb-Douglas Model Relible Computing 4: 63 69, 998. 63 c 998 Kluwer Acdemic Publishers. Printed in the Netherlnds. Relible Optiml Production Control with Cobb-Dougls Model ZHIHUI HUEY HU Texs A&M University, College Sttion,

More information

Electron Correlation Methods

Electron Correlation Methods Electron Correltion Methods HF method: electron-electron interction is replced by n verge interction E HF c E 0 E HF E 0 exct ground stte energy E HF HF energy for given bsis set HF Ec 0 - represents mesure

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient * Interntionl Mthemticl Forum, 4, 9, no., 7-3 Appliction of Exp-Function Method to Huxley Eqution with Vrible Coefficient * Li Yo, Lin Wng nd Xin-Wei Zhou. Deprtment of Mthemtics, Kunming College Kunming,Yunnn,

More information

12 Basic Integration in R

12 Basic Integration in R 14.102, Mt for Economists Fll 2004 Lecture Notes, 10/14/2004 Tese notes re primrily bsed on tose written by Andrei Bremzen for 14.102 in 2002/3, nd by Mrek Pyci for te MIT Mt Cmp in 2003/4. I ve mde only

More information

Module 1. Energy Methods in Structural Analysis

Module 1. Energy Methods in Structural Analysis Module 1 Energy Methods in Structurl Anlysis Lesson 4 Theorem of Lest Work Instructionl Objectives After reding this lesson, the reder will be ble to: 1. Stte nd prove theorem of Lest Work.. Anlyse stticlly

More information

DESIGN OF ANCHORED WALLS: THE INFLUENCE OF DESIGN APPROACHES AND DESIGN METHODS

DESIGN OF ANCHORED WALLS: THE INFLUENCE OF DESIGN APPROACHES AND DESIGN METHODS INTERNATIONAL SCIENTIFIC CONFERENCE CIBv 2010 12 13 November 2010, Brşov DESIGN OF ANCHORED WALLS: THE INFLUENCE OF DESIGN APPROACHES AND DESIGN METHODS Nicolet Mri ILIEŞ *, Trifn Alexndru HULPUȘ *, Augustin

More information

MECHANICS OF MATERIALS

MECHANICS OF MATERIALS 9 The cgrw-hill Compnies, Inc. All rights reserved. Fifth SI Edition CHAPTER 5 ECHANICS OF ATERIALS Ferdinnd P. Beer E. Russell Johnston, Jr. John T. DeWolf Dvid F. zurek Lecture Notes: J. Wlt Oler Texs

More information

Integral equations, eigenvalue, function interpolation

Integral equations, eigenvalue, function interpolation Integrl equtions, eigenvlue, function interpoltion Mrcin Chrząszcz mchrzsz@cernch Monte Crlo methods, 26 My, 2016 1 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion

More information

Comparative Investigation of Trial load and Finite Element Methods in Analysis of Arch Dams

Comparative Investigation of Trial load and Finite Element Methods in Analysis of Arch Dams omprtive Investigtion of ril lod nd Finite Element ethods in Anlysis of Arch Dms Vhid Nourni Abstrct: Becuse of importnt role of dms nd dm construction in humn life, in the present pper the method of nlysis

More information

Construction of Gauss Quadrature Rules

Construction of Gauss Quadrature Rules Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

INVESTIGATION ON THE MODEL OF VORTEX-INDUCED

INVESTIGATION ON THE MODEL OF VORTEX-INDUCED The Seventh Asi-Pcific Conference on Wind Engineering, November 8-1, 9, Tipei, Tiwn ABSTRACT INVESTIGATION ON THE MODEL OF VORTEX-INDUCED VIBRATIONS OF RECTANGULAR SUPER HIGH-RISE BUILDINGS Hi-Yng Wu 1

More information

Module 11.4: nag quad util Numerical Integration Utilities. Contents

Module 11.4: nag quad util Numerical Integration Utilities. Contents Qudrture Module Contents Module 11.4: ng qud util Numericl Integrtion Utilities ng qud util provides utility procedures for computtion involving integrtion of functions, e.g., the computtion of the weights

More information