曾攀老師的《有限元分析基礎教程》第三章有二維桿單元的推導，並結合一個例題進行了解析解和基於Matlab的程序求解。但是我感覺書中的MATLAB代碼有點羅嗦，而且一些實現方法也比較麻煩，比如已經知道了桿單元的起點和終點坐標，仍然需要另外給出單元局部坐標與整體坐標的夾角，這完全沒必要。於是我就用 ...

《有限元編程：菜鳥篇》 一、前言 相信很多做過有限差分之后又想做做有限元的初學者會有和我一樣的困惑，能看懂有限元算法的理論分析，但是真正應用到實際編程當中之前心里發怵，請教學過有限元程序的同學的時候，他們往往會，這個怎么怎么的簡單，這個你怎么能不會？這個不就是什么什么嗎bulabula... ...

對於偏微分方程，使其解成為唯一的輔助條件可分為兩種：一種是表達場的邊界所處的物理情況、稱為邊界條件；另一種是確定場的初始狀態，稱為初始條件。邊界條件和初始條件合稱為定解條件。 2.1 二維電磁場基本理論 電磁場的經典描述是麥克斯韋方程組，電機電磁場分析—般采用位函數表示，位函數 比場量本身更容易 ...

效果圖： ...

不得不說，Mathematica真是個好東西，以前學習有限元的時候，對於書中的方程推導，看到了就看過去了，從沒有想過要自己推導一遍，原因是手工推導太復雜。有了MM，原來很復雜的東西突然變得簡單了。 1.單元幾何描述 上圖是純彎梁單元，長度l，彈模E，面積A，慣性矩I。兩個節點1和2的位移列 ...

域離散為有限個單元，（比如，三角形面片、四面體單元，等。）離散域中的點（網格節點），稱之為 結點。 ...

一、“近似”的兩種分類 一個復雜的函數，可以通過一系列的“基底函數”（base function）的組合來近似，也就是函數逼近，有兩種典型的方法： 基於全域的逼近，如傅立葉級數展開； 基於子域的分段函數組合，如有限元方法。 第一種函數逼近方式，就是力學分析中經 ...

p47.（實習題-李榮華）用線性元求下列邊值問題的數值解 ...