در مجاری باز، ورود هوا باعث افزایش در عمق آب می‌شود.
در سرریزهای شوت، هوای نزدیک کف می‌تواند باعث افزایش سرعت جریان شود.
در مجاری بسته، برای جریان عبوری از یک مقطع مشخص حضور هوا منجر به تغییر در دبی آب یا توزیع فشار سیستم می‌شود.

(( اینجا فقط تکه ای از متن درج شده است. برای خرید متن کامل فایل پایان نامه با فرمت ورد می توانید به سایت feko.ir مراجعه نمایید و کلمه کلیدی مورد نظرتان را جستجو نمایید. ))

حضور هوا بر کارایی سیستم هیدرولیکی تأثیر معکوس دارد.
در جریان­های هیدرولیکی واسطه­ای^[۲۴]، امواج فشاری به صورت قابل ملاحظه توسط هوا مستهلک می­شوند.
تجمع هوا در یک سیستم هیدرولیکی می‌تواند منجر به قطع جریان شود که به این اثرات blow-out یا blow-back گفته می‌شود.
حضور حباب‌های هوا منجر به تشدید انتقال اکسیژن و نیتروژن می‌شود تا آب به حالت اشباع کامل برسد.
ورود هوا به داخل جریان، بسته به شرایط می ­تواند مفید یا مضر باشد. از جمله مزایایی که ورود حباب هوا به داخل جریان دارد می‌توان به پیشگیری از وقوع کاویتاسیون، اکسیژن زایی و کاهش سرعت موج‌های فشاری اشاره کرد و از جمله معایب آن تأثیرات سوء آن بر عملکرد پمپ‌ها، سازه‌های آبگیر و در سیستم‌های انتقال آب تحت فشار می‌باشد. (Wood, 1991)

پیشینه موضوع تحقیق

در پی خرابی‌های ناشی از کاویتاسیون مطالعات آزمایشگاهی و عددی بسیار گسترده‌ای در مورد راه‌های پیشگیری و جلوگیری از صدمات ناشی از کاویتاسیون به سازه‌های هیدرولیکی انجام گرفته که از جمله این تحقیقات می­توان به موارد زیر اشاره کرد:
پترکا^[۲۵] (۱۹۵۳) تحقیقات گسترده­ی آزمایشگاهی جهت تعیین نقش هوادهی به جریان را در کاهش خسارات ناشی از کاویتاسیون انجام داد. این آزمایشات در داخل یک ونتوری در طول دو ساعت و با سرعت جریان در گلوگاه ونتوری بیش از m/s 30 ، در دو حالت با تزریق و بدون تزریق هوا روی نمونه­های بتنی صورت گرفت. به ازای غلظت هوای ۲% افت وزنی نمونه­های بتنی در اثر کاویتاسیون بسیار کاهش یافته و برای غلظت­های ۶% تا ۸% این افت قابل اندازه ­گیری نبود. بر اساس نتایج این آزمایشات مقدار غلظت هوای ۵% تا ۸% برای پیشگیری از کاویتاسیون و کاهش خسارات آن در سطوح بتنی توصیه شد. امروزه نیز در بسیاری از تحقیقات بر روی این موضوع، اغلب مقادیر توصیه شده توسط پترکا ارائه شده و به عنوان معیار هوادهی در نظر گرفته می­ شود. (شکل ۲-۶) (Wood, 1991)
راسموسن^[۲۶] (۱۹۵۶) آزمایشاتی بر روی جریان دو فازه (مخلوط آب و هوا) برای بررسی تاثیر آن در پیشگیری از پدیده کاویتاسیون انجام داده است. در این تحقیق آزمایشگاهی دو روش، یکی دیسک چرخان با سرعت بالا در داخل آب و دیگری جریان مشابه آزمایش پترکا مورد استفاده قرار گرفت. نتایج آزمایش نشان داد که در هر دو نمونه مورد آزمایش مقدار هوای ۸% تا ۱۰% برای از بین بردن کامل پدیده کاویتاسیون مناسب می­باشد. (Krammer, 2004)
شکل ‏۲‑۶- تاثیر ورود هوا در کاهش تخریب نمونه بتنی (Wood, 1991)
روزانوف و کاوشینکف^[۲۷] (۱۹۷۳) مطالعات وسیعی در مورد خسارات ناشی از کاویتاسیون بر روی مستهلک کننده­ های انرژی در حوضچه­های آرامش انجام دادند. اندازه ­گیری­های روی مدل و نمونه اصلی نشان داد که وقوع کاویتاسیون در اثر نوسانات فشار حتی اگر میانگین فشار بزرگتر از فشار بخار مایع باشد، امکان­ پذیر خواهد بود. در نتیجه این آزمایشات مشخص شد که تزریق مقدار هوای ۲۵/۱% تا ۲% به منطقه با فشار کم (ناحیه کاویتاسیون)، خسارت ناشی از کاویتاسیون را تا ۵/۲% کاهش داده و برای مقدار هوای ۶% تا ۷% بلوک­های مستهلک کننده متحمل هیچ­گونه خسارتی نخواهند شد. (Krammer, 2004)
راسل و شیهان^[۲۸] (۱۹۷۴) تحقیقات آزمایشگاهی خود را بر روی نمونه­های بتنی با مقاومت فشاری مختلف جهت تعیین درصد هوای مناسب برای حذف خسارت کاویتاسیون انجام دادند. در این آزمایشات مقاومت نمونه بتنی از Mpa 10 تا Mpa 20 در نظر گرفته شد و تاثیر هوادهی بر جریان مورد بررسی قرار گرفت. نتایج آزمایشات نشان داد که وجود ۸/۲% هوا جهت حذف خسارات ناشی از کاویتاسیون در نمونه­های مورد آزمایش کافی می­باشد. با این­حال آن­ها بیان کردند که مقدار هوای واقعی برای جلوگیری از خسارت پدیده کاویتاسیون در شرایط مختلف به پارامترهای گوناگونی وابسته می­باشد. همچنین قابل ذکر است که در این آزمایشات تنها غلظت هوای متوسط برای کل حجم محاسبه شده و مقدار غلظت موضعی هوا برای منطقه خسارت دیده مشخص نشده است. (Krammer, 2004)
کین^[۲۹] (۱۹۷۸) مطالعاتی بر روی یک مدل واقعی سرریز جهت یافتن تأثیر هوادهی در جلوگیری از تخریب ناشی از کاویتاسیون انجام داد. این تحقیقات بر اساس نظریه‌ی کاهش خسارات ناشی از کاویتاسیون با هوادهی به جریان به علت کاهش سرعت صوت در مخلوط آب و هوا صورت پذیرفت. همانطور که در شکل (۲-۷) مشهود است با افزایش غلظت هوای جریان از ۰% تا ۱۰%، سرعت صوت با شدت زیادی از m/s1460 به m/s50 کاهش یافته است. نتایج بدست آمده از این آزمایشات مشخص کرد که اگر غلظت هوا در نزدیک بستر بیش از ۱% تا ۲% باشد تخریب ناشی از کاویتاسیون کاهش‌یافته و به ازای غلظت هوای ۵% تا ۱۰% تخریبی رخ نمی‌دهد. شکل (۲-۷) (Wood, 1991)
شکل ‏۲‑۷- تاثیر هوادهی به جریان در کاهش سرعت صوت (Wood, 1991)
سمنکوو و اسکلکوو^[۳۰] (۱۹۸۰) بررسی­های تئوریکی و آزمایشگاهی روی هوادهی جریان، انواع هواده‌ها و تاثیر این سازه­ها بر روی جریان انجام دادند. در این آزمایشات ویژگی‌های هندسی مانند ارتفاع و زاویه هواده برای انواع هواده‌ها مثل رمپ، پله و ترکیب آن‌ها مورد بررسی قرار گرفت. نتایج این تحقیقات مشخص نمود که هوادهی به جریان اقتصادی‌ترین راه مقابله با تخریب‌های ناشی از کاویتاسیون می‌باشد ولی انتخاب نوع هواده بسته به شرایط به پارامترهای مختلفی وابسته می‌باشد.
اروین^[۳۱] و همکاران (۱۹۸۰) در طی آزمایشاتی به بررسی تاثیر سرعت و مقدار هوای ورودی توسط جت برخوردی دایروی به سطح آب بر شدت توربولانسی پرداختند. در طی این آزمایشات مشخص شد که ورود هوا توسط جریان جت به آب باعث کاهش شدت توربولانسی شده و با افزایش مقدار هوای ورودی به جریان از شدت توربولانسی کاسته می­ شود. (شکل ۲-۸) (Wood, 1991)
شکل ‏۲‑۸- تغییرات شدت توربولانسی بر اساس سرعت هوای ورودی به جریان (Ervin et al., 1980)
کن و گارود^[۳۲] (۱۹۸۱) مطالعاتی در مورد علل تخریب سرریز تونلی سد تاربلا پاکستان انجام داده­اند. در نتیجه این مطالعات نشان داده شد که خسارات رخ داده در جداره بتنی تونل این سد ناشی از بروز پدیده کاویتاسیون در اثر جریان­های برشی شدید حاصله به همراه سرعت بالای m/s 30 جریان و نوسانات فشار بوده است. تونل انتقال آب سد تاربلا دارای سه دریچه در دهانه ورودی بوده که به هنگام وقوع حادثه دریچه میانی باز و دو دریچه کناری بسته بودند. با نوسانات فشار و همچنین برقراری جریان برشی شدید در مجاورت پایه­ های ورودی تونل و تشکیل جریان­های گردابی، حباب­های کاویتاسیون تولید شده و به همراه جریان به مناطق پایین­دست (دارای فشار بالا) حمل شده و در آن منطقه منفجر شدند و در نتیجه، طی گذشت زمان پوشش بتنی جداره داخلی تونل دچار خوردگی شدید شده است.
چانسون^[۳۳] (۱۹۸۹) تحقیقات آزمایشگاهی خود را جهت تعیین غلظت­های هوای ورودی به جریان در اثر دو نوع هواده انجام داد. در این آزمایشات که روی مدل هیدرولیکی سرریز آزاد سد کلاید با مقیاس ۱:۱۵ انجام گرفت، اعداد فرود جریان در مدل بین ۳ تا ۲۵ و سرعت ورودی m/s 3 تا m/s 14 بوده است. دو نوع هواده اولی شامل ترکیبی از رمپ با زاویه ۷/۵ درجه و پله به ارتفاع mm 30 و دومی تنها شامل پله­ای به ارتفاع mm 30 به طور جداگانه در مدل تعبیه و وضعیت جریان در پایین­دست هواده­ها و جریان هوادهی مورد بررسی قرار گرفتند. در نتیجه این آزمایشات در منطقه با جریان متغیر تدریجی، تحلیل جریان بر اساس معادلات پیوستگی و انرژی منجر به بدست آمدن دو معادله دیفرانسیلی شده که می­توانند برای تعیین ورود هوا به داخل جریان سرریز آزاد در روش­های عددی به کار گرفته شوند. همچنین پیش ­بینی غلظت هوای جریان توسط این روش بسیار قابل اعتماد بوده و می ­تواند در محاسبه فاصله بین هواده­ها مورد استفاده قرار گیرد.
چانسون (۱۹۸۹) با انجام مطالعاتی روند ورود هوا به داخل جریان توسط هواده­ها را مورد بررسی قرار داد. در این تحقیقات که بر روی مدل یک سرریز شوت با سرعت بالا (m/s 4 تا m/s 15) صورت پذیرفت، نحوه­ ورود هوا به داخل جریان به طور مفصل مورد بحث قرار گرفت و روابط جدیدی برای تعیین غلظت هوای ورودی به جریان و اندازه ­گیری سرعت جریان ارائه شد. نتایج این مطالعات مشخص کرد که مقدار هوای ورودی به جریان وابسته به اختلاف فشار داخل فضای کاویتی و خارج آن است و همچنین در مورد کاویتی پایین­دست هواده و موقعیت نقطه برخورد روابطی ارائه شد.
وود^[۳۴] و همکاران (۱۹۹۱) مطالعات گسترده­ای در رابطه با هوادهی سطحی و اجباری به جریان انجام دادند. در این مطالعات تاثیر ورود هوا بر میدان­های فشار و سرعت و همچنین توربولانس جریان مورد بررسی قرار گرفت. نتایج این تحقیق نشان داد که ورود هوا به داخل جریان باعث افزایش سرعت جریان می­ شود. تغییرات فشار جریان در کف کانال به صورت کیفی مشابه شکل (۲-۸) بوده که در اثر ورود هوا به داخل جریان، فشار موضعی خصوصا در ناحیه چرخش که دارای فشار منفی است، افزایش می­یابد.
شکل ‏۲‑۹- نواحی مختلف جریان و دیاگرام توزیع فشار در کف (Wood, 1991)
اروین و همکاران (۱۹۹۵) آزمایشاتی بر روی هواده نصب شده بر روی سرریزها به منظور بررسی تأثیر هندسه‌ی رمپ هواده بر وضعیت هوادهی جریان و همچنین وضعیت توربولانسی جریان انجام دادند. این آزمایشات روی یک سرریز آزاد با شیب طولی بستر ۴۵ درجه و به ازاء سه زاویه‌ی مختلف برای رمپ هواده، انجام شد. نتایج این آزمایشات نشان می‌دهد که در حالت کلی وجود رمپ هواده معمولا ۲۰ تا ۳۰ درصد اغتشاش اضافی در جریان به وجود می‌آورد و همچنین رمپ‌های با زاویه‌ی بزرگ‌تر اغتشاش بیشتری نسبت به رمپ‌های با زاویه‌ی کوچک‌تر ایجاد می‌کنند. (شکل ۲-۹)
شکل ‏۲‑۱۰- تغییرات شدت توربولانسی در طول جت جریان پایین­دست رمپ هواده سرریز آزاد (Ervin et al., 1995)
دانگ^[۳۵] و همکاران (۲۰۰۵) با انجام تحقیقات آزمایشگاهی به مقایسه‌ی وضعیت جریان در دو حالت با هوادهی و بدون هوادهی و تأثیر هوادهی در جلوگیری از خرابی ناشی از کاویتاسیون پرداختند. آزمایشات در یک تونل آب که شامل چهار مقطع متوالی که به ترتیب مقطع هوادهی، مقطع جمع شونده، مقطع مشاهده و مقطع پخش شونده نامگذاری شده بود، انجام گرفت. با ورود هوا به داخل جریان، وضعیت جریان با توجه به پارامترهای تغییرات فشار و اندیس کاویتاسیون در مقطع مشاهده مورد بررسی قرار گرفت. نتیجه آزمایشات حاکی از این است که هوای وارد شده به جریان، فشار در ناحیه‌ی کاویتاسیون یا همان مقطع مشاهده را به طور قابل‌ملاحظه‌ای افزایش می‌دهد و همچنین با رسم نمودار موج فشاری ناشی از انفجار حباب‌های کاویتاسیون در دو حالت با هوادهی و بدون هوادهی جریان نشان داده شد که این نمودار شکلی به صورت پله به خود می‌گیرد که ناشی از کاهش قدرت موج فشاری در اثر ورود هوا می‌باشد.
چانسون و همکاران (۲۰۰۶) آزمایشات گسترده­ای جهت تعیین تاثیر حضور پایاب بر مشخصه­های هیدرولیکی جریان سرریز انجام دادند. با مدل­سازی عددی مدل آزمایشگاهی، مطالعات بر روی هر دو مدل صورت گرفت. بر روی مدل آزمایشگاهی سرریز ۵ حالت و بر روی مدل عددی ۴ حالت برای برای عبور جریان با تغییر در هر ورودی به سرریز و عمق­های متفاوت پایاب صورت گرفت. مقادیر حاصله از مدل­های هیدرولیکی و عددی برای مشخصه­های جریان مانند دبی جریان، فشارهای نسبی و فشارهای مطلق با یکدیگر مقایسه شدند. مقایسه نتایج دبی جریان در مدل هیدرولیکی و عددی بیانگر نتایج مورد قبول مدل عددی در پیش ­بینی وضعیت جریان می­باشد. همچنین فشارهای حاصله از مدل عددی نیز بیانگر توانایی قابل قبول این مدل در سرریزهای دارای پایاب می­باشد. اگرچه مدل عددی خود دارای محدودیت­هایی می­باشد اما قابلیت این را دارد که در نقاطی که امکان برداشت نتایج در مدل هیدرولیکی وجود ندارد تغییرات و گرادیان پارامترها را به دقت محاسبه کند.
شمسایی و سلیمان­زاده (۲۰۰۶) در طی مطالعات عددی به شبیه­سازی جریان آب و هوا در مدل آزمایشگاهی یک تخلیه کننده تحتانی پرداختند. در این تحقیقات از نتایج مدل آزمایشگاهی اسپیرلی و هگر (Speerli and Hager, 2000) استفاده شده و شبیه­سازی عددی جریان دو فازی آب و هوا برای آن صورت گرفت. پارامترهای مختلف حاصله از مدل عددی برای جریان با نتایج مدل آزمایشگاهی مقایسه شدند. نتایج این شبیه­سازی نشان داد که توانایی مدل عددی در تخمین پارامترهای مشخص کننده وضعیت جریان مانند دبی جریان، اندیس کاویتاسیون، غلظت هوا و توزیع هوا در مقاطع مختلف قابل قبول بوده به طوری که مدل عددی این پارامترها را با اختلاف ناچیزی نسبت به نتایج آزمایشگاهی تخمین زده است.
دانگ (۲۰۰۶) تحقیقات آزمایشگاهی خود را در مورد کنترل کاویتاسیون به وسیله هوادهی به جریان در سریزهای آزاد انجام داد. این آزمایشات در سرعت­های متفاوت m/s 20 تا m/s 50 در دو حالت با هوادهی و بدون هوادهی صورت گرفته و در طی آن مقادیر موج فشاری در هر دو حالت اندازه ­گیری شد. نتایج این آزمایشات مشخص نمود که نمودار موج فشاری در حالت با هوادهی و بدون هوادهی در ناحیه کاویتاسیون، دارای پله است و هوادهی به جریان به طور قابل توجهی فشار در این ناحیه را افزایش می­دهد. همچنین در این تحقیقات رابطه­ای بین حداقل غلظت هوا جهت پیشگیری از کاویتاسیون و سرعت جریان ارائه شد.
درگاهی (۲۰۰۶) به صورت سه بعدی و توسط نرم­افزار Fluent به شبیه سازی جریان عبوری از روی سرریز اوجی و تعیین سطح آزاد جریان و پروفیل­های سرعت در ازای استفاده از مدل­های توربولانسی و توابع دیواره مختلف، پرداخته است. پارامترهای اصلی جریان دارای تطابق خوبی با داده ­های آزمایشگاهی هستند با این حال این تطابق وابستگی فراوانی به انتخاب اندازه شبکه، مدل توربولانسی، تابع دیواره و عدد رینولدز دارد. نتایج این تحقیقات ثابت کرد که در سرریزها استفاده از شبکه ­های به نسبت کوچکتر در نزدیکی دیواره از اهمیت بیشتری برخوردار است. مدل توربولانسی RNG k-ε نیز نتایج بهتری نسبت به مدل­های دیگر برای میدان سرعت در نزدیکی دیواره­ها ارائه می­دهد با این حال انتخاب تابع دیواره مناسب برای مرزهای صلب از اهمیت بیشتری در مقایسه با نوع مدل توربولانسی برخوردار است. در کل تعیین سطح آزاد آب با اختلاف ۵/۱% تا ۹/۲% نسبت به داده ­های آزمایشگاهی در هدهای مختلف سرریز به خوبی امکان­ پذیر بود.
روآن^[۳۶] و همکاران (۲۰۰۷) با انجام تحقیقات آزمایشگاهی به مقایسه‌ی عملکرد دو نوع هواده دارای رمپ و همچنین رمپ به همراه گروو در دو طرح مختلف در یک سرریز تونلی با توجه به پارامترهای طول کاویتی، دبی هوای ورودی و غلظت هوا در درون جریان آب پرداختند. در شکل (۲-۱۱) محل قرارگیری هواده‌ها و در شکل (۲-۱۲) طرح هواده‌های به کار رفته نمایش داده شده است.
شکل ‏۲‑۱۱- سرریز تونلی همراه با هواده (Ruan et al., 2007)

شکل ‏۲‑۱۲- طرح شماره ۱ (سمت راست) و طرح شماره ۲ (سمت چپ) (Ruan et al., 2007)
طرح شماره ۱ شامل دو هواده رمپ به همراه گروو و طرح شماره ۲ شامل دو هواده رمپ می‌باشد. نتایج آزمایشات نشان داد که با توجه به پارامترهای جریان مثل طول کاویتی، دبی هوای ورودی به جریان آب و غلظت هوای وارد شده به جریان در پایین‌دست هواده، هر دو هواده در طرح شماره ۲ عملکرد بهتری نسبت به طرح شماره ۱ داشته و همچنین مشخص شد که برای سرریزهای تونلی، شکل و اندازه‌ی هواده‌ها نیازمند طراحی خاص می‌باشد زیرا در پروژه‌های مختلف، شرایط مختلفی برای دبی جریان، سرعت، شیب کف و طول سرریز تونلی وجود دارد. طرح به کار رفته در این آزمایش (طرح شماره ۲) می‌تواند با افزایش فضای کاویتی، طول کاویتاسیون، دبی هوای ورودی و غلظت هوا در نزدیکی کف، از سرریزهای تونلی در برابر کاویتاسیون محافظت نماید.
وو و روآن^[۳۷] (۲۰۰۷) در طی مطالعات خود با ارائه معادلات جدید و حل آن­ها به بررسی زاویه جریان عبوری از روی هواده­ها پرداختند. آن­ها معادلات خود را در رابطه با طول کاویتی پایین­دست هواده و زاویه جریان عبوری از روی هواده برای دو نوع هواده دارای رمپ و پله ارائه دادند. در این معادلات دو پارامتر موثر در زاویه جریان عبوری از روی هواده، یعنی عمق و مؤلفه­ های سرعت جریان، مورد بررسی قرار گرفته و نتایج حاصله با روش­های محققین دیگر در تخمین زاویه جریان و نتایج مدل آزمایشگاهی مقایسه شد. در نتیجه این مطالعات مشخص شد که در هواده رمپ، زاویه جریان عبوری از روی هواده وابستگی زیادی به عمق جریان روی رمپ و مؤلفه­ های سرعت دارد اما در هواده پله این زاویه به عمق جریان وابسته نبوده و فقط تابع مؤلفه­ های سرعت جریان می­باشد.
سو^[۳۸] و همکاران (۲۰۰۹) آزمایشاتی بر روی یک نوع جدید هواده جهت بررسی عملکرد آن در جریان­های با عدد فرود پایین و شیب کف کانال ملایم انجام دادند. این آزمایشات در یک تونل آب بر روی یک هواده جدید که هواده ذوزنقه­ای شیاردار با مقطع شیبدار نامگذاری شده بود، صورت گرفت. در طی این آزمایشات مشخصات جریان هوادهی شده توسط این هواده مانند الگوی جریان، شکل کاویتی، غلظت هوا در نزدیکی کف، توزیع فشار در کف تونل و سرعت هوا مورد بررسی قرار گرفت. نتایج آزمایشات نشان داد که این هواده به علت ایجاد سطح تماس بزرگ بین آب و هوا، قادر است حجم هوای ورودی به جریان را افزایش داده و در نتیجه شرایط هوادهی را بهبود ببخشد. همچنین این هواده توانایی پایدار سازی کاویتی پایین­دست هواده را داشته و طول آن را افزایش می­دهد.
محقق و وو (۲۰۰۹) با انجام آزمایشاتی به بررسی تاثیر پارامترهای هیدرولیکی و هندسی بر روی طول کاویتی پایین­دست دریچه سرویس تونل تخلیه پرداختند. در این آزمایشات رابطه L/d=f(Fr, t_s/d, , α) ارائه شد که حاکی از وابستگی طول کاویتی به عدد فرود جریان، ارتفاع پله، شیب کف تونل و زاویه جداره فوقانی تونل قبل از دریچه سرویس می­باشد. در جدول (۲-۲) مشخصات هندسی تونل در حالت­های مختلف آزمایش آمده است. در نتیجه این آزمایشات مشخص شد که افزایش عدد فرود جریان، ارتفاع پله، شیب کف و زاویه جداره فوقانی باعث افزایش طول کاویتی پایین­دست می­شوند که نقش شیب کف تونل برجسته­تر می­باشد. (شکل ۲-۱۳)
جدول ‏۲‑۲- مشخصات هندسی تونل (Mohaghegh and wu, 2009)

(a) (b)
©
شکل ‏۲‑۱۳- تاثیر پارامترهای هندسی (a) زاویه جداره فوقانی تونل، (b) شیب کف تونل، © ارتفاع پله
بر طول کاویتی پایین­دست (Mohaghegh and Wu, 2009)
فیشر و هگر^[۳۹] (۲۰۱۰) با تحقیقات آزمایشگاهی وسیعی روی مدل هیدرولیکی سرریز آزاد به تعیین خصوصیات جریان هوادهی شده از جمله غلظت هوا در نواحی نزدیک هواده و دورتر از آن پرداختند. این آزمایشات برای دامنه گسترده­ای از هندسه­های مختلف هواده و پارامترهای هیدرولیکی صورت گرفت. در این تحقیق شرایط جریان بعد از هواده به سه منطقه جت پرتاب، ناحیه برخورد جریان به کف و منطقه دور از هواده تقسیم شده است. در نتیجه این آزمایشات مقادیر غلظت­های هوای ورودی از سطح جریان، غلظت هوای ورودی از کاویتی و مقدار میانگین غلظت در مقطع عرضی برای هر سه منطقه نامبرده ارائه شد. مقادیر غلظت هوای جریان در نزدیکی کف دارای حداکثر مقدار خود بوده و تا نقطه برخورد جت ادامه داشته و از نقطه برخورد جت این غلظت کاهش می­یابد. همچنین گرادیان هوادهی بالایی نیز در محل برخورد جت به کف سرریز مشاهده گشت.
ژنگ^[۴۰] و همکاران (۲۰۱۱) در طی تحقیقات گسترده‌ی آزمایشگاهی و عددی با شبیه‌سازی عددی سه بعدی جریان هوادهی شده بعد از هواده در یک تونل، به بررسی و مقایسه‌‌ی روش­های VOF و Mixture پرداختند. در این تحقیقات با تحلیل عددی جریان در نزدیکی هواده‌، مقایسه‌ی نتایج بدست آمده با نتایج آزمایشگاهی صورت گرفت که در آن وضعیت جریان با توجه به پارامترهای جریان مانند طول کاویتی و مقادیر متوسط فشار در پایین‌دست هواده، به ازائ هر دو روش عددی مورد بررسی قرار گرفت. نتایج این مطالعات نشان داد که هر دو روش با دقت خوبی سطح آزاد آب و طول کاویتی را محاسبه می‌کنند. غلظت هوای درون آب و مقادیر متوسط فشار پایین‌دست هواده بدست آمده از مدل توربولانسی Mixture خطای کمتری نسبت به نتایج آزمایشگاهی دارد.
فیشر (۲۰۱۱) با انجام تحقیقات آزمایشگاهی به تعیین مشخصه­های هیدرولیکی جریان دو فازی در سرریز آزاد دارای سازه هواده پرداخته است. در این آزمایشات تغییرات مشخصه­های هیدرولیکی جریان دو فازی از جمله ضریب هوادهی و طول کاویتی نسبت به دو پارامتر مهم زاویه دفلکتور و مقدار زیرفشار داخل کاویتی مورد بررسی قرار گرفت. نتایج تحقیق نشان می­دهد که دفلکتورهای با زاویه بزرگتر عملکرد بهتری در وارد کردن درصد هوای بیشتر به جریان دارند و همچنین افزایش زاویه دفلکتور (α) باعث افزایش طول نسبی کاویتی (L/h₀) می­ شود. (شکل ۲-۱۴) با افزایش زیرفشار داخل کاویتی شکل گرفته بلافاصله در پایین­دست رمپ، طول کاویتی و ضریب هوادهی کاهش می­یابند. (شکل ۲-۱۵) در این شکل­ها h₀ ارتفاع جریان روی رمپ، L_M ماکزیمم طول کاویتی و P زیرفشار داخل کاویتی می­باشد.
وو و همکاران (۲۰۱۱) تاثیر پارامترهای هندسی هواده را بر روی مقادیر آب برگشتی به داخل کاویتی تشکیل شده بلافاصله در پایین­دست رمپ هواده و طول کاویتی را به صورت آزمایشگاهی در سرریز آزاد مورد مطالعه قرار دادند. نتایج تحقیق مشخص نمود که در اعداد فرود پایین و شیب طولی کم شوت، آب برگشتی ناشی از برخورد جت عبوری از روی رمپ هواده به کف شوت می ­تواند کاویتی پایین­دست هواده را پر نموده و موجب بروز کاویتاسیون گردد. نتایج این تحقیق مشخص نمود که مقدار آب برگشتی به داخل کاویتی به طول کاویتی تشکیل شده در پایین­دست رمپ وابسته بوده و برای هواده­های با ارتفاع بزرگتر یا زاویه رمپ بیشتر مقدار طول کاویتی بزرگتر می­باشد.

FCI

۶۷/۰

۲۱/۲

۰۲/۰

Constant

۷۹/۲۰-

۹۸/۲-

۰۰۳/۰

F- لیمر

۰۰۰/۰=, P>F28/12 = (164،۶)F

آزمون هاسمن

۰۰۰/۰=۴۵/۴۸=(۶) =H

منبع: یافته­های پژوهش
همان­طور که در جدول شماره­ (۴-۴) مشاهده می­ شود، ضریب متغیر سرمایه ­گذاری مستقیم خارجی منفی و از نظر آماری معنادار است. به نظر می­رسد، طبق تئوری تجارت، سرمایه ­گذاری مستقیم خارجی در کشور­های در­حال توسعه میزبان، منجر به بهبود تخصیص منابع و از این طریق، افزایش توان اشتغال­زایی و کاهش بیکاری شده است. شاید دلیل این موضوع آن باشد که آثار سرمایه ­گذاری مستقیم خارجی بر اشتغال و بیکاری کشور­های مختلف، به میزان توسعه ­یافتگی و آماده بودن بستر­های مناسب در این کشور­ها بستگی دارد. بنابراین می­توان بیان کرد، منفی بودن این ضریب نشان می­دهد که سرمایه ­گذاری مستقیم خارجی در کشور­های در­حال توسعه منتخب منجر به تسهیل دسترسی بنگاه­ها به سرمایه و کاهش هزینه­ آن­ها شده است. بنابراین فعالیت بنگاه­ها افزایش یافته و با رشد تولید، تقاضا برای نیروی کار افزایش و بیکاری کاهش یافته است، به گونه ­ای که یک درصد افزایش در FDI،بیکاری را به میزان ۲۲/۰ درصد کاهش داده است.

( اینجا فقط تکه ای از متن پایان نامه درج شده است. برای خرید متن کامل فایل پایان نامه با فرمت ورد می توانید به سایت feko.ir مراجعه نمایید و کلمه کلیدی مورد نظرتان را جستجو نمایید. )

درجه­ باز بودن تجاری (نسبت مجموع صادرات و واردات به GDP)، بیکاری را در کشور­های منتخب کاهش داده است و برآورد می­ شود که یک درصد افزایش در درجه­ باز بودن تجاری، بیکاری را به میزان ۵۷/۸ درصد با کاهش رو به رو نموده است. ضریب این متغیر برابر با ۵۷/۸- برآورد شده و در سطح بالایی معنادار است. این مقدار نشان می­دهد که درجه­ باز تجاری در کشور­های منتخب عضو منا توانسته است از طریق افزایش صادرات به ایجاد یک بازار داخلی قدرتمند و سپس ایجاد یک محیط رقابتی بینجامد. این شرایط به بهبود رشد اقتصادی این کشور­ها کمک کرده و افزایش رشد اقتصادی موجب ایجاد و افزایش فرصت­های شغلی جدید و کاهش بیکاری گردیده است. از سوی دیگر، افزایش صادرات باعث افزایش سطح تولید این کشور­ها شده و از آن­جایی که تولید بیشتر، نیازمند استخدام نیروی کار بیشتر است به کاهش بیکاری در این کشور­ها کمک نموده است.
با توجه به نتایج، اثر نرخ ارز واقعی بر بیکاری منفی و معنادار است و یک درصد افزایش در نرخ ارز واقعی، بیکاری را ۰۸۵/۰ درصد کاهش می­دهد. تأثیر منفی نرخ ارز واقعی بر بیکاری را می­توان این­گونه بیان کرد؛ با افزایش نرخ ارز واقعی (کاهش ارزش پول ملی)، تولید کالا­های صادراتی مقرون به صرفه می­ شود و تولیدکنندگان کالا­های صادراتی از این راه، سود سرشاری به دست می­آورند و به منظور حفظ کمیت و کیفیت تولید، مبادرت به افزایش سرمایه ­گذاری می­ کنند که این رشد سرمایه ­گذاری پاسخگوی رشد بالای نرخ بیکاری می­باشد. افزایش سرمایه ­گذاری تولیدکنندگان باعث افزایش تقاضا برای نیروی کار جدید می­ شود، بنابراین با افزایش نرخ واقعی ارز، نرخ بیکاری کاهش می­یابد.
نتایج برآورد معادله، بیانگر اثر منفی و معنادار رشد تولید ناخالص داخلی بر بیکاری است. طبق جدول (۴-۴)، ضریب برآورد شده این متغیر برابر با ۰۵/۰- است که انتظار می­رود افزایش یک درصدی رشد تولید ناخالص داخلی، کاهشی ۰۵/۰ درصدی در بیکاری به همراه داشته باشد. طبق قانون اوکان، رشد ۳ درصدی تولید ناخالص داخلی منجر به یک درصد کاهش بیکاری می­ شود و علامت ضریب متغیر برآورد شده با این قانون مطابقت دارد و افزایش تولید ناخالص داخلی در کشور­های
در­حال توسعه منتخب، کاهش بیکاری را به همراه داشته است. علت این تأثیر منفی را می­توان این­گونه توجیه کرد؛ در کشور­های مورد مطالعه نرخ رشد بهره­وری و نرخ رشد نیروی کار از نرخ رشد تولید ناخالص داخلی پایین­تر بوده است و همین امر باعث شده که با رشد GDP، نرخ بیکاری کاهش یابد.
نرخ مشارکت نیروی کار، دارای تأثیر مثبت بر نرخ بیکاری بوده و بیکاری را در کشور­های مورد مطالعه افزایش داده است؛ به گونه ­ای که افزایش یک درصدی این نرخ به افزایش ۷۶/۰ درصدی در نرخ بیکاری منجر شده است. اگرچه این تأثیر مثبت و معنادار، برخلاف ادبیات و تئوری است اما شاید بتوان علت این امر را این­گونه توجیه کرد:
رشد بی­رویه­ی جمعیت در کشور­های در­حال توسعه منجر به افزایش کل جمعیت و هم­چنین جمعیت فعال در این کشور­ها گردیده است. از سوی دیگر با افزایش جمعیت، تعداد افراد شاغل و در
جست و جوی شغل افزایش یافته که منجر به افزایش نرخ مشارکت نیروی کار در این کشور­ها شده است. کمبود تقاضای نیروی کار در مقایسه با عرضه­ی رو به رشد نیروی کار در این کشور­ها که اغلب از کاهش سطح تقاضای کل، بیشتر بر اثر سرمایه ­گذاری ناکافی یا به کارگیری تکنولوژی سرمایه­بر منشأ می­گیرد به این امر کمک نموده که با وجود افزایش نرخ مشارکت نیروی کار، بیکاری افزایش یابد.
ضریب متغیر شاخص بحران مالی برابر با ۶۷/۰ برآورد شده است که تأثیر مثبت و معناداری بر بیکاری دارد. به عبارتی افزایش شاخص بحران مالی باعث افزایش بیکاری شده که این تأثیر را
می­توان این­طورعنوان کرد؛ در اثر بحران مالی، عواید و درآمد­های صادراتی کشور­های مورد مطالعه
به­ ویژه صادرات مواد خام و اولیه مانند نفت و منابع طبیعی دیگر کاهش یافته و از آن­جایی که در این کشور­ها، صادرات و درآمد ارزی حاصل از آن به عنوان یکی از مهم­ترین منابع رشد اقتصادی
محسوب می­ شود، رشد اقتصادی را کاهش داده که نتیجه­ آن بر اساس قانون اوکان، افزایش بیکاری بوده است.
هم­چنین کاهش تقاضا برای کالا­ها و محصولات تولید شده در کشور­های مورد مطالعه از سوی
کشور­های درگیر بحران، کاهش تولید را به همراه داشته است که نتیجه­ آن کاهش تقاضای نیروی کار، کاهش تعداد ساعات کاری کارگران و بیکاری نیروی کار بوده است. از سوی دیگر، بحران مالی اخیر به دلیل طولانی و عمیق بودن، تعداد نیروی کار بالقوه را از طریق نا­امید کردن کارگران جویای کار کاهش داده که این مسئله افزایش نرخ بیکاری را به همراه داشته است.
۴-۳- تأثیر بحران مالی بر بیکاری کشورهایی با گروه ­های درآمدی متفاوت
در این بخش، با بهره گرفتن از سطوح درآمدی کشورها به بررسی تأثیر بحران مالی بر نرخ بیکاری پرداخته می­ شود.بر اساس طبقه بندی­های بانک جهانی، از ده کشور مورد مطالعه، شش کشور جزء کشورهایی با درآمد بالاتر از متوسط و درآمد بالا هستند که به دلیل محدود بودن جامعه­ آماری مورد مطالعه، این کشورها به عنوان کشورهای در­حال توسعه با درآمد بالا در نظر گرفته شده ­اند و شامل: الجزایر، ایران، اردن،کویت، عربستان سعودی و تونس می­باشند.چهار کشور مصر، مراکش، سودان و سوریه نیز طبق طبقه ­بندی بانک جهانی، کشورهایی با در آمد پایین هستند که در این بخش، تأثیر بحران مالی بر نرخ بیکاری کشورهای مذکور بررسی می­ شود.
۴-۳-۱- تجزیه و تحلیل تأثیر بحران مالی بر نرخ بیکاری کشورهایی با در آمد بالا
به منظور بررسی تأثیر بحران مالی بر بیکاری کشور­های با درآمد بالا، الگوی زیر برآورد می­گردد:
(۴-۴)
به منظور تعیین نوع برآورد به روش داده ­های تلفیقی و داده ­های پانلی، از آزمون F- لیمر و برای تعیین برآورد معادله به روش اثرات ثابت یا اثرات تصادفی، از آزمون هاسمن استفاده می­گردد. بر اساس مدل مورد استفاده و نتایج آزمون F- لیمر که نشان­دهنده بزرگتر بودن F محاسبه شده از F جدول است، فرضیه­ صفر مبنی بر این­که داده ­ها به صورت ترکیبی هستند، به نفع فرضیه­ مقابل یعنی قابلیت استفاده از روش داده ­های پانلی رد شده است. بنابراین از روش داده ­های پانلی برای برآورد مدل استفادهمی­شود.نتایج بدست آمده از برآورد معادله به وسیله­ آزمون هاسمن نیز نشان می­دهد که فرضیه­ صفر مبنی بر این­که باید از روش اثرات تصادفی استفاده شود، در مقابل فرضیه­ یک که مربوط به روش اثرات ثابت است، رد شده و برآورد مدل بر اساس روش اثرات ثابت صورت گرفته است. نتایج برآورد تأثیر بحران مالی بر نرخ بیکاری کشورهایی با درآمد بالا را می­توان در جدول (۴-۵) به صورت زیر نشان داد.
جدول (۴-۵): نتایج برآورد معادله­ (۴-۲) بر مبنای اثرات ثابت

متغیر­های توضیحی

– کاربرد با بازیابی مقید ( استفاده از نیرو ) : به کاربردهایی اطلاق می شود که در آنها یک نیروی خارجی جلوی بازیابی کرنش در آلیاژ را می گیرد . در این حالت هیچ کرنشی بازیابی نمی شود . ولی مقدار زیادی تنش ایجاد می شود ، مثل چفت و بست ها.
– کاربردهایی با بازیابی تحت فشار ( استفاده از کار ) : هم تنش و هم کرنش حین گرم شدن بازیابی می شوند و کار مکانیکی ایجاد می شود.
در واقع نیروی محرک نیروی اعمالی برای انتقال حاصل تفاوت بین انرژی آزاد Gibbs دو فاز است و بستگی دارد به دما و تنش اعمالی.

(( اینجا فقط تکه ای از متن درج شده است. برای خرید متن کامل فایل پایان نامه با فرمت ورد می توانید به سایت feko.ir مراجعه نمایید و کلمه کلیدی مورد نظرتان را جستجو نمایید. ))

در انتقال فاز رفت نیروی محرک با افزایش کرنش الاستیک و انرژی های داخلی و مقاومت ناشی از هرگونه حرکت داخلی متعادل می شود و در انتقال فاز برگشت تنها کرنش الاستیک ذخیره شده و انرژی های داخلی در مقابل نیروهای محرک هستند.
۲–۴: کریستالوگرافی آلیاژ Nitinol
در این قسمت به بررسی فرایند تبدیل فازها از دیدگاه مولکولی می پردازیم:
۲-۴-۱ : تغییر حالت های مارتنزیتی و پدیده حافظه دار شدن
تغییر حالت متالورژیکی جامدات از دو طریقه زیر امکان پذیر است .
الف) حرکت و جابجایی اتم ها وابسته به درجه حرارت و زمان با تغییر در ترکیب شیمیایی فاز جدید نسبت به زمینه قبلی.
ب) تغییر آرایش اتمی به صورت هماهنگ وابسته به دما و بدون وابستگی به زمان و هیچگونه تغییری در ترکیب شیمیایی فاز جدید نسبت به زمینه قبلی .
تغییر حالت های مارتنزیتی به طریقه دوم مرتبط است ودارای مشخصات زیراست:
۱) تغییر مکان به صورت شبه برشی می باشد و در آن اتم ها به صورت هماهنگ و گروهی جابجا می شود.
۲) دیفوزیون اتمی در آن اتفاق نمی افتد.
رفتار حافظه دار شدن کاملاً به مشخصه اول مرتبط بوده و نظم اتم های آلیاژ نباید به هم بخورد.
۲-۴-۲:تغییر حالت تبدیل آستنیت به مارتنزیت
از لحاظ کریستالوگرافی این تغییر در سه مرحله قابل بررسی است .
۱- تغییر فرم شبکه ای
۲- برش ناهمگن
۳- دوران شبکه ای
فرایند تبدیل آستنیت به مارتنزیت در مرحله تغییر فرم شبکه ای در شکل (۲-۱۶) نشان داده شده است . در این مرحله اتم ها با جابه جایی جزئی و هماهنگ، پیشروی فصل مشترک از هر لایه اتمی را موجب می شوند.
شکل۲-۱۶: تغییر فرم شبکه ای آستنیت به مارتنزیت(از راست به چپ)
باید توجه داشت پدیده حافظه داری بدون تغییر حجم و تغییر شکل امکان پذیر بوده و برش ناهمگن توجیه کننده این مطالب می باشد.
برش ناهمگن در مارتنزیت به دو طریق امکان پذیر است:
۱)مکانیزیم لغزش یافتن صفحات اتمی
۲)مکانیزیم تشکیل دوقلویی ها
تصاویر نشان داده شده در شکل ۲-۱۷چگونگی انطباق فاز مارتنزیت بر فاز آستنیت را در هنگام جابجایی جزیی و گروهی اتمها با حفظ شبکه کریستالی نشان می دهد.
باید توجه داشت که لغزش صفحات اتمی به علت شکسته شدن باند های اتمی بعنوان مکانیزیم تغییر فرم پلاستیک دائم محسوب می شود، در صورتی که در مکانیزیم دو قلویی به علت انرژی پایین مرز دوقلویی و برخورداری از تحرک و لغزندگی نسبی تغییر فرم غیر دائم است.
شکل ۲-۱۷: انطباق مارتنزیت بر آستنیت(الف مکانیزم افزایش یافتن صفحات اتمی و ب مکانیزم تشکیل دوقلوئی)
در آلیاژهای حافظه دار ، کرنش های ناشی از تغییر حالت در اثر تشکیل یک جفت از دوقلویی های دو طرف مرز ذخیره سازی می شوند و برای برگشت پذیری از آن استفاده می شود.
توضیح اینکه مارتنزیت ایجاد شده توسط دما دوقلوئی^[۲۰]و ایجاد شده توسط تنش غیر دوقلوئی^[۲۱]نامیده می شود و مارتنزیت غیر دو قلو با اعمال تنش به مارتنزیت دو قلو تبدیل میشود.
پس اگر بارگذاری مکانیکی روی آلیاژ در فاز مارتنزیت دوقلویی انجام شود ، مارتنزیت از حالت دوقلویی خارج شده و تغییر شکل می دهد. به محض برداشتن بار ، مارتنزیت به همان حالت باقی می ماند. با گرم کردن آلیاژ بالا تر از دمای فاز مارتنزیت به آستنیت تغییر می یابد و اگر فقط با بارگذاری مکانیکی اعمال شود نتیجه این عمل مارتنزیت دوقلویی نشده به همراه مقدار زیادی کرنش است. حال اگر دمای آلیژ بالای باشد ؛ به محض عدم بارگذاری ، آلیاژ به شکل اول خود باز می گردد.
شکل (۲-۱۸) مرز دوقلویی را نمایش می دهد و هر یک از دوقلویی های دو طرف مرز دوقلویی یک وا ریانت را شامل می شود. در صورت وارد کردن تنش برشی به مرز دو قلویی باعث حرکت یکی از واریانت ها شده و واریانت دیگری حذف می شود.(شکل ۲-۱۸ ،B) این روند می تواند تا تبدیل تمامی واریانت به یک واریانت واحد ادامه یابد(شکل ۲-۱۸،C) .
شکل ۲-۱۸:تشکیل دوقلوئی
شکل۲-۱۰(تکرار): پروسه انتقال برگشت
۲-۴-۳: بررسی پدیده حافظه داری
بررسی پدیده حافظه داری در تک کریستال آستنیت در شکل (۲-۱۹)نمایش داده شده است:
شکل ۲-۱۹: مکانیسم پدیده حافظه داری(a تک کریستال آستنیت -b تشکیل واریانت های A,B,C,D در اثر سرد شدن نمونه زیر دمای و c رشد واریانت A در اثر حذف واریانتهای دیگر
درمرحله اول همانطور که از شکل پیداست بعد از سرد کردن کریستال در زیر دمای واریانت های A و B و C و D تشکیل می شوند. در مرحله دوم با وارد کردن تنش به کریستال ، واریانتها شروع به حرکت و حذف شدن می کنند تا واریانت واحد A تشکیل گردد. حین تشکیل واریانت واحد A کرنش هایی در جهت واریانتA ذخیره می شود. مرحله سوم مربوط به حرارت دادن کریستال نمونه برای تبدیل مارتنزیت به آستینت می باشد از آنجاییکه کرنش ها تنها در جهت واریانت A ذخیره شده اند، پس تنها مسیر برای برگشت پذیری، واریانت A می باشد و نمونه به شکل اولیه خود باز می گردد.
۲-۵:پیشینه تحقیق
تحقیقات انجام شده درباره کیفیت نقش آلیاژهای حافظه دار در مقایسه با روش های موجودبه شرح زیر می باشد:
۲-۵-۱: مقایسه رفتار لرزه ای سیستم های دارای بادبند ساخته شده از آلیاژهای حافظه دار شکلی و سیستم های دارای بادبندهای ^[۲۲]BRBدر اثر زلزله با آنالیز دینامیکی غیرخطی تحت اثر تاریخچه های زمانی متفاوت [۱۶]
در این تحقیق آلیاژهای حافظه دار شکلی به صورت رشته سیم هائی بر روی مهاربند فولادی نصب می شوند و مشخصات مادی و هندسی آنها از قبیل مدول الاستیسیته و سطح مقطع و طول رشته سیم ها طوری انتخاب می شوند که رفتار مورد نظر از بادبند را تحقق بخشند. سطح مقطع قسمت فولادی مهاربند طوری انتخاب می شود که تغییر شکل های غیر خطی و استهلاک انرژی تنها در قسمت میراگر اتفاق بیفتد و قسمت فولادی دچار تغییر شکل های پلاستیک نشود.
مشخصات مدل:
با فرض تقارن در پلان می توان تنها یک قاب دو بعدی از سازه را مورد تحلیل و بررسی قرار داد.ارتفاع هر طبقه ۳.۹ متر بوده و پلان سازه دارای ابعاد ۹.۱ در ۹.۱ متر می باشد. طبقات دارای سقف کامپوزیت با ارتفاع قسمت فولادی ۷.۶ و اندکی پوشش بتنی می باشد.
مشخصات مادی و هندسی تیرها و ستون ها در هر دو مدل مهار بندی BRB و SMA یکسان است. بادبندهای SMA در هر طبقه به نحوی انتخاب می شوند که سختی محوری و ظرفیت محوری یکسانی با بادبندهای BRB داشته باشند(جدول ۲-۱ و شکل ۲-۲۰)
شکل ۲-۲۰ : مدل به کار گرفته شده در تحقیق[۱۶]
جدول ۲-۱: سایر ویژگی های مدل[۱۶]
جرم هر طبقه به صورت متمرکز بر روی گره های سازه ای در چهار گوشه پلان هر طبقه در نظر گرفته شد. اتصال ستون ها به پی و اتصال تیرها به ستون به صورت مفصلی در نظر گرفته شدند.
از تغییر شکل محوری تیرها صرف نظر شده و سقف طبقات صلب در نظر گرفته شده است.

۸ـ عبدالله بن جبیر.^[۳۲۹]
۹ـ عبیدالله بن عبدالله بن مسعود.^[۳۳۰]
۱۰ـ یزید بن شراحیل انصاری.^[۳۳۱]
۱۱ـ عبدالله بن مسعود.^[۳۳۲]
۱۲ـ ابو ثروان.^[۳۳۳]
۱۳ـ عبدالله بن عباس.^[۳۳۴]
۱۴ـ مسلم، غلام حضرت امیر۷.^[۳۳۵]
۲.۱.۳. گردآوران صدقات (زکات)
۱ـ مخنف بن سلیم : مسئول جمع آوری صدقاتِ^[۳۳۶] بکر بن وائل.^[۳۳۷]
۲ـ عبدالله بن زمعه. (در بعضی از نقلها سمت کارگزارای صدقات برای پدر و یا جدّ وی ذکر شده است).^[۳۳۸]

(( اینجا فقط تکه ای از متن درج شده است. برای خرید متن کامل فایل پایان نامه با فرمت ورد می توانید به سایت feko.ir مراجعه نمایید و کلمه کلیدی مورد نظرتان را جستجو نمایید. ))

۳ـ مسیّب بن نجبه : مسئول جمع آوری زکات و صدقات کوفه.^[۳۳۹]
۴ـ عبدالرحمان کندی : مسئول جمع آوری زکات و صدقات کوفه.^[۳۴۰]
۵ـ شریک بن نمله محاربی.^[۳۴۱]
شایان ذکر است که بنا بر نقل اصابه، شریک بن نمله از جانب عمر، کارگزار صدقات بوده است.^[۳۴۲]
۶ـ جلاس بن عمیر کلبی : کاتب زکات سماوه. ^[۳۴۳]ـ^[۳۴۴]
۷ـ عروه بن عشبه کلبی (عمرو بن مالک بن عشبه)^[۳۴۵] : مسئول جمع آوری زکات سماوه.^[۳۴۶]
۸ـ جعفر بن عبدالله اشجعی : مسئول جمع آوری زکات سماوه.^[۳۴۷]
۹ـ مخنف بن سلیم ازدی (غامدی) : کارگزار صدقات^[۳۴۸] بکر بن وائل.^[۳۴۹]
۲.۱.۴. کارگزاران خراج
۱ـ مصعب بن یزید انصاری : مسئول جمع آوری خراجِ برخی از مناطق مدائن.^[۳۵۰]
۲ـ ضبیعه بن زهیر اسدی (زهیر بن ضبیعه).^[۳۵۱]
۳ـ مردی از ثقیف : مسئول جمع آوری مالیات بانقیا و برخی از مناطق کوفه،^[۳۵۲] و نیز عکبرا (عکبراء ـ عکبری)^[۳۵۳] و بزرج سابور^[۳۵۴] (مدرج سابور).^[۳۵۵]
۴ـ قرظه بن کعب : مسؤول جمع آوری مالیّات منطقه‌ی عین التمر،^[۳۵۶] نهران و فلالیج (فلالیح).^[۳۵۷]
۵ـ زیاد بن ابیه : مسؤول جمع آوری مالیّات بصره.^[۳۵۸]
۲.۱.۵. کارگزاران بیت المال
۱ـ ابوالاسود الدؤلی : مسئول بیت المال بصره.^[۳۵۹]
۲ـ عبدالله بن یحیی : قسّام بیت المال.^[۳۶۰]
۳ـ ابوالهیثم بن تهیان انصاری.^[۳۶۱]
۴ـ عمّار بن یاسر.^[۳۶۲]
۵ـ عبیدالله بن ابی رافع.^[۳۶۳]
۶ـ ابوبکره^[۳۶۴] نفیع بن مسروح (نفیع بن حارث بن کلده).^[۳۶۵]
۷ـ عبدالرحمان بن ابی بکره.^[۳۶۶]
۸ـ أبو رافع.^[۳۶۷] (نام او را ابراهیم^[۳۶۸]، اسلم، سنان و عبدالرحمن^[۳۶۹] ذکر کرده‌اند).
۹ـ زیاد بن ابیه.^[۳۷۰]
۱۰ـ وهب بن عبدالله، ابو جحیفه : مسئول بیت المال کوفه.^[۳۷۱]
۱۱ـ علی بن ابی رافع.^[۳۷۲]
۲.۱.۶. قضات
۱ـ شریح بن حارث : قاضی کوفه.^[۳۷۳]
۲ـ محمد بن یزید (زید)^[۳۷۴] بن خلیده شیبانی : قاضی کوفه.^[۳۷۵]
۳ـ سعید بن نمران : قاضی کوفه.^[۳۷۶]
۴ـ عبیده سلمانی : قاضی کوفه.^[۳۷۷]
۵ـ ابوالاسود الدؤلی : قاضی بصره.^[۳۷۸]
۶ـ عبدالرّحمن بن یزید حُدّانی : قاضی بصره.^[۳۷۹]
۷ـ ضحاک بن عبدالله هلالی : قاضی بصره.^[۳۸۰]
۸ـ عبدالله بن فضاله لیثی : قاضی بصره.^[۳۸۱]
۹ـ رفاعه بن شدّاد بجلی : قاضی اهواز.^[۳۸۲]
۲.۱.۷. بازرسان
۱ـ سعد بن حارث خزاعی
ابن ابی‌الحدید می‌نویسد، امیرالمؤمنین۷ یکی از غلامانش را به نام سعد، نزد زیاد بن ابیه ـ که جانشین ابن عباس در بصره بود ـ فرستاد تا اموال بصره را به کوفه منتقل کند امّا بین سعد و زیاد بگو مگوهایی رخ داد.^[۳۸۳]
بر پایه‌ی این گزارش، کتاب سیمای کارگزاران علی بن ابی‌طالب۷، سعد را جزو بازرسین معرّفی کرده است^[۳۸۴] امّا بعید به نظر می‌رسد که بتوان وی را جزو بازرسین دانست.
۲ـ مالک بن کعب ارحبی
حضرت امیر۷ به وی مأموریت داد تا عملکرد کارگزاران در محدوده‌ی عذیب و دجله را بازرسی کند و نتیجه را به حضرت گزارش دهد.^[۳۸۵]
شایان ذکر است که یعقوبی نامه‌ای را نقل کرده است که حضرت امیر۷ به شخصی به نام «کعب بن مالک» نوشته است و به وی مأموریت داد، عملکرد کارگزاران در محدوده‌ی عذیب و دجله را بازرسی کند و نتیجه را به حضرت گزارش دهد.^[۳۸۶] امّا به نظر می‌رسد در نقل نام اشتباه شده است و نقل درست همان «مالک بن کعب» است که در تاریخ یعقوبی به اشتباه «کعب بن مالک» ذکر شده است؛ مؤید این مطلب نامه‌ای است که بلاذری در انساب الاشراف^[۳۸۷] نقل کرده است که در آن، «مالک بن کعب» ذکر شده است نه «کعب بن مالک» و مضمون نامه‌ی بلاذری با نامه‌ای که یعقوبی نقل کرده است مشابهت دارد ولی نقل یعقوبی مطالب بیشتری را در بردارد و به نظر می‌رسد که هر دو نامه، یکی باشد. علاوه بر این، در کتب رجال، فردی با نام «کعب بن مالک» به عنوان کارگزار امیرالمؤمنین۷ ثبت نشده است.^[۳۸۸] شیخ طوسی در رجال، «کعب بن مالک» را جزو اصحاب پیامبر۶ ذکر کرده است.^[۳۸۹] در برخی از منابع تاریخی نیز آمده است که «کعب بن مالک» جزو چند نفری بوده است که در جنگ تبوک شرکت نکردند و مورد مذمّت پیامبر۶ قرار گرفتند.^[۳۹۰]

شکل۵–۳ پیش ­بینی قیمت تمام شده مرغ در مدل پایه………………………………………………………….. ۷۹
شکل۵–۴ پیش­­بینی میزان تولید مرغ در مدل پایه…………………………………………………………………. ۷۹
شکل۵–۵ پیش ­بینی قیمت مرغ زنده در مدل پایه ………………………………………………………………… ۷۹
شکل۵–۶ پیش ­بینی کل سود تولید گوشت مرغ در مدل پایه…………………………………………………. ۸۰
شکل۵–۷ پیش ­بینی مقدار شاخص تولید در مدل پایه ………………………………………………………….. ۸۰
شکل۵–۸ پیش ­بینی هزینه تولید جوجه در مدل پایه…………………………………………………………….. ۸۰
شکل۵–۹ پیش ­بینی کل سود تولید جوجه در مدل پایه ………………………………………………………… ۸۱
شکل۵-۱۰ پیش ­بینی میزان ضریب تبدیل در مدل پایه…………………………………………………………. ۸۱
شکل۵-۱۱ نتایج سناریو پردازی- قیمت ذرت……………………………………………………………………. ۸۲
شکل۵–۱۲ نتایج سناریو پردازی- قیمت سویا……………………………………………………………………. ۸۲
شکل۵–۱۳ نتایج سناریو پردازی- هزینه تولید جوجه………………………………………………………….. ۸۳
شکل ۵–۱۴ نتایج سناریو پردازی- قیمت جوجه…………………………………………………………………. ۸۳
شکل ۵–۱۵ نتایج سناریو پردازی- کل سود تولید جوجه…………………………………………………….. ۸۴
شکل ۵–۱۶ نتایج سناریو پردازی- قیمت تمام شده تولید مرغ……………………………………………… ۸۴
شکل ۵–۱۷ نتایج سناریو پردازی-شاخص تولید………………………………………………………………… ۸۵
شکل ۵–۱۸ نتایج سناریو پردازی- میزان تولید مرغ آماده…………………………………………………….. ۸۶
شکل ۵–۱۹ نتایج سناریو پردازی- قیمت مرغ زنده…………………………………………………………….. ۸۶
شکل ۵–۲۰ نتایج سناریو پردازی- کل سود مرغ گوشتی……………………………………………………… ۸۷
فصل اول
کلیات پژوهش
۱-۱ مقدمه
مدیریت زنجیره تامین، یکی از استراتژی های رقابتی مهم برای افزایش بهره­وری و سود­دهی در هر سازمان است. مفهوم «مدیریت زنجیره تامین» منجر به تغییر تفکر مدیران در جهت یکپارچه­سازی تمام مراحل تهیه و توزیع کالا و خدمت می­ شود. این نوع نگرش نه تنها مدیران را مجبور می­ کند که کارکرد های سازمان خود را کارا و بهره­ور نمایند، بلکه باید منافع کلیه حلقه­های زنجیره عرضه را نیز مد­نظر قرار­ دهند. چون ناکارآمدی هر یک از آنها روی کل زنجیره اثر منفی دارد. نگرش مدیریت زنجیره عرضه با مدیریت سنتی که قبل از آن متداول بود، کاملا متفاوت است. در مدیریت سنتی هدف بهینه­سازی در درون مرزهای یک حلقه­ی زنجیره است در حالی که در مدیریت زنجیره­ی عرضه، بهینه­سازی کل زنجیره مدنظر است. نتیجه­ این هماهنگی و یکپارچگی در حلقه­های زنجیره عرضه، ارائه کالا و خدمت با هزینه ی کمتر و کیفیت بهتر می­باشد.

در مدیریت سنتی پیکاری همیشگی بین دو حلقه­ی پیاپی زنجیره در میان بود. مثلا تولید­کننده و توزیع­کننده همیشه به دنبال قیمت مطلوب­تر و نگهداشت موجودی کمتر بودند. در حالی که شاید این پیکار به تضعیف یا نابودی هر دوی آنها بیانجامد. پیامد این کار پیدایش یک اختلال در کل زنجیره است که همه حلقه­ها از آن آسیب می­بینند. پس بهینه ­سازی هر یک از زیرسیستم­ها به بهینه ­سازی کل سیستم منجر نمی­ شود. اما در مدیریت زنجیره تامین با دید کل­نگر به سیستم از بروز چنین مشکلاتی جلوگیری می­ شود. این امر نیازمند همکاری و هماهنگی بین تمام حلقه­های زنجیره­ی عرضه و روش­های جدید برنامه­ ریزی می­باشد و دانش مدیریت زنجیره تامین به دنبال گسترش این رویکرد است [۱]. پورتر زنجیره تامین را شامل همه فعالیت­های مورد نیاز برای ارائه یک خدمت یا محصول به مشتری نهایی می­داند. از دیدگاه او به زنجیره تامین، مراحل ساخت و توزیع به عنوان قسمتی از جریان کالا و خدمات افزوده
می­ شود. به عبارت دیگر با این رویکرد، زنجیره تامین در بر گیرنده سه حوزه تدارک، تولید و توزیع
می­گردد [۲]. این رویکرد می ­تواند تمام مراحل ایجاد ارزش در طول زنجیره را شناسایی و جایگاه آن را بررسی نماید. [۱] [۲] [۳]
صنعت مرغداری دومین صنعت بعد از صنعت نفت کشور است. بر اساس آمار فائو صنعت مرغ ایران رتبه­ی هفتم در تولید مرغ در جهان را داراست [۳]. با وجود گستردگی واحدهای تولید و اهمیت آن در سبد غذایی مردم، این صنعت در کشور ما با مشکلات فراوانی روبرو است که در نوسان و افزایش قیمت گوشت مرغ خود را نشان می­دهد و آنها را می­توان در دو دسته کلی برشمرد: (۱) ضعف مدیریت مرغداری­ها و عدم تمایل به آموزش در مدیران سنتی، بالا بودن دوره پرورش، عدم فعالیت و یا نیمه­فعال بودن بسیاری از واحد­های تولیدی، فقدان مدیریت واحد در زنجیره تولید، نوسان تولید در جوجه یک روزه، مستقل بودن واحدهای تولیدی و عدم وجود سیستم قراردادی (۲) هزینه بالای تولید به دلیل بالا بودن ضریب تبدیل، تجهیزات قدیمی، مصرف انرژی بالا، نوسان قیمت نهاده­های دان، عدم استفاده کامل از دان پلت.
نوسانات قیمت گوشت مرغ و نهاده های آن چالش های اساسی این صنعت است که هر ساله موجب عدم تعادل در بازار می­گردد و چاره جویی برای کاهش این نوسانات یکی از نگرانی­های دولت است. این نوسانات می ­تواند تابعی از عوامل عرضه مانند تغییرات قیمت عوامل تولید و یا تقاضای مصرف کنندگان باشد. از این رو با وجود تلاش­هایی که دولت برای تنظیم بازار برداشته، اما هنوز به موفقیت قطعی کامل در تنظیم بازار این محصول دست نیافته است [۴].[۴]
گوشت مرغ به دلایل گوناگون اقتصادی و بهداشتی از جایگاه ویژه­ای در بین اقلام مختلف پروتئین حیوانی برخوردار است. این امر به نوبه خود، رشد و گسترش صنعت مرغداری در ایران و جهان را به همراه داشته است به گونه ­ای که تولید گوشت مرغ در ایران از ۴۲ هزار تن در سال ۱۳۴۴ به رقم ۱۸۷۱ هزار تن در سال ۱۳۹۱ افزایش یافته و از این جهت رتبه هشتم را در بین کشورهای تولید کننده این محصول دارد. همسو با افزایش تولید، سرانه مصرف این محصول نیز از ۱۸ کیلوگرم در سال ۱۳۵۰ به رقم ۲۴٫۴ کیلوگرم در سال ۱۳۹۰ رسیده است و از متوسط سرانه مصرف گوشت مرغ در دنیا (۱۴کیلوگرم) نیز بیشتر می­باشد [۳].[۳]
صنعت مرغداری به دلایل گوناگونی از جمله سرعت بالای رشد طیور در زمان کوتاه نسبت به سایر دامها، ضریب تبدیل غذایی پایین، امکان تولید در تمام شرایط آب­و­هوایی، سرمایه پایین نسبت به سایر صنایع دامپروری دارای اولویت است. افزون بر این، مصرف گوشت مرغ از لحاظ جنبه­ های بهداشتی نیز به دلیل عدم امکان انتقال سریع بیماریهای میکروبی، ویروسی و انگلی به انسان، ارزش غذایی بالا، ترکیب مناسب اسیدهای آمینه ضروری، ضریب هضم بالا و کلسترول پایین­تر نسبت به گوشت قرمز برتری دارد. همچنین، شایان ذکر است که این صنعت در سالهای اخیر پیشرفت های فناوری وسیعی در زمینه اصلاح ژنتیکی طیور داشته که از لحاظ کمی و کیفی تولید را تحت تاثیر قرارداده است. اما به رغم این اهمیت و کوشش های دولت در حمایت از این صنعت، تولیدات این بخش هنوز برای تقاضای داخلی هم کافی نیست.
بررسی­های انجام شده توسط مرکز پشتیبانی طیور شرکت سهامی پشتیبانی امور دام کشور نشان
می­دهد که به طور متوسط حدود ۶۰ درصد هزینه تولید یک کیلوگرم گوشت مرغ اختصاص به هزینه خوراک (دان) دارد. در بین اقلام مختلف تشکیل دهنده غذای مرغ، ذرت، کنجاله سویا به طور متوسط ۶۲ تا ۷۰ درصد از هزینه خوراک مرغ را تشکیل می­دهد. آنچه بر اهمیت ذرت، کنجاله سویا، ذرت به عنوان نهاده های اصلی صنعت مرغداری می افزاید، آن است که این اقلام عمدتاً وارداتی بوده و نوسانات ناشی از قیمت واردات این محصولات همواره باعث نوسانات شدید قیمت محصول نهایی این صنعت (گوشت مرغ) گردیده است. به رغم تلاش های صورت گرفته، صنعت مرغداری کشور هنوز نوسانات شدید قیمت نهاده های تولید و عدم کنترل مناسب قیمت آنها را از یک سو و کنترل شدید قیمت محصول تولیدی (گوشت مرغ) را از سوی دیگر تجربه می­ کند، که این امر به نوبه خود موجب کاهش سود تولید کننده و عدم رغبت به سرمایه گذاری در این بخش می­گردد [۵].[۵]
۱-۲ تعریف مساله و بیان موضوع
در این پژوهش ضمن بررسی عوامل گوناگون موثر بر عملکرد زنجیره تامین و شناسایی متغیرها، تلاش خواهد شد تا با ارائه مدلی شبیه سازی شده از زنجیره ارزش و ارزیابی صحت مدل با اطلاعات در دسترس صنعت مرغ کشور، آثار تغییر عوامل تاثیر­گذار در زنجیره ارزش مرغ بررسی و نتایج به کارگیری سیاست­های گوناگون برای بهبود تولید ارزیابی شود. مساله این طور تعریف می­ شود که با توجه به وجود ظرفیت های فراوان در صنعت مرغ که طی سالهای گذشته ایجاد شده است کشور ما هنوز بخش قابل توجهی از نیازهای داخلی خود را از خارج از کشور وارد می­ کند. نوسان قیمت مرغ باعث شده بسیاری از کسانی که در این صنعت مشغول به کار بوده ­اند زیان ببینند و مانند گذشته با تمام توان خود فعالیت نکرده یا اینکه واحد­های خود را تعطیل کنند. چالشهای پیشروی صنعت مرغ را می­توان این­گونه دسته­بندی کرد.
۱-۲-۱ نوسان قیمت مرغ
از عواملی که موجب نوسان قیمت گوشت مرغ می­ شود، می­توان به موارد زیر اشاره کرد:
۱) نوسانات تولید: کاهش تولید در شش ماهه اول سال و افزایش آن در شش ماهه دوم یکی از دلایل نوسان قیمت است؛ کاهش جوجه­ریزی در فصل گرما به دلیل نبود تجهیزات مناسب در بسیاری از واحدهای تولیدی برای جلوگیری از تلفات بیشتر ناشی از گرمازدگی می­باشد.
۲) نوسان قیمت نهاده­های تولید: همانگونه که در قسمتهای بعد بررسی خواهد شد، قیمت نهاده­ها و بویژه جوجه یک­روزه و دان که بخش اعظم هزینه تولید را تشکیل می­ دهند نوسانات فراوانی دارند که با وقفه­های مشخصی قیمت تمام شده را تحت تأثیر قرار می­دهد. به طوری که سهم هزینه دان و جوجه یکروزه در دوره ۸۷-۹۲ به طور متوسط ۵۷ و ۱۸ درصد است.
۳) نوسانات مقطعی مصرف (تغییرات تقاضا): در برخی ماه­ها )مانند فروردین( و یا بعضی مناسبت­ها (ماه های قمری مانند رمضان، محرم، عید قربان و بازگشت حجاج( مصرف گوشت افزایش
می­یابد. همچنین، در بعضی ماه­ها مانند اردیبهشت از مصرف گوشت مرغ کاسته می­ شود [۶، ۷].[۶] [۷]
۱-۲-۲ افزایش قیمت مرغ
هزینه تمام شده گوشت مرغ در کشورهای مختلف بسته به فن­آوری تولید، شرایط اقلیمی، هزینه دان و مدیریت تولید متفاوت است. در بین کشورهای عمده تولیدکننده، برزیل به دلیل دسترسی ارزان به دان مرغ و برخورداری از فن آوری بالا از هزینه پایین­تر تولید و به تبع آن قیمت تمام شده پایین­تر برخوردار است. هر اندازه که ضریب تبدیل دان به تولید گوشت بالاتر باشد، هزینه تولید نیز بیشتر می­ شود و تلاش برای کاهش این ضریب می ­تواند تاثیر معنی داری در کاهش هزینه های تولید داشته باشد .مقایسه هزینه تولید کشورهای منتخب نشان می­دهد که هزینه تولید در کشورهای چین و برزیل پایین­تر از سایر کشورهای عمده تولیدکننده است. هزینه تمام شده هر کیلوگرم گوشت مرغ در ایران از ۱٫۵۲ دلار از سال ۲۰۰۳ به ۱٫۳ دلار در سال ۲۰۱۳ کاهش ­یافته است که در مقایسه با کشورهای عمده تولیدکننده رقم بسیار بالایی است[۴، ۴۱].[۴] [۸]
دلایل افزایش قیمت گوشت مرغ را می­توان اینگونه نام برد:
۱) افزایش سطح عمومی قیمتها یا نرخ تورم
۲) افزایش نرخ برابری ارز
۳) عدم جانشینی کامل گوشت منجمد با گوشت مرغ تازه
۴) بالا بودن هزینه تولید: بررسی­ها نشان می­دهد که کشورهای تولیدکننده و صادرکننده عمده گوشت مرغ (برزیل و آمریکا و سپس چین و تایلند) دارای کمترین هزینه تولید هستند .واردکنندگان عمده­ی گوشت مرغ (روسیه، عربستان و ژاپن) بالاترین سطح هزینه ها را دارند .
عوامل مهم و اثرگذار بر هزینه تولید عبارتند از: فناوری تولید، ظرفیت تولید، آب­و­هوا و میزان دسترسی به نهاده­های دان. تراکم جوجه­ریزی و در نتیجه میزان تولید نیز در کشورهای مختلف متفاوت است و از عواملی مانند تجهیزات موجود، شرایط آب وهوایی و کیفیت جیره اثر می­پذیرد.
۵) بالا بودن دوره پرورش: با توجه به خواست بازار و پسند مصرف کنندگان، وزن نهایی جوجه (مرغ) وبه تبع آن دوره پرورش در دنیا از ۴۱ روز تا ۶۴ روز متفاوت است. این مسأله بدون شک بر هزینه­ های تولید و میزان تولید نیز اثر می­ گذارد. بررسی­ها نشان می­دهد که سنگین­ترین مرغ ها در ژاپن با متوسط وزن ۲٫۹ کیلوگرم پرورش می یابند. چنانچه از جدول ۱-۱ بر می ­آید وزن متوسط لاشه مرغ در جهان ۱٫۳۸ کیلوگرم است. وزن لاشه در کشور برزیل به عنوان یکی از تولیدکنندگان عمده گوشت مرغ ۱٫۵۴ کیلوگرم است. وزن لاشه در ایران در مقایسه با متوسط جهانی بالاتر است. بالا بودن وزن لاشه به دوره پرورش آن برمی­گردد. طبق استانداردهای جهانی حداکثر زمان پرورش مرغ گوشتی ۴۵ روز و وزن آن حداکثر یک کیلوگرم گوشت آماده طبخ بدون چربی است، ولی تولیدکنندگان مرغ گوشتی در ایران دوران پرورش را گاهی تا ۶۰ روز افزایش می­ دهند. در این حالت ضمن افزایش وزن مرغ تا ۳ کیلوگرم و حتی بیشتر، مرغ چربی دار با پوست ضخیم روانه بازار می­ شود. وزن بالای مرغ در هنگام کشتار، باعث بالارفتن ضریب تبدیل دان به گوشت مرغ و در نتیجه­ نیاز به مصرف بیشتر دان و هزینه­ بیشتر تولید
می­ شود.
شکل ۱-۱ وزن لاشه مرغ در کشور­های مختلف [۴۱] [۸]