دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
اثر بارگذاری سطحی بر لوله های فاضلاب شهری در روش های مختلف اجرا
1
10
FA
پرویز
معارف وند
دانشگاه امیر کبیر
parvizz@aut.ac.ir
مسعود
شمس الدین سعید
masoud_shams90@yahoo.com
از آنجا که طرحهای مربوط به شبکههای آب و فاضلاب شهری از داخل شهر و بعضاً از محلهای پرتردد، تجاری و مسکونی میگذرد، برآورد نشست زمین و آگاهی از میزان بارهای وارده بر لوله، پس از اتمام این پروژهها بسیار حائز اهمیت است. در این مقاله میزان نشست و تنشهای القایی وارد بر لوله، ناشی از بارگذاری سطحی در روشهای مختلف اجرا، سنتی(نقبی، ترانشهباز) و مکانیزه (میکروتونلینگ، لولهرانی) مورد بررسی قرار میگیرد. برای شبیهسازی روشهای مذکور از نرم افزار عددی UDEC استفاده شده است. خاک مورد مطالعه کنگلومرای هزار دره تهران تیپ A است. نتایج شبیهسازی نشان میدهد که بیشترین نشست سطح زمین و مساحت زون پلاستیک تشکیل شده در اطراف لوله مربوط به روش ترانشه باز است که دلیل آن پر کردن ترانشه با خاک دستی و چسبندگی پایین خاک است. بیشترین تنش محوری وارد بر لوله در روش لولهرانی دیده میشود و دلیل آن تماس کامل لوله با محیط اطراف در حین رانش و اصطکاک مابین لوله-خاک است. همچنین تغییرات جابهجایی و تنشهای وارد بر لوله در روشهای میکروتونلینگ و حفر مقنی به دلیل داشتن شکل مقطع تقریباً مشابه این روشها، همخوانی خوبی با یکدیگر دارند.
شبیهسازی,UDEC,ترانشه باز,نقبی,میکروتونلینگ,لولهرانی
https://anm.yazd.ac.ir/article_328.html
https://anm.yazd.ac.ir/article_328_07c50599c014fe9e80d6e2eafa8d5a25.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
تخمین عیار آهن در آنومالی تپه قرمز معدن سنگ آهن سنگان با استفاده از سه الگوریتم فازی عصبی تطبیقی
10
17
FA
سعید
مجدی فر
شهید باهنر کرمان
غلامرضا
کمالی
شهید باهنر کرمان
kamali@uk.ac.ir
تخمین عیار یکی از مراحل حیاتی در بررسیهای فنی و اقتصادی معادن است. بنابراین یافتن روشی که بتواند این مقادیر را با دقت بالا تخمین بزند یک ضرورت به شمار میرود. در این مطالعه، یک الگوریتم فازی عصبی تطبیقی به منظور تخمین عیار آهن کانسار های اسکارنی به کار رفته است. این الگوریتم در حقیقت یک سیستم فازی میباشد که باید توابع عضویت آن توسط دادههای آموزشی تنظیم گردد. قابل ذکر است که تعریف توابع عضویت اولیه در کیفیت نتایج خروجی این الگوریتم تاثیرگذار میباشد که بدین منظور سه الگوریتم تفکیک شبکه، خوشه بندی کاهشی و خوشه بندی فازی برای تنظیم این توابع مرسوم میباشد. در این مقاله سه سیستم فازی عصبی بر پایه ی الگوریتمهای مذکور برای تخمین عیار آهن به کار گرفته شد. نتایج حاصل از دادههای تست، بیانگر توانایی قابل توجه سیستم فازی عصبی بر پایه خوشه بندی کاهشی در تخمین عیار آهن نسبت به دو الگوریتم دیگر میباشد.
تخمین عیار,کانسار آهن اسکارنی,الگوریتم فازی عصبی,الگوریتمهای خوشه بندی
https://anm.yazd.ac.ir/article_329.html
https://anm.yazd.ac.ir/article_329_fc5eb176beb4c30229a095256a33643a.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
مدلسازی ناهمسانگردی سه بعدی عیار در کانسارهای رگهای بر اساس ضخامت و فاصله
18
27
FA
امین حسین
مرشدی
دانشگاه تهران
امیرحسین
کوهساری
دانشگاه یزد
مطالعه بر روی نهشتههای رگهای، به عنوان کوچکترین اهداف اکتشافی، امری مشکل و با عدم قطعیت همراه است. در این پژوهش، با بررسی پارامترهای ضخامت رگه، فاصله از رگه و عیار کانسنگ مورد نظر، یک مدل سه بعدی ارائه شده است. به طور کلی، دو رویکرد مدلسازی سه بعدی برای نهشتههای رگهای متصور است؛ در حالت اول، سیال کانی ساز از رگهها به سمت سنگ دیواره حرکت کرده که در این صورت هرچه رگه دارای ضخامت بیشتری باشد سنگ دیواره در برگیرنده رگه از نظر عیار فلز کانی ساز غنیتر است و دیدگاه دوم حالتی است که بر اساس حرکت سیال کانی ساز از سنگ دیواره به سمت رگه است که در این حالت هرچه رگه دارای ضخامت بیشتری باشد سنگ دیواره میزبان رگه از فلز کانی ساز تهیتر است. <br /> در کانسار رگهای عمارت، با افزایش ضخامت رگه، سنگ میزبان از فلز کانی ساز تهیتر میشود. در این مطالعه، باتوجه به رابطه لگاریتمی بین عیار و فاصله از رگه، حداکثر شعاع تأثیر و پارامترهای بیضوی جستجوی ناهمسانگردی محاسبه شد. با استفاده از مدلسازی سه بعدی توزیع عیار، ضخامت رگه و فاصله از رگه برای کمر بالا و پایین رگه، عیار در هر موقعیت فضایی با توجه به فاصله داری و ضخامت رگه تخمین زده شده است.
مدلسازی سه بعدی,ضخامت – فاصله رگه,شعاع تأثیر,بیضوی جستجوی ناهمسانگردی,معدن عمارت,ایران
https://anm.yazd.ac.ir/article_330.html
https://anm.yazd.ac.ir/article_330_0c892e7c8dcd1206e7da09d652ad55be.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
ارایه روش تکمیلی بهبود کیفی الگوریتم جداسازی آماری سینکلر بر اساس نمودارهای احتمال
28
37
FA
سیدسعید
قنادپور
دانشگاه صنعتی امیرکبیر
احمدرضا
مختاری
دانشگاه صنعتی اصفهان
ar.mokhtari@cc.iut.ac.ir
اردشیر
هزارخانی
دانشگاه صنعتی امیرکبیر
ardehez@aut.ac.ir
نادر
فتحیان پور
دانشگاه صنعتی اصفهان
fathian@cc.iut.ac.ir
تاکنون روشهای گوناگونی برای جداسازی جوامع مختلط مورد بررسی قرار گرفتهاند؛ که یکی از مهمترین و پرکاربردترین آنها روش پیشنهادی سینکلر بر پایه یافتن نقطه عطف نمودار ترکیبی و جداسازی جوامع با استفاده از آن است. در تحقیق حاضر با استفاده از مفهوم توسعه یافته این روش، در ابتدا سعی شده است که به طراحی یک نرمافزار موثر و کاربردی پرداخته شود؛ که در آن با استفاده از نمودارهای احتمال به تحـلیل و جداسازی جوامـع پرداختـه و سپس بررسی میشود که آیا میتوان از مماسهای رسم شده در اولین نقطه جدا شده به روش سینکلر، به عنوان یک روش جداسازی استفاده نمود یا خیر. از مزایای عمده این نرمافزار دقت و سرعت بالای انجام تحلیلها روی نمونههای مورد بررسی است. در حالی که انجام این فرآیند به صورت دستی معمولاً کاری زمانبر و کم دقت است. این نرمافزار دارای انعطافپذیری بالایی در محاسبه خطای فرآیند بوده و میتواند با تحلیل سناریوهای متفاوت، به جداسازی جوامع آماری با کمترین خطای محاسباتی اقدام نماید. در نهایت نیز به منظور آشنایی و بررسی عملکرد کار نرمافزار، عیار عنصر کروم جنوب رشته کوه ایرانکوه به عنوان داده در اختیار نرمافزار قرار گرفته و مشاهده شده است که نرمافزار به خوبی توانسته جوامع را با دقت بالا و در زمان بسیار کوتاه جداسازی نماید. همچنین روش پیشنهادی جداسازی به کمک مماسها، فرآیند جداسازی را نسبت به روش سینکلر با سرعت بیشتری به پایان میرساند.
جداسازی جوامع آماری,نمودار احتمال,روش سینکلر,مماس,نرمافزار جداسازی
https://anm.yazd.ac.ir/article_331.html
https://anm.yazd.ac.ir/article_331_73539c8e8387a16bbe80e5ff5a2f74b5.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
استفاده از روش ناپیوستگی جابهجایی در تحلیل مکانیسم انتشار ترکها در مواد شبه سنگی
38
49
FA
هادی
حائری
کوروش
شهریار
محمد
فاتحی مرجی
0000-0003-2917-6264
mohammad.fatehi@gmail.com
پرویز
معارف وند
parvizz@aut.ac.ir
در طبیعت، مواد ترد سنگی ترک دار (ترکهای از پیش موجود) معمولاً تحت بارگذاری فشاری قرار میگیرند. مطالعه مکانیسم انتشار ترکها در سنگ برای طراحی سازههای سنگی از اهمیت ویژهای برخوردار است. در این تحقیق، مکانیسم شروع، گسترش و بهم پیوستگی ترکهای از پیش موجود در مدلهای عددی و نمونههای شبه سنگی تحلیل شده است . همچنین، تأثیر زاویه ترک بر تنش شروع رشد ترک، تنش بهم پیوستگی و زاویه شروع رشد ترک مورد بحث قرار گرفته است. روش عددی بکارگرفته شده ناپیوستگی جابهجایی به همراه المانهای مرتبه بالا (المانهای کوادراتیک) است که یکی از زیر مجموعههای روش المان مرزی غیر مستقیم میباشد. نمونههای آزمایشگاهی از مواد شبه سنگی ( ترکیبی از سیمان، ماسه و آب) تهیه شدهاند و در آنها ترکهای از پیش موجود با زوایای مختلف و فاصله داری یکسان نسبت به همدیگر ایجاد شده است. با بارگذاری محوری این نمونهها تنش شروع ترک، تنش شکست و مسیر انتشار و بهم پیوستگی ترکها تعیین گردیده است.. اثر زاویه ترکهای اولیه نسبت به همدیگر باعث تغییر مسیر بهم پیوستگی و مسیر نهایی شکست میگردد. بر اساس نتایج بدست آمده از آزمایشات و مدلسازی عددی، ترکهای بالهای ابتدا شروع به انتشار میکنند و نقش اصلی را در شکست نهایی نمونههای شبه سنگی دارند، اما امکان تشکیل شدن ترکهای ثانویه از نوک ترکهای از پیش موجود وجود دارد. مقایسه نتایج عددی و آزمایشگاهی نشان میدهد که روند انتشار و بهم پیوستگی ترکها از تطابق بسیار خوبی برخوردار هستند.
نمونههای شبه سنگی,مسیر بهم پیوستگی ترکها,تحلیلهای عددی و آزمایشگاهی
https://anm.yazd.ac.ir/article_332.html
https://anm.yazd.ac.ir/article_332_a3c89f8591c7516b8bb713218c6edc35.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
مروری بر روشهای جدایش آنومالی های ناحیهای از باقیمانده در دادههای پتانسیل با مطالعه موردی دادههای ژئوفیزیک هوایی چهارگوش آباده به همراه تفسیر زمین شناسی
50
58
FA
کمال
علمدار
kamal.alamdar@gmail.com
احمد
معدنچی
آنومالیهای میدان پتانسیل حاصل از برداشتهای مغناطیسسنجی و گرانیسنجی از آنومالیهای ناحیهای ناشی از ساختارهای بزرگ مقیاس و آنومالیهای باقیمانده ناشی از ساختارهای محلی کوچک مقیاس تشکیل شده است. بسته به هدف و مقیاس برداشتهای میدان پتانسیل از آنومالی های ناحیهای و یا باقیمانده استفاده میشود. در نتیجه جدا کردن آنومالی های باقیمانده و ناحیهای امری ضروری است. جداسازی آنومالیها به دو شیوه ترسیمی و تحلیلی انجام میگیرد. در روشهای ترسیمی گرچه ممکن است نتایج با واقعیتهای زمین شناسی انطباق بیشتری داشته باشد ولی روشهای تحلیلی در مواقعی که با مجموعه دادههای بزرگ مقیاس مواجه هستیم ارجحیت دارند. روشهای مختلفی برای جداسازی آنومالی به روش تحلیلی وجود دارند که مهمترین آنها عبارتند از روش آنالیز روند سطح، روش فیلتر کردن طول موجی، روش گسترش میدان و روش مشتق قائم که در این نوشته این روشها بر روی دادههای مغناطیس هوایی چهارگوش1:250000 آباده به کاربرده شده است. از روی نقشه آنومالی ناحیهای به دست آمده از چهارگوش آباده گرانیت شیرکوه در شمال شرق، گسل دهشیر بافت با روند شمال غرب – جنوب شرق و کفه طاقستان در بخشهای مرکزی مشخص شده است.
آنومالیهای میدان پتانسیل,آنومالی باقیمانده,فیلتر کردن طول موجی,گسترش میدان,مشتق قائم
https://anm.yazd.ac.ir/article_333.html
https://anm.yazd.ac.ir/article_333_b993d4a4f2bf9fc44af5116f93a966ad.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
تعیین پتانسیل فشارندگی به کمک روشهای آماری و ارائه راهکار بهینه جهت عبور از زون گسله پورکان ـ وردیج واقع در مسیر تونل انتقال آب از سد امیرکبیر به تهران
59
65
FA
محمدرضا
جاری
mreza.jari.46@gmail.com
سیدیحیی
روحانی
پرویز
رفعتی
یکی از عوامل مهم در حفاری تونلهای مکانیزه انتخاب مسیر مناسب است. در این تحقیق به بررسی چگونگی انتخاب مسیر قطعه دوم تونل انتقال آب از سد امیرکبیر به تهران پرداخته شده است. با توجه به وجود گسلهای متعدد در طول مسیر این تونل از جمله گسل فعال پورکان – وردیج و احتمال وقوع پدیده فشارندگی و تبعات ناشی از آن، پیش بینی احتمال وقوع این پدیده بر اساس روشهای موجود شامل روش سینگ، جتوا، گوئل و هوک و مارینوس در مسیرهای پیشنهادی شامل مسیر مستقیم و مسیر قوسی( مسیر با طول بیشتر به منظور دور شدن از هسته مرکزی گسل) در دستور کار قرار گرفت و با تخصیص درصد ریسک احتمال وقوع پدیده فشارندگی بر اساس طبقه بندیهای پیشنهادی توسط محققین فوقالذکر، محاسبه درصد ریسک میانگین وقوع پدیده فشارندگی در دو مسیر مستقیم و قوسی با استفاده از تحلیلهای آماری انجام شد. نتایج گویای آن است که درصد ریسک میانگین احتمال وقوع پدیده مذکور در حین حفاری مسیر اجرا شده فعلی(مسیر قوسی) در قیاس با مسیر مستقیم در حدود 69/42% کمتر میباشد. همچنین با بررسی نقشههای زمین شناسی منطقه میتوان دریافت که طول زون گسل پورکان ـ وردیج در حین حفاری در مسیر مستقیم 125 متر است که تقریباً 2/5 برابر مسیر قوسی بوده است. بنابراین با توجه به افزایش احتمال وقوع پدیده فشارندگی در حین حفاری در مسیر مستقیم نسبت به مسیر قوسی و نیز با لحاظ طول محدوده گسله، محاسبه نسبت احتمال وقوع پدیده فشارندگی در حین حفاری مسیر مستقیم نسبت به مسیر قوسی انجام و مشخص گردید که نسبت مذکور 4/17 بوده که گویای میزان درصد ریسک بسیار بالا در حین حفاری مسیر مستقیم نسبت به مسیر قوسی است. بنابراین با در نظر گرفتن ریسک حفاری مسیر مستقیم و احتمال گیر کردن و یا مدفون شدن دستگاه حفاری TBM در حین حفاری مسیر مذکور، اجرای عملیات حفاری مسیر قوسی به منظور جلوگیری از بروز مخاطرات احتمالی در دستور کار قرار گرفت.
فشارندگی,درصد ریسک احتمال وقوع فشارندگی,نسبت احتمال وقوع پدیده فشارندگی,مسیر مستقیم,مسیر قوسی
https://anm.yazd.ac.ir/article_334.html
https://anm.yazd.ac.ir/article_334_f7da936b77aecdaa8318dc8199d722e7.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
ارائه مدلی جدید جهت تخمین هزینههای سنگشکنهای فکی با بازوی مضاعف
66
76
FA
مهسا
خوشفرمان برجی
احمدرضا
صیادی
sayadi@modares.ac.ir
محمدرضا
خالصی
عملیات خردایش از جمله کلیدیترین مراحل فرایند فرآوری مواد معدنی محسوب شده و سنگشکنهای فکی یکی از پرکاربردترین تجهیزات اصلی در این بخش هستند. انتخاب این تجهیزات تابعی از متغیرهای متعددی به ویژه قیمت تمام شده عملیات سنگشکنی است. در این تحقیق مدلی جدید جهت تخمین هزینههای سرمایهای و عملیاتی سنگشکنهای فکی با بازوی مضاعف در قالب توابع رگرسیون تک متغیره و چند متغیره ارائه شده است. در تحلیل چند متغیره، به منظور افزایش دقت مدل، دادههای جمعآوری شده مورد بررسی دقیق قرار گرفته و همبستگی بین آنها با استفاده از روش آنالیز مؤلفههای اصلی بر طرف شده است. در تابع رگرسیون تک متغیره، سطح دهانه ورودی به عنوان متغیر مستقل در نظر گرفته شده است در حالیکه در تابع هزینه چند متغیره سطح دهانه ورودی به همراه توان سنگشکن به عنوان متغیرهای مستقل جهت ساخت مدل تخمین هزینه مورد استفاده قرار گرفتهاند. کارایی هر یک از توابع هزینه با استفاده از میانگین خطای مطلق و بررسی باقیماندهها سنجیده شده است. به طور میانگین، میزان خطای 93/14 و 91/10 درصد به ترتیب برای رگرسیون تک و چند متغیره محاسبه شده که میزان قابل قبولی در مطالعات پیشامکانسنجی و امکانسنجی است. همچنین توزیع نرمال باقیماندهها گواهی بر اعتبار مدلهای ارائه شده است. به منظور نشان دادن چگونگی کارکرد مدل، هزینههای سرمایهای و عملیاتی واحد سنگ شکنی اولیه یک معدنی فرضی طلا به کمک این مدل برآورده شده است.
تخمین هزینه,سنگشکن فکی,رگرسیون,تحلیل مؤلفههای اصلی
https://anm.yazd.ac.ir/article_335.html
https://anm.yazd.ac.ir/article_335_6aab30e1eccbd5367fbcac5e4e52df39.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
ارائه روشی برای بررسی ساختار سیستم اندازهگیری مدارهای پیچیده کارخانههای فرآوری مواد معدنی ـ استخراج بخش قابل مشاهده و انجام موازنه جرم
77
88
FA
فردیس
نخعی
علی اکبر
عبداله زاده
مهدی
ایران نژاد
0000-0002-5469-084X
iranajad@aut.ac.ir
برای انجام محاسبات متالورژیکی، اقتصادی و کنترل روند عملکرد کارخانه، در اختیار داشتن اطلاعات کافی و صحیح از بخشهای مختلف و به خصوص جریانهای مواد ضروری است. معمولاً به دلایل فنی یا اقتصادی امکان اندازهگیری یا نمونهبرداری از تمام جریانهای موجود در مدار وجود ندارد و از طرفی به علت طبیعت کار ابزار اندازهگیری در کارخانه و یا فرآیندهای نمونهبرداری و آنالیز، دادههای جمع آوری شده با خود خطای غیر قابل اجتناب را به همراه دارد، لذا به کارگیری روشهای مختلف تعدیل خطاها، تحت مدلهای دقیق از جمله موازنه جرم و دستیابی به دادههای صحیح الزامی است. در این مقاله با توجه به سیستم اندازهگیری یا نمونهبرداری در مدار واحد فرآوری مورد مطالعه، ضمن معرفی روشی نوین برای انجام موازنه جرم مدار در حالت چند متغیره، نتایج بهکارگیری روش ماتریس سیکلهای پایه برای بررسی سیستم اندازهگیری مورد استفاده از نظر میزان دسترسی به اطلاعات مورد نیاز و تعیین بخشهای قابل مشاهده و غیر قابل مشاهده مدار ارائه شده است. به منظور اعتبار سنجی روش پیشنهادی کاربرد آن در مدار سولفورزدایی از کنسانتره شرکت سنگ آهن گل گهر تشریح شده است. نتایج نشان میدهد که، به کارگیری این آزمون از توانایی بالایی در بررسی صحت و تعدیل دادههای جمع آوری شده برخوردار است و به دلیل ساختار مناسب آن در نرم افزار متلب، امکان مطالعه مدارهای با حجم محاسبات زیاد نیز فراهم است.
موازنه جرم,ماتریس سیکلهای پایه,اصلاح دادهها,قابل مشاهده,مدار سولفورزدایی گلگهر
https://anm.yazd.ac.ir/article_336.html
https://anm.yazd.ac.ir/article_336_0307eea3447454e017d4abed7d8cca07.pdf
دانشگاه یزد
روش های تحلیلی و عددی در مهندسی معدن
2251-6565
2676-6795
3
5
2013
06
22
مقایسه روشهای عددی مختلف در تحلیل پایداری مغار نیروگاه زیر زمینی مسجد سلیمان
89
98
FA
سیف الدین
موسی زاده
کامران
گشتاسبی
goshtasb@modares.ac.ir
عبدالهادی
قزوینیان
هدف از این تحقیق کاربرد روشهای عددی مختلف جهت بررسی رفتار مغار نیروگاه مسجد سلیمان است. برای نیل به این هدف از روشهای عددی پیوسته مانند المان محدود( نرم افزار Phase2) و تفاضل محدود( نرمافزارFLAC) استفاده شد. همچنین از روش عددی المان مجزا( نرمافزار دو بعدی UDEC) جهت مدلسازی محیط به صورت ناپیوسته و مقایسه با دیگر روشها استفاده شد. در این تحقیق مغارهای نیروگاه با اعمال یک سیستم نگهداری بهینه با استفاده از روشهای مختلف عددی مدلسازی شدند. در نهایت نتایج حاصل از روشهای مختلف با نتایج ابزار دقیق مقایسه شدند. نتایج بدست آمده نشان میدهد که مدل ساخته شده توسط نرم افزار UDECدر مقایسه با دو نرم افزار دیگر، به شرایط واقعی زمین نزدیکتر است. با توجه به نتایج بدست آمده از مقایسه مقدار جابجاییهای رخ داده در اطراف فضای حفاری شده میتوان گفت که این محیط دارای یک رفتار شبه پیوسته است. ولی وجود لایههای ضعیف (به خصوص لایه سیلتستون در سقف مغار اصلی) باعث به وجود آمدن درزههایی با پارامترهای مقاومتی پایین شده است که این عامل باعث باز شدگی درزهها و ناپایداری در این نواحی شدهاند. مدلسازی با UDECاین موضوع را تأیید میکند.
روشهای عددی,المان محدود,تفاضل محدود,المان مجزا,تحلیل پایداری,مسجد سلیمان,Phase2,FLAC و UDEC
https://anm.yazd.ac.ir/article_337.html
https://anm.yazd.ac.ir/article_337_23813d48239fd033ef462c2ee0c771db.pdf