材料的彈性模量有時候隨坐標是變化的，例如梯度涂層等，這個時候就需要借助子程序來實現了，下面是成功的模型結果

模量隨y坐標變化的模型的狀態變量

無梯度代表著均勻模量模型，有梯度代表模量從下到上隨y坐標變化（如果想讓它隨x坐標也變化，可以修改程序，很簡單）

模量隨y坐標成拋物線變化，底端固定，上面施加拉伸載荷

結果分析：

1 狀態變量值的大小代表了模量的大小（程序中設置E與狀態變量是線性關系），從狀態變量的云圖結果可以看出，底端模量最小，SDV2值最小，隨著y增加，開始增加很慢，然后增加速率增大，頂端的SDV2值最大，由此可知，模量隨y拋物線變化程序是沒有問題的；

2 從y反向的彈性應變也可以看出：對于均勻的彈性模量模型，因為總體模量都是200GPa，所以相同的拉伸載荷之后，y方向的彈性應變更大，這也印證了程序的正確性；

以下打包文件中包含源程序和例子：

溫馨提示：

大家要學習這種思想

ABAQUS斷裂模擬收徒 ，保證快速學會各種ABAQUS斷裂模擬方法  1200/人（將享有各種插件以及程序，價值3000+、專門定制視頻、全程親自教學、各種模型調試及解答問題等等，傾囊相教）