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.
Accessing Stiffness Matrices and Internal Force Vectors for Incremental Solution in COMSOL
Posted 2 sept. 2024, 05:47 UTC−4 API, Studies & Solvers, Structural Mechanics Version 6.2 5 Replies
Please login with a confirmed email address before reporting spam
Hello COMSOL Community,
I'm currently working on a phase-field fracture problem and trying to implement a custom solver that follows a Newton-Raphson iteration scheme with arc-length control, similar to what is shown in the attached image (equation 17a and 17b).
I need to access the stiffness matrices (K_uu, K_φφ, etc.) and the internal force vectors (f_int,u, f_int,φ) in COMSOL in order to calculate the solution increments (Δu, Δφ, Δλ) at each iteration. I am aware that quantities like displacement (solid.uf) and damage (solid.dmg) can be accessed, but I'm unsure how to retrieve the stiffness matrices and internal force vectors.
I would like to use them in a custom script, potentially via the mph package or COMSOL's Java API.
Thank you!
Attachments: