(完整版)計(jì)算力學(xué)復(fù)習(xí)題答案

計(jì)算力學(xué)試題答案1.有限單元法和經(jīng)典Ritz法的主要區(qū)別是什么？答：經(jīng)典Ritz法是在整個(gè)區(qū)域內(nèi)假設(shè)未知函數(shù)，適用于邊界幾何形狀簡單的情形；有限單元法是將整個(gè)區(qū)域離散，分散成若干個(gè)單元，在單元上假設(shè)未知函數(shù)。有限單元法是單元一級的Ritz法。2、單元?jiǎng)偠染仃嚭驼w剛度矩陣各有什么特征？剛度矩陣[K]奇異有何物理意義？在求解問題時(shí)如何消除奇異性？答：單元?jiǎng)偠染仃嚨奶卣鳎孩艑ΨQ性⑵奇異性⑶主元恒正⑷平面圖形相似、彈性矩陣D、厚度t相同的單元，相同⑸的分塊子矩陣按結(jié)點(diǎn)號排列，每一子矩陣代表一個(gè)結(jié)點(diǎn)，占兩行兩列，其位置與結(jié)點(diǎn)位置對應(yīng)。整體剛度矩陣的特征：⑴對稱性⑵奇異性⑶主元恒正⑷稀疏性⑸非零元素呈帶狀分布。的物理意義是任意給定結(jié)構(gòu)的結(jié)點(diǎn)位移所得到的結(jié)構(gòu)結(jié)點(diǎn)力總體上滿足力和力矩的平衡。為消除的奇異性，需要引入邊界條件，至少需給出能限制剛體位移的約束條件。3.列式說明乘大數(shù)法引入給定位移邊界條件的原理？答：設(shè)：，則將即：修改后的第個(gè)方程為由于得所以對于多個(gè)給定位移時(shí)，則按序?qū)⒚總€(gè)給定位移都作上述修正，得到全部進(jìn)行修正后的K和P，然后解方程即可得到包括給定位移在內(nèi)的全部結(jié)點(diǎn)位移值。4.何為等參數(shù)單元？為什么要引入等參數(shù)單元？答：等參變換是對單元的幾何形狀和單元內(nèi)的場函數(shù)采用相同數(shù)目的結(jié)點(diǎn)參數(shù)及相同的插值函數(shù)進(jìn)行變換，采用等參變換的單元稱之為等參數(shù)單元。借助于等參數(shù)單元可以對于一般的任意幾何形狀的工程問題和物理問題方便地進(jìn)行有限元離散，其優(yōu)點(diǎn)有：對單元形狀的適應(yīng)性強(qiáng)；單元特性矩陣的積分求解方便（積分限標(biāo)準(zhǔn)化）；便于編制通用化程序。5、對于平面4節(jié)點(diǎn)（線性）和8節(jié)點(diǎn)（二次）矩形單元，為了得到精確的剛度矩陣，需要多少個(gè)Gauss積分點(diǎn)？說明理由。答：對于平面4節(jié)點(diǎn)（線性）矩形單元：所以因而積分點(diǎn)數(shù)為：矩陣對于平面8節(jié)點(diǎn)（二次）矩形單元：所以因而積分點(diǎn)數(shù)為：矩陣⑴矩形、正方形、平行四邊形1、有限單元法的解題步驟如何？它與經(jīng)典Ritz法有何區(qū)別？答：⑴劃分單元，輸入結(jié)點(diǎn)和單元信息；⑵單元分析：⑶整體分析：⑷求解方程得到解引入位移邊界條件得到：⑸對位移結(jié)果進(jìn)行有關(guān)整理、計(jì)算單元或結(jié)點(diǎn)的應(yīng)力、應(yīng)變2、總剛度矩陣[K]的任一元素kij的物理意義是什么？如何解釋總剛度矩陣的奇異性和帶狀稀疏性？答：K中元素的物理意義：當(dāng)結(jié)構(gòu)的第個(gè)結(jié)點(diǎn)位移方向上發(fā)生單位位移，而其它結(jié)點(diǎn)位移方向上位移為零時(shí)，需在第個(gè)結(jié)點(diǎn)位移方向上施加的結(jié)點(diǎn)力大小。奇異性：=0,力學(xué)意義是對任意給定的結(jié)點(diǎn)位移所得到的結(jié)構(gòu)結(jié)點(diǎn)力總體上是滿足力和力矩的平衡。反之，給定任意滿足力和力矩平衡的結(jié)點(diǎn)載荷P，由于K的奇異性卻不能解得結(jié)構(gòu)的位移,因而結(jié)構(gòu)仍可能發(fā)生任意的剛體位移。為消除的奇異性，結(jié)構(gòu)至少需給出能限制剛體位移的約束條件。帶狀稀疏性：由于連續(xù)體離散為有限個(gè)單元體時(shí)，每個(gè)結(jié)點(diǎn)的相關(guān)單元只是圍繞在該結(jié)點(diǎn)周圍為數(shù)甚少的幾個(gè)，一個(gè)結(jié)點(diǎn)通過相關(guān)單元與之發(fā)生關(guān)系的相關(guān)結(jié)點(diǎn)也只是它周圍的少數(shù)幾個(gè)，因此雖然總體單元數(shù)和結(jié)點(diǎn)數(shù)很多，結(jié)構(gòu)剛度矩陣的階數(shù)很高，但剛度系數(shù)中非零系數(shù)卻很少，即為總剛度矩陣的稀疏性。另外，只要結(jié)點(diǎn)編號是合理的，這些稀疏的非零元素將集中在以主對角線為中心的一條帶狀區(qū)域內(nèi)，即為總剛度矩陣的帶狀分布特性。3、以3節(jié)點(diǎn)三角形單元為例證明插值函數(shù)特性答：圖形見課本P105圖3.6，n為節(jié)點(diǎn)數(shù)。由面積坐標(biāo)：插值函數(shù)：所以4、什么是等參單元？等參單元的收斂性如何？答：等參變換是對單元的幾何形狀和單元內(nèi)的場函數(shù)采用相同數(shù)目的結(jié)點(diǎn)參數(shù)及相同的插值函數(shù)進(jìn)行變換，采用等參變換的單元稱之為等參元。等參單元滿足收斂性需滿足兩個(gè)條件：即單元必須是協(xié)調(diào)的和完備的。完備性條件：要求插值函數(shù)中包含完全的線性項(xiàng)（包含常數(shù)項(xiàng)和一次項(xiàng)）。協(xié)調(diào)性條件：單元邊界上位移連續(xù)，相鄰單元邊界具有相同的結(jié)點(diǎn)，每一單元沿邊界的坐標(biāo)和未知函數(shù)采用相同的插值函數(shù)。5、對于空間8節(jié)點(diǎn)（線性）和20節(jié)點(diǎn)（二次）六面體單元，為了得到精確的剛度矩陣，需要多少個(gè)Gauss積分點(diǎn)？說明理由。答：對于空間8節(jié)點(diǎn)（線性）六面體單元：所以因而積分點(diǎn)數(shù)為：矩陣對于空間20節(jié)點(diǎn)（二次）六面體單元：所以因而積分點(diǎn)數(shù)為：矩陣1、為什么說3節(jié)點(diǎn)三角形單元是常應(yīng)變單元？答：常應(yīng)變單元指的是在一個(gè)單元內(nèi)的應(yīng)變?yōu)槌?shù)，有限元中的常應(yīng)變單元指的是線性三角形單元，線性三角形單元的位移場為線性的，應(yīng)變?yōu)槲灰频囊浑A導(dǎo)數(shù)，故為常數(shù)，因此稱為常應(yīng)變單元。2.以平面4節(jié)點(diǎn)雙線性單元為例，說明形函數(shù)的兩個(gè)重要特性。答：圖形見課件2.5矩形單元插值函數(shù)（形函數(shù)）的性質(zhì)進(jìn)而驗(yàn)證插值函數(shù)的性質(zhì)：3、何為等參變換？等參元有那些優(yōu)點(diǎn)？答：等參變換是對單元的幾何形狀和單元內(nèi)的場函數(shù)采用相同數(shù)目的結(jié)點(diǎn)參數(shù)及相同的插值函數(shù)進(jìn)行變換，采用等參變換的單元稱之為等參元。借助于等參元可以對于一般的任意幾何形狀的工程問題和物理問題方便地進(jìn)行有限元離散，其優(yōu)點(diǎn)有：對單元形狀的適應(yīng)性強(qiáng)；單元特性矩陣的積分求解方便（積分限標(biāo)準(zhǔn)化）；便于編制通用化程序。二、如圖所示平面問題有限元網(wǎng)格，每個(gè)單元4個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)2個(gè)自由度，1.給出適當(dāng)?shù)墓?jié)點(diǎn)編號，使總的系數(shù)矩陣的半帶寬最小，并計(jì)算半帶寬的值；2.在您的節(jié)點(diǎn)編號下，圖中節(jié)點(diǎn)A的主對角線上的元素在總系數(shù)矩陣中的行號和列號如何？3.哪幾個(gè)單元對節(jié)點(diǎn)A的主對角線上的系數(shù)有非零貢獻(xiàn)？1、答：沿短邊回頭編號，存儲(chǔ)量最小。帶寬的計(jì)算：2、答：由得節(jié)點(diǎn)A的主對角線上的元素、在總系數(shù)矩陣中的行號為11和12,列號為11,12。3、答：2、3、4單元對A的主對角線上的系數(shù)有非零貢獻(xiàn)。注意：桿件單元在每個(gè)節(jié)點(diǎn)上有1個(gè)自由度，帶寬不用乘以2。三、圖示6節(jié)點(diǎn)三角形單元的142邊作用有均布側(cè)壓力q，單元厚度為t，求單元的等效節(jié)點(diǎn)載荷。見課件2.5.3—5采用面積坐標(biāo)時(shí)，單元矩陣的計(jì)算中：（2）均布側(cè)壓力q和（4）x方向三角形分布力中的形函數(shù)。四、見課件3.3.3五、答：⑴觀察插值函數(shù)，包含一次完全多項(xiàng)式，滿足完備性要求，但其含有三次項(xiàng)與，此點(diǎn)與八點(diǎn)矩形單元四個(gè)邊上均僅有三個(gè)點(diǎn)，最多滿足二次項(xiàng)矛盾，不滿足協(xié)調(diào)性，因而不滿足收斂性。⑵觀察差值函數(shù)，包含一次完全多項(xiàng)式，多項(xiàng)式最高次為2，明顯既滿足完備性和協(xié)調(diào)性，故收斂。第三頁P(yáng)DF：六、考慮等截面軸力桿單元，題圖中分別示出2節(jié)點(diǎn)和3節(jié)點(diǎn)單元體，1.寫出它們的位移插值函數(shù)；2.推導(dǎo)這兩種單元體的剛度矩陣；3.對3節(jié)點(diǎn)單元體用靜力凝聚法消去中間節(jié)點(diǎn)自由度，建立單元體剛度矩陣表達(dá)式。解：⑴圖a，令，則有。故有，，。圖b，令，則有。故有，⑵圖a：有應(yīng)變矩陣。。圖b：有應(yīng)變矩陣。。⑶，有代入與，消去，得：。二、用最小勢能原理推導(dǎo)平面彈性問題的有限元方程。解：利用最小勢能（位能）原理推導(dǎo)平面彈性問題有限元方程，見課本P61。第二章課件也有。三、如圖所示的平面內(nèi)部三角形單元網(wǎng)格每節(jié)點(diǎn)2個(gè)自由度，1.用圖中所給節(jié)點(diǎn)編號計(jì)算總剛矩陣的半帶寬；2.對節(jié)點(diǎn)重新編號，使結(jié)構(gòu)總剛矩陣的半帶寬最小，并說明此時(shí)中心節(jié)點(diǎn)主對角線上的元素在總剛矩陣中的行號和列號？3.是否所有單元對中心節(jié)點(diǎn)主對角線上的元素都有非零貢獻(xiàn)？4.兩種編號下，結(jié)構(gòu)總剛矩陣中的非零元素是否相等？為什么？解：⑴半帶寬（相鄰結(jié)點(diǎn)號碼的最大差值+1）自由度。⑵中心節(jié)點(diǎn)元素編號為5，其余元素編號順時(shí)針依次編寫為1-4，6-9。此時(shí)中心節(jié)點(diǎn)對應(yīng)主對角線上的元在總剛度矩陣中的行號、列號分別為第9行和第9列，第10行和第10列素、⑶是。由于8個(gè)三角型單元均有一個(gè)節(jié)點(diǎn)是中心節(jié)點(diǎn)主對角線上元素，有均不為零，故所有單元對主對角線上元素都有非零貢獻(xiàn)。⑷兩種編號方式，非零元素相等。編號的合理化只是將非零元素的位置集中在以主對角線為中心的一條帶狀區(qū)域內(nèi)，但并不改變非零元素的個(gè)數(shù)。四、3節(jié)點(diǎn)三角形單元的jm邊作用有分布側(cè)壓力，如圖示，單元厚度為t，求單元的等效節(jié)點(diǎn)載荷。解：如圖梯形分布力可以分解成為一個(gè)三角形分布力和一個(gè)均勻分布力，然后疊加。參加課本P68頁或者課件，易知：均勻分布部分，三角型分布部分：兩者相加，第四頁P(yáng)DF五、塊體棱柱單元的上下三角形表面?=±1，利用面積坐標(biāo)L(i＝1,2,3)i和等參坐標(biāo)?寫出該線性單元的形函數(shù)。解：；；；。六、對圖示的矩形單元采用如下的插值函數(shù)1.請分析該單元的協(xié)調(diào)性。2.若將插值多項(xiàng)式寫成解：⑴如圖所示的矩形單元若滿足協(xié)調(diào)性，因?yàn)樵趛方向有三個(gè)節(jié)點(diǎn)，x方向有兩個(gè)節(jié)點(diǎn)，故y方向上，差值函數(shù)視x為常數(shù)y的次數(shù)至少為3次，x方向上，插值函數(shù)視y為常數(shù)，x的次數(shù)至少為2次。對比給出的方程知，該單元若采用該插值函數(shù)滿足協(xié)調(diào)性。⑵插值函數(shù)改變后，增加了一個(gè)性。，而由圖中矩形單元知道x方向多項(xiàng)式最高次為2次，顯然不滿足協(xié)調(diào)七、證明4節(jié)點(diǎn)平行四邊形二維單元的雅可比矩陣是常數(shù)矩陣.對于二維平行四邊形有，有（7.1）等參變換下取，，而在自然坐標(biāo)下不妨設(shè)該平行四邊形的四邊的方程分別為：①②③④從而知道：；；。其中，，，，解得四個(gè)節(jié)點(diǎn)的坐標(biāo)，一并代入方程7.1得：的四個(gè)元素均為常數(shù)，故有4節(jié)點(diǎn)平行四邊形二維單元的雅可比矩陣是常數(shù)矩陣。第五頁pdf該單元在結(jié)構(gòu)中的位置使得總體節(jié)點(diǎn)編號分別為19、20、30、31，回答：1.在未引入邊界條件以前，j單元?jiǎng)偠染仃嘯K]e的系數(shù)將貢獻(xiàn)給總體剛度矩陣[K]中的哪些行、列？2.具體寫出剛度矩陣[K]e中的哪些元素對總體剛度矩陣[K]中的下列行和列有貢獻(xiàn)，（1）59行61列；（2）38行39列；（3）59行59列；（4）37行37列。（題解：⑴由節(jié)點(diǎn)位移矢量知該結(jié)構(gòu)節(jié)點(diǎn)自由度為2，設(shè)某一個(gè)節(jié)點(diǎn)的編號為，則該個(gè)節(jié)點(diǎn)將分別對總體剛度矩陣的第、行與列，產(chǎn)生貢獻(xiàn)。所以，總體編號為19、20、30、31的總體節(jié)點(diǎn)將相應(yīng)對總體剛度矩陣的37、38行和列，39、40行