JP3488848B2 - Dynamic response analysis method of earthquake motion - Google Patents
Dynamic response analysis method of earthquake motionInfo
- Publication number
- JP3488848B2 JP3488848B2 JP29835399A JP29835399A JP3488848B2 JP 3488848 B2 JP3488848 B2 JP 3488848B2 JP 29835399 A JP29835399 A JP 29835399A JP 29835399 A JP29835399 A JP 29835399A JP 3488848 B2 JP3488848 B2 JP 3488848B2
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- analysis method
- phenomenon
- physical properties
- liquefaction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
Description
【０００１】[0001]
【発明の属する技術分野】本発明は、地震動の動的応答
解析法に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a dynamic response analysis method for earthquake motion.
【０００２】[0002]
【従来の技術】地震動の動的応答解析法としては、振動
領域で解析を行う周波数応答解析法と、時刻歴領域で解
析を行う時刻歴応答解析法とが知られている。前者の解
析法では、解析手法に、複素応答解析法を採用する方法
であり、後者の解析法では、解析手法に、直接積分法を
採用する方法である。2. Description of the Related Art As a dynamic response analysis method for earthquake motion, a frequency response analysis method for performing analysis in a vibration region and a time history response analysis method for performing analysis in a time history region are known. The former analysis method uses a complex response analysis method as the analysis method, and the latter analysis method uses a direct integration method as the analysis method.
【０００３】図１４は、複素応答解析法による動的応答
特性の解析法の手順の一例を示している。同図は、地盤
中に構築されたカルバートに複素応答解析法を適用した
場合であって、左端の上部に示した地震動が入力した場
合に、カルバートにどのような応答振動が発生するのか
を解析している。FIG. 14 shows an example of the procedure of the dynamic response characteristic analysis method by the complex response analysis method. This figure shows a case in which the complex response analysis method is applied to a culvert built in the ground, and what kind of response vibration is generated in the culvert when the earthquake motion shown at the upper left part is input is analyzed. is doing.
【０００４】解析手順は、まず、．複素フーリエ変換
により、入力動Ｆ(ｔ)のフーリエスペクトルＦ(ｆ)を求
める。The analysis procedure is as follows. The Fourier spectrum F (f) of the input motion F (t) is obtained by the complex Fourier transform.
【０００５】一方、解析モデルに、有限要素法によるメ
ッシュを設定し、．解析モデルの剛性などを評価した
運動方程式を解き、解析モデルの各節点に置ける伝達関
数Ａ(ｆ)を求める。On the other hand, a mesh based on the finite element method is set in the analytical model ,. The equation of motion that evaluates the rigidity of the analytical model is solved, and the transfer function A (f) placed at each node of the analytical model is obtained.
【０００６】入力動Ｆ(ｔ)のフーリエスペクトルＦ(ｆ)
と伝達関数Ａ(ｆ)とが求められると、．両者を乗算し
て、応答のフーリエスペクトルＲ(ｆ)を求める。応答の
フーリエスペクトルＲ(ｆ)が求められると、．複素フ
ーリエ逆変換を行い、応答のフーリエスペクトルＲ(ｆ)
から応答波形Ｒ(ｔ)の時刻歴を求める。Fourier spectrum F (f) of the input motion F (t)
And the transfer function A (f) are obtained ,. The two are multiplied to obtain the Fourier spectrum R (f) of the response. When the Fourier spectrum R (f) of the response is obtained ,. Inverse Fourier transform is performed and Fourier spectrum of response R (f)
Then, the time history of the response waveform R (t) is obtained.
【０００７】他方、時刻歴応答解析法においては、微分
方程式の初期値問題をステップバイステップに順次解く
直接積分法が行われ、例えば、運動方程式が、Ｍｕ”＋
Ｃｕ'＋Ｋｕ＝ｆである場合には、
ｕ(ｔ＋Δｔ)＝ｕ(ｔ)＋Δｔｕ'(ｔ)
ｕ'(ｔ＋Δｔ)＝ｕ'(ｔ)＋Δｔｕ"(ｔ)
ｕ"(ｔ＋Δｔ)＝Ｍ-1{ｆ(ｔ＋Δｔ)−Ｋｕ(ｔ＋Δｔ)−
Ｃｕ'(ｔ＋Δｔ)}
という形で計算を進めて、解を求める。On the other hand, in the time history response analysis method, a direct integration method for sequentially solving the initial value problem of the differential equation is performed step by step. For example, the equation of motion is Mu ″ +
When Cu ′ + Ku = f, u (t + Δt) = u (t) + Δtu ′ (t) u ′ (t + Δt) = u ′ (t) + Δtu ″ (t) u ″ (t + Δt) = M −1 {f (t + Δt) -Ku (t + Δt)-
The solution is obtained by advancing the calculation in the form of Cu '(t + Δt)}.
【０００８】ところで、このような２つの動的応答解析
法においては、前者の場合には、解析が無条件で安定に
なるが、後者の場合には、時刻刻みの大きさなどによ
り、解析が不安定になる場合がある。By the way, in the two dynamic response analysis methods as described above, the analysis is unconditionally stable in the former case, but in the latter case, the analysis is performed due to the size of the time step. It may become unstable.
【０００９】また、地震動の動的応答解析では、変化が
緩やかな範囲が多く、前者の解析法では、解析振動数の
定義や、振動数領域での補間手段を採用することによ
り、演算処理時間の低減を図ることができるが、後者の
場合には、解の安定性を考えると、時刻刻みをあまり大
きくすることが難しいので、演算処理に時間がかかると
いう問題があって、地震動の解析には、比較的多く前者
の解析法が採用されている。In the dynamic response analysis of seismic motion, there are many ranges in which the change is gradual. In the former analysis method, the definition of the analysis frequency and the interpolating means in the frequency domain are used to reduce the calculation processing time. However, in the latter case, considering the stability of the solution, it is difficult to make the time step too large, so there is a problem that the calculation process takes time, and it is difficult to analyze the seismic motion. For, the former analysis method is adopted relatively often.
【００１０】しかし、複素応答解析法は、フーリエ変換
という重ね合わせの原理を採用しているので、線形問題
にしか適用することができない。ところが、地盤の歪な
どの解析は、非線形問題であって、複素応答解析法をそ
のまま適用することができない。However, since the complex response analysis method adopts the principle of superposition called Fourier transform, it can be applied only to a linear problem. However, analysis of ground strain is a non-linear problem, and the complex response analysis method cannot be applied as it is.
【００１１】そこで、このような非線形問題に複素応答
解析法を適用するための手法として、線形解析を複数回
繰り返すことにより、非線形問題の解を求める等価線形
解析法が案出され、非線形問題に広く用いられている。Therefore, as a method for applying the complex response analysis method to such a nonlinear problem, an equivalent linear analysis method for solving the nonlinear problem by repeating the linear analysis a plurality of times has been devised. Widely used.
【００１２】図１５は、この等価線形解析法により、非
線形関係にある地盤の歪と剛性との関係を解く場合の解
析フロー図１５(Ａ)と、物性値の歪依存曲線図１５(Ｂ)
とを示している。FIG. 15 is an analysis flow diagram (A) for solving the relation between the strain and the stiffness of the ground having a non-linear relationship by this equivalent linear analysis method, and a strain dependence curve diagram (B) of the physical property value.
Is shown.
【００１３】物性値の歪依存曲線は、実際の地盤の調査
結果に基づいて作成され、同曲線においては、Ｇが地盤
の剛性を、ｈが歪をそれぞれ示している。The strain-dependent curve of the physical property value is prepared based on the actual ground survey result, and in the curve, G represents the rigidity of the ground and h represents the strain.
【００１４】等価線形解析を行う際には、まず、ステッ
プ１で、初期物性値の定義が行われる。図１５(Ｂ)に示
した歪依存曲線においては、初期値は、Ｇ0およびｈ0と
定義される。In carrying out the equivalent linear analysis, first, in step 1, initial physical property values are defined. In the strain dependence curve shown in FIG. 15B, initial values are defined as G 0 and h 0 .
【００１５】続くステップ２では、剛性の評価が行わ
れ、続くステップ３で、周波数応答計算が行われる。こ
の周波数応答計算は、図１４に示した複素応答解析法と
同じ手順により実行される。In the following step 2, rigidity is evaluated, and in the following step 3, frequency response calculation is performed. This frequency response calculation is executed by the same procedure as the complex response analysis method shown in FIG.
【００１６】そして、周波数応答計算の結果が得られる
と、次ぎのステップ４で、有効歪eγの計算が行われ、
続くステップ５で、計算に用いた物性値(歪値)と、ステ
ップ４で計算された物性値(有効歪値)とが適合するか否
かが判断され、適合していれば、解析フローを終了す
る。When the result of the frequency response calculation is obtained, the effective distortion eγ is calculated in the next step 4,
In the following step 5, it is judged whether or not the physical property value (strain value) used for the calculation and the physical property value (effective strain value) calculated in the step 4 are matched, and if they are matched, the analysis flow is executed. finish.
【００１７】一方、ステップ５で適合していないと判断
された場合には、ステップ６に移行し、歪依存曲線によ
り計算された歪に合う物性値を新たに定義して、ステッ
プ２に戻る。On the other hand, if it is determined in step 5 that they do not match, the process proceeds to step 6, where a physical property value that matches the strain calculated by the strain dependence curve is newly defined, and the process returns to step 2.
【００１８】ステップ４から６までをより具体的に説明
すると、今ここで、ステップ４の計算によって得られた
物性値Ｇ0,ｈ0に対する有効歪がeγ1であったとする
と、ステップ５では、初期歪値ｈ0とこの有効歪eγ1と
が適合しているか否かが判断される。The steps 4 to 6 will be described more specifically. Now, assuming that the effective strain for the physical property values G 0 , h 0 obtained by the calculation in step 4 is eγ 1 , then in step 5, It is determined whether the initial strain value h 0 and this effective strain e γ 1 match.
【００１９】この場合には、図１５(Ｂ)に示した歪依存
曲線から見ると明らかなように、適合していない。そこ
で、ステップ６では、有効歪eγ1に対する物性値Ｇ1,ｈ
1を新たに定義して、ステップ２以降の手順を実行し、
このような手順を複数回繰り返すことにより、収束させ
て解を得る。In this case, as is apparent from the strain dependence curve shown in FIG. 15 (B), it does not match. Therefore, in step 6, the physical property values G 1 , h for the effective strain e γ 1
1 is newly defined, the procedure from step 2 is executed,
By repeating such a procedure a plurality of times, the solution is converged and the solution is obtained.
【００２０】そして、このような解析フローを全ての点
において実行することにより、地盤に対する応力解析が
行われるが、このような従来の等価線形解析法には、以
下に説明する技術的な課題があった。The stress analysis on the ground is performed by executing such an analysis flow at all points. However, such a conventional equivalent linear analysis method has the following technical problems. there were.
【００２１】[0021]
【発明が解決しようとする課題】すなわち、従来の等価
線形解析法は、各要素に対して、動的変形特性から１つ
の代表的な物性値を仮定する線形解析であるので、ある
程度の歪レベルに対しては、信頼性が高いものの、地盤
の液状化現象のように、経時的な関係において、物性が
急激に変化する現象においては、１つの代表的な物性値
でこれを表現することが難しいので、解析結果の信頼性
に欠けると言う問題があった。That is, since the conventional equivalent linear analysis method is a linear analysis in which one representative physical property value is assumed from the dynamic deformation characteristics for each element, the strain level to a certain extent is obtained. In contrast, although it is highly reliable, one representative physical property value can be used to represent a phenomenon in which the physical properties change rapidly over time, such as the liquefaction phenomenon of the ground. Since it is difficult, there is a problem that the reliability of the analysis result is lacking.
【００２２】より具体的に説明すると、図１６は、液状
化が発生した前後の地盤の歪依存特性を示しており、同
図には、実線で液状化前の剛性低下の状態が、また、点
線で液状化後の剛性低下の状態がそれぞれ示されてい
る。More specifically, FIG. 16 shows the strain-dependent characteristics of the ground before and after liquefaction. In the same figure, the solid line shows the state of decreased rigidity before liquefaction. The dotted lines indicate the state of decreased rigidity after liquefaction.
【００２３】このような歪依存特性を有する地盤におい
て従来の等価線形解析法により解析を行う際には、物性
値の選択として、実線か点線のいずれか１つしか選べな
いので、地盤の液状化を評価することが困難で、地盤の
液状化現象などのように、急激に物性が変化する場合つ
いて、等価線形解析法では、厳密に評価できないという
考え方は、当業者の常識的な考え方であった。When performing analysis by the conventional equivalent linear analysis method on the ground having such strain-dependent characteristics, since only one of the solid line and the dotted line can be selected as the physical property value, the ground is liquefied. It is difficult to evaluate, and the idea that the equivalent linear analysis method cannot be used for rigorous evaluation when the physical properties change rapidly, such as the liquefaction phenomenon of the ground, is a common sense of those skilled in the art. It was
【００２４】このような問題の解決に際しては、等価線
形解析法に時間的な概念を導入することも考えられる。In solving such a problem, it is possible to introduce a temporal concept into the equivalent linear analysis method.
【００２５】しかしながら、等価線形解析法では、複素
応答解析法を基礎としており、振動数領域での動的応答
解析を行い、振動数領域での応答解析を考えた場合、当
然のことながら、入力地震動の全てを対象としてフーリ
エ変換を行う。However, the equivalent linear analysis method is based on the complex response analysis method, and when the dynamic response analysis in the frequency domain is performed and the response analysis in the frequency domain is considered, the input Fourier transform is applied to all seismic motions.
【００２６】 ところが、フーリエ変換された入力地震
動は、この時点で、既に時間的概念がなくなっており、
解析自体は、０秒から継続時間までの地震が入力された
場合の応答解析となる。[0026] However, the input ground motion that has been Fourier transform, at this point, has become not already time concept,
The analysis itself is a response analysis when an earthquake from 0 seconds to the duration is input.
【００２７】従って、従来の等価線形解析法では、時間
の概念が全く存在しない状態で解析を行っており、解析
者は、時間的な概念を全く思考の外においていて、入力
地震動をある時間で分割するといった発想をする必然性
がなく、従来の等価線形解析法から、物性の急激変化時
の前後で時間的に分割するという技術思想を簡単に導き
出すことはできなかった。Therefore, in the conventional equivalent linear analysis method, the analysis is performed in the state where the concept of time does not exist at all, and the analyst does not think of the time concept at all, and inputs the seismic motion at a certain time. There is no inevitability to take the idea of division, and it was not possible to easily derive the technical idea of temporal division before and after the sudden change of physical properties from the conventional equivalent linear analysis method.
【００２８】本発明は、このような従来の問題点に鑑み
てなされたものであって、その目的とするところは、物
性が急激に変化いる現象において、信頼性の高い解析が
可能な地震動の動的応答解析法を提供することにある。The present invention has been made in view of the above conventional problems, and an object of the present invention is to obtain a seismic motion capable of highly reliable analysis in a phenomenon in which physical properties are rapidly changed. It is to provide a dynamic response analysis method.
【００２９】[0029]
【課題を解決するための手段】上記目的を達成するた
め、本発明は、線形解析を複数回繰り返すことにより、
非線形問題に適用される等価線形解析法に基づく地震動
の動的応答解析法において、入力地震動を経時的に複数
に分割し、分割した部分にそれぞれ前記等価線形解析法
による解析を実施し、しかる後に、各解析によって得ら
れた結果を重ね合わせる動的応答解析法であって、前記
動的応答解析法を経時的に物性が急激に変化する現象に
適用する際に、前記経時的な分割個所を前記物性の急激
変化点に設定するとともに、分割した部分にそれぞれ実
施される前記等価線形解析法による解析の物性を別にす
るようにした。このように構成した地震動の動的応答解
析法によれば、入力地震動を経時的に複数に分割し、分
割した部分にそれぞれ等価線形解析法による解析を実施
し、しかる後に、各解析によって得られた結果を重ね合
わせる。このことは、分割した部分の前後で、それぞれ
独立して等価線形解析を行うことを意味し、それぞれの
解析で使用する物性に関しては、別の物性を用いる。従
って、分割した部分にそれぞれ別の物性を用いることが
可能になり、従来の等価線形解析法の課題を解決して、
信頼性の高い解析結果が得られ、例えば、経時的に物性
が急激に変化する地盤の液状化に適用すると、それを正
当に評価することができる。この場合、本発明では、物
性が急激に変化する現象として、地盤の液状化現象以外
に、地盤に破砕帯が存在する場合の滑り現象、構造物の
浮き上がり現象,構造物の塑性化現象から選択すること
ができる。また、本発明の解析法は、前記物性が急激に
変化する現象が、地盤の液状化現象であって、前記物性
の急激変化時の前後で複数に分割し、分割した部分にそ
れぞれ前記等価線形解析法による解析を実施する前に、
分割しない状態で等価線形解析法による予備解析を実施
して、この予備解析に基づいて、前記地盤の液状化時期
を求めることができる。この構成を採用すると、地盤の
液状化時期の判断と動的応答解析とを一連の手順で処理
することができる。以上のような本発明の技術的思想の
根幹は、経時的に物性が急激に変化する現象に、等価線
形解析法を適用するに際して、物性の急激変化時の前後
で時間的に分割して、分割した部分にそれぞれ等価線形
解析法を適用し、その後に解析結果を重ね合わせるの
で、複数の等価線形解析法を重合する多重等価線形解析
法であるともいる。In order to achieve the above-mentioned object, the present invention provides linear analysis by repeating a plurality of times.
In the dynamic response analysis method of the earthquake motion based on the equivalent linear analysis method applied to the non-linear problem, the input earthquake motion is divided into multiple over time, and the divided parts are analyzed by the equivalent linear analysis method, respectively. , A dynamic response analysis method for superposing the results obtained by each analysis ,
Dynamic response analysis method is used for phenomena in which physical properties change rapidly with time.
When applying it, make sure that the physical properties of the
In addition to setting the change point,
Apart from the physical properties of the analysis by the equivalent linear analysis method
It was to so. According to the dynamic response analysis method for seismic motions configured in this way, the input seismic motion is divided into a plurality of parts over time, and the divided parts are each analyzed by the equivalent linear analysis method. Stack the results. This means that the equivalent linear analysis is performed independently before and after the divided parts , and different physical properties are used for the physical properties used in each analysis . Therefore, it becomes possible to use different physical properties for the divided parts, solving the problems of the conventional equivalent linear analysis method,
A highly reliable analysis result can be obtained, and when it is applied to, for example, liquefaction of ground whose physical properties change rapidly with time, it can be properly evaluated. In this case, in the present invention, as a phenomenon in which the physical properties change abruptly, in addition to the liquefaction phenomenon of the ground, a sliding phenomenon when a crush zone exists in the ground, a lifting phenomenon of the structure, a plasticization phenomenon of the structure is selected. can do. In the analysis method of the present invention, the phenomenon in which the physical properties change rapidly is a liquefaction phenomenon of the ground, and is divided into a plurality of parts before and after the rapid change in the physical properties, and each of the divided parts has the equivalent linear shape. Before carrying out the analysis by the analysis method,
It is possible to carry out a preliminary analysis by the equivalent linear analysis method in a state without division, and to determine the liquefaction time of the ground based on this preliminary analysis. If this configuration is adopted, the determination of the liquefaction time of the ground and the dynamic response analysis can be processed by a series of procedures. The basis of the technical idea of the present invention as described above, the phenomenon in which the physical properties change rapidly over time, when applying the equivalent linear analysis method, divided in time before and after the rapid change in physical properties, The equivalent linear analysis method is applied to each of the divided parts, and the analysis results are superposed after that, so it is also said to be a multiple equivalent linear analysis method that superimposes a plurality of equivalent linear analysis methods.
【００３０】[0030]
【発明の実施の形態】以下、本発明の好ましい形態につ
いて添付図面を参照して詳細に説明する。図1および図
２は、本発明にかかる地震動の動的応答解析法の第1実
施例を示している。Preferred embodiments of the present invention will be described in detail below with reference to the accompanying drawings. 1 and 2 show a first embodiment of a dynamic response analysis method for earthquake motion according to the present invention.
【００３１】これらの図に示した地震動の動的応答解析
法は、経時的に物性が急激に変化する現象として、地盤
の液状化現象に本発明を適用した場合を示している。図
１は、本実施例の解析法の解析手順のフロー図である。The dynamic response analysis method of seismic motions shown in these figures shows the case where the present invention is applied to the liquefaction phenomenon of the ground as a phenomenon in which the physical properties change rapidly with time. FIG. 1 is a flow chart of the analysis procedure of the analysis method of this embodiment.
【００３２】解析手順がスタートすると、まず、ステッ
プ１０で地盤の予備解析が行われる。この予備解析は、
液状化の判定と、液状化発生時刻の推定とを行うため
に、実施される手順であって、通常の等価線形解析法に
基づいて、全応力解析が実行される。When the analysis procedure starts, first, in step 10, preliminary analysis of the ground is performed. This preliminary analysis
A total stress analysis is executed based on a normal equivalent linear analysis method, which is a procedure that is performed to determine liquefaction and estimate the liquefaction occurrence time.
【００３３】なお、この全応力解析は、液状化前の動的
変形特性を用い、図１４,１５に示した手順と同じ方法
で行われる。この予備解析が終了すると、ステップ１１
で、液状化判定および液状化発生時刻の推定が行われ
る。This total stress analysis is performed by the same method as the procedure shown in FIGS. 14 and 15 using the dynamic deformation characteristics before liquefaction. When this preliminary analysis is completed, step 11
Then, the liquefaction determination and the liquefaction occurrence time are estimated.
【００３４】この判定および推定では、水位,繰り返し
三軸試験による液状化強度曲線,判定領域を入力し、こ
れらのデータとステップ１０で実施した全応力解析によ
るせん断応力時刻歴とを用い、累積損傷度の概念を適用
することで行われる。In this judgment and estimation, the water level, the liquefaction strength curve by the repeated triaxial test, and the judgment area are input, and by using these data and the shear stress time history by the total stress analysis carried out in step 10, the cumulative damage is calculated. It is done by applying the concept of degrees.
【００３５】続くステップ１２では、ステップ１１で推
定した液状化発生時刻ｔの前後で、入力地震動の分割が
行われ、液状化前の第１グループと、液状化後の第２グ
ループとにグループ分けされる。In the following step 12, the input seismic motion is divided before and after the liquefaction occurrence time t estimated in step 11, and is divided into a first group before liquefaction and a second group after liquefaction. To be done.
【００３６】次ぎに、液状化前の第１グループは、ステ
ップ１３で等価線形解析が行われ、この等価線形解析
は、ステップ１０の予備解析と同様に、液状化前の動的
変形特性を用いた通常の等価線形解析となる。Next, the first group before liquefaction is subjected to equivalent linear analysis in step 13, and this equivalent linear analysis uses the dynamic deformation characteristics before liquefaction as in the preliminary analysis in step 10. It becomes the usual equivalent linear analysis.
【００３７】ただし、ここで行われる等価線形解析で
は、入力地震動は、予備解析と異なり、０秒から液状化
発生時刻ｔまでとされる。However, in the equivalent linear analysis performed here, the input seismic motion is from 0 second to the liquefaction occurrence time t, unlike the preliminary analysis.
【００３８】一方、液状化後の第２グループは、ステッ
プ１４で等価線形解析が行われ、この等価線形解析で
は、液状化後の動的変形特性を用い、入力地震動が、液
状化発生時刻ｔから最後までとされる。On the other hand, the second group after liquefaction is subjected to the equivalent linear analysis in step 14, and in this equivalent linear analysis, the dynamic deformation characteristics after the liquefaction are used, and the input seismic motion is liquefaction occurrence time t. From the end to the end.
【００３９】このようにして、第１および第２グループ
のそれぞれに異なった条件の動的変形特性で等価線形解
析を行い、その解が得られると、ステップ１５で解析結
果(応答波形)を重ね合わせて手順が終了する。In this way, the equivalent linear analysis is performed on the first and second groups with different dynamic deformation characteristics under different conditions, and when the solution is obtained, the analysis results (response waveforms) are overlaid in step 15. Then, the procedure ends.
【００４０】図２は、本発明の理解を容易にするため
に、所要部分を模式的に示したものであって、図１に示
した手順中のステツプ１２からステップ１５に対応して
いる。図２において、最上段に示した図形が入力地震動
の波形であり、図１４のＦ(ｔ)に相当する。FIG. 2 is a schematic view of a required portion for facilitating the understanding of the present invention, and corresponds to steps 12 to 15 in the procedure shown in FIG. In FIG. 2, the figure shown at the top is the waveform of the input earthquake motion, and corresponds to F (t) in FIG.
【００４１】この入力地震動の波形は、本実施例の場合
には、ステップ１１で推定された液状化発生時刻ｔに基
づいて、続くステップ１２で、第１および第２グループ
に分けられる。In the case of this embodiment, the waveform of the input seismic motion is divided into the first and second groups in the following step 12 based on the liquefaction occurrence time t estimated in step 11.
【００４２】このようにして、分割された２つのグルー
プが図２の上から２段目に示されている。そして、本実
施例の場合には、ステップ１３および１４で分割された
グループ毎に等価線形解析が実行される。The two groups thus divided are shown in the second row from the top of FIG. Then, in the case of the present embodiment, the equivalent linear analysis is executed for each of the groups divided in steps 13 and 14.
【００４３】この時の解析の条件が図２の四角内に示さ
れていて、この四角の下方に、第１および第２グループ
のそれぞれに異なった条件で等価線形解析を行つて得ら
れた応答結果がそれぞれ示されている。The conditions of the analysis at this time are shown in the square of FIG. 2, and the response obtained by performing the equivalent linear analysis under the different conditions of the first and second groups below the square. The results are shown respectively.
【００４４】そして、最後に、得られた解析結果(応答
波形)を重ね合わせて手順が終了することになる。Finally, the obtained analysis results (response waveforms) are overlaid to complete the procedure.
【００４５】さて、以上のように構成した地震動の動的
応答解析法によれば、物性の急激変化時の前後で複数に
分割し、分割した部分にそれぞれ等価線形解析法による
解析を実施し、しかる後に、各解析によって得られた結
果を重ね合わせるので、以下の作用効果が得られる。According to the dynamic response analysis method for seismic motions configured as described above, a plurality of parts are divided before and after a sudden change in physical properties, and the divided parts are each analyzed by the equivalent linear analysis method. After that, the results obtained by each analysis are overlaid, so that the following operational effects can be obtained.
【００４６】．解析上、物性の急激変化時を境にし
て、全く動的変形特性が異なる場合を考慮することがで
きる。
．分割に伴なうグループ化により、地盤の時間的な相
互作用を評価することができる。
．地震力の時間的な変化に対応する物性の変化を表現
することができる。.. From the viewpoint of analysis, it is possible to consider the case where the dynamic deformation characteristics are completely different from each other when the physical property changes rapidly. ． The grouping that accompanies the division makes it possible to evaluate the temporal interaction of the ground. ． It is possible to express changes in physical properties corresponding to changes in seismic force over time.
【００４７】図３〜図１０は、本発明の作用効果を確認
するため、実際に発生した地震観測記録を基に行ったシ
ミュレーション解析の条件および結果を示している。こ
のシミュレーション解析では、１９８７年１１月に米国
カリフォルニア州で得られた鉛直アレー観測記録(図３,
４)を用いた。FIGS. 3 to 10 show conditions and results of a simulation analysis performed on the basis of an earthquake observation record that actually occurred in order to confirm the action and effect of the present invention. In this simulation analysis, vertical array observation records obtained in November 1987 in California, USA (Fig. 3,
4) was used.
【００４８】図３に示した地震１は、マグニチュードが
５．８で、最大加速度１２６ｇａｌを記録したが、地盤
の液状化は、発生しなかった。図４に示した地震２は、
マグニチュードが６．１で、最大加速度２０１ｇａｌを
記録し、図５に示した過剰間隙水圧の記録から地盤は、
ほぼ液状化したものと思われる。In the earthquake 1 shown in FIG. 3, the magnitude was 5.8 and the maximum acceleration was 126 gal, but the liquefaction of the ground did not occur. Earthquake 2 shown in Figure 4,
With a magnitude of 6.1, a maximum acceleration of 201 gal was recorded, and from the record of excess pore water pressure shown in FIG.
It seems to be almost liquefied.
【００４９】そこで、地震１を用いた同定解析を実施
し、地盤物性の妥当性や不明確である地盤の歪依存性を
算定し、これによって得られた地盤の特性を踏まえて、
地震２に基づく、本発明の解析法を実施した。Therefore, an identification analysis using Earthquake 1 was carried out to calculate the appropriateness of the ground physical properties and the strain dependency of the ground, which is unclear. Based on the characteristics of the ground obtained by this,
The analysis method of the present invention based on earthquake 2 was performed.
【００５０】全ての解析、ＧＬ−７．５ｍで得られた観
測記録を入力とし、地表の実記録と解析結果とを比較す
ることで、解析の適否を判断した。表１に地盤構造およ
びその物性諸元を示している。The suitability of the analysis was judged by comparing all the analyzes and the observation records obtained at GL-7.5 m with the actual records on the ground surface and the analysis results. Table 1 shows the ground structure and its physical properties.
【００５１】[0051]
【表１】 [Table 1]
【００５２】地盤物性の妥当性や不明確である地盤の動
的変形特性を推定する目的で実施した地震１に基づく同
定解析の結果を図６に示しているが、実記録と解析結果
とは、概ね同様に傾向を示していて、同定解析の妥当性
が認められる。FIG. 6 shows the result of the identification analysis based on earthquake 1 which was carried out for the purpose of estimating the dynamic deformation characteristics of the ground where the physical properties of the ground are unclear and unclear. The actual record and the analysis result are shown in FIG. , And the tendency is almost the same, and the validity of the identification analysis is confirmed.
【００５３】図７は、この同定解析で使用した歪依存曲
線を示している。この同定解析によって得られた解析結
果を踏まえて、地震２に関して、本発明の解析法を実施
した。本発明の解析法では、地盤の液状化の発生時刻で
入力地震動を分割する。そこで、まず、以下の表２に示
すように、各地層毎の液状化時刻を判定した。FIG. 7 shows the strain dependence curve used in this identification analysis. Based on the analysis results obtained by this identification analysis, the analysis method of the present invention was carried out for earthquake 2. In the analysis method of the present invention, the input seismic motion is divided according to the time of occurrence of ground liquefaction. Therefore, first, as shown in Table 2 below, the liquefaction time for each layer was determined.
【００５４】[0054]
【表２】 [Table 2]
【００５５】この判定では、各層で若干異なるものの、
約１３．６ｓｅｃで液状化が発生する結果となった。図
５に示した過剰間隙水圧の記録では、約１３．０ｓｅｃ
から過剰間隙水圧の上昇が認められるので、この判定の
妥当性が認められる。In this judgment, although each layer is slightly different,
As a result, liquefaction occurred in about 13.6 seconds. In the record of excess pore water pressure shown in FIG. 5, it is about 13.0 sec.
Since the increase in excess pore water pressure is observed, the validity of this judgment is confirmed.
【００５６】これらの結果から、本シミュレーションで
は、１３．６ｓｅｃを液状化発生時刻として、この時刻
で２分割し、その前後で等価線形解析を行った。From these results, in this simulation, 13.6 sec was set as the liquefaction occurrence time, and this time was divided into two, and the equivalent linear analysis was performed before and after that.
【００５７】図８に、ＧＬ−７．５ｍに対する地表の伝
達関数、図９に、地表の加速度スペクトル、図１０に、
地表の加速度時刻歴をそれぞれ示している。これらの図
中には、比較のために、本発明の解析法とともに、従来
の等価線形解析法で実施した解析結果も合わせて表示し
ている。FIG. 8 shows the ground surface transfer function for GL-7.5 m, FIG. 9 shows the ground surface acceleration spectrum, and FIG.
The acceleration time history of the ground surface is shown. In these figures, for comparison, the analysis results of the conventional equivalent linear analysis method are also shown together with the analysis method of the present invention.
【００５８】これらの結果から明らかなように、従来の
等価線形解析法の解析結果は、実記録と比較して差異が
顕著に認められるが、本発明の解析法では、概ね実記録
と同じ解析結果が得られ、本発明の有効性が確認され
た。As is clear from these results, the analysis results of the conventional equivalent linear analysis method show significant differences as compared with the actual recording, but in the analysis method of the present invention, the analysis is almost the same as the actual recording. The results were obtained and the effectiveness of the present invention was confirmed.
【００５９】図１１は、本発明にかかる地震動の動的応
答解析法の第２実施例を示しており、以下にその特徴点
についてのみ説明する。同図に示した実施例では、本発
明にかかる解析法を地盤に破砕帯がある場合に適用した
場合を示している。FIG. 11 shows a second embodiment of the dynamic response analysis method for earthquake motion according to the present invention, and only the characteristic points will be described below. The embodiment shown in the figure shows a case where the analysis method according to the present invention is applied to the case where the ground has a shatter zone.
【００６０】地盤中に破砕帯があると、地震動が地盤に
加わった際に、破砕帯にある地震応力が生じると、抵抗
がなくなり、すべり現象が発生する。このような地盤の
解析では、時刻歴応答解析法では、逐次解析(ステップ
バイステップ解析)なので、破砕帯の地震応力の増加関
係を時々刻々と把握し、破砕帯の許容応力以上になる
と、抵抗をなくす要素、すなわち、ジョイント要素を用
いることで評価している。If there is a crush zone in the ground, when the earthquake motion is applied to the ground, if the seismic stress in the crush zone occurs, the resistance disappears and the slip phenomenon occurs. In such ground analysis, the time history response analysis method is a sequential analysis (step-by-step analysis), so the increase relationship of seismic stress in the fracture zone is grasped momentarily, and when the stress exceeds the allowable stress in the fracture zone, the resistance is increased. It is evaluated by using the element that eliminates, that is, the joint element.
【００６１】一方、複素応答解析を基本とする従来の等
価線形解析法では、時間的概念がないことから、破砕帯
の地震応力の増加関係を時々刻々に把握することができ
ないので、このようなモデルの解析が困難とされてい
た。On the other hand, in the conventional equivalent linear analysis method based on the complex response analysis, since there is no concept of time, it is not possible to grasp the increasing relation of the seismic stress in the fracture zone moment by moment. It was considered difficult to analyze the model.
【００６２】ところが、本発明にかかる地震動の動的応
答解析法では、経時的に物性が急激に変化する現象に適
用する際に、物性の急激変化時の前後で複数に分割し、
分割した部分にそれぞれ等価線形解析法による解析を実
施し、しかる後に、各解析によって得られた結果を重ね
合わせる。However, in the dynamic response analysis method for earthquake motions according to the present invention, when applied to a phenomenon in which physical properties change rapidly with time, it is divided into a plurality of parts before and after the rapid change in physical properties,
Each of the divided parts is analyzed by the equivalent linear analysis method, and then the results obtained by each analysis are overlaid.
【００６３】従って、破砕帯の地震応力に基づいて、こ
れが許容応力以内か、以上かを判定し、許容応力以上に
なった時刻で、入力地震動を分割し、許容応力以内の場
合には、通常の剛性を破砕帯に持たせ、許容応力以上の
場合には、無抵抗として解析すると、時刻歴応答解析方
と同様に評価することができる。Therefore, based on the seismic stress of the crush zone, it is judged whether the stress is within the allowable stress or not. At the time when the stress exceeds the allowable stress, the input seismic motion is divided. When the fracture zone is provided with the rigidity of and the stress is equal to or more than the allowable stress, it can be evaluated in the same manner as the time history response analysis method by analyzing as no resistance.
【００６４】この場合、本実施例の解析法では、時刻歴
応答解析のように解析が不安定になることがなく、ま
た、解析時間も大幅に短縮することができる。In this case, according to the analysis method of this embodiment, the analysis does not become unstable unlike the time history response analysis, and the analysis time can be greatly shortened.
【００６５】図１２は、本発明にかかる地震動の動的応
答解析法の第３実施例を示しており、以下にその特徴点
についてのみ説明する。同図に示した実施例では、構造
物の浮き上がり現象に本発明を適用した場合である。FIG. 12 shows a third embodiment of the dynamic response analysis method for seismic motion according to the present invention. Only the characteristic points will be described below. The embodiment shown in the figure is a case where the present invention is applied to the floating phenomenon of a structure.
【００６６】ビルなどの建築構造物が直接基礎を採用し
ていると、地震動が加わった際に、地震力によりロッキ
ング挙動を起こし、構造物の一部が地盤から離れる浮き
上がり現象が発生することがある。When a building structure such as a building directly employs a foundation, when an earthquake motion is applied, rocking behavior is caused by seismic force, and a part of the structure is lifted from the ground. is there.
【００６７】構造物が地盤から離れると、その前後で地
盤の物性が急激に変化し、このような浮き上がり現象
は、時刻歴応答解析では、第２実施例と同様に評価でき
るが、従来の等価線形解析法では、評価できないとされ
ていた。When the structure separates from the ground, the physical properties of the ground change abruptly before and after the structure, and such a lifting phenomenon can be evaluated in the time history response analysis in the same manner as in the second embodiment. It was said that the linear analysis method could not be evaluated.
【００６８】ところが、構造物の浮き上がりを判定し、
構造物に浮き上がり現象が発生した時刻の前後で、入力
地震動の分割を行い、分割した部分にそれぞれ等価線形
解析法による解析を実施し、しかる後に、各解析によっ
て得られた結果を重ね合わせることにより、浮き上がり
現象を評価することが可能になり、本実施例の場合に
も、上記実施例と同等の作用効果が得られる。However, it is judged whether the structure is floating,
By dividing the input seismic motion around the time when the floating phenomenon occurred in the structure, performing the analysis by the equivalent linear analysis method on each of the divided parts, and then superimposing the results obtained by each analysis. As a result, it becomes possible to evaluate the floating phenomenon, and in the case of this embodiment, the same operational effect as the above embodiment can be obtained.
【００６９】図１３は、本発明にかかる地震動の動的応
答解析法の第４実施例を示しており、以下にその特徴点
についてのみ説明する。同図に示した実施例では、構造
物の塑性化現象に本発明を適用した場合である。FIG. 13 shows a fourth embodiment of the dynamic response analysis method for seismic motion according to the present invention, and only the characteristic points will be described below. In the embodiment shown in the figure, the present invention is applied to the plasticization phenomenon of a structure.
【００７０】地中に構築されているボックスカルバート
などの構造物は、使用している構造材料の特性から、大
地震の発生時に、塑性化現象を起こすことが知られてい
る。例えば、鉄筋コンクリート構造物においては、ひび
割れ、降伏状態がそれである。It is known that a structure such as a box culvert constructed in the ground causes a plasticizing phenomenon when a large earthquake occurs due to the characteristics of the structural materials used. For example, in a reinforced concrete structure, it is a cracked or yielded state.
【００７１】このように塑性化現象が発生すると、構造
物には、図１３(Ｂ)に示したような塑性ヒンジが発生す
る。When the plasticizing phenomenon occurs in this way, a plastic hinge as shown in FIG. 13B occurs in the structure.
【００７２】 この場合、時刻歴応答解析では、この現
象を評価することができるが、従来の等価線形解析法で
は、時間的な概念がないので評価することができない。In this case, this phenomenon can be evaluated by the time history response analysis, but cannot be evaluated by the conventional equivalent linear analysis method because there is no concept of time.
【００７３】ところが、本実施例の場合には、塑性ヒン
ジの発生時刻を判定し、この発生時刻において、入力地
震動を分割して、それぞれに等価線形解析を行い、得ら
れた解析結果を重ね合わせることで、その評価が可能に
なる。However, in the case of this embodiment, the generation time of the plastic hinge is determined, the input seismic motion is divided at this generation time, the equivalent linear analysis is performed for each, and the obtained analysis results are superposed. This allows the evaluation.
【００７４】[0074]
【発明の効果】以上実施例で詳細に説明したように、本
発明にかかる地震動の動的応答解析法によれば、従来等
価線形解析法での解析が不可能とされていた物性が急激
に変化する現象などの解析が可能になり、しかも、解析
の信頼性も高められる。As described in detail in the above embodiments, according to the dynamic response analysis method of the seismic motion according to the present invention, the physical properties which were previously impossible to be analyzed by the equivalent linear analysis method are drastically changed. Analysis of changing phenomena is possible, and the reliability of analysis is improved.
【図１】本発明にかかる地震動の動的応答解析法の第１
実施例を示す解析手順のフロー図である。FIG. 1 is a first method of analyzing a dynamic response of earthquake motion according to the present invention.
It is a flowchart of the analysis procedure which shows an Example.
【図２】図１のフロー図の要部を模式的示した説明図で
ある。FIG. 2 is an explanatory view schematically showing a main part of the flow chart of FIG.
【図３】図１に示した解析法のシミュレーションに用い
た地震の観測記録図である。3 is an observation record diagram of an earthquake used in the simulation of the analysis method shown in FIG.
【図４】図１に示した解析法のシミュレーションに用い
た地震の観測記録図である。FIG. 4 is an observation record diagram of an earthquake used in the simulation of the analysis method shown in FIG.
【図５】図４の地震の際の過剰間隙水圧の記録図であ
る。5 is a record diagram of excess pore water pressure during the earthquake of FIG.
【図６】図３に示した地震を入力とする等価線形解析法
による解析結果(加速度応答スペクトル)と、観測記録と
を示すグラフである。6 is a graph showing an analysis result (acceleration response spectrum) by an equivalent linear analysis method using the earthquake shown in FIG. 3 as an input, and an observation record.
【図７】図３に示した地震を入力とする等価線形解析法
で用いた歪依存曲線である。FIG. 7 is a strain-dependent curve used in the equivalent linear analysis method using the earthquake shown in FIG. 3 as an input.
【図８】図４に示した地震を入力とする本発明にかかる
解析法による解析結果(伝達関数)と、観測記録とを示す
グラフである。FIG. 8 is a graph showing an analysis result (transfer function) by the analysis method according to the present invention, which uses the earthquake shown in FIG. 4 as an input, and an observation record.
【図９】図４に示した地震を入力とする本発明にかかる
解析法による解析結果(加速度応答スペクトル)と、観測
記録とを示すグラフである。FIG. 9 is a graph showing an analysis result (acceleration response spectrum) by an analysis method according to the present invention using the earthquake shown in FIG. 4 as an input, and an observation record.
【図１０】図４に示した地震を入力とする本発明にかか
る解析法による解析結果(地表の加速度時刻歴)と、観測
記録とを示す波形図である。10 is a waveform chart showing an analysis result (acceleration time history of the ground surface) and an observation record by the analysis method according to the present invention in which the earthquake shown in FIG. 4 is input.
【図１１】本発明にかかる地震動の動的応答解析法の第
２実施例を示す解析モデルの説明図である。FIG. 11 is an explanatory diagram of an analysis model showing a second embodiment of the dynamic response analysis method for earthquake motion according to the present invention.
【図１２】本発明にかかる地震動の動的応答解析法の第
３実施例を示す解析モデルの説明図である。FIG. 12 is an explanatory diagram of an analysis model showing a third embodiment of the dynamic response analysis method for earthquake motion according to the present invention.
【図１３】本発明にかかる地震動の動的応答解析法の第
４実施例を示す解析モデルの説明図である。FIG. 13 is an explanatory diagram of an analysis model showing a fourth embodiment of the dynamic response analysis method for earthquake motion according to the present invention.
【図１４】複素応答解析法の解析手順の説明図である。FIG. 14 is an explanatory diagram of an analysis procedure of a complex response analysis method.
【図１５】等価線形解析の解析手順の説明図である。FIG. 15 is an explanatory diagram of an analysis procedure of equivalent linear analysis.
【図１６】地盤の液状化前後の歪依存特性を示すグラフ
である。FIG. 16 is a graph showing strain-dependent characteristics before and after liquefaction of ground.
───────────────────────────────────────────────────── フロントページの続き (72)発明者 田中 典明 神奈川県横浜市鶴見区江ヶ崎町４番１号 東京電力株式会社 電力技術研究所内 (72)発明者 佐藤 博 神奈川県横浜市鶴見区江ヶ崎町４番１号 東京電力株式会社 電力技術研究所内 (56)参考文献 特開 平５−157655（ＪＰ，Ａ) 特開 平６−330518（ＪＰ，Ａ) 特開 平10−90051（ＪＰ，Ａ) 特開 平10−307859（ＪＰ，Ａ) 特開 平11−84017（ＪＰ，Ａ) 特開 平11−160144（ＪＰ，Ａ) 特開2000−257054（ＪＰ，Ａ) 柳下文雄、土方勝一郎、西村功、久保 賀也、成川匡文、大島豊、内川祐一郎、 塩見忠彦，ＳＲモデルによる杭基礎建築 物の地震応答解析における液状化の考慮 方法−その９ 自由地盤の有効応力解析 に基，1999年度日本建築学会大会（中 国） 学術講演梗概集 Ｂ−２，1999年 ７月30日，ｐｐ．377−378 塩見忠彦、貫井泰、柳下文雄、土方勝 一郎、横山治男、小山桂介，累積損傷度 を考慮した簡易液状化解析手法−その１ 簡易液状化解析手法の提案，1999年度 日本建築学会大会（中国） 学術講演梗 概集 Ｂ−２，1999年 ７月30日，ｐ ｐ．379−380 吉田望，地震応答解析に用いる地盤物 性をどう評価するか，建築基礎の設計施 工に関する研究資料４、液状化地盤にお ける基礎設計の考え方，日本建築学会構 造委員会基礎構造運営委員会編，1998 年，ｐｐ． 29−45、1998，ＵＲＬ，ｈ ｔｔｐ：／／ｂｏｈ0709．ｉｎｆｏｓｅ ｅｋ．ｌｉｖｅｄｏｏｒ．ｃｏｍ／Ｔｅ ｃｈｎｏｔｅ／ｐ03−ｍａｔ．ｐｄｆ 吉田望，液状化に伴う流動のメカニズ ム，液状化に伴う流動のメカニズム，地 震時の地盤・土構造物の流動性と永久変 形に関するシンポジウム発表論文集，地 盤工学会，1998年12月31日，ｐｐ． 53 −70，ＵＲＬ，ｈｔｔｐ：／／ｂｏｈ 0709．ｉｎｆｏｓｅｅｋ．ｌｉｖｅｄｏ ｏｒ．ｃｏｍ／Ｔｅｃｈｎｏｔｅ／ｐ01 −ｆｌｏｗ．ｐｄｆ (58)調査した分野(Int.Cl.7，ＤＢ名) G01M 7/00 - 7/02 G01D 21/00 ─────────────────────────────────────────────────── ─── Continued front page (72) Inventor Noriaki Tanaka 4-1, Egasaki-cho, Tsurumi-ku, Yokohama-shi, Kanagawa, Tokyo Electric Power Co., Inc. (72) Inventor Hiroshi Sato E, Tsurumi-ku, Yokohama-shi, Kanagawa No. 4 No. 1 Kagasaki Town, Electric Power Research Laboratory, Tokyo Electric Power Co., Inc. (56) Reference JP-A-5-157655 (JP, A) JP-A-6-330518 (JP, A) JP-A-10-90051 (JP , A) JP 10-307859 (JP, A) JP 11-84017 (JP, A) JP 11-160144 (JP, A) JP 2000-257054 (JP, A) Fumio Yanashita, Hijikata Shoichiro, Isao Nishimura, Kaya Kubo, Masafumi Narukawa, Yutaka Oshima, Yuichiro Uchikawa, Tadahiko Shiomi, Method of Considering Liquefaction in Seismic Response Analysis of Pile Foundation Building by SR Model-Part 9 Effective Stress Solution of Free Ground Based on the analysis, the 1999 Architectural Institute of Japan Conference (China), Academic Lecture Summary B-2, July 30, 1999, pp. 377-378 Tadahiko Shiomi, Yasushi Nakai, Fumio Yanagita, Shoichiro Hijikata, Haruo Yokoyama, Keisuke Koyama, Simplified Liquefaction Analysis Method Considering Cumulative Damage-Part 1 Proposal of Simple Liquefaction Analysis Method, 1999 Architectural Institute of Japan Conference (China) Summary of Academic Lectures B-2, July 30, 1999, p. 379-380 Nozomi Yoshida, How to evaluate the physical properties of ground used for seismic response analysis, Research Material 4 on the design and construction of building foundations, Concept of basic design on liquefied ground, Foundation of the Japan Institute of Architecture Structure Steering Committee, 1998, pp. 29-45, 1998, URL, http: // boh0709. info ek. livedoor. com / Technote / p03-mat. pdf Nozomi Yoshida, Mechanism of flow due to liquefaction, Mechanism of flow due to liquefaction, Proceedings of symposium on fluidity and permanent deformation of ground and soil structures during earthquake, Geotechnical Society, 1998 December 31, pp. 53-70, URL, http: // boh 0709. infoseek. livedo or. com / Technote / p01-flow. pdf (58) Fields investigated (Int.Cl. 7 , DB name) G01M 7 /00-7/02 G01D 21/00
Claims (3)
非線形問題に適用される等価線形解析法に基づく地震動
の動的応答解析法において、 入力地震動を経時的に複数に分割し、分割した部分にそ
れぞれ前記等価線形解析法による解析を実施し、 しかる後に、各解析によって得られた結果を重ね合わせ
る動的応答解析法であって、前記動的応答解析法を経時
的に物性が急激に変化する現象に適用する際に、前記経
時的な分割個所を前記物性の急激変化点に設定するとと
もに、分割した部分にそれぞれ実施される前記等価線形
解析法による解析の物性を別にすることを特徴とする地
震動の動的応答解析法。1. By repeating the linear analysis a plurality of times,
In the dynamic response analysis method for seismic motions based on the equivalent linear analysis method applied to nonlinear problems, the input seismic motion is divided into multiple parts over time, and the divided parts are each analyzed by the equivalent linear analysis method. , A dynamic response analysis method that superimposes the results obtained by each analysis,
When applied to a phenomenon in which physical properties suddenly change,
If you set the temporal division point to the point where the physical properties change drastically,
In the same way, the equivalent linear
Dynamic response analysis method of seismic motion characterized by different physical properties of analysis by analysis method.
の液状化現象 , 地盤に破砕帯が存在する場合の滑り現
象、構造物の浮き上がり現象 , 構造物の塑性化現象から
選択されることを特徴とする請求項１記載の地震動の動
的応答解析法。 2. A phenomenon in which the physical properties change drastically is
Liquefaction Phenomenon of Soil, Slip Behavior in the Presence of Fracture Zone
From the elephant and the floating phenomenon of the structure, the plasticization phenomenon of the structure
The motion of the seismic motion according to claim 1, which is selected.
Response analysis method.
の液状化現象であって、前記物性の急激変化時の前後で
複数に分割し、分割した部分にそれぞれ前記等価線形解
析法による解析を実施する前に、分割しない状態で等価
線形解析法による予備解析を実施して、この予備解析に
基づいて前記地盤の液状化時期を求めることを特徴とす
る請求項２記載の地震動の動的応答解析法。 3. The phenomenon that the physical properties change rapidly is
Liquefaction phenomenon of before and after the sudden change of physical properties
It is divided into multiple parts, and the equivalent linear solution is
Equivalent without splitting before analysis by analysis method
Perform a preliminary analysis using the linear analysis method
The liquefaction time of the ground is obtained based on
The method for analyzing a dynamic response of earthquake motion according to claim 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29835399A JP3488848B2 (en) | 1999-10-20 | 1999-10-20 | Dynamic response analysis method of earthquake motion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29835399A JP3488848B2 (en) | 1999-10-20 | 1999-10-20 | Dynamic response analysis method of earthquake motion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001116651A JP2001116651A (en) | 2001-04-27 |
| JP3488848B2 true JP3488848B2 (en) | 2004-01-19 |
Family
ID=17858596
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP29835399A Expired - Fee Related JP3488848B2 (en) | 1999-10-20 | 1999-10-20 | Dynamic response analysis method of earthquake motion |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3488848B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4616101B2 (en) * | 2005-07-15 | 2011-01-19 | 中部電力株式会社 | Dynamic response analysis |
| JP4987503B2 (en) * | 2007-02-06 | 2012-07-25 | 株式会社東芝 | Ultrasonic diagnostic apparatus, image data display apparatus, and image data display control program |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH05157655A (en) * | 1991-12-04 | 1993-06-25 | Hitachi Ltd | On-line earthquake response load experiment device |
| JP3153385B2 (en) * | 1993-05-21 | 2001-04-09 | 株式会社竹中工務店 | Ground analysis method to capture soil behavior |
| JP3567407B2 (en) * | 1996-09-13 | 2004-09-22 | 清水建設株式会社 | Ground-structure coupled nonlinear seismic response analysis system |
| JPH10307859A (en) * | 1997-05-07 | 1998-11-17 | Hironori Nagai | Computer readable record medium recording analysis program of building structure and elasto-plastic vibration response analysis method for building structure by computer |
| JPH1184017A (en) * | 1997-09-03 | 1999-03-26 | Hitachi Ltd | Estimating device of earthquake damage and medium storing execution program of the same |
| JP3938993B2 (en) * | 1997-11-26 | 2007-06-27 | 関西電力株式会社 | 3D seismic response analysis method |
| JP2000257054A (en) * | 1999-03-11 | 2000-09-19 | Sekisui Chem Co Ltd | Earthquake response analyzing method for ground and liquefaction judging method for ground |
-
1999
- 1999-10-20 JP JP29835399A patent/JP3488848B2/en not_active Expired - Fee Related
Non-Patent Citations (4)
| Title |
|---|
| 吉田望，地震応答解析に用いる地盤物性をどう評価するか，建築基礎の設計施工に関する研究資料４、液状化地盤における基礎設計の考え方，日本建築学会構造委員会基礎構造運営委員会編，１９９８年，ｐｐ． ２９−４５、１９９８，ＵＲＬ，ｈｔｔｐ：／／ｂｏｈ０７０９．ｉｎｆｏｓｅｅｋ．ｌｉｖｅｄｏｏｒ．ｃｏｍ／Ｔｅｃｈｎｏｔｅ／ｐ０３−ｍａｔ．ｐｄｆ |
| 吉田望，液状化に伴う流動のメカニズム，液状化に伴う流動のメカニズム，地震時の地盤・土構造物の流動性と永久変形に関するシンポジウム発表論文集，地盤工学会，１９９８年１２月３１日，ｐｐ． ５３−７０，ＵＲＬ，ｈｔｔｐ：／／ｂｏｈ０７０９．ｉｎｆｏｓｅｅｋ．ｌｉｖｅｄｏｏｒ．ｃｏｍ／Ｔｅｃｈｎｏｔｅ／ｐ０１−ｆｌｏｗ．ｐｄｆ |
| 塩見忠彦、貫井泰、柳下文雄、土方勝一郎、横山治男、小山桂介，累積損傷度を考慮した簡易液状化解析手法−その１ 簡易液状化解析手法の提案，１９９９年度日本建築学会大会（中国） 学術講演梗概集 Ｂ−２，１９９９年 ７月３０日，ｐｐ．３７９−３８０ |
| 柳下文雄、土方勝一郎、西村功、久保賀也、成川匡文、大島豊、内川祐一郎、塩見忠彦，ＳＲモデルによる杭基礎建築物の地震応答解析における液状化の考慮方法−その９ 自由地盤の有効応力解析に基，１９９９年度日本建築学会大会（中国） 学術講演梗概集 Ｂ−２，１９９９年 ７月３０日，ｐｐ．３７７−３７８ |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2001116651A (en) | 2001-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kampitsis et al. | Seismic soil–pile–structure kinematic and inertial interaction—A new beam approach | |
| Militano et al. | Dynamic response of a pile in a multi-layered soil to transient torsional and axial loading | |
| TW202101253A (en) | Analysis method, computer product and device for discontinuous structure | |
| Lu et al. | Computer simulation of the dynamic layered soil pile structure interaction system | |
| CN104462641A (en) | Bridge pile foundation anti-earthquake analysis simplified method considering whole soil liquefaction process | |
| JP3488848B2 (en) | Dynamic response analysis method of earthquake motion | |
| Huh et al. | A novel risk assessment for complex structural systems | |
| JP4616101B2 (en) | Dynamic response analysis | |
| CN111310321B (en) | Layered soil single pile horizontal vibration analysis method based on Pasternak foundation model | |
| Tsiapas et al. | Selective filtering of numerical noise in liquefiable site response analyses | |
| Smyrou et al. | Wavelet analysis for relating soil amplification and liquefaction effects with seismic performance of precast structures | |
| Pender et al. | Rocking controlled design of shallow foundations | |
| JP2004069302A (en) | Dynamic response analyzing method for earthquake motion | |
| Besseling | Soil-structure interaction modelling in performance-based Seismic Jetty Design | |
| Takewaki | Response spectrum method for nonlinear surface ground analysis | |
| JP2006266940A (en) | Method of estimating input earthquake motion | |
| Folić et al. | Analysis of Seismic Interactions Soil-Foundation—Bridge Structures for Different Foundations | |
| Ahmadi et al. | Dynamic analysis of piles in sand based on soil-pile interaction | |
| Riahi et al. | Comparison of the Cosserat Continuum approach with finite element interface models in a simulation of layered materials | |
| Subhani et al. | A study of guided wave propagation in timber pole using spectral finite element method | |
| Biswas et al. | Caisson foundation response during liquefaction induced lateral spreading | |
| Yoshida et al. | Practical three-dimensional effective stress analysis considering cyclic mobility behavior | |
| CN116305480A (en) | Single pile horizontal dynamic response analysis method under seismic S shear wave effect | |
| JP3878626B2 (en) | Impulse response calculation method, apparatus and program | |
| Mandal et al. | Response of rectangular raft foundations under transient loading |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071031 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081031 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081031 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091031 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091031 Year of fee payment: 6 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091031 Year of fee payment: 6 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101031 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111031 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111031 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121031 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131031 Year of fee payment: 10 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |