編計算力學程序犯過的錯誤

來自專欄 無力不起早編過些IGA的程序，由於本身就是等參元的結構，除了單元影響域和控制點外其與等參有限元並無其他區別，犯過最嚴重的錯誤就是在求解形函數對坐標導數的鏈式法則出了問題。待續吧，但願少犯錯誤，多寫經驗少寫教訓。

超彈性材料本構關係多用應變能來給出而且是全量形式，一般用三個不變數來表示，文獻中寫的是G的三個不變數（具體形式是兩個偏不變數加第一不變數），G means Cauchy Green Deformation Tensor 誤以為是格林應變，之後發現問題請教大神才知道原來是右柯西格林應變張量C，看來統一符號很重要O__O "…

切應力應變張量與矢量之間的轉換關係中間有個二倍 需要注意

本構在更新構型之後需要進行一步變換，原因在於坐標變換時第二pk應力和格林應變張量施加的變換是不同的

在更換構型後關聯節點信息肯定全變了 剛度陣也要進行相應轉換

在有限元類方法中，形函數的光滑性與邊界條件施加難度成正比，例如有限元與無網格，而等幾何分析同時做到形函數光滑性與邊界條件的方便施加，但形函數卻被局限於張量積（大單元）中，這大概就是這一類數值方法的魚與熊掌吧。

推薦閱讀：