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

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

इतिहास

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

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

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

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

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

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

एप्लिकेशन के दृष्टिकोण से, सिस्टम को इस तरह से मॉडल करना महत्वपूर्ण है:

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

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

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

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

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

${\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}$ = समतुल्य नोडल बलों के वेक्टर, नोडल बलों के अलावा अन्य सभी बाहरी प्रभावों का प्रतिनिधित्व करते हैं जो पहले से ही पूर्ववर्ती नोडल बल वेक्टर आर में शामिल हैं। इन बाहरी प्रभावों में वितरित या केंद्रित सतह बल, शरीर बल, थर्मल प्रभाव, प्रारंभिक तनाव और तनाव शामिल हो सकते हैं।

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

$\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)

फिर से, संख्यात्मक एकीकरण उनके मूल्यांकन के लिए सुविधाजनक है। क्यू का एक समान प्रतिस्थापन (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 का उपयोग करके मोडल विश्लेषण
• परिमित तत्व सॉफ्टवेयर पैकेजों की सूची
• संरचनात्मक विश्लेषण
• आभासी कार्य
• अंतराल परिमित तत्व

संदर्भ