تعیین تحریک بحرانی سیستم‌های چند درجه‌آزاد بر اساس تابع هدف جابه‌جایی و قید انرژی

نوع مقاله : Articles

نویسندگان

دانشگاه تربیت مدرس، تهران، ایران

چکیده

در مقاله‌ی حاضر کاربرد روش تحریک بحرانی احتمالاتی برای تابع هدف جابه‌جایی در سیستم‌های چند ‌درجه‌آزاد مورد بررسی قرار گرفته است. نشان داده می‌شود که حد بالای انرژی ورودی کل زلزله در واحد جرم می‌تواند یک معیار منطقی برای تخمین محدوده‌ی مجاز انرژی زلزله‌های مشابه احتمالی باشد؛ بنابراین می‌توان این حد را ثابت نگه‌داشته و در مقابل با تغییر در دامنه و محتوای فرکانسی تحریک، رکورد دیگری که تابع هدف انتخابی را بیشینه نماید، تولید نمود. با معرفی این حد به‌عنوان قید مسئله، روش تحریک بحرانی برای تولید شتاب‌نگاشت مصنوعی سازگار با سه مدل قاب برشی 8، 14 و 20 طبقه به‌کار گرفته می‌شود. به‌منظور نشان دادن کارایی روش، با استفاده از تحریک‌های بحرانی یافته شده و سه رکورد واقعی متناظر با سه مدل، به‌عنوان تحریک‌های بحرانی فرضی، تحلیل دینامیکی خطی انجام گرفته و کمیت‌های مورد علاقه، شامل بیشینه‌ی جابه‌جایی، دریفت و شتاب طبقات مورد مقایسه قرار می‌گیرند. نتایج نشان می‌دهد که شتاب‌نگاشت‌های مصنوعی ساخته‌شده می‌توانند تخمین قابل‌قبولی از رفتار سازه به دست دهند.

کلیدواژه‌ها


عنوان مقاله [English]

Determining Critical Excitation of MDOF Systems for Displacement Objective Function and Energy Constraint

نویسندگان [English]

  • Arsalan Bazrafshan
  • Naser Khaji
Tarbiat Modares University, Tehran, Iran
چکیده [English]

Considering the increasing use of dynamic analysis methods in structural design, the selection of appropriate design earthquake has been an important part of the design procedure. For intermediate and low importance structures, the design earthquake is typically provided by the seismic codes as a design spectrum. However, for important structures for which time-history analysis should be performed, the use of recorded ground motions as input is inevitable. On the other hand, the existing ground motions only show a small part of the reality. Experiences of past earthquakes indicate that sole reliance on existing data will never resolve all issues, and new damage problems have occurred recently. In order to overcome this problem, a new paradigm has to be used. The concept of “critical excitation” and the structural design based on this concept can become one of such new paradigms.
In the present paper, the probabilistic critical excitation method is used to determine the critical excitations for three shear building models, which are modeled as MDF system. Selection of appropriate constraints is the main problem when using this method. It is shown that the upper bound of earthquake input energy per unit mass can be considered as suitable constraint for the critical excitation. This bound can be a reasonable benchmark to estimate the allowable range of possible energy in similar earthquake. Considering this bound as constraint, the method is used to determine the critical PSD functions and generating synthetic accelerograms. In order to demonstrate the effectiveness of the method, three real accelerograms are selected as benchmarks and linear dynamic analysis was conducted using these accelerograms and the generated critical excitations, and the key parameters of response including maximum displacement, drift and acceleration of stories are compared.
Input base acceleration is defined as the product of an envelope function and a stationary Gaussian process with zero mean. Considering the sum of the mean-square interstory drift of the system as the objective function, the critical excitation problem is defined as follows:
Given the mass, stiffness and the viscous damping matrix of a linear elastic MDOF system, as well as the envelope function, find the critical PSDF, so that the objective function is maximized under specific constraints.
In order to solve this problem, first the constraints have to be selected. According to the concept of the constraints, there is no straight way to recommend a specific value of them. Thus that is very difficult to estimate that what value of one constraint is sufficient to make the final nonstationary input critical. Nevertheless, at least a reasonable assumption can be made for the ratio of constraints. Therefore, the method can be used to locate the rectangular function.
Takewaki [1] by setting constraints on the acceleration and velocity time history of the ground motion determined an upper bound of total input energy per unit mass of a record for a damped linear elastic system. This bound is well defined by two curves that perfectly bound the actual input energy curve in the range of short and long periods. Investigations on the time-history of various ground motions indicate that even for the same level of energy bound for acceleration constraint, the maximum amount of actual energy is not constant. Moreover, this maximum value may occur in different periods. The upper bound of the total input energy for acceleration constraint can be used to estimate the possible range of energy in similar ground motions. This bound is only related to the area of the PSDF of excitation and the maximum value of it. These two parameters can define a class of ground motions that existing record is just one realization of them. As a result it is desirable to fix the upper bound of input energy and let the excitation to change in amplitude and frequency content in such a way that the objective function is maximized. By this mean, using a sample of ground motion, a synthetic accelerogram can be generated with the same upper bound but different energy content.
In order to examine the effect of the critical excitations on the structures, a numerical simulation was carried out using the three typical shear building of 8, 14 and 20 stories. These buildings were designed as steel moment frames, using the ETABS commercial software, and then the mass and stiffness of stories were determined and the corresponding MDOF models were used for linear dynamic analysis.
This study shows that the upper bound of the total input energy for acceleration constraint can be a reasonable benchmark to estimate the possible range of energy in similar earthquakes. This bound is only related to the area of the PSDF of excitation and the maximum value of it. By controlling these two parameters, the upper bound of input energy can be fixed and then by changing the amplitude and frequency content using the critical excitation method required accelerograms are determine so that the objective function is maximized. Comparison of linear dynamic analysis of three designed models under critical excitations and the corresponding actual records, which have been selected in such a way that the energy content of the record is the maximum in the natural frequency of the fundamental mode of structures, shows that the synthetic accelerograms can reasonably estimate the behavior of structures, including the maximum displacement, interstory drift and absolute acceleration of stories.
References

Takewaki, I. (2004) Bound of earthquake input energy. J. Struct. Eng., 130(9), 1289-1297.

کلیدواژه‌ها [English]

  • Critical Excitation
  • nonstationary random process
  • Total Input Energy per Unit Mass
  • Synthetic Accelerograms
  1. Drenick, R.F. (1970) Model-free design of aseismic structures. Journal of Engineering Mechanics, 96(4), 483-493.
  2. Shinozuka, M. (1970) Maximum structural response to seismic excitations. Journal of Engineering Mechanics, 96(5), 729-738.
  3. Srinivasan, M., Ellingwood, B., and Corotis, R. (1991) Critical base excitations of structural systems. Journal of Engineering Mechanics, 117(6), 1403-1422.
  4. Drenick, R.F. (1973) Aseismic design by way of critical excitation. Journal of Engineering Mechanics, 99(4), 649-667.
  5. Wang, P.C., Wang, W., Drenick, R., and Vellozzi, J. (1976) Critical excitation and response of free standing chimneys. Proc. of the Int. Symposium on Earthquake Struct. Engrg, vol. I. Mo., St. Louis. Aug., 269-284.
  6. Drenick, R.F., Wang, P.C., Yun, C.B., Philippacopoulos, A.J. (1980) Critical seismic response of nuclear reactors. J. Nuclear Eng. and Design, 59, 425-435.
  7. Ahmadi, G. (1979) On the application of the critical excitation method to aseismic design. J. Struct. Mech., 7(1), 55-63.
  8. Takewaki, I. (2001b) Nonstationary random critical excitation for acceleration response. Journal of Engineering Mechanics, 127(6), 544-556.
  9. Westermo, B.D. (1985) The critical excitation and response of simple dynamic systems. J. Sound and Vibration, 100(2), 233-242.
  10. Takewaki, I. (2004) Bound of earthquake input energy. J. Struct. Eng., 130(9), 1289-1297.
  11. Takewaki, I. (2006) Probabilistic critical excitation method for earthquake energy input rate. Journal of Engineering Mechanics, 132(9), 990-1000.
  12. Takewaki, I. (2004) Critical envelope functions for non-stationary random earthquake input. Computers & Structures, 82(20-21), 1671-1683.
  13. Iyengar, R.N. (1972) Worst inputs and a bound on the highest peak statistics of a class of non-linear systems. J. Sound and Vibration, 25(1), 29-37.
  14. Drenick, R.F. (1977) The critical excitation of nonlinear systems. J. Appl. Mech., 44(2), 333-336.
  15. Philippacopoulos, A. and Wang, P. (1984) Seismic inputs for nonlinear structures. Journal of Engineering Mechanics, 110(5), 828-836.
  16. Takewaki, I. (2002) Critical excitation for elastic–plastic structures via statistical equivalent linearization. Probabilistic Engineering Mechanic, 17(1), 73-84.
  17. Iyengar, R.N. and Manohar, C.S. (1985) System dependent critical stochastic seismic excitations. In: M15/6, Proc. of the 8th Int. Conf. on SMiRT. Belgium, Brussels.
  18. Iyengar, R. and Manohar, C. (1987) Nonstationary random critical seismic excitations. Journal of Engineering Mechanics, 113(4), 529-541.
  19. Srinivasan, M., Corotis, R., and Ellingwood, B. (1992) Generation of critical stochastic earthquakes. Earthquake Engineering and Structural Dynamics, 21(4), 275-288.
  20. Manohar, C.S. and Sarkar, A. (1995) Critical earthquake input power spectral density function models for engineering structures. Earthquake Engineering and Structural Dynamics, 24(12), 1549-1566.
  21. Takewaki, I. (2001) A new method for non-stationary random critical excitation. Earthquake Engineering and Structural Dynamics, 30(4), 519-535.