کد متلب محاسبه ی شبکه ی خطوط لوله

این تابع برای تحلیل شبکه لوله ای در حالت جریان پایدار (steady state) به کار می رود. ورودی ها : m: بردار دبی حجمی در هر لوله (واحد: m^3/s) – طول آن تعداد لوله ها را نشان می دهد.rho: چگالی سیال (واحد: kg/m^3)mu: ویسکوزیته دینامیکی سیال (واحد: Pa.s)D: بردار قطر لوله ها (واحد: m) – طول آن تعداد لوله ها را نشان می دهد.L: بردار طول لوله ها (واحد: m) – طول آن تعداد لوله ها را نشان می دهد.rf: بردار زبری سطح داخلی لوله ها (واحد: m) – طول آن تعداد لوله ها را نشان می دهد.خروجی : fm: ماتریس معادلات حاکم بر شبکه لوله ایعملکرد کد : محاسبه سرعت و عدد رینولدز : سرعت (v) در هر لوله با استفاده از دبی حجمی و قطر لوله محاسبه می شود.عدد رینولدز (Nre) برای هر لوله با استفاده از سرعت، قطر، چگالی و ویسکوزیته محاسبه می شود.محاسبه ضریب اصطکاک : پارامتر C با توجه به عدد رینولدز و زبری سطح لوله محاسبه می شود.پارامترهای A و B با استفاده از لگاریتم و عدد رینولدز محاسبه می شوند.نهایتا، ضریب اصطکاک (f) با در نظر گرفتن عدد رینولدز و پارامترهای A و B محاسبه می شود.محاسبه مشتقات : مشتقات پارامتر C نسبت به عدد رینولدز (d …