بخشی از مقاله
ﭼﻜﻴﺪه:
در اﻳﻦ ﻣﻘﺎﻟﻪ ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻳﺎ ﺗﻮزﻳﻊ ﻏﻴﺮ ﻳﻜﻨﻮاﺧﺖ ﺗﻨﺸﻬﺎي ﻧﺮﻣﺎل ﺧﻤﺸﻲ در ﻋﺮض داﻟﻬﺎي ﻓﻮﻗﺎﻧﻲ وﺗﺤﺘـﺎﻧﻲ ﻋﺮﺷـﻪ ﻫـﺎي ﺻـﻨﺪوﻗﻪ اي ﻣﻮرد ﻣﻄﺎﻟﻌﻪ ﻗﺮار ﮔﺮﻓﺘﻪ اﺳﺖ. اﺛﺮ ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻣﺜﺒﺖ و ﻣﻨﻔﻲ ﺗﺸﺮﻳﺢ ﺷﺪه و ﻧﺮم اﻓﺰار اﺟﺰا ﻣﺤﺪود Ansys 6.1 ﺟﻬﺖ ﺗﺤﻠﻴﻞ ﺗﺎﺧﻴﺮ ﺑﺮﺷـﻲ ﻋﺮﺷﻪ ﻫﺎي ﺻﻨﺪوﻗﻪ اي ﻣﻮرد اﺳﺘﻔﺎده ﻗﺮار ﮔﺮﻓﺘﻪ اﺳﺖ. ﻣﺪل اﺟﺰا ﻣﺤﺪود ﭘﻞ ﺻﻨﺪوﻗﻪ اي ﺑﺘﻨﻲ ﭘﻴﺶ ﺗﻨﻴﺪة ﺷﻴﺮﻳﻦ ﺷـﻬﺮ اﻫـﻮاز، ﺑـﻪ ﻃـﻮر ﻛﺎﻣﻞ ﺑﺎ ﻛﻤﻚ اﻳﻦ ﻧﺮم اﻓﺰار ﺳﺎﺧﺘﻪ ﺷﺪه اﺳﺖ. ﻋﺮﺷﻪ ﻣﺬﻛﻮر ﺑﺼﻮرت ﭘﺲ ﻛﺸﻴﺪه و ﺑﺎ روش ﻃﺮه اي آزاد اﺟـﺮا ﺷـﺪه اﺳـﺖ. اﺛـﺮ ﺗـﺎﺧﻴﺮ ﺑﺮﺷﻲ ﺗﺤﺖ ﺑﺎرﻫﺎي ﻣﺘﻤﺮﻛﺰ، وزن و ﻧﻴﺮوي ﻛﺎﺑﻞ ﻫﺎي ﭘﺲ ﻛﺸﻴﺪﮔﻲ ﻃﺮه اي، در ﺣﻴﻦ ﺳﺎﺧﺖ ﻣﻮرد ﻣﻼﺣﻈﻪ ﻗﺮار ﮔﺮﻓـﺖ. ﺗـﺎﺧﻴﺮ ﺑﺮﺷـﻲ ﻣﺜﺒﺖ وﻣﻨﻔﻲ ﻗﺎﺑﻞ ﻣﺸﺎﻫﺪه ﺑﻮد. ﻫﻤﭽﻨﻴﻦ ﻧﺸﺎن داده ﺷﺪ ﻛﻪ ﻧﻴﺮوي ﺣﺎﺻﻞ از ﻛﺎﺑﻞ ﻫﺎي ﭘﺲ ﻛﺸﻴﺪﮔﻲ ﻃﺮه اي ﺑﺎﻋﺚ ﺗﺸـﺪﻳﺪ اﺛـﺮ ﺗـﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻣﻲ ﺷﻮد و اﺛﺮ دﻳﺎﻓﺮاﮔﻢ ﺑﺮ روي ﭘﺪﻳﺪة ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻧﻴﺰ ﻗﺎﺑﻞ ﺗﻮﺟﻪ ﻣﻲ ﺑﺎﺷﺪ.
ﻛﻠﻴﺪ واژه ﻫﺎ: ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ، ﻋﺮﺷﻪ ﺻﻨﺪوﻗﻪ اي، اﺟﺰا ﻣﺤﺪود، ﭘﻴﺶ ﺗﻨﻴﺪﮔﻲ، دﻳﺎﻓﺮاﮔﻢ
-1 ﻣﻘﺪﻣﻪ:
وﻗﺘﻲ ﻛﻪ ﺗﻴﺮ ﺻﻨﺪوﻗﻪ اي ﺗﺤﺖ ﺑﺎر ﻗﺎﺋﻢ ﻗﺮار ﻣـﻲ ﮔﻴـﺮد ﺗـﻨﺶ ﻧﺮﻣـﺎل در ﻋـﺮض داﻟﻬـﺎي ﻓﻮﻗـﺎﻧﻲ و ﺗﺤﺘـﺎﻧﻲ داراي ﺗﻮزﻳﻊ ﻏﻴﺮﻳﻜﻨﻮاﺧﺘﻲ ﻣﻲ ﺑﺎﺷﺪ ﻛﻪ ﺑﻪ اﻳﻦ ﭘﺪﻳﺪه ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ ﮔﻔﺘﻪ ﻣﻲ ﺷﻮد و ﻣﻄﺎﺑﻖ ﺷﻜﻞ1 ﺑﻪ دو ﺻﻮرت رخ ﻣﻲ دﻫـﺪ:
در ﺣﺎﻟﺖ اول، ﺗﻨﺶ ﻋﻤﻮدي در ﻣﺤﻞ ﺗﻘﺎﻃﻊ ﺟﺎن و ﺑﺎل ﺑﻴﺸﺘﺮ از ﻣﻘﺪار آن در وﺳـﻂ ﺑـﺎل ﻣـﻲ ﺑﺎﺷـﺪ ﻛـﻪ ﺑـﻪ آن ﺗـﺄﺧﻴﺮ ﺑﺮﺷﻲ ﻣﺜﺒﺖ ﮔﻔﺘﻪ ﻣﻲ ﺷﻮد. در ﺣﺎﻟﺖ دوم، ﻣﻴﺰان ﺗﻨﺶ در وﺳﻂ ﺑﺎل ﺑﻴﺸﺘﺮ ازﻣﺤﻞ اﺗﺼﺎل ﺑﺎل و ﺟﺎن ﻣـﻲ ﺑﺎﺷـﺪ ﻛـﻪ ﺑـﻪ اﻳﻦ ﺣﺎﻟﺖ ﻧﻴﺰ ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ ﻣﻨﻔﻲ ﮔﻔﺘﻪ ﻣﻲ ﺷﻮد. ﭘﺪﻳﺪة ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ در ﻣﻘﺎﻃﻊ ﺻﻨﺪوﻗﻪ اي ﺑﺘﻨﻲ در 20 ﺳﺎل اﺧﻴـﺮ ﻣـﻮرد ﺑﺮرﺳﻲ ﻗﺮار ﮔﺮﻓﺘﻪ اﺳﺖ. ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ در ﻣﺮاﺣﻞ اوﻟﻴﺔ ﺗﺤﻠﻴﻞ ﺑﻪ ﺻـﻮرت ﺗـﺎﺧﻴﺮ ﺑﺮﺷـﻲ ﻣﺜﺒـﺖ ﻣـﻮرد ﺗﺤﻠﻴـﻞ ﻗـﺮار ﻣـﻲ ﮔﺮﻓﺖ.[1] ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻣﻨﻔﻲ ﺑﺮاي اوﻟﻴﻦ ﺑﺎر در ﺗﻴﺮﻫﺎي ﺻﻨﺪوﻗﻪ اي ﻃﺮه اي ﻣﺸﺎﻫﺪه ﮔﺮدﻳﺪ.[2] ﺗـﺎﺧﻴﺮ ﺑﺮﺷـﻲ ﺗﻴﺮﻫـﺎي
1
ﺻﻨﺪوﻗﻪ اي ﺑﺎ ﻣﻘﻄﻊ ﻣﺘﻐﻴﺮ ﻧﻴﺰ ﺟﺰو ﻣﻄﺎﻟﻌﺎﺗﻲ اﺳﺖ ﻛـﻪ در اﻳـﻦ زﻣﻴﻨـﻪ اﻧﺠـﺎم ﮔﺮﻓﺘـﻪ اﺳـﺖ.[3] در ﻫﻤـﻴﻦ زﻣﻴﻨـﻪ ﻣﻄﺎﻟﻌـﺔ ﭘﺎراﻣﺘﺮي ﺟﺎﻣﻊ وﻣﻔﻴﺪي ﺻﻮرت ﮔﺮﻓﺘﻪ ﻛﻪ اﺛﺮ ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻣﺜﺒﺖ وﻣﻨﻔﻲ را در ﺗﻴﺮﻫﺎي ﺑﺎ ﻣﻘﻄـﻊ ﻣﺘﻐﻴـﺮ وﻳﻜﻨﻮاﺧـﺖ ﺑـﻪ ازاي ﺷﺮاﻳﻂ ﺗﻜﻴﻪ ﮔﺎﻫﻲ و ﺑﺎرﮔﺬاري ﻣﺨﺘﻠﻒ ﻣﻮرد ارزﻳﺎﺑﻲ ﻗﺮار داده اﺳﺖ .[4] ﻣﻄﺎﻟﻌﺔ دﻳﮕـﺮي در اﻳـﻦ زﻣﻴﻨـﻪ اﻧﺠـﺎم ﮔﺮﻓﺘﻪ اﺳﺖ ﻛﻪ دﻟﻴﻞ وﻗﻮع ﺗﻨﺶ ﻫﺎي ﻛﺸﺸﻲ در ﻧﻮاﺣﻲ ﻓﺸﺎري و ﺑﺮﻋﻜﺲ را ﻣـﻮرد ﺗﺤﻠﻴـﻞ ﻗـﺮار داده و آن را ﺗﺤـﺖ ﻋﻨﻮان اﺛﺮafter effect ﻣﻌﺮﻓﻲ ﻧﻤﻮده اﺳﺖ.[5] روﺷـﻬﺎي ﺟـﺎﻣﻌﻲ ازﻗﺒﻴـﻞ روش ﺗﺌـﻮري ﺻـﻔﺤﺎت ﺗﺎﺷـﺪه، روش ﻧـﻮار ﻣﺤﺪود و روش اﺟﺰا ﻣﺤﺪود ﻧﻴﺰ ﺟﻬﺖ ﺗﺤﻠﻴﻞ اﻳﻦ ﭘﺪﻳﺪه ﻣﻮرد اﺳﺘﻔﺎده ﻗﺮار ﮔﺮﻓﺘﻪ اﺳﺖ.[6]
ﺷﻜﻞ:1 ﺗﻮزﻳﻊ ﺗﻨﺶ ﻧﺮﻣﺎل در ﺣﺎﻟﺖ ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻣﺜﺒﺖ وﻣﻨﻔﻲ[4]
-2 روش ﺗﺤﻠﻴﻞ
ﺑﺮاي ﺗﺤﻠﻴﻞ ﭘﺪﻳﺪة ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﺑﺎ اﺳﺘﻔﺎده از روش اﻧﺮژي، ﺗﺎﺑﻌﻲ ﻓﺮﺿﻲ ﺑﺮاي ﺗﻐﻴﻴﺮ ﺷـﻜﻞ ﺑﺎﻟﻬـﺎ در اﻣﺘـﺪاد ﻃـﻮل ﻋﺮﺷﻪ در ﻧﻈﺮ ﮔﺮﻓﺘﻪ ﻣﻲ ﺷﻮد ﻛﻪ ﻣﻌﻤﻮﻻ ﺑﻪ ﺻﻮرت ﺳﻬﻤﻲ ﻣﻲ ﺑﺎﺷﺪ ﺳﭙﺲ راﺑﻄﺔ اﻧﺮژي ﺑﺮاي ﺗﻐﻴﻴـﺮ ﺷـﻜﻞ ﻓـﺮض ﺷـﺪه ﻧﻮﺷﺘﻪ ﻣﻲ ﺷﻮد و ﺑﺎ ﻣﺸﺘﻖ ﮔﻴﺮي از راﺑﻄﺔ اﻧﺮژي و ﻣﻴﻨﻴﻤﻢ ﻛﺮدن آن، ﭘﺎراﻣﺘﺮﻫﺎي ﺗﺎﺑﻊ ﻓﺮﺿـﻲ، ﻣﺸـﺨﺺ ﺷـﺪه و ﺗﻮزﻳـﻊ
ﺗﻨﺶ ﺑﺮاي ﺗﻐﻴﻴﺮ ﺷﻜﻞ ﺑﺪﺳﺖ آﻣﺪه، ﺗﻌﻴﻴﻦ ﻣﻲ ﺷﻮد. راﺑﻄﺔ ﻧﻬﺎﺋﻲ ﺑﻪ ﺻﻮرت زﻳﺮ ﻣﻲ ﺑﺎﺷﺪ:
(1) (M ( x) M f ( x)) C t σ e x ±
I
ﻛﻪ در آن Mf ﻟﻨﮕﺮاﺿﺎﻓﻲ،σxe ﺗﻨﺶ ﻧﺮﻣﺎل در ﻣﺤﻞ ﺗﻘﺎﻃﻊ ﺟﺎن و ﺑﺎل و Ct ﻓﺎﺻﻠﺔ دال ﻓﻮﻗﺎﻧﻲ از ﻣﺤـﻮر ﺧﻨﺜـﻲ اﺳـﺖ.
در ﺻﻮرﺗﻴﻜﻪ Mf و M ﻫﻢ ﻋﻼﻣﺖ ﺑﺎﺷﻨﺪ ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ ﻣﺜﺒﺖ، و در ﻏﻴﺮ اﻳﻦ ﺻﻮرت ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ ﻣﻨﻔﻲ رخ ﻣﻲ دﻫﺪ.
اﺛﺮ ﺗﺎﺧﻴﺮ ﺑﺮﺷﻲ ﻣﺜﺒﺖ وﻣﻨﻔﻲ در ﺗﻴﺮﻫﺎي ﺑﺎ ﻣﻘﻄﻊ ﻣﺘﻐﻴﺮ و ﻳﻜﻨﻮاﺧﺖ ﺑﻪ ازاي ﺷﺮاﻳﻂ ﺗﻜﻴﻪ ﮔﺎﻫﻲ وﺑﺎرﮔﺬاري ﻣﺨﺘﻠﻒ ﻣﻮرد ارزﻳﺎﺑﻲ ﻗﺮارﮔﺮﻓﺘﻪ اﺳﺖ.[4] در اﻳﻦ ﻣﻄﺎﻟﻌﻪ اﻧﻮاع ﺑﺎرﮔﺬاري ﻫﺎ و ﺷﺮاﻳﻂ ﺗﻜﻴﻪ ﮔﺎﻫﻲ ﻣﺨﺘﻠﻒ ﻣﻮرد ﻣﻄﺎﻟﻌﻪ ﻗﺮارﮔﺮﻓﺘﻪ اﺳﺖ. ﭘﺎراﻣﺘﺮﻫﺎي ﻣﻮرد اﺳﺘﻔﺎده، ﻧﺴﺒﺖ ﻫﺎي h/H وb/L ﻣﻴﺒﺎﺷﺪ ﻛﻪ ﺑﻪ ﺗﺮﺗﻴﺐ ﻧﺴﺒﺖ ارﺗﻔﺎع وﺳﻂ دﻫﺎﻧﻪ ﺑﻪ ارﺗﻔﺎع ﺗﻜﻴﻪ ﮔﺎه و ﻧﺴﺒﺖ ﻋﺮض ﺗﻴﺮ ﺑﻪ دﻫﺎﻧﺔ آن ﻣﻲ ﺑﺎﺷﺪ.
2
ﺿﺮﻳﺐ ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ λ ﺑﺼﻮرت ﻧﺴﺒﺖ ﺗﻨﺶ ﻧﺮﻣﺎل، ﺑﺎ در ﻧﻈﺮﮔﺮﻓﺘﻦ ﭘﺪﻳﺪة ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ، ﺑﻪ ﺗﻨﺶ ﻧﺮﻣﺎل ﻣﺤﺎﺳﺒﻪ ﺷﺪه ﺑﺎ ﺗﺌﻮري ﺧﻤﺶ ﺳﺎده ﺗﻌﺮﻳﻒ ﻣﻴﺸﻮد. λe وλc ﺑﺘﺮﺗﻴﺐ ﺿﺮا ﺋﺐ ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ در ﻣﺤﻞ اﺗﺼﺎل ﺑﺎل و ﺟﺎن و در وﺳﻂ ﺑﺎل ﻣﻲ ﺑﺎﺷﻨﺪ. ﺗﻐﻴﻴﺮات ﻓﺎﻛﺘﻮر ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲ در ﻃﻮل ﺗﻴﺮ، ﺑﺮاي دﻫﺎﻧﺔ ﺳﺎده ﺑﺎ b/L=0/3 ،ﻛﻪ در ﻣﻌﺮض ﺑﺎر ﻣﺘﻤﺮﻛﺰ ﻗﺮارﮔﺮﻓﺘﻪ ﺑﺎﺷﺪ، در ﺷﻜﻞ2 ﻧﺸﺎن داده اﺳﺖ اﻳﻦ ﺗﻮزﻳﻊ ﺑﺮﺣﺴﺐ h/H ﻫﺎي ﻣﺨﺘﻠﻒ، ﻣﺘﻔﺎوت ﻣﻲ ﺑﺎﺷﺪ ﺑﺮاي ﺗﻴﺮ ﺻﻨﺪوﻗﻪ اي ﻛﻨﺴﻮﻟﻲ ﺑﺎ b/L=./2 ﻛﻪ در ﻣﻌﺮض ﺑﺎر ذوزﻧﻘﻪ اي ﻗﺮار ﮔﺮﻓﺘﻪ ﺑﺎﺷﺪ، ﺗﻐﻴﻴﺮ ﻃﻮﻟﻲ ﺿﺮﻳﺐ ﺗﺄﺧﻴﺮ ﺑﺮﺷﻲλc ، در ﺷﻜﻞ3
ﻧﺸﺎن داده ﺷﺪه اﺳﺖ، ﻫﻤﺎﻧﻄﻮر ﻛﻪ در اﻳﻦ ﺷﻜﻞ ﻧﺸﺎن داده ﺷﺪه اﺳﺖ اﺛﺮﺗﺄﺧﻴﺮﺑﺮﺷﻲ ﻣﻨﻔﻲ ﺑﺎ اﻓﺰاﻳﺶh/H ﻣﺸﻬﻮدﺗﺮ ﻣﻲ ﺑﺎﺷﺪ.
ﺷﻜﻞ:2 ﺗﻐﻴﻴﺮات ﻃﻮﻟﻲ ﻓﺎﻛﺘﻮر λe در ﺗﻴﺮ ﺑﺎ دﻫﺎﻧﺔ ﺳﺎده ﺗﺤﺖ ﺑﺎر ﻣﺘﻤﺮﻛﺰ [4 ]( b/L=./3)
ﺷﻜﻞ:3 ﺗﻐﻴﻴﺮات ﻃﻮﻟﻲ ﻓﺎﻛﺘﻮر λc در ﺗﻴﺮﻛﻨﺴﻮﻟﻲ ﺗﺤﺖ ﺑﺎر ذوزﻧﻘﻪ اي[4]
3
-3 ﻣﻌﺮﻓﻲ ﭘﻞ ﺻﻨﺪوﻗﻪ اي ﺷﻴﺮﻳﻦ ﺷﻬﺮ
ﭘﻞ ﺻﻨﺪوﻗﻪاي ﺑﺘﻨﻲ ﺷﻴﺮﻳﻦ ﺷﻬﺮ ﺑﺮ روي رودﺧﺎﻧﺔ ﻛﺎرون ﺳﺎﺧﺘﻪ ﺷﺪه اﺳﺖ. اﺟﺮاي اﻳﻦ ﭘﻞ ﺑﺼﻮرت ﻃﺮهاي آزاد ﺑﺎ ﺑﺘﻦرﻳﺰي درﺟﺎ اﻧﺠﺎم ﺷﺪه و ﺑﻪ ﺻﻮرت ﭘﺲ ﻛﺸﻴﺪه ﻣﻮرد ﭘﻴﺶ ﺗﻨﻴﺪﮔﻲ ﻗﺮار ﮔﺮﻓﺘﻪ اﺳﺖ. ﻋﺮﺷﻪ اﻳﻦ ﭘﻞ داراي ﺳﻪ دﻫﺎﻧﻪ اﺳﺖ ﻛﻪ ﻃﻮل دﻫﺎﻧﻪ ﻣﻴﺎﻧﻲ146 ﻣﺘﺮ و ﻃﻮل دﻫﺎﻧﻪﻫﺎي ﻛﻨﺎري 97 ﻣﺘﺮ ﻣﻲﺑﺎﺷﺪ. ﻋﻼوه ﺑﺮ ﺗﻐﻴﻴﺮ ارﺗﻔﺎع ﻣﻘﻄﻊ در ﻃﻮل ﭘﻞ، ﻣﺸﺨﺼﺎت ﻣﻘﻄﻊ ﻋﺮﺿﻲ از ﻗﺒﻴﻞ ﺿﺨﺎﻣﺖ دال ﺗﺤﺘﺎﻧﻲ و ﺟﺎن، در ﻃﻮل ﭘﻞ، ﻣﺘﻐﻴﺮ ﻣﻲﺑﺎﺷﺪ. ﺑﺮاي ﻣﺪﻟﺴﺎزي ﻋﺮﺷﺔ ﺻﻨﺪوﻗﻪ اي ﭘﻞ ﺷﻴﺮﻳﻦ ﺷﻬﺮ از ﻧﺮم اﻓﺰار اﺟﺰا ﻣﺤﺪودAnsys 6.1 اﺳﺘﻔﺎده ﺷﺪ. اﻟﻤﺎن ﭘﻮﺳﺘﻪ اي Shell 63 ﺟﻬﺖ اﻟﻤﺎن ﺑﻨﺪي داﻟﻬﺎي ﻋﺮﺷﻪ اﻧﺘﺨﺎب ﺷﺪ، اﻳﻦ اﻟﻤﺎن داراي ﭼﻬﺎرﮔﺮه ﻣﻲ ﺑﺎﺷﺪﻛﻪ ﻫﺮﻛﺪام از ﮔﺮه ﻫﺎ داراي ﺳﻪ درﺟﺔ آزادي اﻧﺘﻘﺎﻟﻲ وﺳﻪ درﺟﺔ آزادي ﺧﻤﺸﻲ ﻣﻲ ﺑﺎﺷﻨﺪ.[7] ﺑﺮاي ﻣﺪل ﻧﻤﻮدن ﻛﺎﺑﻞ ﻫﺎ از اﻟﻤﺎن Link8 اﺳﺘﻔﺎده ﺷﺪ. از ﻣﺸﺨﺼﺎت ورودي اﻳﻦ اﻟﻤﺎن، ﺳﻄﺢ ﻣﻘﻄﻊ و ﻛﺮﻧﺶ اوﻟﻴﻪ اﺳﺖ ﻛﻪ ﺑﺎ اﻋﻤﺎل ﻛﺮﻧﺶ اوﻟﻴﻪ ﻧﻴﺮوي ﭘﺲ ﻛﺸﻴﺪﮔﻲ ﻗﺎﺑﻞ ﻣﺪﻟﺴﺎزي ﻣﻲ ﺑﺎﺷﺪ. ﻣﺪل ﻣﺮﺑﻮط ﺑﻪ دﻫﺎﻧﺔ ﻣﻴﺎﻧﻲ ﭘﻞ ﺷﻴﺮﻳﻦ ﺷﻬﺮ در ﺷﻜﻞ4 ﻧﺸﺎن داده ﺷﺪه اﺳﺖ.
