Thank you for your comments.
Ok,
Here is my problem:
I am doing material optimization, combining Abaqus with Matlab.
I use Matlab intrinsic function to update my variable(E modulus).
Matlab uses Quasi-Newton method to update variable, for that reason it needs to calculate gradient at each...