Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Jacobian matrix for External material model

Please login with a confirmed email address before reporting spam

Dear all:

I developed an external material model in which the stress was updated incrementally, which is differnt from the stress update algorithm in the BLOG which is defined by total value. Now I have a problem about the calculation of Jacobian matrix. According to the HELP document, the Jacobian is defined as: , where ds and de are the derivatives with respect to total value. However, for the incremental stress update algorithm, this Jacobian definition seems difficult to be determined. It may be calculated as: d(s_pre+delta_s)=Jacd(e_pre+delta_e), where _pre means the value at previous converged step and delta means the incremental values. As the s_pre and e_pre are constants at current step, this defination may be simplified as: d(delta_s)=Jacd(delta_e) , which is similar to the first one, i.e. “delta_s = Jac * delta_e”. I have tried this defination, but it does not work well. So I think it may be better to use the Jacobian defined with respect to total value.

Could you please share me some ideas how to determine the Jacobian matrix?

Thanks a lot.

-------------------
Yu Zhang

0 Replies Last Post Jan 29, 2020, 10:39 a.m. EST
COMSOL Moderator

Hello Yu Zhang

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.