سابروتینی برای تغییر مدول الاستیسیته نسبت به کرنش کل

سلام،

من یک مدل توی آباکوس ساختم و می خواهم ماده ای براش تعریف کنم که مدول الاستیسیته اش نسبت به کرنش کل (engineering strain) تغییر کنه.

این رابطه بین مدول الاستیسیته و کرنش کل است:

E [Gpa]      Engineering Strain
203050                   0
186750                0.01
174750                0.02
168450                0.03
163900                0.04
159800                0.05
156100                0.06
152850                0.07
148950                0.08

من یک سابروتین VUSDFLD نوشتم (به این دلیل که حلگر من Abaqus Explicit است و چون هیچ پیش زمینه ای در کد نویسی ندارم سعی کردم از کدهای پیچیده ای مانند VUMAT استفاده نکنم). سابروتینم را براتون به این پست پیوست کردم. حالا سوالم راجع بهش این است:

1- مدول الاستیسیته ی مدل من با کرنش توتال (که هم شامل الاستیک و پلاستیک میشه) تغییر می کنه. آیا سابروتین من هر دو قسمت کرنش را پوشش می دهد؟ و اینکه آیا اصلا کدی که نوشتم صحیحه؟

2- آیا من این رابطه ی بین مدول الاستیسیته و کرنش کل را که پیشتر بهش اشاره کردم رو در قسمت متغییر میدانی (field variable) قسمت الاستیکه خواص ماده در آباکوس (Material Property/Elastic) وارد کنم کافیه؟ یعنی آیا آباکوس تشخیص می دهد که این کرنش ها قسمت پلاستیک را هم شامل می شود؟
با تشکر