مهندسی عمران Civil Engineering | آموزش و دانلود PDF | AutoCAD , Etabs ,SAP2000

آموزش حامع مباحث مهندسی عمران و سازه | آموزش نرم افزارهای AutoCAD , Etabs ,SAP2000

مهندسی عمران Civil Engineering | آموزش و دانلود PDF | AutoCAD , Etabs ,SAP2000

آموزش حامع مباحث مهندسی عمران و سازه | آموزش نرم افزارهای AutoCAD , Etabs ,SAP2000

۱۲ مطلب با موضوع «مباحث پیشرفته :: تحلیل طیفی» ثبت شده است

در حین زلزله سه مولفه انتقالی آن به ثبت می‌رسد که دو مولفه افقی و یک مولفه قائم می‌باشد. با استفاده از هر یک از این مولفه‌ها، می‌توان طیف آن را ترسیم نمود. طیف پاسخ، نموداری از حداکثر پاسخ‌های (شتاب، سرعت و یا جابجایی) یک سری سیستم یک درجه آزادی با دوره‌های تناوب مختلف، و میرایی مشخص در برابر یک رکورد زلزله می‌باشد. در واقع یک رکورد موجود به یک سیستم یک درجه آزادی زده شده و مقدار حداکثر پاسخ سیستم محاسبه می‌شود. دوباره با ثابت نگه داشتن مقدار میرایی و تغییر دوره تناوب سیستم، این آزمایش تکرار می‌شود تا طیفی از حداکثر پاسخ‌ها بدست آید. از بهم پیوستن این نقاط، نمودار طیف پاسخ حاصل می‌شود. برای بدست آوردن طیف پاسخ ارتجاعی، معمولاً از انتگرال دیوهامل استفاده می‌شود. انتگرال دیوهامل روش عمومی برای محاسبه پاسخ سیستم یک درجه آزادی خطی به نیروی دینامیکی دلخواه است.
[ تحلیل دینامیکی طیفی[/caption]
 
حال آیین‌نامه‌ها برای راحتی طیفی را ارائه می‌دهند که بدون حالت دندانه بوده و به طیف طرح مشهور است. این طیف، طیف ارتجاعی است. زیرا از حداکثر پاسخ سیستم‌های یک درجه آزادی بدست آمده که هیچ‌گاه غیر ارتجاعی نشده‌اند، بنابراین پاسخ آنها نیز ارتجاعی است. این طیف طبق 2800، همان AB است. مقدار نیروی بدست آمده از این طیف برای طراحی خیلی زیاد است. از طرف دیگر، دوره بازگشت زمین لرزه هایی که آنقدر قوی باشند تا بتوانند چنین نیرویی را در یک سازه ارتجاعی ایجاد کنند، بسیار طولانی است. از این رو، طرح سازه ها برای چنین نیرویی انجام نمی شود، مگر در سازه های بسیار مهم مانند نیروگاه های هسته ای. در حالت معمول فرض می شود که سازه بتواند جابجایی را فراتر از حد الاستیک خود تحمل کند، بدون اینکه فرو بریزد. در چنین سازه-ای، به دنبال زلزله های قوی، تغییر شکل های ماندگاری به وجود می آید که می تواند چندین برابر تغییر شکل تسلیم آن باشد، یعنی با آنکه سازه در طول زلزله پایدار مانده است، نیاز به تعمیرات اساسی یا حتی بازسازی دارد. با فرض اینکه سیستم در طول زلزله می تواند وارد حالت غیرارتجاعی شود، دیگر لازم نیست برای برش پایه الاستیک طراحی شود. اگر این طیف را کاهش دهیم، سازه برای نیروی کمتری طرح شده و در زلزله واقعی انتظار رفتار غیر ارتجاعی از آن داریم. به این طیف کاهش داده شده، طیف غیر ارتجاعی می‌گویند که طبق 2800 خودمان، همان AB/R است.
@AlirezaeiChannel

  • مهندس علیرضا خویه

اطلاعیه دفتر فنی سازمان نظام مهندسی درخصوص تحلیل طیفی سازه ها

 

  • مهندس علیرضا خویه

پارامتر T که دوره تناوب سازه بوده و به مشخصات دینامیکی آن بستگی دارد. برای تعیین آن از روش‌های تجربی، موضوع بند 3-3-3 یا روش‌های تحلیلی استفاده می‌شود. در روش تحلیلی نیاز به ایجاد یک مدل کامپیوتری است.
پارامتر Ts دوره تناوب خاک است. در صورتی که خاک از یک لایه تشکیل شده باشد، با استفاده از رابطه زیر می‌توان مقدار دوره تناوب خاک را تعیین نمود:
Ts=(4H)/vs
که در رابطه فوق، TS دوره تناوب خاک بر حسب ثانیه، H ارتفاع لایه خاک بر حسب متر و vS سرعت موج برشی بر حسب متر بر ثانیه می‌باشد. مقدار سرعت موج برشی در خاک تابعی از نوع خاک آن می‌باشد. در عمل لایه‌های زیادی از خاک وجود داشته و رابطه خیلی پیچیده‌تر است لیکن 2800 در جدول 2-2 مقادیری برای Ts ارائه می‌دهد.
پارامتر T0 تقریبا 20% مقدار Ts بوده و شاخه اول و دوم طیف را از هم تمیز می‌دهد.
پارامتر S شکل طیف را ایجاد می‌کند. شکل طیف را برای خاک‌های مختلف بالا یا پایین می‌کند. همچنین برای خاک‌های سست و لرزه‌خیزی زیاد، اثر غیرخطی شدن خاک در کاهش پاسخ را لحاظ می‌کند.
پارامتر S0 اثرات ساختگاهی را نشان می‌دهد. برای خاک‌های سست اثر بزرگنمایی زلزله بر روی لایه‌های بالایی زمین‌های آبرفتی را لحاظ می‌کند. بدین صورت که در حالتی زمین نوع III یا IV باشد، مقدار طیف B1 از یک شروع نشده و از مقداری بیشتر از 1 شروع خواهد شد.

  • مهندس علیرضا خویه

آموزش تحلیل طیفی

مهندس علیرضا خویه | | ۰ نظر

مطلب از کانال دکتر علیرضایی@AlirezaeiChannel اقتباس شده است.


نکته طراحی در ETABS
روند انجام تحلیل طیفی در برنامه ETABS
در برنامه ETABS برای انجام یک تحلیل طیفی بعد از مدلسازی، سه گام زیر بایستی انجام شود:
1- معرفی تعداد مودهای نوسانی مورد نظر
2- معرفی طیف طرح (شبه‌شتاب) به برنامه
3-تعریف حالت بار طیفی در جهات مورد نظر
در شکل زیر مسیر تعریف تعداد مودهای نوسانی سازه نشان داده شده است. همانطور ابتدا از منوی Analyze گزینه Set Analysis Options را انتخاب نمایید تا پنجره وسط ظاهر شود. در این حالت از قسمت Dynamic Analysis گزینه Set Dynamic Parameters… را انتخاب نمایید تا گام سوم (پنجره سمت راست) نمایان شود. در بخش Number of Modes از پنجره Dynamic Analysis Parameters، تعداد مودهای نوسان سازه را تعیین نمایید. طبق استاندارد 2800 داریم:
تعداد مدهای نوسان، در هر یک از دو امتداد متعامد ساختمان باید به حدی باشد که مجموع جرم‌های مؤثر در آنها بیش تر از 90 درصد جرم کل سازه شود، در نظر گرفته شود.
در اکثر حالات در نظر گرفتن سه مود نوسانی برای هر طبقه کفایت می‌کند. در صورتی که برای Frequency Shift (Center) و Cutoff Frequency (Radius) صفر در نظر گرفته شود، تمام موهای نوسان سازه محاسبه شده و در صورتی که مرکز و شعاع قطع فرکانس‌های سازه داده شود، تنها مودهایی که در این حیطه باشند در نظر گرفته می‌شوند. مقدار تلرانس همگرایی نسبی (Relative Tolerance) بطور پیش فرض برابر 1E-7 بوده که مقدار مناسبی می‌باشد

بعد از تعریف مودها، مطابق شکل زیر، بایستی طیف طرح غیرارتجاعی به برنامه معرفی شود. در این حالت از مسیر نشان داده شده در شکل زیر(سمت چپ- گام اول) پنجره Define Response Spectrum Functions ظاهر خواهد شد (پنجره وسط،گام دوم). در صورتی که بخواهید مولفه‌های طیف را بطور دستی و یکی یکی وارد نمایید، از قسمت Choose Function Type to Add گزینه User Spectrum را انتخاب نموده و طیف را معرفی نمایید.

در صورتی که به مانند شکل زیر قبلاً در یک فایل متنی (فایل Notepad) طیف بازتاب را ساخته باشید می‌توانید با استفاده از گزینه Spectrum from File از زیر بخش Choose Function Type to Add، و زدن دکمه Add New Function… بطور مستقیم و سریع طیف بازتاب را تعریف نمایید.
در صورتی که روش دوم را انتخاب نمایید را انتخاب نمایید، پنجره سمت راست بالا ظاهر شده (پنجره Response Spectrum Function Definition) و با استفاده از دکمه Browse… مسیر فایل متنی را انتخاب نمایید. در صورتی که فایل متنی حاوی طیف شبه شتاب بوده که محور افقی آن دوره تناوب باشد، گزینه Period vs Value و در صورتی که محور افقی آن فرکانس (عکس دوره تناوب) باشد، گزینه Frequency vs Value را از زیر بخش Values are: انتخاب نمایید. در صورتی یک سطر یا بیشتر از فایل متنی حاوی متنی برای مشخص نمودن محتویات فایل بوده و جزو طیف به حساب نمی‌آید، در قسمت Header Lines to Skip تعداد این سطرها را برای صرف نظر نمودن برنامه از خواندن آنها وارد کنید. با زدن دکمه View File محتویات فایل توسط برنامه Notepad یا برنامه‌های مشابه آن باز خواهد شود و شما قادر به دیدن محتویات فایل انتخابی هستید. توجه داشته باشید، در صورتی که فایل محاسبات خود را بخواهید به جای دیگری منتقل نمایید حتماً بایستی فایل متنی حاوی طیف پاسخ نیز همراه فایل شما باشد. برای جلوگیری از این مشکل دکمه Convert to User Defined را انتخاب نمایید تا تمام محتوای فایل متنی به فایل EDB منتقل شوند. در این حالت نیازی به همراه بودن فایل متنی نبوده و طیف پاسخ، در فایل شما ماندگار خواهد شد. در بخش Function Graph نیز در صورتی که گزینه Display Graph زده شود، شکل طیف نمایش داده می‌شود. مقدار میرایی نیز 5% در نظر گرفته شود.
 
در شکل زیر مسیر تعریف حالت بار طیفی نشان داده شده است. بعد از ظاهر شدن پنجره Define Response Spectra بر روی گزینه Add New Spectrum… از قسمت Spectra click to: کلیک نمایید تا پنجره Response Spectrum Case Data ظاهر شود. در صورتی بخواهید یک حالت بار طیفی ساخته شده را مشاهده نمایید، بعد از انتخاب آن بر روی گزینه Modify/Show Spectrum… و در صورتی که خواهید یک حالت بار طیفی را پاک نمایید بعد از انتخاب آن بر روی گزینه Delete Spectrum کلیک نمایید. در پنجره Response Spectrum Case Data، و در قسمت Spectrum Case Name نام حالت طیفی و در بخش Structural and Function Damping مقدار میرایی که معولاً 5% در نظر گرفته می‌شود وارد نمایید. در قسمت Modal Combination روش ترکیب آثار مدها را روش CQC قرار دهید.
این روش ترکیب بهترین روش ترکیب آثار مودها برای سازه‌های سه بعدی است. در صورتی که میرایی سازه صفر باشد، روش SRSS و CQC یکسان خواهند بود. در صورت استفاده از روش ABS نتایج تحلیل طیفی محافظه کارانه خواهند شد. در بخش Directional Combination نحوه ترکیب نیروهای حاصل از حالت طیفی را در جهات مختلف سازه را تعیین نمایید. در ارتباط با این گزینه در ادامه توضیحات بیشتر ارائه خواهد شد. در بخش Input Response Spectra طیف مورد نظر را در جهت مورد نظر خود معرفی نمایید. در صورتی که بخواهید طیف را در جهت محور X کلی به سازه اعمال نمایید از جعبه کرکره‌ای U1 (از قسمت Direction) تابع (Function) مورد نظر که همان طیف باشد را معرفی نمایید و در صورتی که بخواهید جهت محور Y کلی به سازه اعمال نمایید از جعبه کرکره‌ای U2 طیف ساخته شده در گام قبل را معرفی نمایید. در بخش Scale Factor بایستی ضرب مناسبی مطابق با الزامات آیین‌نامه طراحی وارد نمایید. این گزینه طبق استاندارد 2800 همان AI/Rمی‌باشد که در واقع مقدار B قبلاً توسط طیف، به برنامه داده شده است و با ضرب B در این ضریب، مقدار ضریب زلزله هر مود نوسانی مطابق با مقدار دوره تناوب آن مود حاصل می‌شود. گزینه Excitation angle زاویه اعمال طیف را نسبت به محور X کلی سازه مشخص می‌کند. در صورتی که در بخش Excitation angle زوایه 90 درجه انتخاب شده باشد و در بخش Direction، گزینه U1 برای اعمال طیف در نظر گرفته شده باشد، طیف تعریف شده با زوایه 90 درجه نسبت به محور X (به عبارتی در جهت Y) بر سازه اعمال می‌شود.

تعریف حالت بار طیفی در جهات مورد نظر

  • مهندس علیرضا خویه

مفاهیم تحلیل طیفی

مهندس علیرضا خویه | | ۰ نظر

- می دانیم یک سازه چند درجه آزاد دارای فرکانس های ارتعاشی متفاوتی است که کمترین آن مربوط به مود اول و بیشترین آن مربوط به مود آخر است . از طرفی زلزله های مختلف نیز دارای محتوای فرکانسی متفاوتی هستند. از دینامیک سازه‏ها می دانیم که اگر فرکانس بارگذاری که به یک سیستم اعمال می شود، نزدیک به فرکانس طبیعی آن سیستم باشد، در آن سیستم حالت تشدید رخ می‏دهد . با توجه به همین واقعیت، اگر محتوای فرکانسی غالب یک زلزله نیز، نزدیک به یکی از فرکانس‏های ارتعاشی سازه مورد نظر باشد‏، باید در آن مود خاص ، تشدید رخ دهد. اثر مشارکت جرمی در مودها در یک سازه به منظمی آن بستگی دارد. هرچه سازه منظم‌تر باشد، درصد مشارکت مود اول بالاتر خواهد بود.
چرا هیچوقت در یک مد ۱۰۰% اثر مشارکت جرمی نداریم ؟
- اینکه چرا هیچوقت در یک مد ۱۰۰% اثر مشارکت جرمی نداریم، باید به اصول تحلیل مودال برگردیم. دو رهیافت کلی برای تحلیل لرزه‌ای سازه‌های چند درجه آزاد وجود دارد. در رهیافت اول، می‌توان پاسخ هر درجه‌ی آزادی را در هر لحظه‌ی زمانی از وقوع زلزله محاسبه نمود و در نهایت تاریخچه‌ی زمانی پاسخ را محاسبه و ترسیم نمود، این روش تحلیل، اصطلاحاً «روش تحلیل تاریخچه‌ی زمانی پاسخ» نام دارد. پاسخ می‌تواند جابجایی، سرعت و یا شتاب باشد. در رهیافت دوم که «روش تحلیل طیف پاسخ» یا «روش تحلیل طیفی» نام دارد، حداکثر پاسخ محتمل سازه بدون داشتن اطلاعات تاریخچه‌ی زمانی تخمین زده می‌شود. یعنی مهندس محاسب بدون اطلاع از کل تاریخچه‌ی زمانی پاسخ، حداکثر پاسخ محتمل را ـ که در طراحی بسیار با اهمیت است ـ به دست می‌آورد. حداکثر پاسخ از آن جهت اهمیت دارد که با داشتن آن می‌تواند بدترین شرایط سازه را بررسی نمود و سازه را برای مقاومت در برار آن طراحی نمود. در عمل، معمولاً مهندسان روش دوم را ترجیح می‌دهند. زیرا در تحلیل تاریخچه‌ی زمانی حجم بسیار زیادی از داده‌ها تولید می‌شود. این داده‌ها پاسخهای مختلف برای کلیه‌ی درجات آزادی سازه است که بررسی آن دشوار و زمان‌بر است. در این روش یک سازه nدرجه آزاد به nتا سیستم یکدرجه آزاد تبدیل میشود. یعنی اینکه ارتعاش پیچیده (چیزی در واقعیت رخ می‌دهد) را میتوان به n تا ارتعاش ساده تبدیل کرد. بنابراین اگر سازه شما یک درجه آزادی داشته باشد، 100% مشارکت جرمی آن سهم آن یک مود ارتعاشی آن است.
- اصولاً معیاری از نظر آیین‌نامه برای درصد قابل قبول جهت مشارکت جرمی در آن مود وجود ندارد و ملاک میزان درصد تجمی جرمی است.
علت اینکه در ۳ مد اول و یا حتی در مد اول ما مشارکت جرمی کمی داریم (درصد مشارکت مثلا زیر ۱۰% برای هر راستا) چه هست ؟ (در مدهای بالاتر اثر مشارکت شکوفا میشود)
- اگر درصد مشارکت جرمی در مود اول پایین باشه، نشان دهنده نامنظمی سازه و اثر قابل ملاحظه در مودهای بالا است. این نامنظمی میتواند در جرم یا در سختی باشد. یعنی در این حالت سازه شما نامنظم بوده که این اتفاق افتاده.
- ملاک همان مجموع 90% است. لزوماً غلط مدل نشده است. البته اگر سازه شما منظم باشد و این اتفاق ایجاد شود باید حساس شوید و مدل را بازبینی کنید.
-اگر در مود اول درصد مشارکت جرمی 40% باشد، بقیه جرم سازه در مودهای دیگر نهفته است. اشکال مودی و فرکانس‌های متناظر آن‌ها (یا زمان‌های تناوب متناظر)، جزو مشخصات بسیار مهم هر سازه‌ی چنددرجه آزاد است که در تعیین پاسخ‌ آن‌ها تحت اثر شرایط ارتعاش اجباری، از جمله تحریک پایه‌ی زلزله، نقش بسیار اساسی ایفا می‌کند. اشکال مودی و فرکانس‌های آن‌ها با روش‌های ریاضی قابل محاسبه است.
برای دیدن جزئیات بیشتر در ارتباط با این سوال پیشنهاد میکنم این مقاله را مطالعه نمایید:
* حسینی، محمود؛ یعقوبی وایقان، فریبرز. "نگاهی دقیق تر به مفهوم ضریب مشارکت جرمی و تعاریف آن در تحلیل لرزه ای سازه ها"، پژوهشنامه پژوهشگاه بین المللی زلزله شناسی و مهندسی زلزله، شماره 3، پاییز 1378
 
منبع: کانال دکتر علیرضایی

  • مهندس علیرضا خویه

نکاتی در مورد مودهای مختلط

مهندس علیرضا خویه | | ۰ نظر
  • نکاتی در مورد مودهای مختلط
  • در حالت اول، ماتریس های جرم و سختی معمولاً متقارن می باشند. در صورت وجود نیروهای ژیروسکوپی ناشی از دوران، یک ماتریس میرایی پادمتقارن پدیدار می گردد. در این حالت حتی اگر هیچ ماتریس میرایی ویسکوز و یا سازه ای وجود نداشته باشد، مساله مقدار ویژه منجر به شکل مودهایی مختلط خواهد شد.
  • در حالت دوم، اگر میرایی غیرتناسبی باشد بردارهای ویژه سیستم، مختلط خواهند بود. اگر معادلات در فضای حالت نوشته شوند، مشاهده خواهد شد بردارهای ویژه بدست آمده که می توانند دو ماتریس سیستم را قطری نمایند، نمی توانند ماتریس میرایی را قطری نمایند.
  • بر خلاف مودهای حقیقی، المان های شکل مود دارای زوایای فاز مختلفی بین 0 تا 180 درجه می باشند.
  • این بدان معناست که نقاط واقع بر سازه، با وجود آنکه با فرکانس یکسانی ارتعاش می نمایند، بطور همزمان از موقعیت تعادل عبور نخواهند کرد. در نتیجه، نقاط گره ای سازه ثابت نخواهند بود. این نقاط با فرکانسی که توسط میزان مختلط بودن شکل مود تعیین می شود، حرکت می کنند.
  • مشکل واقعی مودهای مختلط در آنالیز مودال بررسی و شبیه سازی آنها نمی باشد، بلکه مشکل در شناسایی آنها به کمک مفاهیم و تئوری استخراج داده های مودال است.
  • در چنین حالتی معمولاً پاسخ های صحیح برای مقایسه وجود ندارند. در نتیجه ممکن است از یک سری اطلاعات FRF نتایج و تفسیرهای مختلفی حاصل شود.
  • سوال اصلی این است که چگونه می توان از صحت مودهای مختلط حاصل، اطمینان یافت. عوامل مختلفی در آنالیز مودال وجود دارند که می توانند به استخراج مودهای مختلط محاسباتی منجر شوند. بعنوان مثال، اگر دو مود حقیقی نزدیک به هم به درستی و با دقت لازم تحلیل نشوند، ممکن است به دو مود مختلط تبدیل شوند. وجود انحراف در فاز ثابت های مودال نیز ممکن است به مودهای مختلط غیر واقعی منجر شود. وجود خواص غیرخطی نیز ممکن است باعث استخراج مودهای مختلط نادرست گردد.
  • تعیین کمّی میزان مختلط بودن یک مود مختلط مشکل است اما به پارامتری برای این کار جهت تحلیل مودها نیاز است. روشن ترین معیار مختلط بودن یک شکل مود، مقدار زوایای فاز المان های آن است.
  • هر چقدر این زوایا از 0 و 180 درجه دورتر باشند، بخش موهومی بزرگتری خواهند داشت. حداکثر این انحراف زاویه، 90 درجه می باشد. با این وجود، می بینیم که انحراف فاز تنها چیزی نیست که حائز اهمیت است. اگر تمام المان ها در یک مود دارای فاز 90 درجه باشند آن مود اصولاً مختلط نخواهد بود.
  • بنابراین، می توانیم پارامتری را بعنوان معیاری برای مختلط بودن مودال تعریف کرده و آن را ضریب مختلط بودن(؟!) مودال بنامیم. این ضریب، معیاری کمّی از انحراف نرمال فاز شکل مود مختلط با محاسبه مقدار متوسط اختلاف فاز با فاز میانگین می باشد. درصد نسبت بین این مقدار با ماکزیمم انحراف 90 درجه بیانگر میزان مختلط بودن مودال می باشد.
  • مهندس علیرضا خویه

فایل سمینار (به صورت پاورپوینت) که در آن به  نحوه ی به دست آوردن طیف­ های خطر یکنواخت ، طیف­ های میانگین شرطی و طیف طرح مقیاس شده  می پردازد که برای تحلیل و طراحی سازه­ ها می­توانند مورد استفاده قرار گیرند.

همچنین در این سمینار به بررسی و مقایسه طیف های یکنواخت با طیف های ایین نامه 2800 پرداخته شده است.

UNIFORM HAZARD SPECTRUM & CONDITIONAL MEAN SPECTRUM

تعداد اسلایدها: 30

به صورت مصور و همراه با گراف های متعدد می باشد.

هزینه: 50 هزار تومان

 

برای سفارش می توانید به شماره: 09382904800 پیامک دهید یا در تلگرام پیغام ارسال کنید.

  • مهندس علیرضا خویه

نکته ای کوچک در تحلیل طیفی

مهندس علیرضا خویه | | ۰ نظر

در تحلیل طیفی همیشه در نظر داشته باشید که طیف آیین نامه ی 2800 را می بایست در ضریب AI ضرب نمایید

  • مهندس علیرضا خویه

دریافت فایل
عنوان: استاتیکی معادل مودال(طیفی) و تاریخچه زمانی
حجم: 954 کیلوبایت
توضیحات: استاتیکی معادل مودال(طیفی) و تاریخچه زمانی

  • مهندس علیرضا خویه

عوامل مختلفی بر روی طیف های زلزله موثر می باشد و هدف ما،
بررسی نحوه ی تاثیر این عوامل بر روی طیف های زلزله می باشد.

مهمترین این موارد به شرح زیر می باشد:

1- نوع خاک منطقه

2- بزرگی زلزله (ریشتر)

3- مکانیزم نوع گسل

4- فاصله

5- مدت زمان

در این پروژه تاثیر چهار عامل اول را بر روی طیف های
پاسخ بررسی میکنیم

دریافت فایل

عنوان: بررسی پارامتر های مختلف بر طیف پاسخ زلزله
حجم: 1.16 مگابایت
توضیحات: بررسی پارامتر های مختلف بر طیف پاسخ زلزله

  • مهندس علیرضا خویه