本文原創，因為是自己的想法，難免有考慮不周的地方，如發現錯誤或者需要討論請發送郵件至only_xiaofei@hotmail.com

前段時間一直在思考C3D20單元的存在意義，因為如果將C3D20單元進行三個方向的剖分（也就是將單元長度減小一半）。它與8個C3D8單元幾乎是一致的，現在就來討論下這兩個單元。

C3D20的空間結構如下

它的形函數是這樣的
8個頂點上的形函數為
（這部分和C3D8一致）
其它中點處的形函數為
相比之下，C3D8單元就顯得很簡單了

它的形函數也要簡單很多
可以看出，C3D20單元在描述位移場的時候是二次的，而C3D8單元是一次，這會導致C3D20單元的邊界會出現曲面（線），而C3D8只能為直線。
這樣問題就變得簡單了，如果單純考慮兩個單元的計算成本的話，那么一個C3D20單元和8個C3D8單元就是等效的，因為兩個的剛度矩陣都是20X20階的，兩者在計算成本上沒有任何區別，如果考慮積分點引起單元剛度矩陣的精度問題，8個C3D8單元擁有32個積分點，其精度一定是大于擁有27個積分點的C3D20單元的。此外，影響精度的還有邊界的表達方式，實際模型的邊界形狀是類似于二次曲線的話，C3D20單元更有利于對位移場的描述。
如果考慮彎曲的話，C3D8單元在求解彎曲問題是有很大的剪切自鎖誤差，而C3D20單元由于擁有眾多的積分點，可以大幅度削弱剪切鎖定，但是由于C3D8單元的尺寸只是C3D20單元尺寸的一半，因此8個C3D8單元的誤差與一個C3D20單元的誤差相差多少，期待有人可以推導一下。