جستجوی مقالات مرتبط با کلیدواژه "gas pipeline" در نشریات گروه "مکانیک"
تکرار جستجوی کلیدواژه «gas pipeline» در نشریات گروه «فنی و مهندسی»-
فولادهای ترمومکانیکال به دلیل چقرمگی و مقامت بالا در برابر رشد ترک، به طور گسترده در خطوط انتقال نفت و گاز استفاده می شوند. بخش وسیعی از خطوط لوله فولادی به کار رفته در صنعت نفت و گاز جمهوری اسلامی ایران از جنس فولاد API X65 است. نوسان فشار داخلی گاز در لوله های فولادی می تواند باعث شکست خستگی و انفجار شود. به همین دلیل، بررسی آسیب و یکپارچگی سازه ای این لوله ها از اهمیت بالایی برخوردار است. در این تحقیق منحنی تنش- عمر و استحکام خستگی فلز پایه لوله فولادی API X65، با انجام آزمایش خستگی تخمین زده شده است. به این منظور، تعداد 24 و 25 نمونه آزمایشگاهی به ترتیب در راستای درز جوش (جهت طولی یا غلتک کاری) و عمود بر درز جوش (جهت عرضی کلاف اولیه) طبق استاندارد، از لوله در مقیاس صنعتی با قطر خارجی 1219 میلی متر و ضخامت 3/14 میلی متر تهیه شد. نمونه های تهیه شده تحت آزمایش خستگی خمشی- چرخشی کاملا معکوس شونده قرار گرفت و با در نظر گرفتن توزیع نرمال لگاریتمی، تحلیل آماری نتایج انجام شد. منحنی میانگین، منحنی مشخصه و بازه اطمینان نتایج آزمایش در ناحیه های عمر خستگی محدود و استحکام خستگی براساس استانداردهای ISO 12107 و ASTM-E739 به دست آمد. مقدار میانگین حد دوام فلز پایه در راستای درز جوش و عمود بر درز جوش به ترتیب برابر 291 و 305 مگاپاسکال به دست آمد. این مقادیر به خوبی در محدوده پیش بینی 4/0 تا 6/0 استحکام نهایی فولاد آزمایش شده و بالاتر از حد دوام فلز جوش (258 مگاپاسکال) این لوله قرار دارد.
کلید واژگان: حد دوام, منحنی تنش- عمر, شکست خستگی, لوله انتقال گاز, فولاد API X65Thermomechanical steels are widely used in oil and gas pipelines due to their high toughness and high resitance against crak growth. A large part of the steel pipelines used in the oil and gas industry in Iran is made of API X65 steel. The fluctuations of internal gas pressure in steel pipes can cause fatigue failure and lead to gas leakage and explosion. So, the control of damage initiation and structural integrity of gas pipelines is of great importance. In this study, the S-N curve and the fatigue strength of the base metal of the API X65 steel were estimated by performing fatigue tests. For this purpose, 24 and 25 test specimens along the seam weld in the coil transverse direction, and perpendicular to the seam weld along the coil rolling direction were prepared according to ISO 1143 standard, respectively. All test samples were cut from an spirally welded pipe with 1219mm outside diameter and 14.3mm wall thickness and were tested on a completely reverse rotating-bending fatigue machine. Statistical analysis of the results was performed by considering the normal logarithmic distribution. The mean curve, characteristic curve, and confidence interval of the results were obtained both in the finite fatigue life range and in the fatigue resistance. The mean endurance limit of the base metal perpendicular to and parallel to the seam seam were 305 and 291 MPa, respectively which were in the range of 0.4 to 0.6 of material tensile strength and above the seam weld endurance limit (258 MPa).
Keywords: Endurance Limit, S-N Curve, Fatigue Fracture, Gas Pipeline, API X65 Steel -
در این مقاله موقعیت نشتی ایجاد شده در یک لوله گاز تحت فشار تشخیص داده شده است. برای این منظور در بخش آزمایشگاهی این تحقیق، نشتی ایجاد شده در یک لوله فولادی به وسیله نازل های مختلف شبیه سازی شده است. به دلایل زیست محیطی و کاهش خطرات احتمالی در حین انجام آزمایش، از هوا به عنوان جایگزین گاز طبیعی استفاده شده است. پیش از این اکثر پژوهش های منتشر شده متکی به استفاده از حسگرهایی با فرکانس بالا در دو سمت نشتی ایجاد شده، بوده اند. در پژوهش حاضر روش نوینی برای تعیین موقعیت محل نشتی در یک لوله ارایه شده است. در این روش با بهره گیری از آنالیز میرایی امواج منتشر شده در لوله، از دو حسگر نصب شده در یک سمت لوله استفاده شده است. در روش پیشنهادی به منظور کاهش اثرگذاری عوامل محیطی در پاسخ نهایی، داده های ذخیره شده در تنها دریک بازه فرکانسی مشخص مورد استفاده قرار گرفته اند. بررسی نتایج حاصل از آزمایش های تجربی، نشان می دهد با افزایش فشار داخل لوله خطای تشخیص موقعیت کاهش می یابد. در نهایت با پردازش سیگنال های ذخیره شده، میانگین مقدار خطا در تشخیص موقعیت نشتی 77/13 سانتی متر حاصل شده است.
کلید واژگان: تشخیص موقعیت صدا, لوله گاز, تشخیص نشتی, آنالیز میراییIn this paper, the position of the leakage in a pressurized gas pipe is localized. As an experimental example, the leakage in a steel pipe is simulated by different orifices. To decrease the risk of explosion and environmental contamination during experiments, pressurized air was used instead of natural gas. Most of the published papers are based on using high sampling rate sensors placed on both sides of the leakage. In this study, a new method is extended due to leakage localization in a gas pipe. To achieve this goal, attenuation analysis is implementing by installing two sensors on one side of the pipe. In the proposed technique, the dominant frequency of recorded signals is used, to decrease the effect of background noises on final results. Experimental results reveal by increasing the pressure the leakage localization error is decreased. Finally, by implementing signal processing into the data the average value of leakage localization error is reported 13.77 cm.
Keywords: Sound Localization, gas pipeline, Leakage detection, Attenuation analysis -
جوشکاری قوس فلزی الکترود پوشش دار بر روی فولادهای کم آلیاژ استحکام بالا در خطوط لوله انتقال گاز طبیعی ایران اهمیت ویژه ای دارد. در این مقاله خواص مکانیکی و متالورژیکی موقعیت های مختلف لوله (با قطر خارجی 36 اینچ و در ساعت های 1:30 - 12، 3 - 1:30، 4:30 - 3 و 4:30 - 6) برای تعیین منطقه بحرانی جوش سر به سر چند پاسه، ناحیه متاثر از حرارت و فلز پایه به وسیله آزمایشهای کوانتومتری، متالوگرافی، کشش، ضربه و سختیسنجی بررسی شد. نتایج آزمون کشش، کمترین میزان استحکام نهایی را (در راستای عمود بر جوش و در موقعیت 3 - 4:30) برابر 463 مگاپاسکال و کمترین انرژی شکست را (در راستای عمود بر جوش و در موقعیت 1:30 - 3) برابر 96 ژول ارزیابی میکند که 61 درصد کاهش نسبت به فلز پایه را نشان می دهد. میزان تغییرات سختی در نواحی و موقعیت های گوناگون ناچیز (کمتر از 5 درصد) است. تصاویر حاصل از آزمون متالوگرافی با میکروسکوپ های نوری و الکترونی نشان می دهد که مقدار پرلیت ناحیه پاس سطح جوش و منطقه تحت تاثیر حرارت نزدیک فلز پایه نسبت به فلز اصلی به ترتیب 12. 5 و 25 درصد کاهش داشته است. مقدار کاهش نیوبیوم در موقعیت 4:30 - 6 نسبت به فلز 0. 025 می باشد.کلید واژگان: فولاد کم آلیاژ استحکام بالا, خطوط انتقال گاز, جوش سر به سرShield metal arc welding on the high strength low alloy steels in pipelines to transport natural gas from Iran is of great importance. In this paper, mechanical and metallurgical properties of multi-pass weld in weld zone, heat affected zone (HAZ), and base metal (with 36 inch outside diameter) for determination of critical area is investigated. Also Chemical analysis, metallography, tensile, impact and hardness testes in different positions of pipe (hours 12-1:30, 1:30-3, 3-4:30 and 4:30-6) were carried out. Tensile test results showed the lowest ultimate strength (Vertical to weld and in position 3-4:30) equal to 463mpa, and the lowest energy impact (Vertical to weld and in position 1:30-3) equal to 96J that determine 61 percent reduction compared to the base metal. The amount of hardness variation in different areas and positions is negligible (less than 5 percent). Images of metallographic test made by light and electron microscopes demonstrated that the amount of perlite cap pass weld and heat affected zone near the base metal were decreased compared to main metal to respectively 12.5 and 25 percent. The amount of Nb reduction in position (4:30-6) relative to the base metal is 0.025.Keywords: High strength low alloy steel, Gas pipeline, Girth weld
-
یکی از مسائل مطرح و مورد توجه در صنعت گازرسانی، تحلیل جریان گذرای گاز طبیعی درون خطوط لوله انتقال می باشد. با وجود مطالعات صورت گرفته، هنوز دقت و زمان محاسبات به عنوان دو چالش مهم در این زمینه مطرح هستند. در این مقاله، طراحی یک الگوریتم موازی جهت شبیه سازی عددی جریان هم دما و غیرهم دمای گاز مطرح شده است. تحلیل عددی جریان با استفاده از روش اختلاف محدود ضمنی تجزیه بردار شار استیگر- وارمینگ صورت پذیرفته است. موازی سازی مسئله نیز به روش انتقال پیام با کتابخانه MPI انجام شده است. به منظور نشان دادن قابلیت های برنامه کامپیوتری توسعه داده شده، دو مسئله شامل جریان درون دو خط لوله با شرایط متفاوت تحلیل و صحت سنجی شده است. در انتها، پس از صحت سنجی نتایج، معیارهای عملکرد روش پیشنهادی مثل زمان محاسبات، کاهش زمان به دست آمده و میزان افزایش سرعت از نتایج شبیه سازی استخراج شده اند. بررسی نتایج به دست آمده نشان می دهد که با استفاده از پردازش موازی می توان تا حد زیادی زمان محاسبات جریان درون خطوط لوله انتقال گاز طویل را کاهش داد. علاوه بر این، با به کارگیری این روش در شبکه های محاسباتی ریز نسبت به شبکه های درشت، کارایی الگوریتم موازی افزایش می یابد.کلید واژگان: خط لوله گاز, جریان گذرا, پردازش موازی, کتابخانه MPIAnalysis of natural gas transient flow in transmission pipelines is one of the most important issues in the gas industry. Despite the previous studies, the accuracy and the computational time have yet considered as two important challenges in this field. In this paper, a parallel algorithm for numerical simulation of isothermal and non-isothermal gas flows is presented. Numerical analysis of the flow is performed using the implicit Steger-Warming flux vector splitting method. For parallelization, the computer program has been parallelized using Message Passing Interface library. In order to demonstrate the capabilities of the developed computer program, the flow inside two pipelines with different conditions is solved, and the results are validated. Then, some factors such as the computational time, reduction of the time, and the speed up criteria are obtained to demonstrate the computational efficiency of the proposed method. The results show that parallel processing method can significantly reduce computational time of natural gas flow in long transmission pipelines. Moreover, it is shown that application of this approach on the fine computational grids is more efficient than on the coarse grids.Keywords: Gas Pipeline, Transient Flow, Parallel Processing, Message Passing Interface Library
-
وجود تنش های پسماند باعث کاهش مقدار تنش طراحی مجاز در خطوط انتقال گاز می شود. در این مقاله ابتدا دو قطعه لوله فولادی X70 با قطر 56اینچ به روش دستی جوشکاری شد. سپس آزمون تجربی کرنش سنجی سوراخ برای ارزیابی توزیع تنش های پسماند در این اتصال انجام شد. در ادامه، شبیه سازی المان محدود فرآیند جوشکاری برای ارزیابی دقیق توزیع تنش های پسماند در این اتصال انجام شد که نتایج آن با مقایسه با نتایج تجربی صحه گذاری شد. نتایج کیفی به دست آمده به میزان مناسبی با نتایج آزمایشگاهی همخوانی دارد. نتایج شبیه سازی (با اختلاف 15% نسبت به نتایج تجربی) حداکثر تنش پسماند را در راستای محیطی درز جوش خارجی لوله ارزیابی نمود. نتایج تجربی نشان داد که حداکثر تنش پسماند کششی در مرکز ناحیه جوش و در راستای محیطی سطح خارجی لوله به وجود آمده است. همچنین حداکثر تنش پسماند فشاری در راستاهای محیطی و محوری سطح داخلی لوله در منطقه متاثر از حرارت جوش قرار دارد. تغییرات تنش های پسماند محیطی در سطوح داخل و خارج لوله رفتارهای مشابه داشته و این تنش ها در مرکز ناحیه جوش به صورت کششی است. اما با فاصله گرفتن از مرکز ناحیه جوش رفتاری متضاد (در سطح خارج لوله به صورت کششی و در سطح داخل آن به صورت فشاری) مشاهده شد.کلید واژگان: شبیه سازی, تنش های پسماند, جوش سر به سر لوله, کرنش سنجی سوراخ, لوله انتقال گازThe weld residual stresses decrease the design stress in gas transportation pipelines. In this paper, two X70 steel pipes of 56 inch outside diameter were firstly girth welded. Experimental hole drilling test was conducted to evaluate the residual stress distribution in this joint. Then, the finite element simulation of the welding process was performed to evaluate the residual stress distribution precisely. The numerical results were verified by comparison with the obtained experimental measurements. The qualitative results achieved match properly with the experimental results. Simulation results (with a difference about 15% compared to experimental results) evaluated the maximum residual stress in hoop direction of pipe’s external weld metal. The experimental data showed that the maximum tensile residual stress was located on the center line of the weld gap on the pipe outer surface alongside with the pipe hoop direction. Moreover, the maximum compressive (hoop and axial) residual stresses occurred on the pipe inner surface in heat affected zone. The variations of the hoop residual stresses on the inner and outer surfaces of the pipe had similar trend with tensile distribution at the center line of the weld gap. However, these stresses showed different trends (tensile stress on the outer surface, and compressive stress on the inner surface) with distancing from the weld center line.Keywords: Simulation , Residual stress , Hole Drilling , Multi pass girth weld , Gas pipeline
-
When a natural gas pipeline ruptures, the adjacent automatic line control valves (ALCVs) should close quickly to prevent leakage or explosion. The differential pressure set point (DPS) at each valve location is the main criteria for value setting in ALCV action. If the DPS is not properly adjusted, the ALCV may mistakenly close or it may not take any action at proper time. This study focused on the DPS values prediction for setting ALCV installed on a gas pipeline with 1mm orifice diameter. The effect of characteristic parameters such as pipeline operational pressure (POP) and pipeline pressure drop rate (ROD) due to rupture or major leak was experimentally investigated on DPS. Twenty-five different conditions with double set of typical mentioned characteristic parameters were chosen. For each condition, the differential pressure (DP) was measured over 180 seconds by analyzing the experimental values. Therefore, 25 maximum DP values (DPSs) were obtained. The DPS increases by increase in ROD or decreasing POP parameters. Because of using nitrogen gas instead of natural gas due to safety reasons, the DPS results can be practically applied by adding a safety factor of 15%. The diagram of DPS with respect to ROD and non-dimensional DPS (DOP) versus non-dimensional ROD (RTP) was provided for different POPs.Keywords: Automatic Control Valve, Operating Pressure, Adjusting Values, Pressure Drop Rate, Gas Pipeline
-
اندازه گیری چقرمگی شکست KIC بعنوان معیار شروع ترک (در نوع اول تنش در حالت کرنش صفحه ای) برای رسم منحنی های ارزیابی آسیب (FAD) برای فولادهای پراستحکام خطوط انتقال گاز اهمیت ویژه ای دارد. در تحقیق حاضر از آزمون خمش سه نقطه ای بر روی شش نمونه آزمایشگاهی استاندارد برای محاسبه چقرمگی شکست لوله فولادی API X65 (با قطر خارجی 1219 میلی متر و ضخامت جداره 3/14 میلی متر) استفاده شده است. از نتایج حاصله، مقدار چقرمگی شکست برابر MPa√m 308 بدست آمد. سپس حل المان محدود نمونه خمش سه نقطه ای با استفاده از تئوری اصلاح شده شکست گرسون انجام شد. از طرح المان محدود، میزان گشودگی دهانه ترک در حالت بحرانی (δ) محاسبه شد. از این طرح با استفاده از روابط مکانیک شکست ارتجاعی خطی مقدار چقرمگی شکست MPa√m 297 متناظر با مقدار δ بدست آمد. مقایسه نتایج تجربی و عددی نشان دهنده برابری مناسب مقادیر بدست آمده از روش های مرسوم تعیین چقرمگی شکست در دمای محیط است.کلید واژگان: چقرمگی شکست, لوله انتقال گاز, رشد ترک نرم, مدل گرسون, آزمون خمش سه نقطه ای, API X65Plain strain fracture toughness is extremely important for failure assessment in high strength steels used in gas transmission pipelines. In the current research an experimental method based on three point bend test specimens was used to calculate fracture toughness of steel pipes of grade API X65 (with outer diameter of 1219 mm and wall thickness of 14.3 mm). A value of 308 MPa√m was found for fracture toughness of test material. Then a finite element solution for the test specimens was conducted using modified Gursons damage theory. In this model, the critical crack tip opening displacement was calculated from damaged elements in the early stage of crack growth. This model resulted in 297 MPa√m of fracture toughness. A comparison between the experimental and numerical results illustrated the fitness of common methods for determining fracture toughness of tested steel in ambient temperature. Due to the lack of experimental data for this steel, the obtained results can be used for safe performance of domestic gas pipelines made from API X65 steel.Keywords: Fracture toughness, Gas Pipeline, Ductile Crack Growth, Gurson Model, Three Point Bend, API X65
-
اندازه گیری چقرمگی شکست KIC بعنوان معیار شروع ترک (در نوع اول تنش در حالت کرنش صفحه ای) برای رسم منحنی های ارزیابی آسیب (FAD) برای فولادهای پراستحکام خطوط انتقال گاز اهمیت ویژه ای دارد. در تحقیق حاضر از آزمون خمش سه نقطه ای بر روی شش نمونه آزمایشگاهی استاندارد برای محاسبه چقرمگی شکست لوله فولادی API X65 (با قطر خارجی 1219 میلی متر و ضخامت جداره 3/14 میلی متر) استفاده شده است. از نتایج حاصله، مقدار چقرمگی شکست برابر MPa√m 308 بدست آمد. سپس حل المان محدود نمونه خمش سه نقطه ای با استفاده از تئوری اصلاح شده شکست گرسون انجام شد. از طرح المان محدود، میزان گشودگی دهانه ترک در حالت بحرانی (δ) محاسبه شد. از این طرح با استفاده از روابط مکانیک شکست ارتجاعی خطی مقدار چقرمگی شکست MPa√m 297 متناظر با مقدار δ بدست آمد. مقایسه نتایج تجربی و عددی نشان دهنده برابری مناسب مقادیر بدست آمده از روش های مرسوم تعیین چقرمگی شکست در دمای محیط است.کلید واژگان: چقرمگی شکست, لوله انتقال گاز, رشد ترک نرم, مدل گرسون, آزمون خمش سه نقطه ای, API X65Plain strain fracture toughness is extremely important for failure assessment in high strength steels used in gas transmission pipelines. In the current research an experimental method based on three point bend test specimens was used to calculate fracture toughness of steel pipes of grade API X65 (with outer diameter of 1219 mm and wall thickness of 14.3 mm). A value of 308 MPa√m was found for fracture toughness of test material. Then a finite element solution for the test specimens was conducted using modified Gursons damage theory. In this model¡ the critical crack tip opening displacement was calculated from damaged elements in the early stage of crack growth. This model resulted in 297 MPa√m of fracture toughness. A comparison between the experimental and numerical results illustrated the fitness of common methods for determining fracture toughness of tested steel in ambient temperature. Due to the lack of experimental data for this steel¡ the obtained results can be used for safe performance of domestic gas pipelines made from API X65 steel.Keywords: Fracture toughness, Gas Pipeline, Ductile Crack Growth, Gurson Model, Three Point Bend, API X65
-
روش تفاضل محدود ضمنی استیگر- وارمینگ، برای شبیه سازی عددی جریان گذرای غیرهم دمای گاز مورد استفاده قرار گرفته است. با توجه به غیرخطی بودن معادلات، در ابتدا خطی سازی حول هر گام زمانی انجام شده است که این امر ضمن کاهش تلاش محاسباتی جریان مورد نظر را نیز با سهولت بیشتری تحلیل می نماید. به منظور اعتبارسنجی و ارزیابی دقت روش حل عددی، ابتدا جریان فانو و جریان درون یک لوله شوک مورد بررسی قرار گرفته است. سپس جریان گذرا درون یک خط لوله گاز که فشار در ورودی آن به صورت تابعی از زمان تغییر می کند، شبیه سازی شده است. نتایج شبیه سازی عددی حاضر نشان می دهد که روش عددی تفاضل محدود خطی شده مورد استفاده در این مطالعه، تغییرات ناگهانی در شرایط جریان را به خوبی روش حل غیرخطی نمایان می سازد. همچنین روش عددی حاضر قادر است جریان های گذرای گاز را با دقت مناسب در مقایسه روش حل غیرخطی با تلاش محاسباتی کمتر تحلیل نماید.
کلیدواژگان جریان غیرهم دما؛ جریان گذرا؛ خطوط لوله گاز؛ روش های خطی شده؛ گاز حقیقیNumerical simulation of non-isothermal transient gas flow is performed using implicit Steger-Warming finite difference method. Because of nonlinearity of the governing equations، they are linearized at each time step. The linearization either reduces computational effort or analyzes the flowfield more conveniently. In order to validate and evaluate the accuracy of current numerical method، Fanno and shock tube flows are investigated first. Then، transient flow in a gas pipeline that its inlet pressure changes with time is simulated. The results of present study show that Steger-Warming finite difference scheme can well captured the sudden changes in the flowfield. Moreover، the present method is able to analyze transient gas flows as nearly accurate as the nonlinear one with less computational effort.Keywords: non, isothermal flow, transient flow, gas pipeline, linearized schemes, real gas -
در بررسی و مدل سازی شبکه خطوط انتقال گاز دو نوع جریان پایا و گذرا وجود دارد. مدل سازی جریان گذرا نسبت به جریان پایا پیچیده تر و مشکلتر است ولی به خوبی می تواند جوابگوی پیش بینی شبکه انتقال در حالات مختلف باشد. شبکه خطوط انتقال گاز شامل خط لوله، ایستگاه تقویت فشار، ایستگاه تقلیل فشار، شیرهای مسدودکننده و میادین ذخیره است که اثر هر یک در رفتار گذرای شبکه نسبت به دیگری کاملا متفاوت است. در این تحقیق خط لوله بعنوان مهمترین پارامتر تاثیرگذار در رفتار گذرای شبکه انتقال گاز برای حالت گذرای همدما مدل سازی و بررسی شده است. مدل سازی خط لوله با استفاده از معادلات بقای جرم، ممنتوم و معادله حالت، منتهی به معادلات دیفرانسیلی جزئی می شود که باید نسبت به زمان و مکان حل شوند. در این تحقیق سه گزینه از شرایط عملیاتی خط لوله انتقال گاز در حالت گذرا مورد بررسی قرار گرفته است که نتایج مدل سازی هر یک بصورت منحنی های تغییرات دبی جریان و فشار در طول خط لوله برای زمانهای مختلف رسم و تفسیر شده اند.
کلید واژگان: خط لوله, جریان گذرا, مدل سازیIn case of studying and modeling gas flow in the pipeline network, there are two types of flow: steady flow and transient flow. The natural gas pipeline network consists of pipelines, pressure compressor stations, pressure reduction stations, line break valves and storage fields. The pipeline is the key element in determining the overall network dynamic behavior. Dynamic behavior of a gas pipeline due to changes of gas flow rate for isothermal transient status has been modeled and investigated in this study. The modeling of isothermal transient flow is resulted in differential equations using mass conservation equations, momentum and an equation of state. By solving these equations according to time and place, the behavior of gas pipeline will be obtained. Three cases of the operating conditions for gas transmission pipelines in transient status have been studied. The modeling results of each case have been shown and interpreted in curves of flow rate and pressure variations along pipeline in different time spots.Keywords: gas pipeline, transient flow, modeling
- نتایج بر اساس تاریخ انتشار مرتب شدهاند.
- کلیدواژه مورد نظر شما تنها در فیلد کلیدواژگان مقالات جستجو شدهاست. به منظور حذف نتایج غیر مرتبط، جستجو تنها در مقالات مجلاتی انجام شده که با مجله ماخذ هم موضوع هستند.
- در صورتی که میخواهید جستجو را در همه موضوعات و با شرایط دیگر تکرار کنید به صفحه جستجوی پیشرفته مجلات مراجعه کنید.