清明節作業

核動力三輪車

+ h$ }$ M* N6 t; D* b! R* c/ x/ `
" X: P3 w- O# i& d) E& T3 y6 o清明節沒有回老家，窗外雨下個不停。自己就在家里宅到現在。清靜的時候最好學習，做了一份作業，鞏固下基礎，有些收獲和大家分享。

題目來自大學課本，劉鴻文教授寫的《材料力學》第二冊。原題用的是有限差分解法。有限差分法數學意義很明確，解法也比較簡單。但用過都知道，有限差分法編程不方便。當然，對于簡單的問題，也不一定就要編程。基于假設位移場和最小勢能原理的瑞利里茲法就是一種很好用的方法。在有限元法出現前，和有限差分法一同是很好用的方法。第一張圖里我用的就是瑞利里茲法。這種方法優點就是物理意義清楚。但缺點也很明顯，如圖里所示。我只能取自己假設的位移場v(x)的第一項。如果要取第二項，那么計算就會很繁瑣。幸運的是，即便只取第一項，位移的結果也是足夠精確的。
第二張圖里我用的是有限元法。有限元法可以看作傳統的瑞利里茲法改進。如圖一所示，假設了一個函數v(x)描述軸的每一點的位移。對于高維問題、復雜結構和復雜邊界條件，就可能需要構建極為復雜的場<u(x,y,z,t,f),v(x,y,z,t,f),w(x,y,z,t,f)>描述整個的構建。而更要命的是，這個過程必須人工干預，所以是很難實現的。如果能把瑞利里茲法應用在一個簡單的結構中，而想辦法把一個極為復雜的結構分解為若干簡單的結構組合那不就可以保留這種方法的優點嗎？這在數學上的提法叫做分片插值，這就是有限元法的本質。有的書說有限元法來自瑞利里茲法，有的書說來自伽遼金法，其實都不妥帖。發明有限法的目的是為了在計算機上繼續應用上述兩種方法，這不是要求把算法復雜化，反而是要求把算法簡化，流程化。計算機適合做簡單枯燥的循環，所以必須把太靈活的傳統方法“簡化”為計算機上可以實現的算法。圖三是最終形成的剛度矩陣。求解過程中要對矩陣進行重排、分塊、求逆、求乘以及后處理。對于這根軸的計算，有限元法大幅增加了計算量，但把人從枯燥的計算和針對復雜對象的不可能任務中解放了出來。
( H3 K% p3 [3 I( ~& Y9 [在這道題中，采用的是歐拉粱單元。既然歐拉梁是在xy坐標系中定義的，為什么要用截面轉交和撓度描述其變形呢？一個很重要的原因就是這樣做更簡潔。。。由于theta=dy/dx，連續性條件要求假設的位移場存在一階導數，而應變能中包含位移場的二階導數，這要求位移場存在二階導數。那么，如圖二所示，合適的位移場是三次多項式。當然，也幸好只是三次的，高次的很難處理。
這道題還有其他的解法，比如采用常應變三角形——我非常想嘗試。這可以練習使用形函數、等參單元，并且可以用到雅可比矩陣。但可以預見的是，即便網格更細密的三角形單元，也很難得到粱單元的精度。這提示我們有限元分析時一定要采用合理的單元類型。還有就是直接積分，大學時這樣干過，寫了好幾頁。當然，伽遼金法也可以。
( {0 U0 X7 q, U$ q說了這些，給大家出道題：圖二的底部是節點位移，怎么從節點廣義位移得到節點位移呢？對于這道題，劃分更細密的粱單元，可以得到更精確的結果嗎？

qinghuap123

不明覺厲

hkslljf

大學學的材料力學好久都沒有用過了。難道是我們公司太低端了？