متلب به عنوان یک زبان برنامهنویسی سطح بالا و محیط محاسبات عددی، ابزاری بسیار مفید برای مهندسان است. این نرمافزار به دلیل قدرت محاسباتی بالا، قابلیت تحلیل دادههای پیچیده، و توانایی شبیهسازی مسائل مهندسی، به یکی از ابزارهای ضروری در مهندسی تبدیل شده است که در رشتههای مختلف مهندسی، از جمله مهندسی عمران، کاربردهای گستردهای دارد. از کاربردهای متلب در رشته مهندسی عمران میتوان به تحلیل ماتریسی قابهای سازهای (Frame Analysis) اشاره کرد که از روشهای متداول در مهندسی سازه است. این روش برای تحلیل نیروهای داخلی و تغییر مکانها در قابهای دوبعدی و سهبعدی استفاده میشود. تحلیل تیرهای تحت بارگذاری (Beam Analysis) نیز در مهندسی عمران، یکی از مسائل مهمی است که تحلیل رفتار تیرهای ساده یا مرکب تحت بارگذاریهای مختلف قابلیت پیادهسازی در متلب دارد. شبیهسازی جریان آب در کانال باز (Open Channel Flow) نیز از مسائل مهم در مهندسی هیدرولیک و مهندسی عمران است. تحلیل پایداری شیبها (Slope Stability Analysis) از مسائل مهندسی ژئوتکنیک است که به منظور جلوگیری از لغزش زمین یا مصالح ساختمانی انجام میشود. طراحی دال بتونی تحت بار متمرکز و یکنواخت در مهندسی سازه، نیز میتواند در نرم افزار متلب شبیهسازی شود. تحلیل تنش و کرنش در ستون بتونی، محاسبه عمق بحرانی جریان در سرریز، شبیهسازی انتقال حرارت در دیوارهای ساختمانی، تحلیل تغییر مکان پل تحت بارگذاری ترافیکی و طراحی سیستم زهکشی برای سازههای هیدرولیکی نیز قابلیت برنامهنویسی در نرم افزار متلب دارند.
همچنین تحلیل نفوذ آب در خاک (Seepage Analysis in Soil) یکی از مسائل مهم در مهندسی ژئوتکنیک و محیطزیست است. این مسئله به بررسی جریان آب در خاک و تأثیر آن بر پایداری سازهها، تحلیل پایداری شیبها، و طراحی سازههای هیدرولیکی مربوط میشود. یکی از مدلهای پایه برای تحلیل نفوذ آب، معادله دارسی است که بیان میکند جریان آب در خاک متناسب با گرادیان هیدرولیکی و نفوذپذیری خاک است که در نرمافزار متلب شبیهسازی میشود.
موسسه پروژستان آماده خدمات رسانی به دانشجویان عزیز رشته مهندسی عمران در تمامی مقاطع برای انجام پروژه متلب می باشد.
مثال 1: تحلیل ماتریسی قابها (Frame Analysis)
clc; clear;close all
% تعداد گرهها و اعضای قاب
n_nodes = 3;
n_elements = 2;
% مختصات گرهها
nodes = [0, 0; 3, 0; 6, 0];
% مشخصات اعضا (طول، مدول الاستیسیته و سطح مقطع)
E = 210e9; % مدول الاستیسیته (پاسکال)
A = 0.02; % سطح مقطع (متر مربع)
% تعریف ماتریس سختی
K = zeros(2*n_nodes);
% کد برای محاسبه ماتریس سختی برای هر عضو
for i = 1:n_elements
% محاسبه طول عضو
node1 = nodes(i, :);
node2 = nodes(i+1, :);
L = norm(node2 – node1);
% محاسبه ماتریس سختی عضو
k = E*A/L * [1 -1; -1 1];
% اضافه کردن به ماتریس سختی کلی
K(i:i+1, i:i+1) = K(i:i+1, i:i+1) + k;
end
% نمایش ماتریس سختی نهایی
disp(K);
خروجی:
1.0e+09 *
1.4000 -1.4000 0 0 0 0
-1.4000 2.8000 -1.4000 0 0 0
0 -1.4000 1.4000 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
توضیح- مرحله بندی کد به شرح زیر است:
- ابتدا تعداد گرهها و اعضا مشخص میشود.
- مختصات گرهها در یک ماتریس ذخیره میشود.
- مدول الاستیسیته (E) و سطح مقطع (A) تعریف میشوند.
- ماتریس سختی کل قاب به صورت یک ماتریس صفر آغاز میشود.
- طول هر عضو قاب با استفاده از فاصله اقلیدسی گرهها محاسبه شده و ماتریس سختی هر عضو تعیین و به ماتریس کلی اضافه میشود.
مثال 2: تحلیل تیرهای تحت بارگذاری (Beam Analysis)
clc; clear;close all
% طول تیر و بارگذاری
L = 5; % متر
q = 10; % کیلو نیوتون بر متر
% معادله برش و لنگر
x = linspace(0, L, 100);
V = q * (L – x); % معادله برش
M = q * x .* (L – x) / 2; % معادله لنگر
% رسم نمودار برش و لنگر
figure;
subplot(2,1,1);
plot(x, V);
title(‘نمودار برش’);
xlabel(‘مکان (متر)’);
ylabel(‘نیرو (کیلو نیوتون)’);
subplot(2,1,2);
plot(x, M);
title(‘نمودار لنگر’);
xlabel(‘مکان (متر)’);
ylabel(‘لنگر (کیلو نیوتون-متر)’);
خروجی:
توضیح:
- طول تیر (L) و مقدار بارگذاری یکنواخت (q) تعریف میشوند.
- معادلات برش و لنگر بر اساس روابط کلاسیک تیرها محاسبه میشوند.
- با استفاده از دستورplot نمودارهای برش و لنگر رسم میشوند.
مثال 3: شبیهسازی جریان آب در کانال باز (Open Channel Flow)
clc; clear;close all
% مشخصات کانال
b = 5; % عرض کانال (متر)
h = 2; % عمق آب (متر)
g = 9.81; % شتاب گرانش (متر بر ثانیه مربع)
% محاسبه عدد فرود
Fr = sqrt(g * h) / (h * b);
% نمایش نتیجه
if Fr < 1
disp(‘جریان زیر بحرانی است.’);
else
disp(‘جریان فوق بحرانی است.’);
end
خروجی:
جریان زیر بحرانی است.
توضیح:
- عرض و عمق کانال و شتاب گرانش تعریف میشوند.
- عدد فرود (Fr) محاسبه میشود تا نوع جریان (زیر بحرانی یا فوق بحرانی) تعیین شود.
- با استفاده از شرط if نوع جریان نمایش داده میشود.
مثال 4: تحلیل پایداری شیبها (Slope Stability Analysis)
clc; clear;close all
% زاویه شیب و ضریب اصطکاک
theta = 30; % درجه
phi = 35; % ضریب اصطکاک
% محاسبه ضریب اطمینان
FS = tan(deg2rad(phi)) / tan(deg2rad(theta));
% نمایش نتیجه
disp([‘ضریب اطمینان: ‘, num2str(FS)]);
خروجی:
ضریب اطمینان: 1.2128
توضیح:
- زاویه شیب (theta) و ضریب اصطکاک داخلی (phi) تعریف میشوند.
- ضریب اطمینان پایداری شیب (FS) محاسبه شده و نمایش داده میشود.
مثال 5: طراحی دال بتنی تحت بار متمرکز
clc; clear;close all
% مشخصات دال بتنی
b = 2; % عرض دال (متر)
h = 0.15; % ضخامت دال (متر)
P = 5000; % بار متمرکز (نیوتون)
% محاسبه تنش خمشی
M = P * b / 8;
sigma = M / (b * h^2 / 6);
% نمایش نتیجه
disp([‘تنش خمشی: ‘, num2str(sigma), ‘ پاسکال’]);
خروجی:
تنش خمشی: 166666.6667 پاسکال
توضیح:
- مشخصات دال بتنی و بار متمرکز تعریف میشوند.
- تنش خمشی در دال محاسبه شده و نمایش داده میشود.
مثال 6: تحلیل نفوذ آب در خاک (Seepage Analysis in Soil)
clc; clear;close all
% پارامترهای ورودی
k = 1e-6; % ضریب نفوذپذیری خاک (متر بر ثانیه)
h1 = 10; % ارتفاع آب در بالادست (متر)
h2 = 3; % ارتفاع آب در پایین دست (متر)
L = 20; % طول خاکریز (متر)
% محاسبه جریان نفوذ
Q = k * (h1 – h2) / L;
% نمایش نتیجه
disp([‘دبی نفوذ آب: ‘, num2str(Q), ‘ متر مکعب بر ثانیه’]);
خروجی:
دبی نفوذ آب: 3.5e-07 متر مکعب بر ثانیه
توضیح:
- تعریف پارامترها:
k = 1e-6 این خط ضریب نفوذپذیری خاک را تعریف میکند که نشاندهنده سرعت نفوذ آب در خاک است. واحد آن متر بر ثانیه است.
h1 = 10 ارتفاع آب در بالادست را تعریف میکند. این مقدار نشاندهنده سطح آب در سمت بالای خاکریز است.
h2 = 3 ارتفاع آب در پایین دست را تعریف میکند. این مقدار سطح آب در سمت پایینی خاکریز را نشان میدهد.
L = 20 طول خاکریز را مشخص میکند که فاصله بین دو نقطه (بالادست و پایین دست) است.
- محاسبه جریان نفوذ:
Q = k * (h1 – h2) / L در این خط، دبی نفوذ آب محاسبه میشود. بر اساس معادله دارسی، جریان آب از طریق خاک به صورت زیر محاسبه میشود:
(h1 – h2) این عبارت نشاندهنده اختلاف ارتفاع آب در بالادست و پایین دست است که به عنوان گرادیان هیدرولیکی در نظر گرفته میشود.
k ضریب نفوذپذیری خاک است که نشاندهنده قابلیت خاک در عبور آب میباشد.
L طول خاکریز است که جریان آب از آن عبور میکند. با تقسیم بر طول، دبی به دست میآید.
اگر علاقه مند به یادگیری نرم افزار متلب هستید، مطالعه مطالب زیر به شما پیشنهاد می شود:
