‫ﺗﺤﻠﻴﻞ و ﺑﺮرﺳﻲ ﭘﺪﻳﺪه ﺿﺮﺑﻪ ﻗﻮچ ﺑﺎ اﺳﺘﻔﺎده از روش ﻣﺸﺨﺼﻪ‪،‬ﻣﻄﺎﻟﻌﻪ‬

‫ﻣﻮردي اﻳﺴﺘﮕﺎه ﭘﻤﭙﺎژ ﺑﺨﺶ اﻳﺮي‬

‫‪ -1‬دﻛﺘﺮ اﺑﻮاﻟﻔﻀﻞ ﺷﻤﺴﺎﻳﻲ )اﺳﺘﺎد داﻧﺸﮕﺎه ﺻﻨﻌﺘﻲ ﺷﺮﻳﻒ ﮔﺮوه ﻋﻤﺮان(‬

‫‪021-66014828‬‬

‫‪ -2‬ﺳﻴﻨﺎ ﻓﺮد ﻣﺮادي ﻧﻴﺎ ) ﻋﻀﻮ ﻫﻴﺌﺖ ﻋﻠﻤﻲ ﮔﺮوه ﻋﻤﺮان داﻧﺸﮕﺎه آزاداﺳﻼﻣﻲ واﺣﺪ ﺗﺒﺮﻳﺰ(‬
‫‪[email protected]‬‬
‫‪09143142619---0411-4799365‬‬

‫ﭼﻜﻴﺪه‪:‬‬
‫ﭘﺪﻳﺪه ﺿﺮﺑﻪ ﻗﻮچ از ﻣﻬﻤﺘﺮﻳﻦ ﭘﺪﻳﺪه ﻫﺎﻳﻲ اﺳﺖ ﻛﻪ در اﻛﺜﺮ اﻳﺴﺘﮕﺎه ﻫﺎي ﭘﻤﭙﺎژ و ﺧﻄﻮط اﻧﺘﻘﺎل اﺗﻔﺎق ﻣﻲ اﻓﺘﺪ‪.‬از‬

‫آﻧﺠﺎﻳﻴﻜﻪ اﻳﻦ ﭘﺪﻳﺪه ﻋﺎﻣﻞ اﺻﻠﻲ ﺗﺨﺮﻳﺐ در اﻳﺴﺘﮕﺎﻫﻬﺎي ﭘﻤﭙﺎژ وﺧﻄﻮط اﻧﺘﻘﺎل ﻣﻲ ﺑﺎﺷﺪ‪ ،‬ﺑﺮرﺳﻲ آن ﺟﺰء ﻣﺮاﺣﻞ اوﻟﻴﻪ‬

‫ﻃﺮاﺣﻲ ﻣﺤﺴﻮب ﻣﻲ ﮔﺮدد‪ .‬در اﻳﻦ ﻣﻘﺎﻟﻪ ﺑﺮاي ﺗﺠﺰﻳﻪ و ﺗﺤﻠﻴﻞ ﺟﺮﻳﺎﻧﻬﺎي ﻧﺎﭘﺎﻳﺪاراز روش ﺿﺮﺑﺖ ﻗﻮچ ﻳﺎ ﺗﺌﻮري اﻻﺳﺘﻴﺴﻴﺘﻪ‬

‫اﺳﺘﻔﺎده ﺷﺪه اﺳﺖ ﻛﻪ در آن ﺧﻮاص اﻻﺳﺘﻴﻜﻲ ﺧﻂ ﻟﻮﻟﻪ و ﺳﻴﺎل ﻫﺮ دو وارد ﻣﺤﺎﺳﺒﺎت ﻣﻲ ﮔﺮدﻧﺪ ‪.‬ﻣﻌﺎدﻻت ﺗﺸﺮﻳﺢ ﻛﻨﻨﺪه‬

‫ﺟﺮﻳﺎن ﺳﻴﺎل در اﻳﻦ روش دو ﻣﻌﺎدﻟﺔ دﻳﻔﺮاﻧﺴﻴﻞ ﻏﻴﺮﺧﻄﻲ ﻣﻲ ﺑﺎﺷﻨﺪ ﻛﻪ ﺣﻞ آﻧﻬﺎ ﺑﺴﻴﺎر ﻣﺸﻜﻞ و ﭘﻴﭽﻴﺪه اﺳﺖ ‪.‬ﻟﻴﻜﻦ در ﻋﻮض‬

‫ﺟﻮاﺑﻬﺎي ﺣﺎﺻﻠﻪ ﺑﺴﻴﺎر دﻗﻴﻘﺘﺮ و ﻣﻄﻤﺌﻦ ﺗﺮ ﺧﻮاﻫﺪ ﺑﻮد‪.‬ﺑﺎ اﺳﺘﻔﺎده از روش ﻣﺸﺨﺼﻪ و ﺑﻜﺎرﮔﻴﺮي زﺑﺎن ﺑﺮﻧﺎﻣﻪ ﻧﻮﻳﺴﻲ ‪TURBO‬‬

‫‪C‬ﻣﻌﺎدﻻت دﻳﻔﺮاﻧﺴﻴﻞ ﻏﻴﺮﺧﻄﻲ ﻣﺮﺑﻮﻃﻪ ﺣﻞ ﮔﺮدﻳﺪه و ﺟﻮاﺑﻬﺎي ﻣﻨﺎﺳﺐ ﺑﺪﺳﺖ آﻣﺪه اﺳﺖ‪ .‬ﺑﺎ اﺳﺘﻔﺎده از روش ﻣﺸﺨﺼﻪ اﺑﺘﺪا‬

‫ﻣﻌﺎدﻻت دﻳﻔﺮاﻧﺴﻴﻞ ﺟﺰﺋﻲ را ﺑﻪ ﻣﻌﺎدﻻت دﻳﻔﺮاﻧﺴﻴﻞ ﻣﻌﻤﻮﻟﻲ ﺗﺒﺪﻳﻞ ‪ ،‬و ﺳﭙﺲ ﺑﺎ روش ﺗﻔﺎوت ﻣﺤﺪود ﮔﻮﻳﺎﺣﻞ ﮔﺮدﻳﺪه اﺳﺖ‪.‬‬

‫ﺑﻌﻠﺖ اﻳﻨﻜﻪ ﻫﺮ ﻳﻚ از ﺷﺮاﻳﻂ ﻣﺮزي و ﻫﺮ ﻳﻚ از ﻣﻘﺎﻃﻊ ﻟﻮﻟﻪ در ﻫﺮ ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ﺑﻄﻮر ﺟﺪاﮔﺎﻧﻪ ﻣﻮرد ﺗﺤﻠﻴﻞ ﻗﺮار ﻣﻲ ﮔﻴﺮﻧﺪ ‪.‬‬

‫اﻳﻦ روش ﺑﺮاي ﺳﻴﺴﺘﻤﻬﺎي داراي ﺷﺮاﻳﻂ ﻣﺮزي ﭘﻴﭽﻴﺪه ﻣﻨﺎﺳﺒﺘﺮ اﺳﺖ‪.‬دراﻳﻦ ﻣﻘﺎﻟﻪ ﻫﻤﭽﻨﻴﻦ ﭘﺪﻳﺪه ﺿﺮﺑﻪ ﻗﻮچ در اﻳﺴﺘﮕﺎه‬

‫ﭘﻤﭙﺎژ ﺑﺨﺶ اﻳﺮي ﻣﻮرد ﻣﻄﺎﻟﻌﻪ ﻗﺮار ﮔﺮﻓﺘﻪ و ﻣﻨﺤﻨﻲ ﻫﺎي ﻣﺮﺑﻮﻃﻪ اراﺋﻪ ﮔﺮدﻳﺪه اﺳﺖ‪.‬‬

‫ﻫﺎ‪ :‬ﺿﺮﺑﻪ ﻗﻮچ‪-‬روش ﺗﻔﺎوت ﻣﺤﺪود ﮔﻮﻳﺎ‪ -‬روش ﻣﺸﺨﺼﻪ‪-‬ﺟﺮﻳﺎن ﻣﻴﺮا‪-‬ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭗ‬ ‫ﻛﻠﻴﺪ واژه‬

‫ﺑﺮاي ﺑﺪﺳﺖ آوردن ﻣﻌﺎدﻟﻪ ﺟﺮﻳﺎن ﻏﻴﺮ داﺋﻤﻲ در داﺧﻞ ﻣﺠﺎري ﺑﺴﺘﻪ ﻧﻴﺎزﻣﻨﺪ ﻣﻌﺎدﻻت ﭘﻴﻮﺳﺘﮕﻲ و‬

‫دﻳﻨﺎﻣﻴﻚ ﻫﺴﺘﻴﻢ‪.‬ﺑﺮاي ﺗﻌﻴﻴﻦ ﻣﻌﺎدﻻت ﺟﺮﻳﺎن ﻫﺎي ﻣﻴﺮا در ﺧﻄﻮط اﻧﺘﻘﺎل ﻧﻘﻠﻲ‪ ،‬ﺗﺤﺖ ﭘﻤﭙﺎژ ﻣﻔﺮوﺿﺎت‬

‫زﻳﺮﺿﺮوري ﻣﻲ ﺑﺎﺷﻨﺪ‪[١]:‬‬

‫‪١‬‬
 ‫‪1‬ـ ﺟﺪاره ﻣﺠﺮا و ﺧﻮد ﺳﻴﺎل اﻻﺳﺘﻴﻚ ﺧﻄﻲ ﻣﻲ ﺑﺎﺷﻨﺪ‪.‬‬

‫‪2‬ـ رواﺑﻂ ﻣﺮﺑﻮط ﺑﻪ ﻣﺤﺎﺳﺒﻪ اﻓﺘﻬﺎي اﺻﻄﻜﺎﻛﻲ در داﺧﻞ ﻟﻮﻟﻪ در ﺣﺎﻟﺖ ﺟﺮﻳﺎن داﺋﻤﻲ ﺑﺮاي ﺟﺮﻳﺎن‬

‫ﻣﻴﺮا ﻧﻴﺰ ﻣﻮرد ﻗﺒﻮل ﻫﺴﺘﻨﺪ‪.‬‬

‫‪3‬ـ ﺟﺮﻳﺎن ﺳﻴﺎل در ﻣﺠﺮاي ﻳﻚ ﺑﻌﺪي ﻣﻲ ﺑﺎﺷﺪ‪.‬‬

‫‪ -1‬ﻣﻌﺎدﻟﻪ دﻳﻨﺎﻣﻴﻚ‪:‬‬

‫ﻣﻌﺎدﻟﻪ دﻳﻨﺎﻣﻴﻚ ﺣﺎﻛﻢ ﺑﺮ ﭘﺪﻳﺪه راﺑﻄﻪ ‪ L1‬ﺑﺎدو ﻣﺘﻐﻴﺮ ‪ x‬و ‪ t‬ﻣﻲ ﺑﺎﺷﺪ‪. .‬‬

‫‪∂v 1 ∂p‬‬ ‫‪Fv‬‬

‫= ‪L1‬‬ ‫‪+ . + g sin θ +‬‬ ‫‪=0‬‬
‫‪∂t ρ ∂x‬‬ ‫‪2D‬‬

‫‪ -2‬ﻣﻌﺎدﻟﻪ ﭘﻴﻮﺳﺘﮕﻲ‪:‬‬

‫دوﻣﻴﻦ ﻣﻌﺎدﻟﻪ اي ﻛﻪ ‪) H‬ارﺗﻔﺎع ﻧﻈﻴﺮ اﻧﺮژي( و ‪) V‬ﺳﺮﻋﺖ ﺟﺮﻳﺎن( را ﺑﻪ ﻫﻤﺪﻳﮕﺮ ﻣﺮﺑﻮط ﻣﻲ ﺳﺎزد‬

‫ﻣﻌﺎدﻟﻪ ﭘﻴﻮﺳﺘﮕﻲ ﻣﻲ ﺑﺎﺷﺪ ﻛﻪ ﺑﺎ ‪ L2‬ﻧﺸﺎن داده ﻣﻲ ﺷﻮد‪.‬‬

‫‪∂p‬‬ ‫‪∂v‬‬
‫= ‪L2‬‬ ‫‪+ ρ a2 + = 0‬‬
‫‪∂t‬‬ ‫‪∂x‬‬

‫دو ﻣﻌﺎدﻟﻪ ‪ L1‬و ‪ L2‬ﺟﺰ ﻣﻌﺎدﻻت دﻳﻔﺮاﻧﺴﻴﻞ ﺟﺰﺋﻲ ﻣﺮﺗﺒﻪ اول ﻫﺴﺘﻨﺪو در اﻳﻦ ﻣﻌﺎدﻻت دو ﻣﺘﻐﻴﺮ‬

‫ﻫﺴﺘﻨﺪ ﻛﻪ‬
‫‪ρ‬‬ ‫ﻣﺴﺘﻘﻞ ‪ t‬و ‪ x‬و دو ﻣﺘﻐﻴﺮ ﻏﻴﺮ ﻣﺴﺘﻘﻞ ‪ P‬و ‪ V‬وﺟﻮد دارد‪ .‬دﻳﮕﺮ ﻣﺘﻐﻴﺮ ﻫﺎ ‪ D‬و ‪ A‬و‬

‫ﻣﺸﺨﺼﺎت ﻣﺠﺮاي ﻋﺒﻮر ﻣﺎﻳﻊ و ﻣﺸﺨﺼﻪ ﺧﻮد ﻣﺎﻳﻊ اﺳﺖ و ﺑﺴﺘﮕﻲ ﺑﻪ زﻣﺎن ﻧﺪارﻧﺪ‪ ،‬ﻟﻜﻦ ﻣﻤﻜﻦ اﺳﺖ ﺗﺎﺑﻊ ‪x‬‬

‫ﺑﺎﺷﻨﺪ‪ .‬ﺿﺮﻳﺐ اﺻﻄﻜﺎك ‪ f‬ﺗﺎﺑﻊ ﻋﺪد رﻳﻨﻮﻟﺪز اﺳﺖ‪ .‬ﺑﺎ وﺟﻮد اﻳﻦ در ﻣﺴﺎﻳﻞ ﻣﻘﺪار ‪ f‬را ﺛﺎﺑﺖ ﻓﺮض ﻣﻲ ﻛﻨﻴﻢ‬

‫زﻳﺮا اﺛﺮ ﺗﻐﻴﻴﺮات ‪ f‬در ﺟﺮﻳﺎن ﻫﺎي ﻣﻴﺮا ﻗﺎﺑﻞ اﻏﻤﺎض اﺳﺖ‪ .‬زﻣﺎﻧﻲ ﻛﻪ ﺷﺮاﻳﻂ ﺟﺮﻳﺎن از ﻳﻚ ﺣﺎﻟﺖ داﻳﻤﻲ ﺑﻪ‬

‫ﺣﺎﻟﺖ داﻳﻤﻲ دﻳﮕﺮ ﺗﻐﻴﻴﺮ ﻛﻨﺪ‪،‬ﺷﺮاﻳﻂ ﺟﺮﻳﺎن ﻣﺎ ﺑﻴﻦ اﻳﻦ دو ﺣﺎﻟﺖ را ﺟﺮﻳﺎن ﻣﻴﺮا ‪(TRANSIENT‬‬

‫)‪FLOW‬ﻣﻲ ﻧﺎﻣﻨﺪ‪.‬‬

‫ﺑﺮاي ﻧﻤﺎﻳﺶ رﻳﺎﺿﻲ ﭘﻤﭗ ﺑﺎﻳﺴﺘﻲ ﭼﻬﺎر ﻣﺘﻐﻴﺮ دﺑﻲ)‪ ،(Q‬ارﺗﻔﺎع ﭘﻤﭙﺎژ )‪ ،(H‬ﺳﺮﻋﺖ دوراﻧﻲ )‪ (N‬و‬

‫ﮔﺸﺘﺎور ﭼﺮﺧﺸﻲ )‪ (T‬ﻣﺸﺨﺺ ﺷﻮﻧﺪ‪ .‬ﻣﻨﺤﻨﻲ ﻫﺎﻳﻲ ﻛﻪ ﻧﺸﺎن دﻫﻨﺪة راﺑﻄﻪ اي ﻣﺎ ﺑﻴﻦ اﻳﻦ ﻣﺘﻐﻴﺮﻫﺎ ﺑﺎﺷﺪ ﺑﻪ‬

‫ﻣﻨﺤﻨﻲ ﻫﺎي ﻣﺸﺨﺼﻪ ﭘﻤﭗ )‪ (Pump characteristic curve‬ﻣﻮﺳﻮم ﻫﺴﺘﻨﺪ‪.‬‬

‫‪٢‬‬
 ‫اﮔﺮ ﭼﻨﺎن ﭼﻪ اﻃﻼﻋﺎت ﻣﺮﺑﻮط ﺑﻪ ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭗ در ﻳﻚ ﻳﺎ ﭼﻨﺪ ﻧﺎﺣﻴﻪ از ﻧﻮاﺣﻲ ﺳﻪ ﮔﺎﻧﻪ‬

‫ﻛﺎرﻛﺮد‪ .‬ﭘﻤﭙﻲ‪ ،‬اﺳﺘﻬﻼك اﻧﺮژي و ﻋﻤﻠﻜﺮد ﺗﻮرﺑﻴﻨﻲ در دﺳﺘﺮس ﻧﺒﺎﺷﺪ ﻣﻲ ﺗﻮان از ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭙﻲ‬

‫ﻛﻪ ﺗﻘﺮﻳﺒﺎً داراي ﺳﺮﻋﺖ ﻣﺨﺼﻮص ﻣﺴﺎوي ﺑﺎ ﭘﻤﭗ ﻣﻮرد ﻧﻈﺮ اﺳﺖ اﺳﺘﻔﺎده ﻛﺮد‪ .‬اﻃﻼﻋﺎت در ﻣﻮرد ﻣﻨﺤﻨﻲ‬

‫ﻣﺸﺨﺼﺔ ﭘﻤﭗ اﺻﻠﻲ را ﻣﻲ ﺗﻮان ﺑﺎ اﺳﺘﻔﺎده از ﻧﺘﺎﻳﺞ ﺣﺎﺻﻞ از آزﻣﺎﻳﺶ ﺑﺮ روي ﻣﺪﻟﻬﺎ و ﺑﺎ اﺳﺘﻔﺎده از رواﺑﻂ‬

‫ﻫﻢ ﺧﺎﻧﻮاده و ﺗﺸﺎﺑﻪ ﭘﻴﺪا ﻛﺮد‪[٢].‬‬

‫اﮔﺮ اﻧﺪﻳﺲ ‪ R‬ﺷﺮاﻳﻂ ﺟﺮﻳﺎن داﻳﻤﻲ ﻳﺎ ﭘﺎﻳﺪار را ﻧﺸﺎن دﻫﺪ ﻣﻲ ﺗﻮان رواﺑﻂ زﻳﺮ را ﺗﻌﺮﻳﻒ ﻛﺮد‪.‬‬

‫)‪(1‬‬
‫‪N‬‬ ‫‪H‬‬
‫=‪α‬‬ ‫‪,‬‬ ‫=‪h‬‬
‫‪NR‬‬ ‫‪HR‬‬
‫‪Q‬‬ ‫‪T‬‬
‫=‪v‬‬ ‫‪,‬‬ ‫=‪β‬‬
‫‪QR‬‬ ‫‪TR‬‬

‫ﻛﻪ رواﺑﻂ اﺧﻴﺮ ﻧﺴﺒﺘﻬﺎي ﺳﺮﻋﺖ دوراﻧﻲ‪ ،‬دﺑﻲ‪ ،‬ارﺗﻔﺎع ﭘﻤﭙﺎژ‪ ،‬ﮔﺸﺘﺎور ﭘﻤﭗ ﺑﻪ ﻣﻘﺎدﻳﺮ ﻧﻈﻴﺮﺷﺎن در‬

‫ﺣﺎﻟﺖ ﺟﺮﻳﺎن داﻳﻤﻲ را ﻧﺸﺎن ﻣﻲ دﻫﻨﺪ‪.‬‬

‫‪H‬‬ ‫‪N‬‬
‫‪N2‬‬ ‫ﺛﺎﺑﺖ =‬ ‫‪Q‬‬ ‫ﺛﺎﺑﺖ =‬

‫‪h‬‬ ‫‪H N R2‬‬

‫=‬ ‫‪.‬‬
‫‪α2‬‬ ‫‪N 2 HR‬‬ ‫ﻣﻘﺪار ﺛﺎﺑﺖ =‬ ‫)‪(2‬‬
‫‪α N QR‬‬
‫‪= .‬‬
‫ﻣﻘﺪار ﺛﺎﺑﺖ = ‪v Q N R‬‬ ‫)‪(3‬‬

‫ﭘﺎراﻣﺘﺮ ‪ θ‬را ﭼﻨﻴﻦ ﺗﻌﺮﻳﻒ ﻣﻲ ﻛﻨﻴﻢ‪.‬‬

‫)‪(4‬‬
‫‪α‬‬ ‫‪α‬‬
‫= ‪tan θ‬‬ ‫‪Q = Arc tan‬‬
‫‪v‬‬ ‫‪v‬‬

‫‪٣‬‬
 ‫ﺑﻪ دﻟﻴﻞ اﻳﻨﻜﻪ درﺟﺮﻳﺎن ﻣﻴﺮا اﺣﺘﻤﺎل ﺻﻔﺮ ﺷﺪن ‪ α‬ﻳﺎ ‪ v‬ﺑﻪ ﺗﻨﻬﺎﻳﻲ وﺟﻮد دارد ﺑﻨﺎﺑﺮاﻳﻦ ﺑﻪ ﺟﺎي‬
‫‪β‬‬ ‫‪h‬‬ ‫‪h‬‬ ‫‪h‬‬
‫‪ α + v‬اﺳﺘﻔﺎده ﺧﻮاﻫﺪ ﺷﺪ‪. .‬ﻫﻤﺎﻧﻨﺪ ﭘﺎراﻣﺘﺮ ‪ ، α + v‬ﭘﺎراﻣﺘﺮ ‪ α + v‬ﻧﻴﺰ در ﻣﺤﺎﺳﺒﺎت‬ ‫‪α2‬‬
‫‪2‬‬ ‫‪2‬‬ ‫‪2‬‬ ‫‪2‬‬ ‫‪2‬‬ ‫‪2‬‬

‫‪h‬‬
‫= ‪Fh‬‬
‫و‬ ‫ﺑﻜﺎر ﺧﻮاﻫﺪ رﻓﺖ ﺑﺎ اﻳﻦ اوﺻﺎف ﻣﻲ ﺗﻮان ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﻣﻘﺎدﻳﺮ ‪ θ‬را ﺑﺮ ﺣﺴﺐ ‪α + v2‬‬
‫‪2‬‬

‫‪β‬‬
‫= ‪Fβ‬‬
‫ﻫﻤﺮاه ﺑﺎ ﺗﻐﻴﻴﺮات ﺳﺮﻋﺖ ﻣﺨﺼﻮص ‪ N S‬رﺳﻢ ﻧﻤﻮد‪.‬ﺑﺮاي ﺗﻌﺮﻳﻒ‬ ‫ﻫﻤﭽﻨﻴﻦ ‪ θ‬ﺑﺮ ﺣﺴﺐ ‪α + v2‬‬
‫‪2‬‬

‫‪N R QR‬‬
‫= ‪NS‬‬
‫ﺳﺮﻋﺖ ﻣﺨﺼﻮص ﻣﻲ ﺗﻮان ﻧﻮﺷﺖ‪( HR ) 3/ 4 :‬‬

‫‪ -3‬ﺷﺮاﻳﻂ ﻣﺮزي ﺑﺮاي از ﺣﺮﻛﺖ اﻓﺘﺎدن ﭘﻤﭗ‬

‫‪ -1-3‬ﻣﻌﺎدﻻت ﺷﺮاﻳﻂ ﻣﺮزي اﻋﻤﺎل ﺷﺪه از ﻃﺮف ﭘﻤﭗ‬

‫ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭗ را ﻣﻲ ﺗﻮان ﺑﻄﻮر رﻳﺎﺿﻲ ﺑﻮﺳﻴﻠﺔ ﻣﻨﺤﻨﻲ ﻫﺎي ‪ Fβ , Fh‬ﺑﺮ ﺣﺴﺐ ‪ θ‬ﻧﺸﺎن داد‪.‬‬

‫ﺑﺮاي اﺳﺘﻔﺎده از اﻳﻦ ﻣﻨﺤﻨﻲ ﻫﺎ‪ ،‬ﻣﻘﺎدﻳﺮ ‪ Fβ , Fh‬در ﻓﻮاﺻﻞ ﻣﺴﺎوي ‪ ) θ‬ﭘﻨﺞ درﺟﻪ( در ﻛﺎﻣﭙﻴﻮﺗﺮ ذﺧﻴﺮه ﻣﻲ‬

‫ﺷﻮﻧﺪ‪ .‬ﻫﺮ ﻗﺴﻤﺖ از اﻳﻦ ﻣﻨﺤﻨﻲ ﻫﺎ ﻣﺎ ﺑﻴﻦ ﻧﻘﺎﻃﻲ ﻛﻪ در ﻛﺎﻣﭙﻴﻮﺗﺮ ذﺧﻴﺮه ﺷﺪه اﻧﺪ را ﻣﻲ ﺗﻮان ﺧﻂ راﺳﺖ‬

‫ﻓﺮض ﻛﺮد ﺑﺸﺮﻃﻲ ﻛﻪ ﺗﻌﺪاد ﻧﻘﺎط ﻛﺎﻓﻲ و ﺑﻪ ﺣﺪ ﻻزم ﻣﺪ ﻧﻈﺮ ﻗﺮا ر ﮔﻴﺮﻧﺪ ﻛﻪ ﺑﺎ اﻳﻦ ﻓﺮض ﺧﻄﺎي ﻧﺎﺷﻲ از‬

‫ﻓﺮض ﻛﺮدن ﻗﻄﻌﺎت ﻣﻨﺤﻨﻲ ﺑﺼﻮرت ﺧﻂ ﻣﺴﺘﻘﻴﻢ ﻗﺎﺑﻞ ﭼﻤﺸﭙﻮﺷﻲ ﻣﻲ ﺑﺎﺷﺪ‪[٣].‬‬

‫⎞‪⎛α‬‬
‫⎟ ⎜‪θ = Arc tan‬‬
‫ﻣﻘﺪار ‪ θ‬را ﻣﺎ ﺑﻴﻦ ) ‪ ( 0,π‬و ) ‪ ( −π ,0‬ﻣﺤﺎﺳﺒﻪ ﻣﻲ ﻛﻨﺪ در ﺻﻮرﺗﻴﻜﻪ‬ ‫ﻋﺒﺎرت ⎠ ‪⎝ v‬‬

‫ﻓﺎﺻﻠﻪ ﻣﻮرد ﻧﻈﺮ ﻣﺎ ﺑﻴﻦ ‪ 0,2π‬اﺳﺖ‪ .‬اﻳﻦ ﻣﺤﺪودﻳﺖ را ﭼﻨﻴﻦ ﻣﻲ ﺗﻮان از ﺑﻴﻦ ﺑﺮد ﻛﻪ در ﻫﻨﮕﺎم ﻣﻨﻔﻲ ﺷﺪن‬

‫‪ θ‬ﺑﻪ اﻧﺪازة ‪ 2π‬ﺑﻪ ﻣﻘﺪار ﻣﺤﺎﺳﺒﻪ ﺷﺪه ‪ θ‬ﻣﻲ اﻓﺰاﻳﻴﻢ‪ .‬اﻛﻨﻮن ﻓﺮض ﻣﻲ ﻛﻨﻴﻢ ﻛﻪ ﻣﺤﺎﺳﺒﺎت ﻣﺎ ﺑﻪ‪ I‬اﻣﻴﻦ‬

‫ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ رﺳﻴﺪه اﺳﺖ و ﻣﻘﺎدﻳﺮ ‪ α , v, β , h‬در آﻏﺎز اﻳﻦ ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ﻣﻌﻠﻮم اﺳﺖ و ﻫﺪف ﻣﺎ ﻳﺎﻓﺘﻦ‬

‫ﻣﻘﺎدﻳﺮ ﻣﺠﻬﻮل ‪ α P , vP , β P , hP‬در اﻧﺘﻬﺎي اﻳﻦ ﻣﺮﺣﻠﻪ زﻣﺎﻧﻲ ﻣﻲ ﺑﺎﺷﺪ‪.‬‬

‫‪٤‬‬
 ‫ﺑﺮاي ﺗﻌﻴﻴﻦ ﻣﺠﻬﻮﻻت ﻣﺮﺑﻮﻃﻪ اﺑﺘﺪا ﻣﻲ ﺑﺎﻳﺴﺖ ﻣﻌﺎدﻟﻪ ﻗﻄﻌﻪ اي از ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭗ را ﻛﻪ ﺑﻪ‬

‫‪ α P , vP‬ﻣﺮﺑﻮط ﻣﻲ ﮔﺮدد ﺣﺴﺎب ﻛﻨﻴﻢ‪ .‬اﻣﺎ ﭼﻮن ﻣﻘﺎدﻳﺮ اﻳﻦ ﻣﺘﻐﻴﺮﻫﺎ در اﺑﺘﺪا ﻣﺠﻬﻮﻟﻨﺪ ﻣﻲ ﺗﻮاﻧﻴﻢ ﺑﻌﻨﻮان‬

‫اوﻟﻴﻦ ﺗﻘﺮﻳﺐ ﻣﻘﺎدﻳﺮ آﻧﻬﺎ را از ﻃﺮﻳﻖ ﺑﺮوﻧﻴﺎﺑﻲ ﺑﺎ اﺳﺘﻔﺎده از ﻣﻘﺎدﻳﺮ ﻣﻌﻠﻮم زﻣﺎن ﻗﺒﻠﻲ ﺣﺴﺎب ﻧﻤﺎﻳﻴﻢ‪.‬‬

‫اﮔﺮ ‪ α e , ve‬ﻣﻘﺎدﻳﺮ ﺗﻘﺮﻳﺒﻲ در اﻧﺘﻬﺎي ‪ I‬اﻣﻴﻦ ﻓﺎﺻﻠﺔ زﻣﺎﻧﻲ ﺑﺎﺷﺪ و ﻫﻤﭽﻨﻴﻦ ‪ α I vI‬ﻣﻘﺎدﻳﺮ ﻣﻌﻠﻮم در اﺑﺘﺪاي ‪I‬‬

‫اﻣﻴﻦ ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ‪ ∆α I −1 , ∆VI −1‬ﻧﻴﺰ ﺑﻴﺎﻧﮕﺮ ﻧﻤﻮ اﻳﻦ ﻣﺘﻐﻴﺮﻫﺎ در )‪ ( I-1‬اُﻣﻴﻦ ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ﺑﺎﺷﺪ دارﻳﻢ‪:‬‬
‫‪α e = α I + ∆α I −1‬‬ ‫)‪(5‬‬
‫‪ve = vI + ∆vI −1‬‬ ‫)‪(6‬‬

‫اﮔﺮ ‪ ∆t‬ﻣﻘﺪار ﻛﻮﭼﻜﻲ ﺑﺎﺷﺪ اﻳﻦ ﺑﺮوﻧﻴﺎﺑﻲ‪،‬ﺗﺨﻤﻴﻨﻬﺎي ﻧﺰدﻳﻚ ﺑﻪ واﻗﻌﻴﺖ اراﻳﻪ ﺧﻮاﻫﺪ داد‪.‬‬
‫‪αe‬‬
‫‪Arc tan‬‬ ‫‪, hP‬‬
‫( را ﺑﺎ ﺗﻮﺟﻪ ﺑﻪ ﺷﻜﻞ)‪ (1‬ﻣﻲ ﻧﻮﻳﺴﻴﻢ‪ .‬ﺑﻪ‬ ‫‪ve‬‬ ‫ﻣﻌﺎدﻟﺔ ﻗﻄﻌﻪ ﺧﻂ ﻣﺴﺘﻘﻴﻢ ﮔﺬرﻧﺪه از ﻧﻘﻄﺔ )‬

‫⎛‬ ‫‪α‬‬ ‫⎞‬

‫⎟ ‪⎜ Arc tan P , β P‬‬
‫⎝‬ ‫‪vP‬‬ ‫⎠‬
‫را ﻧﻮﺷﺖ‪:‬‬ ‫ﻫﻤﻴﻦ ﺗﺮﺗﻴﺐ ﻣﻲ ﺗﻮان ﻣﻌﺎدﻟﺔ ﻗﻄﻌﻪ ﺧﻂ ﻣﺴﺘﻘﻴﻢ ﮔﺬرﻧﺪه از‬

‫ﺷﻜﻞ‪ -1-‬ﺗﻘﺮﻳﺐ زﻧﻲ ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭙﺒﻮﺳﻴﻠﻪ ﺧﻄﻮط ﻣﺴﺘﻘﻴﻢ‬

‫‪hP‬‬ ‫⎞‪⎛α p‬‬

‫⎟ ⎜ ‪= a1 + a 2 tan −1‬‬
‫‪α P + vP‬‬
‫‪2‬‬ ‫‪2‬‬
‫⎠ ‪⎝ vP‬‬
‫)‪(7‬‬

‫‪βP‬‬ ‫⎞ ‪⎛α‬‬
‫⎟ ‪= a 3 + a 4 tan1− ⎜ P‬‬
‫‪α P + vP‬‬
‫‪2‬‬ ‫‪2‬‬
‫⎠ ‪⎝ vP‬‬
‫)‪(8‬‬

‫‪٥‬‬
 ‫‪ a1 , a2 , a3 , a4‬ﻣﻘﺎدﻳﺮ ﺛﺎﺑﺖ ﻗﻄﻌﻪ ﺧﻄﻮط ﻣﺴﺘﻘﻴﻢ ﻧﻮﺷﺘﻪ ﺷﺪه ﻣﻲ ﺑﺎﺷﻨﺪ‪ .‬اﮔﺮ ﭼﻨﺎﻧﭽﻪ ‪ H P‬ﺑﻴﺎﻧﮕﺮ‬

‫ارﺗﻔﺎع ﻧﻈﻴﺮ اﻧﺮژي ﭘﻤﭙﺎژ در اﻧﺘﻬﺎي ﻓﺎﺻﻠﺔ زﻣﺎﻧﻲ و ‪ Hsuction‬ﻧﺸﺎﻧﮕﺮ ارﺗﻔﺎع ﺳﻄﺢ آزاد ﻣﺎﻳﻊ ﻣﺨﺰن ﻧﺴﺒﺖ ﺑﻪ‬

‫ﺧﻂ ﻣﺒﻨﺎ و ‪ ∆H PV‬ﻧﺸﺎن دﻫﻨﺪة اﻓﺖ ارﺗﻔﺎع ﻧﻈﻴﺮ اﻧﺮژي در ﻣﺤﻞ ﺷﻴﺮ ﺧﺮوﺟﻲ ﺑﺎﺷﺪ ﺑﺎ ﺗﻮﺟﻪ ﺑﻪ ﺷﻜﻞ )‪(2‬‬

‫ﻣﻲ ﺗﻮان ﻧﻮﺷﺖ‪.‬‬

‫‪H PI ,1 = HSUCTION + H P − ∆H PV‬‬
‫)‪(9‬‬

‫اﻓﺖ ارﺗﻔﺎع ﻧﻈﻴﺮ اﻧﺮژي در ﻣﺤﻞ ﺷﻴﺮ از ﻣﻌﺎدﻟﺔ زﻳﺮ ﺑﺪﺳﺖ ﻣﻲ آﻳﺪ‬

‫‪∆H PV = CV . QPI ,1 . QPI ,1‬‬

‫)‪(10‬‬

‫در ﻣﻌﺎدﻟﺔ )‪ CV (10‬ﺿﺮﻳﺐ اﻓﺖ ارﺗﻔﺎع ﻧﻈﻴﺮ اﻧﺮژي ﺷﻴﺮ اﺳﺖ‪ .‬در ﺿﻤﻦ ﺑﻪ ﻋﻠﺖ وﺟﻮد ﺟﺮﻳﺎن‬

‫‪QPI ,1 . QPI ,1‬‬

‫اﺳﺘﻔﺎده ﺷﺪه اﺳﺖ‪.‬‬ ‫ﻣﻌﻜﻮس از‬

‫ﺷﻜﻞ‪ -2-‬ﺷﺮاﻳﻂ ﻣﺮزي ﭘﻤﭗ‬

‫‪ -4‬ﻣﻌﺎدﻟﻪ اﺟﺮام دوراﻧﻲ‬

‫ﮔﺸﺘﺎور ﭼﺮﺧﺸﻲ ﻳﻚ ﺳﻴﺴﺘﻢ در ﺣﺎل دوران ﺑﺮاﺑﺮ ﺣﺎﺻﻠﻀﺮب ﻣﻤﺎن اﻳﻨﺮﺳﻲ ﻗﻄﺒﻲ ﺳﻴﺴﺘﻢ و ﺷﺘﺎب‬

‫زاوﻳﻪ اي ﻣﻲ ﺑﺎﺷﺪ‪ .‬ﺑﻌﺪ از ﻗﻄﻊ ﻧﻴﺮوي ﻣﺤﺮﻛﻪ ﭘﻤﭗ ﮔﺸﺘﺎور ﺧﺎرﺟﻲ وارد ﺑﺮ ﭘﻤﭗ وﺟﻮد ﻧﺪارد و ﮔﺸﺘﺎور‬

‫ﻛﻨﺪ ﻛﻨﻨﺪه ﺳﺮﻋﺖ ﭘﻤﭗ‪ ،‬ﮔﺸﺘﺎور ﭘﻤﭗ ﺧﻮاﻫﺪ ﺑﻮد‪[۴].‬‬

‫‪dω‬‬
‫‪T = −I‬‬
‫‪dt‬‬ ‫)‪(11‬‬

‫از ﻃﺮﻓﻲ‪ ،‬ﺗﻮان دوم ﻓﺎﺻﻠﻪ × ﺟﺮم = ﻣﻤﺎن اﻳﻨﺮﺳﻲ ﻗﻄﺒﻲ‬

‫‪٦‬‬
 ‫‪dω‬‬
‫‪T = −WR 2‬‬
‫‪dt‬‬ ‫)‪(12‬‬
‫‪2π dN‬‬
‫‪T = −WR 2 .‬‬ ‫‪.‬‬
‫‪60 dt‬‬ ‫)‪(13‬‬

‫در رواﺑﻂ ﻓﻮق ‪ I=WR2‬ﻣﻤﺎن اﻳﻨﺮﺳﻲ ﭘﻤﭗ‪ ،‬ﻣﻮﺗﻮر ﭘﻤﭗ‪ ،‬ﻗﻄﻌﺎت ﮔﺮدﻧﺪه و ﻣﺎﻳﻊ داﺧﻞ آن اﺳﺖ‪.‬‬

‫‪ ω‬ﺳﺮﻋﺖ زاوﻳﻪ اي ﭘﺮواﻧﻪ ﭘﻤﭗ ﺑﺮ ﺣﺴﺐ رادﻳﺎن ﺑﺮ ﺛﺎﻧﻴﻪ و ‪ N‬ﺗﻌﺪاد دور ﭘﻤﭗ در دﻗﻴﻘﻪ ﻣﻲ ﺑﺎﺷﺪ‪.‬‬

‫ﺑﺎ ﺗﻮﺟﻪ ﺑﻪ رواﺑﻂ )‪ (10‬و )‪ (11‬و )‪ (12‬و )‪(13‬ﻣﻌﺎدﻟﺔ زﻳﺮ ﺑﺪﺳﺖ ﻣﻲ آﻳﺪ‪.‬‬
‫‪2π . N R . dα‬‬
‫‪β = −WR 2 .‬‬
‫‪60TR . dt‬‬ ‫)‪(14‬‬

‫ﻛﻪ در ﻣﻌﺎدﻟﻪ )‪ (14‬دارﻳﻢ‪.‬‬

‫‪60γHR QR‬‬
‫= ‪TR‬‬
‫‪2πη R . N R‬‬ ‫)‪(15‬‬

‫در ﻣﻌﺎدﻟﺔ )‪ γ (15‬وزن ﻣﺨﺼﻮص ﺳﻴﺎل‪ ηá ،‬ﺑﺎزده ﭘﻤﭗ‪ HR ،‬ارﺗﻔﺎع ﭘﻤﭙﺎژ‪ N R ،‬ﺳﺮﻋﺖ دوراﻧﻲ‬

‫ﭘﻤﭗ‪ QR ،‬دﺑﻲ ﭘﻤﭙﺎژ در ﺷﺮاﻳﻂ ﺟﺮﻳﺎن داﻳﻤﻲ ﻫﺴﺘﻨﺪ‪.‬‬

‫ﺑﺎ اﺳﺘﻔﺎده از ﻣﻘﺪار ﻣﺘﻮﺳﻂ ‪ β‬در ﻃﻲ ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ‪ I‬ام ﻣﻌﺎدﻟﺔ )‪ (14‬ﭼﻨﻴﻦ ﻧﺘﻴﺠﻪ ﻣﻲ ﮔﺮدد‪.‬‬
‫‪αP −α‬‬ ‫‪60TR‬‬ ‫‪β +β‬‬
‫=‬ ‫‪. P‬‬
‫‪∆t‬‬ ‫‪2πWR . N R‬‬
‫‪2‬‬
‫‪2‬‬ ‫)‪(16‬‬
‫‪−30TR . ∆t‬‬
‫= ‪C6‬‬
‫ﻣﻌﺎدﻟﺔ )‪ (16‬ﺗﺒﺪﻳﻞ ﺑﻪ راﺑﻄﻪ زﻳﺮ ﻣﻲ ﮔﺮدد‪.‬‬ ‫‪2πWR2 . NR‬‬ ‫اﮔﺮ ‪ C6‬را اﻳﻦ ﭼﻨﻴﻦ ﺗﻌﺮﻳﻒ ﻛﻨﻴﻢ‬

‫⎧‬
‫‪⎪α − C β = α + C β‬‬ ‫)‪(١٧‬‬
‫‪⎪⎪ P‬‬ ‫‪6 P‬‬ ‫‪6‬‬

‫⎨‬
‫⎪‬ ‫‪−30TR∆t‬‬ ‫)‪(١٨‬‬
‫= ‪⎪C6‬‬
‫⎩⎪‬ ‫‪2πWR2 NR‬‬

‫‪ -5‬ﻣﻌﺎدﻻت ﻣﺸﺨﺼﻪ ﻣﺮﺑﻮط ﺑﻪ ﻟﻮﻟﻪ راﻧﺶ‬

‫‪٧‬‬
 ‫ﺑﺎ ﻧﻮﺷﺘﻦ ﻣﻌﺎدﻟﺔ ﻣﺸﺨﺼﻪ ﻣﻨﻔﻲ ﺑﺮاﺳﺎس راﺑﻄﺔ )‪ (98‬در ﻧﻘﻄﻪ ﻣﺠﺎور ﭘﻤﭗ دارﻳﻢ‪.‬‬
‫‪QPI ,1 = CN + Ca H PI ,1‬‬
‫)‪(19‬‬

‫ﻣﻌﺎدﻟﻪ ﭘﻴﻮﺳﺘﮕﻲ ﻫﻢ ﻧﻮﺷﺘﻪ ﻣﻲ ﺷﻮد‪.‬‬

‫‪QPI ,1 = QP‬‬
‫)‪(20‬‬

‫در راﺑﻄﻪ اﺧﻴﺮ ‪ QP‬دﺑﻲ ﭘﻤﭙﺎژ در اﻧﺘﻬﺎي ﻓﺎﺻﻠﺔ زﻣﺎﻧﻲ ‪ I‬ﻣﻲ ﺑﺎﺷﺪ‪.‬‬

‫ﺣﻞ ﻣﻌﺎدﻻت‬

‫ﻣﻌﺎدﻟﻪ)‪ ( 19‬را ﻣﺠﺪداً ﻣﻲ ﻧﻮﻳﺴﻴﻢ‪.‬‬

‫‪H PI ,1 = Hsuction + H P − ∆H PV‬‬

‫‪H PI ,1‬‬
‫ﻣﺤﺎﺳﺒﻪ و در ﻣﻌﺎدﻟﺔ ) ‪ ( 19‬ﺟﺎ ﮔﺬاري ﻣﻲ ﺷﻮد‪.‬‬ ‫از ﻣﻌﺎدﻟﻪ ﻓﻮق‪،‬‬
‫‪QPI ,1 − CN‬‬
‫‪= Hsuction + H P − CV QPI ,1 . QPI ,1‬‬
‫‪Ca‬‬

‫اﮔﺮ ﻣﻌﺎدﻟﺔ اﺧﻴﺮ ﺧﻼﺻﻪ ﺷﻮد‪.‬‬

‫‪QPI ,1 = CN + Ca Hsuction + Ca H P − Ca CV QPI ,1 . QPI ,1‬‬

‫از ﻣﻌﺎدﻟﻪ ) ‪( 11‬ﻣﻲ ﺗﻮان ﺑﻪ ﺟﺎي ‪ P‬ﻣﻘﺪار ‪ vP . QP‬و ﺑﻪ ﺟﺎي ‪ H P‬ﻣﻘﺪار ‪ HR . hP‬را ﺟﺎﻳﮕﺰﻳﻦ‬
‫‪Q‬‬ ‫‪I ,1‬‬

‫ﻛﺮد‪.‬‬
‫‪vP . QR = CN + Ca Hsuction + Ca . HR . hP − Ca CV QR2 . vP . vP‬‬ ‫)‪( 21‬‬

‫‪٨‬‬
 ‫اﻛﻨﻮن ﻣﻌﺎدﻻت ) ‪(87‬و ) ‪( 88‬و ) ‪( 97‬و ) ‪( 98‬را ﺑﺎزﻧﻮﻳﺴﻲ ﻣﻲ ﻛﻨﻴﻢ‪.‬‬

‫⎧‬ ‫‪hP‬‬ ‫⎞ ‪−1 ⎛ α P‬‬

‫‪⎪ 2‬‬ ‫=‬ ‫‪a‬‬ ‫‪+‬‬ ‫‪a‬‬ ‫‪t‬‬ ‫‪a‬‬ ‫‪n‬‬ ‫⎜‬ ‫⎟‬
‫‪⎪ α P + vP‬‬
‫‪2‬‬ ‫‪1‬‬ ‫‪2‬‬
‫⎠ ‪⎝ vP‬‬
‫⎪‬
‫⎪‬
‫⎪‬ ‫‪βP‬‬ ‫⎞ ‪−1 ⎛ α P‬‬
‫⎟ ‪⎪ α 2 + v 2 = a 3 + a 4 ta n ⎜ v‬‬
‫‪⎪ P‬‬ ‫‪P‬‬ ‫⎠ ‪⎝ P‬‬
‫⎪‬
‫⎨‬
‫‪⎪α − C . β = α + C . β‬‬
‫‪⎪ P‬‬ ‫‪6‬‬ ‫‪P‬‬ ‫‪6‬‬

‫⎪‬
‫⎪‬
‫‪⎪ v P . Q R = C N + C a H s u c tio n + C a H R h P − C a . C V . Q R . v P . v P‬‬
‫‪2‬‬

‫⎪‬
‫⎪‬
‫⎩⎪‬

‫از ﭼﻬﺎر ﻣﻌﺎدﻟﻪ ﭼﻬﺎر ﻣﺠﻬﻮل اﺧﻴﺮ ‪ β P , hP‬ﺣﺬف ﻣﻲ ﺷﻮﻧﺪ و دارﻳﻢ‪:‬‬

‫⎞ ‪⎛α P‬‬
‫⎜ ‪F1 ≡ − vP . QR + CN + Ca Hsuction + Ca . HR . a1 (α 2P + vP2 ) + Ca . HR . a 2 .(α 2P + vP2 ).tan −1‬‬ ‫⎟‬
‫⎠ ‪⎝ vP‬‬
‫‪− Ca . CV .. QR2 . vP . vP ≡ 0‬‬ ‫)‪(٢٢‬‬

‫)‪(23‬‬

‫⎞ ‪⎛α P‬‬
‫⎜ ‪F2 ≡ α P − a 3 . C6 (α 2P + vP2 ) − a 4 . C6 (α 2P + vP2 ).tan −1‬‬ ‫‪⎟ − α − C6 β ≡ 0‬‬
‫⎠ ‪⎝ vP‬‬

‫در دو ﻣﻌﺎدﻟﻪ اﺧﻴﺮ ) ‪ ( F1 , F2‬ﻓﻘﻂ دو ﻣﺠﻬﻮل ‪ α P , vP‬را دارﻳﻢ‪ .‬ﺟﻬﺖ ﺣﻞ اﻳﻦ دو ﻣﻌﺎدﻟﻪ دو ﻣﺠﻬﻮل‬

‫)ﻣﻌﺎدﻻت ﻏﻴﺮ ﺧﻄﻲ( روش ﻧﻴﻮﺗﻮن را ﻓﺴﻮن را ﺑﻜﺎر ﻣﻲ ﺑﺮﻳﻢ‪.‬‬

‫در اﻳﻦ روش اﺑﺘﺪا ﺟﻮاب ﻣﻌﺎدﻻت ﺣﺪس زده ﻣﻲ ﺷﻮد و ﺑﻌﺪ ﺑﻪ روش ﺗﻜﺮار ﺗﺎ ﺣﺪ دﻗﺖ ﻣﻮرد ﻧﻈﺮ‬

‫ﺑﻪ ﺟﻮاب واﻗﻌﻲ ﻧﺰدﻳﻚ ﻣﻲ ﺷﻮﻳﻢ‪ .‬اﮔﺮ ﭼﻨﺎﻧﭽﻪ ﺣﺪﺳﻬﺎي اوﻟﻴﻪ را ﺑﺎ ‪ ve ,α e‬ﻧﺸﺎن دﻫﻴﻢ‪.‬‬

‫‪⎧α e = α I + ∆α I −1‬‬
‫⎨‬
‫‪⎩ve = vI + ∆vI −1‬‬
‫اﮔﺮ اﻧﺪﻳﺲ )‪ (1‬ﻣﻌﺮف اوﻟﻴﻦ ﺣﺪس و اﻧﺪﻳﺲ )‪ (2‬ﺟﻮاب ﺑﺪﺳﺖ آﻣﺪه از ﺣﺪس اول ﺑﺎﺷﺪ‪ .‬ﻣﻲ ﺗﻮان‬

‫ﻧﻮﺷﺖ‪:‬‬

‫‪٩‬‬
 ‫‪⎧⎪α (P2 ) = α (P1) + δα P‬‬
‫)‪⎨ (2‬‬
‫‪⎪⎩vP = v(P1) + δvP‬‬

‫ﻃﺒﻖ روش ﻧﻴﻮﺗﻮن راﻓﺴﻮن ﻣﻘﺎدﻳﺮ ‪ δvP ,δα P‬ﻋﺒﺎرﺗﻨﺪاز‪:‬‬

‫‪∂F1‬‬ ‫‪∂F‬‬
‫‪F2 .‬‬ ‫‪− F1. 2‬‬
‫‪∂vP‬‬ ‫‪∂vP‬‬
‫= ‪δα P‬‬
‫‪∂F1 ∂F2 ∂F1 ∂F2‬‬
‫‪.‬‬ ‫‪−‬‬ ‫‪.‬‬
‫‪∂α P ∂vP ∂vP ∂α P‬‬ ‫) ‪( 24‬‬
‫‪∂F1‬‬ ‫‪∂F‬‬
‫‪F2 .‬‬ ‫‪− F1. 2‬‬
‫‪∂α P‬‬ ‫‪∂α P‬‬
‫= ‪δvP‬‬
‫‪∂F1 ∂F2‬‬ ‫‪∂ F ∂F‬‬
‫‪.‬‬ ‫‪− 1 . 2‬‬
‫‪∂vP ∂α P ∂α P ∂vP‬‬ ‫)‪( 25‬‬

‫در ﻣﻌﺎدﻻت )‪ [(25) (24‬ﻣﻘﺎدﻳﺮ ﺗﻮاﺑﻊ ‪ F2 , F1‬و ﻣﺸﺘﻘﺎت آﻧﻬﺎ ﻧﺴﺒﺖ ﺑﻪ ‪ α P , vP‬ﺑﻪ ازاري )‪α (P1) , v(P1‬‬

‫ﻣﺤﺎﺳﺒﻪ ﺧﻮاﻫﻨﺪ ﺷﺪ‪ .‬ﻣﺸﺘﻘﮕﻴﺮي از ﻣﻌﺎدﻻت )‪ (23‬و )‪ (22‬رواﺑﻂ زﻳﺮ را ﺑﺮاي اﻳﻦ ﻣﺸﺘﻘﻬﺎ ﻧﺘﻴﺠﻪ ﻣﻲ دﻫﺪ‪.‬‬

‫‪∂F1‬‬ ‫⎛‬ ‫⎞⎞ ‪⎛ α‬‬

‫⎟ ⎟ ‪= Ca . H R ⎜ 2a1.α P + a 2 . vP + 2a 2 .α P .tan −1 ⎜ P‬‬
‫‪∂α P‬‬ ‫⎝‬ ‫⎠ ⎠ ‪⎝ vP‬‬
‫)‪( 26‬‬

‫‪∂F1‬‬ ‫⎛‬ ‫⎞⎞ ‪⎛ α‬‬

‫‪= Ca . HR⎜ 2a1. vP − a2 .α P + 2a2 . vP .tan−1⎜ P ⎟⎟ − QR − 2Ca .CV.QR2 . vP‬‬
‫‪∂vP‬‬ ‫⎝‬ ‫⎠⎠ ‪⎝ vP‬‬
‫) ‪( 27‬‬

‫‪∂F2‬‬ ‫⎛‬ ‫⎞⎞ ‪⎛ α‬‬

‫⎟ ⎟ ‪= 1 − C6 ⎜ 2a 3 .α P + a 4 . vP + 2a 4 .α P .tan −1 ⎜ P‬‬
‫‪∂α P‬‬ ‫⎝‬ ‫⎠ ⎠ ‪⎝ vP‬‬
‫)‪( 28‬‬

‫‪∂F2‬‬ ‫⎛‬ ‫⎞⎞ ‪⎛ α‬‬

‫⎟ ⎟ ‪= C6 ⎜ −2a 3 . vP + a 4 .α P − 2a 4 . vP .tan −1 ⎜ P‬‬
‫‪∂vP‬‬ ‫⎝‬ ‫⎠ ⎠ ‪⎝ vP‬‬
‫) ‪( 29‬‬
‫)‪(2) (2‬‬
‫اﮔﺮ ‪ δvP , δα P‬از ﺣﺪ ﻣﺠﺎز ﺧﻄﺎي ﻣﻮرد ﻧﻈﺮ )ﻣﺜﻼُ ‪ ( 0,001‬ﻛﻤﺘﺮ ﺷﻮﻧﺪ ‪ α P , vP‬ﺟﻮاﺑﻬﺎي‬
‫)‪(1) (1‬‬
‫ﻣﻌﺎدﻻت )‪ (22‬و )‪ (23‬ﺧﻮاﻫﻨﺪ ﺑﻮد‪ .‬در ﻏﻴﺮ اﻳﻨﺼﻮرت ‪ α P , vP‬ﻣﺴﺎوي ﺑﺎ‬
‫) ‪(2‬‬
‫‪ α‬ﻗﺮار داده ﺧﻮاﻫﻨﺪ‬ ‫‪P‬‬ ‫) ‪, v (P2‬‬

‫ﺷﺪ و ﻣﺤﺎﺳﺒﺎت ﺗﻜﺮار ﻣﻲ ﮔﺮدد ﺗﺎ ﺟﻮاﺑﻬﺎي ﻣﻄﻠﻮب ﺣﺎﺻﻞ ﮔﺮدﻧﺪ‪ .‬ﺑﻌﺪ از ﺗﻌﻴﻴﻦ ‪ α P , vP‬ﺑﺎﻳﺪ ﺑﺮرﺳﻲ‬

‫ﺷﻮد ﻛﻪ آن ﻗﻄﻌﻪ از ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﭘﻤﭗ ﻛﻪ در ﻣﺤﺎﺳﺒﺎت ﺑﻜﺎر رﻓﺘﻪ ﻣﺮﺑﻮط ﺑﻪ ‪ vP‬و ‪ α P‬ﻣﻲ ﺷﻮد ﻳﺎ ﻧﻪ‪.‬‬

‫‪١٠‬‬
‫ﭼﻨﺎﻧﭽﻪ اﻳﻨﻄﻮر ﻧﺒﺎﺷﺪ ‪ ve‬و ‪ α e‬ﻣﺴﺎوي ﺑﺎ ‪ vP‬و ‪ α P‬ﻓﺮض ﺧﻮاﻫﻨﺪ ﺷﺪ و ﻫﻤﺎن روﻧﺪ ﺣﻞ ﺗﻜﺮار ﻣﻲ ﮔﺮدد‪.‬‬

‫وﻟﻲ ﭼﻨﺎﻧﭽﻪ آن ﻗﺴﻤﺖ از ﻣﻨﺤﻨﻲ ﻣﺸﺨﺼﻪ ﻛﻪ در ﻣﺤﺎﺳﺒﺎت ﺑﻜﺎر رﻓﺘﻪ درﺳﺖ ﺑﺎﺷﺪ ‪ hP‬و ‪ β P‬از ﻣﻌﺎدﻻت‬
‫‪H PI‬‬ ‫‪QPI‬‬
‫از ﻣﻌﺎدﻻت ) ‪( 9‬و ) ‪( 20‬‬ ‫‪,1‬‬
‫و‬ ‫‪,1‬‬
‫) ‪ ( 7‬و)‪ ( 8‬ﻣﺤﺎﺳﺒﻪ ﺷﺪه و ‪ H P‬و ‪ QP‬از ﻣﻌﺎدﻻت )‪ ( 1‬و‬

‫ﺣﺴﺎب ﻣﻲ ﮔﺮدﻧﺪ‪.‬‬

‫ﻣﻘﺎدﻳﺮ ‪ v‬و ‪ α‬ﺑﺮاي ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ﺑﻌﺪي ﺑﻌﻨﻮان ﻣﻘﺎدﻳﺮ اوﻟﻴﻪ در ﻧﻈﺮ ﮔﺮﻓﺘﻪ ﺷﺪه و ﺣﻞ ﻣﻌﺎدﻟﻪ در‬

‫ﻓﺎﺻﻠﻪ زﻣﺎﻧﻲ ﺑﻌﺪي اﻧﺠﺎم ﺧﻮاﻫﺪ ﺷﺪ‪ .‬ﺑﺮاي ﭘﻴﺸﮕﻴﺮي از ﻧﺎﻣﺤﺪود ﺷﺪن ﺗﻌﺪاد ﺗﻜﺮارﻫﺎ ﻣﻲ ﺗﻮان از ﻳﻚ ﻋﺪد‬

‫ﺷﻤﺎرﻧﺪه اﺳﺘﻔﺎده ﻛﺮد ﺗﺎ اﻳﻨﻜﻪ ﻣﺜﻼ" وﻗﺘﻲ ﺗﻌﺪاد ﺗﻜﺮارﻫﺎ از ﺣﺪ ﻣﻌﻴﻨﻲ ﻣﺜﻼ" ‪ 40‬ﺗﺠﺎوز ﻛﺮد ﻣﺤﺎﺳﺒﺎت‬

‫ﻣﺘﻮﻗﻒ ﺷﻮد‪.‬‬

‫‪ -6‬ﻧﺘﻴﺠﻪ ﮔﻴﺮي‪:‬‬

‫ﺑﺮﻧﺎﻣﻪ ﻛﺎﻣﭙﻴﻮﺗﺮي در زﺑﺎن ﺑﺮﻧﺎﻣﻪ ﻧﻮﻳﺴﻲ ‪ Turbo c‬ﻧﻮﺷﺘﻪ ﺷﺪه اﺳﺖ ‪.‬اﻳﻦ ﺑﺮﻧﺎﻣﻪ اﻣﻮاج ﻓﺸﺎري ﺣﺎﺻﻠﻪ ﺑﻌـﺪ از ﻗﻄـﻊ‬
‫ﻧﺎﮔﻬﺎﻧﻲ ﺑﺮق ﭘﻤﭗ ) ﻓﺸﺎر ﻣﺎﻛﺰﻳﻤﻢ و ﻣﻴﻨﻴﻤﻢ ( را در ﻧﺰدﻳﻜﻲ ﭘﻤﭗ ودر ﻧﻘﺎط ﻣﺨﺘﻠﻒ در ﻃﻮل ﺧﻂ اﻧﺘﻘﺎل را ﺣﺴﺎب‬
‫ﻣﻲ ﻛﻨﺪ‪ ).‬ﻛﻪ در اﻳﻦ ﺑﺮﻧﺎﻣﻪ ‪ 5‬ﻧﻘﻄﻪ را در ﻧﻈﺮ ﮔﺮﻓﺘﻪ اﻳﻢ( اﮔﺮ ﭼﻨﺎﻧﭽﻪ ﻓﺸﺎر ﻣﺎﻛﺰﻳﻤﻢ از ﻓﺸﺎر ﻗﺎﺑـﻞ ﺗﺤﻤـﻞ ﻟﻮﻟـﻪ و ﻳـﺎ‬
‫ﻓﺸﺎر ﻣﻴﻨﻴﻤﻢ ﺑﻪ ﺻﻮرت ﻓﺸﺎر ﻣﻨﻔﻲ ﻇﺎﻫﺮ ﺷﻮد در اﻳﻦ ﺣﺎﻟﺖ ﻧﻴﺎز ﺑﻪ راﻫﻬﺎي ﭘﻴﺸﮕﻴﺮي ﻛﻨﻨﺪه دارﻳﻢ ﺗـﺎ از اﻳـﻦ ﭘﺪﻳـﺪه‬
‫ﻣﺨﺮب ﺟﻠﻮﮔﻴﺮي ﮔﺮدد‪ .‬در اﻳﻦ ﻣﺮﺣﻠﻪ‪ ،‬ﺑﺮﻧﺎﻣﻪ ﻣﻘﺎدﻳﺮدﺑﻲ و ﻓﺸﺎرﻫﺎي ﺣﺎﺻﻠﻪ از ﺿﺮﺑﻪ ﻗﻮچ در اﻳـﺴﺘﮕﺎه را ﻣﺤﺎﺳـﺒﻪ‬
‫ﻣﻲ ﻛﻨﺪ‪.‬ﺑﻌﺪ از اﺟﺮاي ﺑﺮﻧﺎﻣﻪ در ﻣﺤﻴﻂ ﺑﺮﻧﺎﻣﻪ ﻧﻮﻳﺴﻲ ازﻣﺎ ﻣﻘﺎدﻳﺮي را ﺑﻪ ﺻـﻮرت زﻳـﺮ ﺧﻮاﻫـﺪ ﺧﻮاﺳـﺖ ﻛـﻪ ﺑﻌـﺪ از‬
‫وارد ﻛﺮدن اﻳﻦ ﻣﻘﺎدﻳﺮ‪ ،‬ﻧﺘﺎﻳﺞ ﺑﻪ ﺻﻮرت ﻣﻨﺤﻨﻴﻬﺎي زﻳﺮ اراﺋﻪ ﻣﻴﺸﻮد‪.‬‬

‫‪1- Number of parallel pump? 2‬‬

‫? ‪2- moment intertia‬‬ ‫‪16.85‬‬
‫‪3- number of rotation per minute ? 1100‬‬
‫‪4- diameter of pipe ? 0.75‬‬
‫‪5- one pump debit ? 0.25‬‬
‫‪6- efficiency value ? 0.84‬‬
‫‪7- energy dissipation coefficient at valve ? 0.01‬‬

‫‪١١‬‬
 ‫ﻧﻤﻮدار دﺑﻲ ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن ﺑﻼﻓﺎﺻﻠﻪ ﺑﻌﺪ از ﭘﻤﭗ‬

‫‪0.6‬‬
‫دﺑﻲ ) ﻣﺘﺮ ﻣﻜﻌﺐ ﺑﺮ ﺛﺎﻧﻴﻪ(‬ ‫‪0.4‬‬
‫‪0.2‬‬
‫‪0‬‬
‫‪Series1‬‬
‫‪-0.2 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬

‫‪-0.4‬‬
‫‪-0.6‬‬
‫‪-0.8‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار دﺑﻲ ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪1-‬‬

‫‪0.6‬‬
‫‪0.4‬‬
‫دﺑﻲ ) ﻣﺘﺮ ﻣﻜﻌﺐ ﺑﺮ ﺛﺎﻧﻴﻪ(‬

‫‪0.2‬‬
‫‪0‬‬
‫‪Series1‬‬
‫‪-0.2 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬

‫‪-0.4‬‬
‫‪-0.6‬‬
‫‪-0.8‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار ﻓﺸﺎر ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪1-‬‬

‫‪250‬‬

‫‪200‬‬
‫‪150‬‬
‫ﻓﺸﺎر ) ﻣﺘﺮ (‬

‫‪100‬‬ ‫‪Series1‬‬
‫‪50‬‬

‫‪0‬‬
‫‪-50 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬

‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫‪١٢‬‬
 ‫ﻧﻤﻮدار دﺑﻲ ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪2-‬‬

‫‪0.6‬‬
‫دﺑﻲ ) ﻣﺘﺮ ﻣﻜﻌﺐ ﺑﺮ ﺛﺎﻧﻴﻪ(‬ ‫‪0.4‬‬
‫‪0.2‬‬
‫‪0‬‬
‫‪Series1‬‬
‫‪-0.2 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬
‫‪-0.4‬‬
‫‪-0.6‬‬
‫‪-0.8‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار ﻓﺸﺎر ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪2-‬‬

‫‪250‬‬

‫‪200‬‬
‫‪150‬‬
‫ﻓﺸﺎر ) ﻣﺘﺮ (‬

‫‪100‬‬ ‫‪Series1‬‬
‫‪50‬‬

‫‪0‬‬
‫‪-50 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬

‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار دﺑﻲ ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪3-‬‬

‫‪0.6‬‬
‫‪0.4‬‬
‫دﺑﻲ ) ﻣﺘﺮﻣﻜﻌﺐ ﺑﺮﺛﺎﻧﻴﻪ (‬

‫‪0.2‬‬
‫‪0‬‬
‫‪Series1‬‬
‫‪-0.2 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬
‫‪-0.4‬‬
‫‪-0.6‬‬
‫‪-0.8‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫‪١٣‬‬
 ‫ﻧﻤﻮدار ﻓﺸﺎر ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪3-‬‬

‫‪250‬‬

‫‪200‬‬
‫‪150‬‬
‫ﻓﺸﺎر ) ﻣﺘﺮ(‬

‫‪100‬‬ ‫‪Series1‬‬
‫‪50‬‬

‫‪0‬‬
‫‪-50 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬

‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار دﺑﻲ ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪4-‬‬

‫‪0.6‬‬
‫‪0.4‬‬
‫دﺑﻲ ) ﻣﺘﺮﻣﻜﻌﺐ ﺑﺮ ﺛﺎﻧﻴﻪ(‬

‫‪0.2‬‬
‫‪0‬‬
‫‪-0.2 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬ ‫‪Series1‬‬
‫‪-0.4‬‬
‫‪-0.6‬‬
‫‪-0.8‬‬
‫‪-1‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار ﻓﺸﺎر ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪4-‬‬

‫‪250‬‬

‫‪200‬‬
‫‪150‬‬
‫ﻓﺸﺎر) ﻣﺘﺮ(‬

‫‪100‬‬ ‫‪Series1‬‬
‫‪50‬‬

‫‪0‬‬
‫‪-50 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬

‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫‪١٤‬‬
 ‫ﻧﻤﻮدار دﺑﻲ ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪5-‬‬

‫‪0.6‬‬
‫دﺑﻲ) ﻣﺘﺮﻣﻜﻌﺐ ﺑﺮ ﺛﺎﻧﻴﻪ(‬ ‫‪0.4‬‬
‫‪0.2‬‬
‫‪0‬‬
‫‪-0.2 0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬ ‫‪Series1‬‬
‫‪-0.4‬‬
‫‪-0.6‬‬
‫‪-0.8‬‬
‫‪-1‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻧﻤﻮدار ﻓﺸﺎر ﻧﺴﺒﺖ ﺑﻪ زﻣﺎن در ﻣﻘﻄﻊ‪5-‬‬

‫‪200‬‬

‫‪150‬‬
‫‪100‬‬
‫ﻓﺸﺎر) ﻣﺘﺮ(‬

‫‪50‬‬ ‫‪Series1‬‬
‫‪0‬‬
‫‪0‬‬ ‫‪50‬‬ ‫‪100‬‬ ‫‪150‬‬ ‫‪200‬‬
‫‪-50‬‬
‫‪-100‬‬
‫زﻣﺎن )ﺛﺎﻧﻴﻪ(‬

‫ﻣﻨﺎﺑﻊ و ﻣﺎﺧﺬ‪:‬‬
‫‪1- Axworthy, DH. and karney, B.W.journal of hydraulic engineering.‬‬
‫‪Vol.124.no.3 march1998‬‬
‫‪2- Streeter,v.l. and Wylie,e.b hydraulics transient",mcgraw-hill book co.‬‬
‫‪newyork1967‬‬
‫‪ -٣‬ﻣﺤﻤﺪ ﻧﺠﻤﺎﻳﻲ داﻧﺸﻴﺎر داﻧﺸﮕﺎه ﻋﻠﻢ و ﺻﻨﻌﺖ اﻳﺮان ‪-‬ﺿﺮﺑﻪ ﻗﻮچ‪1374 -‬‬
‫‪ -۴‬آﺷﻔﺘﻪ ﺟﻼل‪ -‬ﭘﺰﺷﻜﻲ ﻧﮋاد ﻋﻠﻴﺮﺿﺎ‪.‬ﻫﻴﺪروﻟﻴﻚ ﻛﺎرﺑﺮدي ﺟﺮﻳﺎن ﻫﺎي ﻣﻴﺮا در ﺧﻄﻮط اﻧﺘﻘﺎل‪ ،‬ﺷـﺒﻜﻪ ﻫـﺎي ﺗﻮزﻳـﻊ‪،‬‬
‫ﻛﺎﻧﺎﻟﻬﺎ و ﻧﻴﺮوﮔﺎﻫﻬﺎ ) ﺿﺮﺑﻪ ﻗﻮچ( ﺟﻠﺪ اول ‪1369‬‬

‫‪١٥‬‬