Application Study on Synthesis Method of Earthquake Motion

  • 格式:pdf
  • 大小:492.18 KB
  • 文档页数:10

Journal ofMechanics Engineering andAutomation 1(201 1)181・190 眷Ll瓣 lN憾 

Application Study on Synthesis Method of Earthquake Motion 

Kahori Iiyama ,Fumio Sasaki ,Masahiko Nakamura ,Akira Tanabe ,Tetsuo Tamaoki ,Wataru Mizuma ̄i and MichioYamada ̄ i.ArchitecturalDepartment,FaculCyofEngineering,Tokyo Universi ofScience.tokyo 102-0073,Japan 2.GeneralProductionCenter,TokyoBuildingDept.,ShimizuCorporation,Tokyo105—8007,Japan 3.Social&EnvironmentalDivision,AtomicEnergySocietyofJapan。Tokyo105-0004。Japan 4.Toshiba Nuclear Engineering Services Corporation,Yokohama 235-8523,Japan 5.Japan NuclearEnergySafetyOrganization, ̄okyo 105—0001,Japan Research InstituteforMathematicalSciences。Kyoto University,Kyoto 606-8502,Japan 

Received:May 17,201 1/Accepted:June 10,201 1/Published:August 25,201 1 Abstract:For seimnicdesignofstructureandmachinery,itisimportantto reproduceinputearthquakemotionsthat arelikelytooccur at a target site.Among the various methods used for generating artificia1 earthquake motions the Synthesis Method of Trigonometric Function is used widely.In this method,artificial waves are reproduced by superimposing sine waves and then adding information about amplitude and phase in the frequency domain.In the Japanese architectural design code,the amplitude is standardized as the target re ̄onse spectrum,andthephasecallbedefinedby randomnumbersorbythephaseofoneobservedwave.However,arandom phaseisdistinctly diferentfromthephase ofanactual earthquake.Further,thephase ofone observedwaveis confinedtothephase characteristic ofthe artificial wave ofonly one specific earthquake motion.In this paper,the authors introduce a new convenient method to reproduce artificia1 waves that not only satisfy the standardized spectrum property but also have the time一 quency characteristics ofmultiple observedwaves.11le authors showthefeature ofthe artificialwaves.discussthemerits ofthemethodby comparingwith existingmethods.andreportthetendencies ofthenon-linear responsebyusing simplemodel_ 

Key words:Artificial earthquake motion.wavelet transform,time-frequency characteristics,non-stationarity 1.IntrOductiOn In recent years,we have experienced many large-scale earthquakes around the country and finally, on March 11 in 2011,the massive earthquake of magnitude of 9.0 occurred off the Pacific coast of the northeasternpart ofthe Japanesemainlandand caused devastating damages in the wide area.Under these circumstances,development of the seismic design for structures and machinery is becoming increasingly important.Practical seismic designs inevitably require 

Fumio Sasaki,professor,Ph.D.,research field:applied mathematics Corresponding author:Kahori Iiyama,research associate, research fields:seismic engineering,seismic performance of buildings.E—mail:kiiyama@rs.kagu tus.ac.jp. 

adequate inputs of earthquake ground motion that mimic possible ground vibrations at all actual site.In order to produce adequate input motions,we have to take into account not only the property of the Fourier amplitude but also that of the Fourier phase that correctly represents the non—stationarity of the frequency.In previous discussions,Tamaoki et a1. proposed a new practical non—parametric method to generate an artificial wave that has the required spectrum property by synthesizing multiple observed waves【1].Recently,Nakamura et a1.showed that all artificial wave includes the phase characteristics of each observed wave,regardless ofthe intensity ofeach wave[2】. 182 Application Study on Synthesis Method of Ea ̄hquake Motion In section 2 ofthis paper,we introduce the algorithm of the method,including a discussion on the general characteristics of the method itself and the artificial waves.In section 3,the merits of the method are demonstrated by comparing with existing methods, focusing especially on the non-stationarity of the frequency.In section 4,a brief comparison of the response performance of the nonlinear structure is reported,and conclusions are summarized in section 5. 

2.Method 2. Algorithm ofMethod The process of this method is shown in Fig.1. Following Nakamura et a1.,the details of each step in Fig.1 are explained as follows【2】: Step 1:Selection of acceleration records First,we select multiple acceleration records from a database.Here,the wave number is expressed by the sub—index f O:l,2 .工),where£is the number of acceleration records needed to reproduce oue artificial wave.In this step,we also adjust the start time of the record.In concrete terms,we set the earthquake occurrence time as the start time,which corresponds to =0.0 s.for each record. Step 2:Wavelet transformation Each of acceleration records (f)is transformed to time-frequency domain by orthonormal wavelet transform and decomposed into particular wavelet angular frequency bands.By using a mother wavelet ( ,orthonormal wavelet transfo咖is defined by i, --E ) , (Oat (1) where{n I 尼∈Z}is the orthonomal wavelet coefficient of the acceleration records ( )and Z is the set of all integers[31.The symbol denotes a complex conjugate. , (f)are so—called‘'wavelets” defined as follows: 哆, (f)=2 (2sf—k) (2) Here,J denotes the scaling parameter;忘shift parameter.Each record (f)can be decomposed into (f)asEq.(3)andthedecomposedwaveseat)to each scaleparanaeterj are expressedbyEq.(4). xi( )=∑ ( ) (3) J =∑o,j, , (f) (4) Step 3:Weight multiplication Multiplying by the weight W (see Eq.(6)),the decomposed time functions of the artificial wave ,( for eachiarcgenerated as ( )= ・ ( ) (5) Here,the weight is defined as =C ・ P j= (_/ ). (6) where is the standard coefficient that indicates the ratio of the average(symbol一)of the target response spectrum ST to that of the response spectrum of each acceleration record sj. w is the pure weight and its initial value is decided at random.In this study,the spectrum intensity is normalized in order to treat all the acceleration records equally,regardless of the magnitude and hypocentral distance. Step 4:Integration By integrating all the members of ( ,the time function of the artificial wave z(t)is calculated as z( )=∑EYS(t) (7) i i Step 5:Comparison of response spectra