संरचनात्मक यांत्रिकी में परिमित तत्व विधि

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

इतिहास

परिमित विधि की उत्पत्ति संरचनाओं के आव्यूह विश्लेषण द्रारा चिन्हित की जा सकती है जहां एक विस्थापन या कठोरता आव्यूह प्रस्ताव की एक अवधारणा प्रस्तुत की गई थी। 1950 के दशक में अभियांत्रिकी विधियों के आधार पर परिमित तत्व अवधारणाएँ विकसित की गईं। परिमित तत्व पद्धति ने 1960 और 1970 के दशक में जॉन आरगाईरिस और सहकर्मियों द्वारा अपनी वास्तविक प्रेरणा प्राप्त की; जैसे रे डब्ल्यू क्लो द्वारा स्टटगार्ट विश्वविद्यालय में; कैलिफोर्निया विश्वविद्यालय, बर्कले में, ओल्गिएर्ड ज़िएनक्यूविज़ द्वारा, और सहकर्मी अर्नेस्ट हिंटन, ब्रूस आयरन्स फिलिप जी सियारलेट द्वारा स्वानसी विश्वविद्यालय में; पियरे-एंड-मैरी-क्यूरी विश्वविद्यालय में; कॉर्नेल विश्वविद्यालय में, रिचर्ड गैलाघेर और सहकर्मियों द्वारा इत्यादि। मूल कृतियाँ जैसे कि आरगाईरिस और क्लो के कार्य आज के परिमित तत्व संरचनात्मक विश्लेषण विधियों का आधार बन गए।

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

• द्वि-आयामी तत्व जो केवल झिल्ली क्रिया जैसे समतल तनाव या समतल विकृति द्वारा अंतस्तल बलों का विरोध करते हैं, और प्लेटें जो अनुप्रस्थ कतरनी और झुकने की क्रिया द्वारा अनुप्रस्थ भार का विरोध करती हैं। तथा उनके पास कई प्रकार के आकार हो सकते हैं जैसे समतल या घूर्णित त्रिकोण और चतुर्भुज। नोड्स को सामान्यतः तत्व के कोनों पर रखा जाता है, और यदि उच्च सटीकता के लिए आवश्यक हो, तो अतिरिक्त बिन्दु को तत्व किनारों के साथ या तत्व के भीतर भी रखा जा सकता है। तत्व वास्तविक परत मोटाई की मध्य-सतह पर स्थित होते हैं।
• झिल्लियों, मोटी प्लेटों, आवरणों और ठोसों जैसी अक्षीय समस्याओं के लिए टोरस के आकार के तत्व होते है। इन तत्वों का अन्तः वर्ग पहले वर्णित प्रकारों के समान है इस प्रकार पतली प्लेटों और गोले के लिए एक आयामी, और ठोस, मोटी प्लेटों और गोले के लिए द्वि-आयामी तत्व होते है ।
• 3-डी ठोस जैसे यंत्र घटकों, बांधों, तटबंध परिवहन या मिट्टी के द्रव्यमान प्रतिरूपण के लिए त्रि-आयामी तत्व का प्रयोग किया जाता है। सरल तत्व आकृतियों में चतुष्फलकीय और षट्फलकीय तत्व सम्मिलित हैं। नोड्स को शीर्ष और संभवतः तत्व के फलकों या तत्व के भीतर रखा जाता है।

तत्व अंतर्संबंध और विस्थापन

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

व्यावहारिक विचार

अनुप्रयोग के दृष्टिकोण से, प्रतिरूप को इस तरह से प्रारूप करना महत्वपूर्ण है:

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

बड़े पैमाने पर वाणिज्यिक सॉफ्टवेयर का संकुल प्रायः जाल उत्पन्न करने और इनपुट और आउटपुट के चित्रमय प्रदर्शन के लिए सुविधाएं प्रदान करते हैं, जो इनपुट डेटा और परिणामों की व्याख्या दोनों के सत्यापन की सुविधा प्रदान करते हैं।

परिमित तत्व विधि-विस्थापन सूत्रीकरण का सैद्धांतिक अवलोकन: तत्वों से, प्रतिरूप समाधान तक

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

संरचनात्मक प्रतिरूप के लिए आभासी कार्य का सिद्धांत बाहरी और आंतरिक आभासी कार्य की गणितीय पहचान को व्यक्त करता है:

${\mbox{External virtual work}}=\int _{V}\delta {\boldsymbol {\epsilon }}^{T}{\boldsymbol {\sigma }}\,dV$

(1)

दूसरे शब्दों में, बाह्य बलों के समुच्चय द्वारा तंत्र पर किए गए कार्य का योग तंत्र को बनाने वाले तत्वों में तनाव ऊर्जा के रूप में संग्रहीत कार्य के बराबर होता है।

उपरोक्त समीकरण के दाईं ओर के आभासी आंतरिक कार्य को अलग-अलग तत्वों पर किए गए आभासी कार्य का योग करके पाया जा सकता है। उत्तरार्द्ध की आवश्यकता है कि बल-विस्थापन कार्यों का उपयोग किया जाए जो प्रत्येक व्यक्तिगत तत्व के लिए प्रतिक्रिया का वर्णन करता है। इसलिए, संरचना के विस्थापन को सामूहिक रूप से असतत तत्वों की प्रतिक्रिया से वर्णित किया गया है। समीकरण मात्र एक समीकरण के अतिरिक्त संरचना के अलग-अलग तत्वों के छोटे डोमेन के लिए लिखे गए हैं जो पूरे प्रतिरूप के रूप में प्रतिक्रिया का वर्णन करता है। उत्तरार्द्ध के परिणामस्वरूप एक जटिल समस्या होगी, इसलिए परिमित तत्व विधि की उपयोगिता है, जैसा कि बाद के अनुभागों में दिखाया गया है, Eq.(1) प्रतिरूप के लिए निम्नलिखित शासी संतुलन समीकरण की ओर जाता है:

$\mathbf {R} =\mathbf {Kr} +\mathbf {R} ^{o}$

(2)

जहाँ

$\mathbf {R}$ = नोडल बलों का सदिश , प्रतिरूप के नोड्स पर लागू बाहरी बलों का प्रतिनिधित्व करता है।

$\mathbf {K}$ = प्रतिरूप कठोरता आव्यूह, जो अलग-अलग तत्वों की कठोरता आव्यूह का सामूहिक प्रभाव है:$\mathbf {k} ^{e}$.

$\mathbf {r}$ = प्रतिरूप के नोडल विस्थापन का सदिश ।

$\mathbf {R} ^{o}$ = समतुल्य नोडल बलों के सदिश , नोडल बलों के अलावा अन्य सभी बाहरी प्रभावों का प्रतिनिधित्व करते हैं जो पहले से ही पूर्ववर्ती नोडल बल सदिश आर में शामिल हैं। इन बाहरी प्रभावों में वितरित या केंद्रित सतह बल, शरीर बल, थर्मल प्रभाव, प्रारंभिक तनाव और तनाव शामिल हो सकते हैं।

प्रतीकात्मक रूप से एक बार समर्थन की बाधाओं के लिए जिम्मेदार होने के बाद, रैखिक समीकरणों की प्रतिरूप को हल करके नोडल विस्थापन पाया जाता है

$\mathbf {r} =\mathbf {K} ^{-1}(\mathbf {R} -\mathbf {R} ^{o})$

(3)

इसके बाद, अलग-अलग तत्वों में तनाव निम्नानुसार पाया जा सकता है:

$\mathbf {\epsilon } =\mathbf {Bq}$

(4)

$\mathbf {\sigma } =\mathbf {E} (\mathbf {\epsilon } -\mathbf {\epsilon } ^{o})+\mathbf {\sigma } ^{o}=\mathbf {E} (\mathbf {Bq} -\mathbf {\epsilon } ^{o})+\mathbf {\sigma } ^{o}$

(5)

जहाँ

$\mathbf {q}$ = एक नोडल विस्थापन का सदिश - प्रतिरूप विस्थापन सदिश आर का एक उपसमुच्चय जो विचाराधीन तत्वों से संबंधित है।

$\mathbf {B}$ = तनाव-विस्थापन आव्यूह जो तत्व में किसी भी बिंदु पर नोडल विस्थापन क्यू को उपभेदों में बदल देता है।

$\mathbf {E}$ = लोच आव्यूह जो प्रभावी उपभेदों को तत्व में किसी भी बिंदु पर तनाव में बदल देता है।

$\mathbf {\epsilon } ^{o}$ = तत्वों में प्रारंभिक उपभेदों का सदिश ।

$\mathbf {\sigma } ^{o}$ = तत्वों में प्रारंभिक तनाव का सदिश ।

आभासी कार्य समीकरण को लागू करने से (1) प्रतिरूप के लिए, हम तत्व आव्यूह स्थापित कर सकते हैं जहां $\mathbf {B}$, $\mathbf {k} ^{e}$ साथ प्रतिरूप मैट्रिसेस को समन्वायोजन करने की तकनीक $\mathbf {R} ^{o}$ और $\mathbf {K}$. अन्य मैट्रिसेस जैसे $\mathbf {\epsilon } ^{o}$, $\mathbf {\sigma } ^{o}$, $\mathbf {R}$ और $\mathbf {E}$ ज्ञात मूल्य हैं और इन्हें सीधे डेटा इनपुट से सेट किया जा सकता है।

प्रक्षेप या आकृति कार्य

मान लीजिए ; $\mathbf {q}$ एक विशिष्ट तत्व के नोडल विस्थापन के सदिश है। तत्व के किसी भी अन्य बिंदु पर विस्थापन प्रक्षेप कार्यों के उपयोग से प्रतीकात्मक रूप से पाया जा सकता है:

$\mathbf {u} =\mathbf {N} \mathbf {q}$

(6)

कहाँ

$\mathbf {u}$ = तत्व के किसी बिंदु {x, y, z} पर विस्थापन का सदिश।

$\mathbf {N}$ = प्रक्षेप कार्यों के रूप में कार्य करने वाले आकृति कार्यों का आव्यूह।

समीकरण (6) बहुत रुचि की अन्य मात्राओं को जन्म देता है:

आभासी विस्थापन जो आभासी नोडल विस्थापन का एक कार्य है:

$\delta \mathbf {u} =\mathbf {N} \delta \mathbf {q}$

(6b)

तत्वों में तनाव जो तत्व के नोड्स के विस्थापन से उत्पन्न होते हैं:

$\mathbf {\epsilon } =\mathbf {Du} =\mathbf {DNq}$

(7)

जहाँ $\mathbf {D}$ = तनाव-विस्थापन संबंधो का आव्यूह जो विस्थापन को रैखिक लोच सिद्धांत का उपयोग करके तनाव में परिवर्तित करता है। समीकरण (7) से पता चलता है कि आव्यूह बी में (4) है

$\mathbf {B} =\mathbf {DN}$

(8)

तत्व के आभासी नोडल विस्थापन के अनुरूप आभासी तनाव:

$\delta {\boldsymbol {\epsilon }}=\mathbf {B} \delta \mathbf {q}$

(9)

मात्रा के एक विशिष्ट तत्व के लिए $V^{e}$, आभासी विस्थापन के कारण आंतरिक आभासी कार्य (5) और (9) में (1) के प्रतिस्थापन द्वारा प्राप्त किया जाता है:

${\mbox{Internal virtual work}}=\int _{V^{e}}\delta {\boldsymbol {\epsilon }}^{T}{\boldsymbol {\sigma }}\,dV^{e}=\delta \ \mathbf {q} ^{T}\int _{V^{e}}\mathbf {B} ^{T}{\big \{}\mathbf {E} (\mathbf {Bq} -\mathbf {\epsilon } ^{o})+\mathbf {\sigma } ^{o}{\big \}}\,dV^{e}$

(10)

एलिमेंट मेट्रिसेस

मुख्य रूप से संदर्भ की सुविधा के लिए, विशिष्ट तत्वों से संबंधित निम्नलिखित आव्यूह को अब परिभाषित किया जा सकता है:

तत्व कठोरता आव्यूह

$\mathbf {K} ^{e}=\int _{V^{e}}\mathbf {B} ^{T}\mathbf {E} \mathbf {B} \,dV^{e}$

(11)

समतुल्य तत्व भार सदिश

$\mathbf {Q} ^{oe}=\int _{V^{e}}-\mathbf {B} ^{T}{\big (}\mathbf {E} \mathbf {\epsilon } ^{o}-\mathbf {\sigma } ^{o}{\big )}\,dV^{e}$

(12)

संख्यात्मक एकीकरण के लिए गॉसियन चतुर्भुज का उपयोग करके सामान्यतः इन मेट्रिसेस का संख्यात्मक रूप से मूल्यांकन किया जाता है। उनका उपयोग निम्नलिखित को सरल करता है (10)

${\mbox{Internal virtual work}}=\delta \ \mathbf {q} ^{T}{\big (}\mathbf {K} ^{e}\mathbf {q} +\mathbf {Q} ^{oe}{\big )}$

(13)

प्रतिरूप नोडल विस्थापन के संदर्भ में तत्व आभासी कार्य

चूंकि नोडल विस्थापन सदिश क्यू प्रतिरूप नोडल विस्थापन आर का एक उपसमुच्चय है, हम नए कॉलम और शून्य की पंक्तियों के साथ तत्व आव्यूह के आकार का विस्तार करके क्यू को आर से बदल सकते हैं:

${\mbox{Internal virtual work}}=\delta \ \mathbf {r} ^{T}\left(\mathbf {K} ^{e}\mathbf {r} +\mathbf {Q} ^{oe}\right)$

(14)

जहां, सरलता के लिए, हम तत्व आव्यूहों के लिए उन्हीं प्रतीकों का उपयोग करते हैं, जिनका आकार अब विस्तारित हो गया है और साथ ही पंक्तियों और स्तंभों को उचित रूप से पुनर्व्यवस्थित किया गया है।

प्रतिरूप आभासी कार्य

सभी तत्वों के लिए आंतरिक आभासी कार्य (14) को समेटने से (1) का दाहिना हाथ मिलता है:

${\mbox{System internal virtual work}}=\sum _{e}\delta \ \mathbf {r} ^{T}\left(\mathbf {k} ^{e}\mathbf {r} +\mathbf {Q} ^{oe}\right)=\delta \ \mathbf {r} ^{T}\left(\sum _{e}\mathbf {k} ^{e}\right)\mathbf {r} +\delta \ \mathbf {r} ^{T}\sum _{e}\mathbf {Q} ^{oe}$

(15)

अब (1) के बायीं ओर को ध्यान में रखते हुए, प्रतिरूप बाहरी आभासी कार्य में निम्न सम्मिलित हैं:

नोडल बलों R द्वारा किया गया कार्य:

$\delta \ \mathbf {r} ^{T}\mathbf {R}$

(16)

बाह्य बलों द्वारा किया गया कार्य $\mathbf {T} ^{e}$ भार पर $\mathbf {S} ^{e}$ तत्वों के किनारों या सतहों और शरीर बलों द्वारा $\mathbf {f} ^{e}$

$\sum _{e}\int _{S^{e}}\delta \ \mathbf {u} ^{T}\mathbf {T} ^{e}\,dS^{e}+\sum _{e}\int _{V^{e}}\delta \ \mathbf {u} ^{T}\mathbf {f} ^{e}\,dV^{e}$

का प्रतिस्थापन (6b) देता है:

$\delta \ \mathbf {q} ^{T}\sum _{e}\int _{S^{e}}\mathbf {N} ^{T}\mathbf {T} ^{e}\,dS^{e}+\delta \ \mathbf {q} ^{T}\sum _{e}\int _{V^{e}}\mathbf {N} ^{T}\mathbf {f} ^{e}\,dV^{e}$

या

$-\delta \ \mathbf {q} ^{T}\sum _{e}\left(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)$

(17a)

जहां हमने नीचे परिभाषित अतिरिक्त तत्व के मैट्रिसेस प्रस्तुत किए हैं:

$\mathbf {Q} ^{te}=-\int _{S^{e}}\mathbf {N} ^{T}\mathbf {T} ^{e}\,dS^{e}$

(18a)

$\mathbf {Q} ^{fe}=-\int _{V^{e}}\mathbf {N} ^{T}\mathbf {f} ^{e}\,dV^{e}$

(18b)

$-\delta \ \mathbf {r} ^{T}\sum _{e}\left(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)$

(17b)

पुनः,संख्यात्मक एकीकरण उनके मूल्यांकन के लिए सुविधाजनक है। क्यू का एक समान प्रतिस्थापन (17a) r के साथ सदिशों को पुनर्व्यवस्थित और विस्तारित करने के बाद देता है $\mathbf {Q} ^{te},\mathbf {Q} ^{fe}$:

$-\delta \ \mathbf {r} ^{T}\sum _{e}\left(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)$

(17b)

प्रतिरूप मैट्रिसेस की समन्वायोजन

जोड़ना (16), (17b) और योग के बराबर (15) देता है: $\delta \ \mathbf {r} ^{T}\mathbf {R} -\delta \ \mathbf {r} ^{T}\sum _{e}\left(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)=\delta \ \mathbf {r} ^{T}\left(\sum _{e}\mathbf {k} ^{e}\right)\mathbf {r} +\delta \ \mathbf {r} ^{T}\sum _{e}\mathbf {Q} ^{oe}$आभासी विस्थापन के बाद से $\delta \ \mathbf {r}$ मनमाने हैं, पूर्ववर्ती समानता कम हो जाती है:

$\mathbf {R} =\left(\sum _{e}\mathbf {k} ^{e}\right)\mathbf {r} +\sum _{e}\left(\mathbf {Q} ^{oe}+\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)$इसके साथ तुलना (2) पता चलता है कि:

• प्रतिरूप कठोरता आव्यूह तत्वों की कठोरता आव्यूह को जोड़कर प्राप्त की जाती है:

  $\mathbf {K} =\sum _{e}\mathbf {k} ^{e}$

• समतुल्य नोडल बलों का सदिश तत्वों के भार को जोड़कर प्राप्त किया जाता है:

  $\mathbf {R} ^{o}=\sum _{e}\left(\mathbf {Q} ^{oe}+\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)$

व्यवहार में, तत्व मैट्रिसेस न तो विस्तारित होते हैं और न ही पुनर्व्यवस्थित होते हैं। इसके अतिरिक्त, प्रतिरूप कठोरता आव्यूह $\mathbf {K}$ अलग-अलग गुणांक जोड़कर इकट्ठा किया जाता है ${k}_{ij}^{e}$ को ${K}_{kl}$ जहां सबस्क्रिप्ट ij, kl का अर्थ है कि तत्व का नोडल विस्थापन ${q}_{i}^{e},{q}_{j}^{e}$ प्रतिरूप के नोडल विस्थापन के साथ क्रमशः मेल खाते हैं ${r}_{k},{r}_{l}$. इसी प्रकार, $\mathbf {R} ^{o}$ अलग-अलग गुणांक जोड़कर इकट्ठा किया जाता है ${Q}_{i}^{e}$ को ${R}_{k}^{o}$ जहाँ ${q}_{i}^{e}$ माचिस ${r}_{k}$. इसका सीधा जोड़ ${k}_{ij}^{e}$ में ${K}_{kl}$ प्रक्रिया को प्रत्यक्ष कठोरता विधि का नाम देता है।

यह भी देखें

• सीमित तत्व विधि
• लचीलापन विधि
• आव्यूह कठोरता विधि
• FEM का उपयोग करके मोडल विश्लेषण
• परिमित तत्व सॉफ्टवेयर पैकेजों की सूची
• संरचनात्मक विश्लेषण
• आभासी कार्य
• अंतराल परिमित तत्व

संदर्भ