संदृढ़ता आव्युह: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
 
(10 intermediate revisions by 4 users not shown)
Line 1: Line 1:
{{Short description|Matrix used in finite element analysis}}
{{Short description|Matrix used in finite element analysis}}
{{For|the stiffness tensor in solid mechanics|Hooke's law#Matrix representation (stiffness tensor)}}
{{For|ठोस यांत्रिकी में कठोरता टेंसर|हुक का नियम मैट्रिक्स प्रतिनिधित्व (कठोरता टेंसर)}}


अण्डाकार [[आंशिक अंतर समीकरण]]ों के संख्यात्मक समाधान के लिए परिमित तत्व विधि में, कठोरता मैट्रिक्स [[मैट्रिक्स (गणित)]] है जो [[रैखिक समीकरणों की प्रणाली]] का प्रतिनिधित्व करता है जिसे अंतर समीकरण के अनुमानित समाधान का पता लगाने के लिए हल किया जाना चाहिए।
अण्डाकार [[आंशिक अंतर समीकरण|आंशिक अंतर समीकरणों]] के संख्यात्मक समाधान के लिए परिमित तत्व विधि में, '''संदृढ़ता  [[मैट्रिक्स (गणित)|आव्युह]]''' [[मैट्रिक्स (गणित)|(गणित)]] है जो [[रैखिक समीकरणों की प्रणाली]] का प्रतिनिधित्व करता है जिसे अंतर समीकरण के अनुमानित समाधान का पता लगाने के लिए हल किया जाना चाहिए।


==[[पॉइसन समस्या]] के लिए कठोरता मैट्रिक्स==
=='''[[पॉइसन समस्या]] के लिए संदृढ़ता  आव्युह'''==
सरलता के लिए, हम पहले पॉइसन समस्या पर विचार करेंगे
सरलता के लिए, हम पहले पॉइसन समस्या पर विचार करेंगे


:<math> -\nabla^2 u = f</math>
:<math> -\nabla^2 u = f</math>
कुछ डोमेन पर {{math|Ω}}, सीमा शर्त के अधीन {{math|1=''u'' = 0}} की सीमा पर {{math|Ω}}. परिमित तत्व विधि द्वारा इस समीकरण को अलग करने के लिए, कोई [[आधार कार्य]] का सेट चुनता है {{math|{''φ''{{sub|1}}, …, ''φ{{sub|n}}''} }} पर परिभाषित किया गया {{math|Ω}} जो सीमा पर लुप्त भी हो जाते हैं। फिर अनुमान लगाता है
कुछ डोमेन पर {{math|Ω}}, सीमा शर्त के अधीन {{math|Ω}} की सीमा पर {{math|1=''u'' = 0}}. परिमित तत्व विधि द्वारा इस समीकरण को भिन्न करने के लिए, व्यक्ति Ω पर परिभाषित [[आधार कार्य|आधार कार्यों]] {φ1, …, φn} का एक समूह चुनता है जो सीमा पर भी लुप्त हो जाता है। फिर एक अनुमान लगाता है


:<math> u \approx u^h = u_1\varphi_1+\cdots+u_n\varphi_n.</math>
:<math> u \approx u^h = u_1\varphi_1+\cdots+u_n\varphi_n.</math>
गुणांक {{math|''u''{{sub|1}}, ''u''{{sub|2}}, …, ''u{{sub|n}}''}} निर्धारित किया जाता है ताकि सन्निकटन में त्रुटि प्रत्येक आधार फ़ंक्शन के लिए ऑर्थोगोनल हो {{mvar|φ{{sub|i}}}}:
गुणांक {{math|''u''{{sub|1}}, ''u''{{sub|2}}, …, ''u{{sub|n}}''}} निर्धारित किया जाता है जिससे कि सन्निकटन में त्रुटि प्रत्येक आधार फलन {{mvar|φ{{sub|i}}}} के लिए ऑर्थोगोनल हो :


:<math> \int_{x \in \Omega} \varphi_i\cdot f \, dx = -\int_{x \in \Omega} \varphi_i\nabla^2u^h \, dx = -\sum_j\left(\int_{x \in \Omega} \varphi_i\nabla^2\varphi_j\,dx\right)\, u_j = \sum_j\left(\int_{x \in \Omega} \nabla\varphi_i\cdot\nabla\varphi_j\, dx\right)u_j.</math>
:<math> \int_{x \in \Omega} \varphi_i\cdot f \, dx = -\int_{x \in \Omega} \varphi_i\nabla^2u^h \, dx = -\sum_j\left(\int_{x \in \Omega} \varphi_i\nabla^2\varphi_j\,dx\right)\, u_j = \sum_j\left(\int_{x \in \Omega} \nabla\varphi_i\cdot\nabla\varphi_j\, dx\right)u_j.</math>
कठोरता मैट्रिक्स है {{mvar|n}}-तत्व [[वर्ग मैट्रिक्स]] {{math|'''A'''}} द्वारा परिभाषित
'''संदृढ़ता  आव्युह''' {{mvar|n}}-तत्व [[वर्ग मैट्रिक्स|वर्ग आव्युह]] {{math|'''A'''}} द्वारा परिभाषित है


:<math> \mathbf A_{ij} = \int_{x \in \Omega} \nabla\varphi_i\cdot\nabla\varphi_j\, dx.</math>
:<math> \mathbf A_{ij} = \int_{x \in \Omega} \nabla\varphi_i\cdot\nabla\varphi_j\, dx.</math>
वेक्टर को परिभाषित करके {{math|'''F'''}}घटकों के साथ <math display="inline">\mathbf F_i = \int_\Omega\varphi_i f\,dx,</math> गुणांक {{mvar|u{{sub|i}}}}रेखीय प्रणाली द्वारा निर्धारित होते हैं {{math|1='''Au''' = '''F'''}}. कठोरता मैट्रिक्स [[सममित मैट्रिक्स]] है, अर्थात। {{math|1='''A'''{{sub|''ij''}} = '''A'''{{sub|''ji''}}}}, इसलिए इसके सभी [[eigenvalue]]s ​​​​वास्तविक हैं। इसके अलावा, यह सख्ती से [[सकारात्मक-निश्चित मैट्रिक्स]] है, ताकि सिस्टम {{math|1='''Au''' = '''F'''}} के पास हमेशा अनोखा समाधान होता है। (अन्य समस्याओं के लिए, ये अच्छी संपत्तियाँ खो जाएँगी।)
सदिश {{math|'''F'''}} घटकों के साथ परिभाषित करके <math display="inline">\mathbf F_i = \int_\Omega\varphi_i f\,dx,</math> गुणांक {{mvar|u{{sub|i}}}} रेखीय प्रणाली {{math|1='''Au''' = '''F'''}} द्वारा निर्धारित किए जाते हैं। संदृढ़ता  आव्युह [[सममित मैट्रिक्स|सममित आव्युह]] है, अर्थात। {{math|1='''A'''{{sub|''ij''}} = '''A'''{{sub|''ji''}}}}, इसलिए इसके सभी स्वदेशी मूल्य ​​​​वास्तविक हैं। इसके अतिरिक्त, यह सख्ती से [[सकारात्मक-निश्चित मैट्रिक्स|धनात्मक-निश्चित आव्युह]] है, जिससे कि प्रणाली {{math|1='''Au''' = '''F'''}} के पास सदैव एक अद्वितीय समाधान होता है। (अन्य समस्याओं के लिए, यह अच्छी संपत्तियाँ खो जाएँगी।)


ध्यान दें कि कठोरता मैट्रिक्स डोमेन के लिए उपयोग किए गए कम्प्यूटेशनल ग्रिड और किस प्रकार के परिमित तत्व का उपयोग किया जाता है, इसके आधार पर भिन्न होगा। उदाहरण के लिए, जब टुकड़ेवार द्विघात परिमित तत्वों का उपयोग किया जाता है तो कठोरता मैट्रिक्स में टुकड़ेवार रैखिक तत्वों की तुलना में स्वतंत्रता की अधिक डिग्री होगी।
ध्यान दें कि संदृढ़ता  आव्युह डोमेन के लिए उपयोग किए गए कम्प्यूटेशनल ग्रिड और किस प्रकार के परिमित तत्व का उपयोग किया जाता है, इसके आधार पर भिन्न होगा। उदाहरण के लिए, जब टुकड़ेवार द्विघात परिमित तत्वों का उपयोग किया जाता है तब संदृढ़ता  आव्युह में टुकड़ेवार रैखिक तत्वों की तुलना में स्वतंत्रता की अधिक डिग्री होगी।


=='''अन्य समस्याओं के लिए कठोरता मैट्रिक्स'''==
=='''अन्य समस्याओं के लिए संदृढ़ता  आव्युह'''==
अन्य पीडीई के लिए कठोरता मैट्रिक्स का निर्धारण अनिवार्य रूप से ही प्रक्रिया का पालन करता है, लेकिन यह सीमा स्थितियों की पसंद से जटिल हो सकता है। अधिक जटिल उदाहरण के रूप में, [[अण्डाकार वक्र]] पर विचार करें
अन्य पीडीई के लिए संदृढ़ता  आव्युह का निर्धारण अनिवार्य रूप से ही प्रक्रिया का पालन करता है, किन्तु यह सीमा स्थितियों की पसंद से समष्टि हो सकता है। अधिक समष्टि उदाहरण के रूप में, [[अण्डाकार वक्र|अण्डाकार समीकरण]] पर विचार करें


:<math> -\sum_{k,l}\frac{\partial}{\partial x_k}\left(a^{kl}\frac{\partial u}{\partial x_l}\right)= f</math>
:<math> -\sum_{k,l}\frac{\partial}{\partial x_k}\left(a^{kl}\frac{\partial u}{\partial x_l}\right)= f</math>
कहाँ <math displaystyle="inline">\mathbf A(x) = a^{kl}(x)</math> प्रत्येक बिंदु के लिए परिभाषित सकारात्मक-निश्चित मैट्रिक्स है {{mvar|x}} डोमेन में. हम रॉबिन सीमा शर्त लागू करते हैं
कहाँ <math displaystyle="inline">\mathbf A(x) = a^{kl}(x)</math> {{mvar|x}} डोमेन में प्रत्येक बिंदु के लिए परिभाषित धनात्मक-निश्चित आव्युह है हम रॉबिन सीमा शर्त क्रियान्वित करते हैं


:<math> -\sum_{k,l}\nu_k a^{kl}\frac{\partial u}{\partial x_l} = c(u-g),</math>
:<math> -\sum_{k,l}\nu_k a^{kl}\frac{\partial u}{\partial x_l} = c(u-g),</math>
कहाँ {{mvar|ν{{sub|k}}}} इकाई जावक [[सामान्य वेक्टर]] का घटक है {{mvar|ν}} में {{mvar|k}}-वीं दिशा. हल करने की प्रणाली है
कहाँ {{mvar|ν{{sub|k}}}} {{mvar|k}}-वीं दिशा में इकाई जावक [[सामान्य वेक्टर|सामान्य सदिश]] {{mvar|ν}} का घटक है हल करने की प्रणाली है


:<math> \sum_j\left(\sum_{k,l}\int_\Omega a^{kl}\frac{\partial\varphi_i}{\partial x_k}\frac{\partial\varphi_j}{\partial x_l}dx+\int_{\partial\Omega}c\varphi_i\varphi_j\, ds\right)u_j = \int_\Omega\varphi_i f\, dx+\int_{\partial\Omega}c\varphi_i g\, ds,</math>
:<math> \sum_j\left(\sum_{k,l}\int_\Omega a^{kl}\frac{\partial\varphi_i}{\partial x_k}\frac{\partial\varphi_j}{\partial x_l}dx+\int_{\partial\Omega}c\varphi_i\varphi_j\, ds\right)u_j = \int_\Omega\varphi_i f\, dx+\int_{\partial\Omega}c\varphi_i g\, ds,</math>
जैसा कि ग्रीन की पहचान के एनालॉग का उपयोग करके दिखाया जा सकता है। गुणांक {{mvar|u{{sub|i}}}} अभी भी रैखिक समीकरणों की प्रणाली को हल करके पाए जाते हैं, लेकिन प्रणाली का प्रतिनिधित्व करने वाला मैट्रिक्स सामान्य पॉइसन समस्या से स्पष्ट रूप से भिन्न है।
जैसा कि ग्रीन की पहचान के एनालॉग का उपयोग करके दिखाया जा सकता है। गुणांक {{mvar|u{{sub|i}}}} अभी भी रैखिक समीकरणों की प्रणाली को हल करके पाए जाते हैं, किन्तु प्रणाली का प्रतिनिधित्व करने वाला आव्युह सामान्य पॉइसन समस्या से स्पष्ट रूप से भिन्न है।


सामान्य तौर पर, प्रत्येक अदिश अण्डाकार ऑपरेटर के लिए {{mvar|L}} आदेश की {{math|2''k''}}, [[द्विरेखीय रूप]] संबद्ध है {{mvar|B}} सोबोलेव क्षेत्र पर {{mvar|H{{sup|k}}}}, ताकि समीकरण का [[कमजोर सूत्रीकरण]] हो सके {{math|1=''Lu'' = ''f''}} है
सामान्यतः, क्रम 2k के प्रत्येक अदिश अण्डाकार ऑपरेटर L के लिए, सोबोलेव स्पेस {{mvar|H{{sup|k}}}} पर एक [[द्विरेखीय रूप]] {{mvar|B}} जुड़ा होता है, ताकि समीकरण {{math|1=''Lu'' = ''f''}} का [[कमजोर सूत्रीकरण|अशक्त सूत्रीकरण]] हो।


:<math> B[u,v] = (f,v)</math>
:<math> B[u,v] = (f,v)</math>
सभी कार्यों के लिए {{mvar|v}} में {{mvar|H{{sup|k}}}}. फिर इस समस्या के लिए कठोरता मैट्रिक्स है
सभी कार्यों के लिए {{mvar|v}} में {{mvar|H{{sup|k}}}}. फिर इस समस्या के लिए संदृढ़ता  आव्युह है


:<math> \mathbf A_{ij} = B[\varphi_j,\varphi_i].</math>
:<math> \mathbf A_{ij} = B[\varphi_j,\varphi_i].</math>
=='''कठोरता मैट्रिक्स की व्यावहारिक असेंबली'''==
=='''संदृढ़ता  आव्युह की व्यावहारिक असेंबली'''==
कंप्यूटर पर परिमित तत्व विधि को लागू करने के लिए, किसी को पहले आधार कार्यों का सेट चुनना होगा और फिर कठोरता मैट्रिक्स को परिभाषित करने वाले इंटीग्रल्स की गणना करनी होगी। आमतौर पर, डोमेन {{math|Ω}} को जाल निर्माण के कुछ रूपों द्वारा विभेदित किया जाता है, जिसमें इसे गैर-अतिव्यापी त्रिभुज जाल या [[जाल के प्रकार]]ों में विभाजित किया जाता है, जिन्हें आम तौर पर तत्वों के रूप में जाना जाता है। फिर आधार कार्यों को प्रत्येक तत्व के भीतर कुछ क्रम के [[बहुपद]] और तत्व सीमाओं के पार निरंतर चुना जाता है। सबसे सरल विकल्प त्रिकोणीय तत्वों के लिए टुकड़ावार रैखिक फ़ंक्शन और आयताकार तत्वों के लिए टुकड़ावार द्विरेखीय हैं।
कंप्यूटर पर परिमित तत्व विधि को क्रियान्वित करने के लिए, किसी को पहले आधार कार्यों का समूह चुनना होगा और फिर संदृढ़ता  आव्युह को परिभाषित करने वाले इंटीग्रल्स की गणना करनी होगी। सामान्यतः, डोमेन {{math|Ω}} को जाल निर्माण के कुछ रूपों द्वारा विभेदित किया जाता है, जिसमें इसे गैर-अतिव्यापी त्रिभुज जाल या [[जाल के प्रकार|जाल के प्रकारों]] में विभाजित किया जाता है, जिन्हें सामान्यतः तत्वों के रूप में जाना जाता है। फिर आधार कार्यों को प्रत्येक तत्व के अंदर कुछ क्रम के [[बहुपद]] और तत्व सीमाओं के पार निरंतर चुना जाता है। सबसे सरल विकल्प त्रिकोणीय तत्वों के लिए टुकड़ावार रैखिक फलन और आयताकार तत्वों के लिए टुकड़ावार द्विरेखीय हैं।


तत्व कठोरता मैट्रिक्स {{math|'''A'''<sup>[''k'']</sup>}}तत्व के लिए {{mvar|T{{sub|k}}}} मैट्रिक्स है
'''तत्व {{mvar|T{{sub|k}}}}''' के लिए तत्व '''संदृढ़ता  आव्युह''' {{math|'''A'''<sup>[''k'']</sup>}}आव्युह है


:<math> \mathbf A^{[k]}_{ij} = \int_{T_k}\nabla\varphi_i\cdot\nabla\varphi_j\, dx.</math>
:<math> \mathbf A^{[k]}_{ij} = \int_{T_k}\nabla\varphi_i\cdot\nabla\varphi_j\, dx.</math>
अधिकांश मानों के लिए तत्व कठोरता मैट्रिक्स शून्य है {{mvar|i}} और {{mvar|j}}, जिसके लिए संबंधित आधार फ़ंक्शन शून्य हैं {{mvar|T{{sub|k}}}}. पूर्ण कठोरता मैट्रिक्स {{math|'''A'''}} तत्व कठोरता मैट्रिक्स का योग है। विशेष रूप से, उन आधार कार्यों के लिए जो केवल स्थानीय रूप से समर्थित हैं, कठोरता मैट्रिक्स [[विरल मैट्रिक्स]] है।
{{mvar|i}} और {{mvar|j}}, अधिकांश मानों के लिए तत्व संदृढ़ता  आव्युह शून्य है जिसके लिए संबंधित आधार फलन {{mvar|T{{sub|k}}}} के भीतर शून्य हैं पूर्ण संदृढ़ता  आव्युह {{math|'''A'''}} तत्व संदृढ़ता  आव्युह का योग है। विशेष रूप से, उन आधार कार्यों के लिए जो केवल स्थानीय रूप से समर्थित हैं, संदृढ़ता  आव्युह [[विरल मैट्रिक्स|विरल]] है।


आधार कार्यों के कई मानक विकल्पों के लिए, यानी त्रिकोणों पर टुकड़े-टुकड़े रैखिक आधार कार्यों के लिए, तत्व कठोरता मैट्रिक्स के लिए सरल सूत्र हैं। उदाहरण के लिए, टुकड़ों में रैखिक तत्वों के लिए, शीर्षों वाले त्रिभुज पर विचार करें {{math|(''x''{{sub|1}}, ''y''{{sub|1}})}}, {{math|(''x''{{sub|2}}, ''y''{{sub|2}})}}, {{math|(''x''{{sub|3}}, ''y''{{sub|3}})}}, और 2×3 मैट्रिक्स को परिभाषित करें
आधार कार्यों के अनेक मानक विकल्पों के लिए, अर्थात त्रिकोणों पर टुकड़े-टुकड़े रैखिक आधार कार्यों के लिए, तत्व संदृढ़ता  आव्युह के लिए सरल सूत्र हैं। उदाहरण के लिए, टुकड़ों में रैखिक तत्वों के लिए, शीर्षों {{math|(''x''{{sub|1}}, ''y''{{sub|1}})}}, {{math|(''x''{{sub|2}}, ''y''{{sub|2}})}}, {{math|(''x''{{sub|3}}, ''y''{{sub|3}})}} वाले त्रिभुज पर विचार करें और 2×3 आव्युह को परिभाषित करें


:<math> \mathbf D = \left[\begin{matrix}x_3 - x_2 & x_1 - x_3 & x_2 - x_1 \\ y_3 - y_2 & y_1 - y_3 & y_2 - y_1\end{matrix}\right].</math>
:<math> \mathbf D = \left[\begin{matrix}x_3 - x_2 & x_1 - x_3 & x_2 - x_1 \\ y_3 - y_2 & y_1 - y_3 & y_2 - y_1\end{matrix}\right].</math>
फिर तत्व कठोरता मैट्रिक्स है
फिर तत्व संदृढ़ता  आव्युह है


:<math> \mathbf A^{[k]} = \frac{\mathbf D^\mathsf{T} \mathbf D}{4\operatorname{area}(T)}.</math>
:<math> \mathbf A^{[k]} = \frac{\mathbf D^\mathsf{T} \mathbf D}{4\operatorname{area}(T)}.</math>
जब अंतर समीकरण अधिक जटिल होता है, मान लीजिए कि अमानवीय प्रसार गुणांक होता है, तो तत्व कठोरता मैट्रिक्स को परिभाषित करने वाले अभिन्न अंग का मूल्यांकन गॉसियन चतुर्भुज द्वारा किया जा सकता है।
जब अंतर समीकरण अधिक समष्टि होता है, मान लीजिए कि अमानवीय प्रसार गुणांक होता है, तब तत्व संदृढ़ता  आव्युह को परिभाषित करने वाले अभिन्न अंग का मूल्यांकन गॉसियन चतुर्भुज द्वारा किया जा सकता है।


कठोरता मैट्रिक्स की स्थिति संख्या संख्यात्मक ग्रिड की गुणवत्ता पर दृढ़ता से निर्भर करती है। विशेष रूप से, परिमित तत्व जाल में छोटे कोण वाले त्रिकोण कठोरता मैट्रिक्स के बड़े eigenvalues ​​​​को प्रेरित करते हैं, जिससे समाधान की गुणवत्ता खराब हो जाती है।
संदृढ़ता  आव्युह की स्थिति संख्या संख्यात्मक ग्रिड की गुणवत्ता पर दृढ़ता से निर्भर करती है। विशेष रूप से, परिमित तत्व जाल में छोटे कोण वाले त्रिकोण संदृढ़ता  आव्युह के बड़े आइगेनवैल्यू ​​​​को प्रेरित करते हैं, जिससे समाधान की गुणवत्ता खराब हो जाती है।


==संदर्भ==
=='''संदर्भ'''==
* {{citation |first1=A. |last1=Ern |first2=J.-L. |last2=Guermond |year=2004 |title=Theory and Practice of Finite Elements |publisher=Springer-Verlag |location=New York, NY |isbn=0387205748 }}
* {{citation |first1=. |last1=Ern |first2=जे.-एल. |last2=गुरमोंड |year=2004 |title=परिमित तत्वों का सिद्धांत और अभ्यास |publisher=स्प्रिंगर-वेरलाग |location=न्यूयॉर्क, एनवाई |isbn=0387205748 }}
* {{citation |first=M.S. |last=Gockenbach |year=2006 |title=Understanding and Implementing the Finite Element Method |publisher=SIAM |location=Philadelphia, PA |isbn=0898716144 }}
* {{citation |first=एम.एस. |last=गोकेनबैक |year=2006 |title=परिमित तत्व विधि को समझना और लागू करना |publisher=एस.आई.ए.एम |location=फिलाडेल्फिया, पीए |isbn=0898716144 }}
* {{citation |first1=C. |last1=Grossmann |first2=H.-G. |last2=Roos |first3=M. |last3=Stynes |year=2007 |title=Numerical Treatment of Partial Differential Equations |publisher=Springer-Verlag |location=Berlin, Germany |isbn=978-3-540-71584-9 }}
* {{citation |first1=सी. |last1=ग्रॉसमैन |first2=एच.-जी. |last2=रूस |first3=एम. |last3=स्टाइन्स |year=2007 |title=आंशिक विभेदक समीकरणों का संख्यात्मक उपचार |publisher=स्प्रिंगर-वेरलाग |location=बर्लिन, जर्मनी |isbn=978-3-540-71584-9 }}
* {{citation |first=C. |last=Johnson |year=2009 |title=Numemerical Solution of Partial Differential Equations by the Finite Element Method |publisher=Dover |isbn=978-0486469003 }}
* {{citation |first=सी. |last=जॉनसन |year=2009 |title=परिमित तत्व विधि द्वारा आंशिक विभेदक समीकरणों का संख्यात्मक समाधान |publisher=डोवर |isbn=978-0486469003 }}
* {{citation |first1=O.C. |last1=Zienkiewicz |author1-link=Olgierd Zienkiewicz |first2=R.L. |last2=Taylor |first3=J.Z. |last3=Zhu |year=2005 |title=The Finite Element Method: Its Basis and Fundamentals |publisher=Elsevier Butterworth-Heinemann |edition=6th |location=Oxford, UK |isbn=978-0750663205 }}
* {{citation |first1=.सी. |last1=ज़िएनकिविज़ |author1-link=ओल्गिएर्ड ज़िएनक्यूविक्ज़ |first2=आर.एल. |last2=टेलर |first3=जे.जेड. |last3=Zhu |year=2005 |title=परिमित तत्व विधि: इसका आधार और बुनियादी बातें |publisher=एल्सेवियर बटरवर्थ-हेनमैन |edition=6th |location=ऑक्सफोर्ड, यूके |isbn=978-0750663205 }}
[[Category: व्यावहारिक गणित]] [[Category: संख्यात्मक विश्लेषण]]


 
[[Category:Articles with hatnote templates targeting a nonexistent page]]
 
[[Category: Machine Translated Page]]
[[Category:Created On 24/07/2023]]
[[Category:Created On 24/07/2023]]
[[Category:Lua-based templates]]
[[Category:Machine Translated Page]]
[[Category:Pages with script errors]]
[[Category:Short description with empty Wikidata description]]
[[Category:Templates Vigyan Ready]]
[[Category:Templates that add a tracking category]]
[[Category:Templates that generate short descriptions]]
[[Category:Templates using TemplateData]]
[[Category:व्यावहारिक गणित]]
[[Category:संख्यात्मक विश्लेषण]]

Latest revision as of 13:58, 14 August 2023

अण्डाकार आंशिक अंतर समीकरणों के संख्यात्मक समाधान के लिए परिमित तत्व विधि में, संदृढ़ता आव्युह (गणित) है जो रैखिक समीकरणों की प्रणाली का प्रतिनिधित्व करता है जिसे अंतर समीकरण के अनुमानित समाधान का पता लगाने के लिए हल किया जाना चाहिए।

पॉइसन समस्या के लिए संदृढ़ता आव्युह

सरलता के लिए, हम पहले पॉइसन समस्या पर विचार करेंगे

कुछ डोमेन पर Ω, सीमा शर्त के अधीन Ω की सीमा पर u = 0. परिमित तत्व विधि द्वारा इस समीकरण को भिन्न करने के लिए, व्यक्ति Ω पर परिभाषित आधार कार्यों {φ1, …, φn} का एक समूह चुनता है जो सीमा पर भी लुप्त हो जाता है। फिर एक अनुमान लगाता है

गुणांक u1, u2, …, un निर्धारित किया जाता है जिससे कि सन्निकटन में त्रुटि प्रत्येक आधार फलन φi के लिए ऑर्थोगोनल हो :

संदृढ़ता आव्युह n-तत्व वर्ग आव्युह A द्वारा परिभाषित है

सदिश F घटकों के साथ परिभाषित करके गुणांक ui रेखीय प्रणाली Au = F द्वारा निर्धारित किए जाते हैं। संदृढ़ता आव्युह सममित आव्युह है, अर्थात। Aij = Aji, इसलिए इसके सभी स्वदेशी मूल्य ​​​​वास्तविक हैं। इसके अतिरिक्त, यह सख्ती से धनात्मक-निश्चित आव्युह है, जिससे कि प्रणाली Au = F के पास सदैव एक अद्वितीय समाधान होता है। (अन्य समस्याओं के लिए, यह अच्छी संपत्तियाँ खो जाएँगी।)

ध्यान दें कि संदृढ़ता आव्युह डोमेन के लिए उपयोग किए गए कम्प्यूटेशनल ग्रिड और किस प्रकार के परिमित तत्व का उपयोग किया जाता है, इसके आधार पर भिन्न होगा। उदाहरण के लिए, जब टुकड़ेवार द्विघात परिमित तत्वों का उपयोग किया जाता है तब संदृढ़ता आव्युह में टुकड़ेवार रैखिक तत्वों की तुलना में स्वतंत्रता की अधिक डिग्री होगी।

अन्य समस्याओं के लिए संदृढ़ता आव्युह

अन्य पीडीई के लिए संदृढ़ता आव्युह का निर्धारण अनिवार्य रूप से ही प्रक्रिया का पालन करता है, किन्तु यह सीमा स्थितियों की पसंद से समष्टि हो सकता है। अधिक समष्टि उदाहरण के रूप में, अण्डाकार समीकरण पर विचार करें

कहाँ x डोमेन में प्रत्येक बिंदु के लिए परिभाषित धनात्मक-निश्चित आव्युह है हम रॉबिन सीमा शर्त क्रियान्वित करते हैं

कहाँ νk k-वीं दिशा में इकाई जावक सामान्य सदिश ν का घटक है हल करने की प्रणाली है

जैसा कि ग्रीन की पहचान के एनालॉग का उपयोग करके दिखाया जा सकता है। गुणांक ui अभी भी रैखिक समीकरणों की प्रणाली को हल करके पाए जाते हैं, किन्तु प्रणाली का प्रतिनिधित्व करने वाला आव्युह सामान्य पॉइसन समस्या से स्पष्ट रूप से भिन्न है।

सामान्यतः, क्रम 2k के प्रत्येक अदिश अण्डाकार ऑपरेटर L के लिए, सोबोलेव स्पेस Hk पर एक द्विरेखीय रूप B जुड़ा होता है, ताकि समीकरण Lu = f का अशक्त सूत्रीकरण हो।

सभी कार्यों के लिए v में Hk. फिर इस समस्या के लिए संदृढ़ता आव्युह है

संदृढ़ता आव्युह की व्यावहारिक असेंबली

कंप्यूटर पर परिमित तत्व विधि को क्रियान्वित करने के लिए, किसी को पहले आधार कार्यों का समूह चुनना होगा और फिर संदृढ़ता आव्युह को परिभाषित करने वाले इंटीग्रल्स की गणना करनी होगी। सामान्यतः, डोमेन Ω को जाल निर्माण के कुछ रूपों द्वारा विभेदित किया जाता है, जिसमें इसे गैर-अतिव्यापी त्रिभुज जाल या जाल के प्रकारों में विभाजित किया जाता है, जिन्हें सामान्यतः तत्वों के रूप में जाना जाता है। फिर आधार कार्यों को प्रत्येक तत्व के अंदर कुछ क्रम के बहुपद और तत्व सीमाओं के पार निरंतर चुना जाता है। सबसे सरल विकल्प त्रिकोणीय तत्वों के लिए टुकड़ावार रैखिक फलन और आयताकार तत्वों के लिए टुकड़ावार द्विरेखीय हैं।

तत्व Tk के लिए तत्व संदृढ़ता आव्युह A[k]आव्युह है

i और j, अधिकांश मानों के लिए तत्व संदृढ़ता आव्युह शून्य है जिसके लिए संबंधित आधार फलन Tk के भीतर शून्य हैं पूर्ण संदृढ़ता आव्युह A तत्व संदृढ़ता आव्युह का योग है। विशेष रूप से, उन आधार कार्यों के लिए जो केवल स्थानीय रूप से समर्थित हैं, संदृढ़ता आव्युह विरल है।

आधार कार्यों के अनेक मानक विकल्पों के लिए, अर्थात त्रिकोणों पर टुकड़े-टुकड़े रैखिक आधार कार्यों के लिए, तत्व संदृढ़ता आव्युह के लिए सरल सूत्र हैं। उदाहरण के लिए, टुकड़ों में रैखिक तत्वों के लिए, शीर्षों (x1, y1), (x2, y2), (x3, y3) वाले त्रिभुज पर विचार करें और 2×3 आव्युह को परिभाषित करें

फिर तत्व संदृढ़ता आव्युह है

जब अंतर समीकरण अधिक समष्टि होता है, मान लीजिए कि अमानवीय प्रसार गुणांक होता है, तब तत्व संदृढ़ता आव्युह को परिभाषित करने वाले अभिन्न अंग का मूल्यांकन गॉसियन चतुर्भुज द्वारा किया जा सकता है।

संदृढ़ता आव्युह की स्थिति संख्या संख्यात्मक ग्रिड की गुणवत्ता पर दृढ़ता से निर्भर करती है। विशेष रूप से, परिमित तत्व जाल में छोटे कोण वाले त्रिकोण संदृढ़ता आव्युह के बड़े आइगेनवैल्यू ​​​​को प्रेरित करते हैं, जिससे समाधान की गुणवत्ता खराब हो जाती है।

संदर्भ

  • Ern, ए.; गुरमोंड, जे.-एल. (2004), परिमित तत्वों का सिद्धांत और अभ्यास, न्यूयॉर्क, एनवाई: स्प्रिंगर-वेरलाग, ISBN 0387205748
  • गोकेनबैक, एम.एस. (2006), परिमित तत्व विधि को समझना और लागू करना, फिलाडेल्फिया, पीए: एस.आई.ए.एम, ISBN 0898716144
  • ग्रॉसमैन, सी.; रूस, एच.-जी.; स्टाइन्स, एम. (2007), आंशिक विभेदक समीकरणों का संख्यात्मक उपचार, बर्लिन, जर्मनी: स्प्रिंगर-वेरलाग, ISBN 978-3-540-71584-9
  • जॉनसन, सी. (2009), परिमित तत्व विधि द्वारा आंशिक विभेदक समीकरणों का संख्यात्मक समाधान, डोवर, ISBN 978-0486469003
  • ज़िएनकिविज़, ओ.सी.; टेलर, आर.एल.; Zhu, जे.जेड. (2005), परिमित तत्व विधि: इसका आधार और बुनियादी बातें (6th ed.), ऑक्सफोर्ड, यूके: एल्सेवियर बटरवर्थ-हेनमैन, ISBN 978-0750663205