مهندسی عمران-زلزله(محمدجواد خسرویانی)
مهندسی عمران-زلزله(محمدجواد خسرویانی)

مهندسی عمران-زلزله(محمدجواد خسرویانی)

محاسبه فرکانس و شکل مودهای سازه

پس از محاسبه سختی و تشکیل ماتریس جرم می توانید با استفاده از کد زیر در متلب شکل مودها و فرکانس سازه را به صورت مرتب شده از کوچک به بزرگ بدست آورید

K=xlsread('stiffness Matrix ETABS.xlsx',1,'A1:L12');
M=xlsread('Mass Matrix dasti.xlsx',1,'A1:L12');
[V,D]=eig(K,M);
[W,k]=sort(diag(D));
V=V(:,k);
Factor=diag(V'*M*V);
Phi=V*inv(sqrt(diag(Factor)));
Omega=diag(sqrt(Phi'*K*Phi));

 

شکل مود های نرمال شده می باشد  Phi  همان فرکانس ها می باشد و  Omega

نظرات 13 + ارسال نظر
محمد سه‌شنبه 11 خرداد 1400 ساعت 16:35

سلام
وقتتون بخیر
ممنون از مطالب مفیدتون
در خصوص کد بالا یک سوال برایم پیش اومده
چرا در محاسبه ی Omega بطور مستقیم از W ( که در اصل مرتب شده ی D است) استفاده نکردید و فرکانس ها را با استفاده از شکل مودهای نرمال شده بدست آوردید؟
به عبارت دیگر بنظر میرسد که همان ابتدا با یک عمل جذر از W میشد به فرکانس ها رسید.
با تشکر فراوان

salam
bale sahih hast
in code baraye yek kare khas tarhi shode bod

مجید پنج‌شنبه 2 اردیبهشت 1400 ساعت 20:28

سلام جناب دکتر
سوالی داشتم آیا این امکان وجود داره ماتریس جرم بر روی درایه های غیر از قطر اصلی غیرصفر شود؟در چه مواقعی این اتفاق می افتد
باتشکر از شما

سلام هر المان در ماتریس، بخشی از قطرهای آن است ، بنابراین در عمل هیچ المان غیر قطری وجود ندارد. المان صفر در ماتریس سختی و (یا ماتریس جرم ) به این دلیل است که شبکه های خاص ( Special grids) هیچ ارتباطی با یکدیگر ندارند و بنابراین هیچ مقداری در ماتریس نشان داده نمی شود. قططر اصلی ماتریس همیشه دارای مقادیر صفر است ، زیرا المان ها ارتباط بین شبکه خود را نشان می دهند. از لحاظ عددی نیز ، این مقدار باید صفر باشد زیرا شما معادله ای بین یک شبکه( grid) و آن المان مورد نظر را حل نمی کنید ، درصورتی که این معدلات و روش حل همیشه بین یک شبکه و شبکه های همسایگی آن هست.
بنابراین وجود کانکشن نشان دهنده عدد و صفر عدم ارتباط می باشد می تونید به ماتریس های جرم لامپ (lump) و همچنین تحلیل های ماتریسی مراجعه کنید. که ماتریس جرم لامپ جرم سازه را در بین نقاط مختلف ( در درجات آزادی تعریف شده ) تقسیم می کند.
اما در ماتریس های جرم سازگار به دلیل این کانکشن ها در همه درایه ها مقدار وجود دارد و صفرنیست نمونه المان های جرم گسترده فصل 8 دینامیک سازه چوپرا موجود است.
در پایان به تعریف ماتریس جرم توجه کنید که هر درایه های ماتریس جرم، بردار نیروی لازم برای شتاب واحد هستند

Parvin سه‌شنبه 9 مرداد 1397 ساعت 17:23

سلام در محاسبه فرکانس بار وارده مهم نیس؟
من میخوام دو تا پارامتر از میرایی رایلی به دست بیارم ... برای محاسبه احتیاج به فرکانس دارم ...باید از همون رابطه زیر حساب کنم؟
چون توی اثبات رابطه فرکانس از ارتعاش ازاد استفاده شده اخه

salam

w be bare vorodi rabti nadare
faghat pasokh ha ba bre vorodi taghir mikonaand

bale bad az mohasebe w az rabete raighly mitond zarayeb va C=algha(m)+beta{(k) ro hesav va mirai bedast biyad

شبنم دوشنبه 21 خرداد 1397 ساعت 13:46

با سلام، ببخشید چطور میشه تابعی بنویسیم ک وقتی سختی وجرم و تعدادطبقات را دادی ماتریس جرم و سختی سازه رابسازه و بعد از ان ماتریس جرم مودال و سختی مودال را بسازه

salam
dar har sorat bayad M va K ro beshesh bedin
M' va k' ro hesab kone
manzoreton chi hast
be linke zir ar bezanid
http://mkhosraviani.blogsky.com/1392/10/03/post-32/

حکم ابادی سه‌شنبه 7 آذر 1396 ساعت 11:56

با سلام.
کدی برای محاسبه ماتریس جرم و فرکانسهای طبیعی یک قاب سه بعدی میخواستم. ممنون میشم اگر کمکی کرده باشید.

salam

matrise jerm ro ke be rhati mishe az barname ha va ya hata dasti bedast avord hata agar saze picheshi bashe
dar nahayat biyarin to hamin code mojod dar blog w bedast miyad

پورعلی دوشنبه 6 آذر 1396 ساعت 03:29

با سلام
بنده کد برنامه متلب برای استخراج ماتریس سختی با روشهای پنالتی و حذفی نوشتم اما نمیدانم چرا در آخر به یک جواب نمیرسم.
لطفا اگر کدی آماده برای استخراج ماتریس سختی به روش حذفی دارید برایم بفرستید.
از کمکتان سپاسگزارم

salam
kheyr vali shayad linke zir mofid bashe
http://mkhosraviani.blogsky.com/1393/11/26/post-51/%D9%85%D8%AD%D8%A7%D8%B3%D8%A8%D9%87-%D9%85%D8%A7%D9%86%D8%B1%DB%8C%D8%B3-%D8%B3%D8%AE%D8%AA%DB%8C-%D9%88-%D8%AC%D8%B1%D9%85-%D8%A7%D8%B2-%D8%B7%D8%B1%DB%8C%D9%82-%D9%86%D8%B1%D9%85-%D8%A7%D9%81%D8%B2%D8%A7%D8%B1-ETABS

خالد پورعلی پنج‌شنبه 18 آبان 1396 ساعت 22:29

با سلام من کد استخراج ماتریس سختی را با استفاده از روش گوس میخوام ؟؟اگر دارید لطفا برام بفرستید..ممنون

salam
yek file ersal kardam baraton

فرشید یکشنبه 14 شهریور 1395 ساعت 21:36

با سلام ... مهندس لطفا در مورد نحو ه ی استفاده از شکل مود وفرکانس در تشخیص سازه با الگوریتم ژنتیک راهنمایی کنید مرسی

علی اکبر سه‌شنبه 11 اسفند 1394 ساعت 09:51

با سلام.....مهندس اگر واستون مقدوره در مورد کد مختصر توضیح داده و هم چنین از چه فرمول هایی استفاده نمودید.....با تشکر

salam

kole ravesh hamon ravebete dynamici hastesh ke det(K-,Mw^2) hast ke be sortae eig value dar matlab neveshte shode
va mod ha ham norma shodan va W va frequency mohasebe mishe

MaliheShahvali شنبه 9 خرداد 1394 ساعت 08:26

سلام ببخشید چه طور ماتریس سختی و جرم از ایتبس به این کد بدیم

MaliheShahvali جمعه 8 خرداد 1394 ساعت 22:37

سلام ببخشید چه جوری ماتریس سختی و جرم از ایتبس به این کد بدیم

ba salam
agar be linke zir morajee konid kamel tozih dadam
onvaght pas az mohasebe K va M mitonid on ro dar matlab baz konid ya dar excel berizid az matlab bekhonidesh
moafagh bashid
http://mkhosraviani.blogsky.com/1393/11/26/post-51/%D9%85%D8%AD%D8%A7%D8%B3%D8%A8%D9%87-%D9%85%D8%A7%D9%86%D8%B1%DB%8C%D8%B3-%D8%B3%D8%AE%D8%AA%DB%8C-%D9%88-%D8%AC%D8%B1%D9%85-%D8%A7%D8%B2-%D8%B7%D8%B1%DB%8C%D9%82-%D9%86%D8%B1%D9%85-%D8%A7%D9%81%D8%B2%D8%A7%D8%B1-ETABS

عباس شنبه 5 اردیبهشت 1394 ساعت 00:33

سلام ببخشید چه جوری ماتریس سختی و جرم از ایتبس به این کد بدیم اگر میشه راهنمایی کنید

salam mitonid be linke zir dar belogam berid kamel tozih dadam
moafagh bashid

amin سه‌شنبه 14 بهمن 1393 ساعت 13:30

سلام ببخشید چه جوری ماتریس سختی و جرم از ایتبس به این کد بدیم

slalam be emailetan ersla shod

برای نمایش آواتار خود در این وبلاگ در سایت Gravatar.com ثبت نام کنید. (راهنما)
ایمیل شما بعد از ثبت نمایش داده نخواهد شد