Moving temperature boundary condition subroutine

Dear All,

I am trying to model an hydrogen embrittlement problem with Cohesive Elements in Abaqus.

I am using a standard Abaqus Cohesive Elements with a specific TSL.
The specified TSL is also linked to temperature in the model (which is the hydrogen concentration).

The process should be as per follow:
• Abaqus will need to perform the diffusion analysis with cohesive elements.
• Once the cohesive element is degraded (SDEG=1), Abaqus should move the temperature boundary condition to the next element (to simulate the hydrogen penetration)


I have read the following thread ( and tried with DISP and UTEMP subroutine without success.

The questions are:
• The mentioned thread proposed use of MPC as a suitable solution to solve the problem. Is there any further suggestion on how to implement MPC Subroutine to solve the above problem?
• Any other suggestion on how to tackle the mentioned issue?

Thanks again
FEA_1981