COMSOL with stiffness or constitutive matrix extraction

Hello, all:

    Can anyone please help me how to deal with the following problem in COMSOL? Thanks a lot.

    How to extract the global stiffness matrix [K] or material constitutive matrix [D_ep] for each load increment after convergence using MATLAB?

   Thanks again.

Regards,

Xu