You are here
Possibility of modifying the tangential stifness matrix using external subroutines in abaqus
Mon, 2009-06-15 12:54 - unsrohith
Hi,
I am working on implementation of non local damage mechanics model in abaqus using external subroutines, But in that case I have to modify the stifness matrix by adding the newly formed nonlocal stifness, Out of the available external subroutine options in abaqus documentation if any one gives the option to modify the tangential stifness matrix Please help me in the same.
Thanks and Regards,
Rohith
»
- unsrohith's blog
- Log in or register to post comments
- 8250 reads
![Subscribe to Comments for "Possibility of modifying the tangential stifness matrix using external subroutines in abaqus"](https://imechanica.org/misc/feed.png)
Comments
stifness matrix modification in ABAQUS through user subroutines
Hi,
I am working on implementation of non local damage mechanics model in
abaqus using external subroutines, But in that case I have to modify
the stifness matrix by adding the newly formed nonlocal stifness, Out
of the available external subroutine options in abaqus documentation if
any one gives the option to modify the tangential stifness matrix.
Rohith
UEL subroutine
Greetings,
For your situation, UEL subroutine can be used to define a new element for non-local damage model. In SR/ UEL you should defined the RHS(residual force vector) and AMATRX(stiffness matrix) for your constitutive model.
Best regards,
Sang Ha
Sang Ha
Hi thank you even I thought the same after going through the documentation..
Rohith
nonlocal damage (integral formulation)
Hi, what you are doing is also interesting for me, since I am trying to implement the nonlocal damage based on integral formulations in ABAQUS (UEL). Definitely, UEL gives us the possibility to modify the element stiffness matrix for the newly fomed part of it as you mentioned. However, what I am struggling about is how to stop the simulation after it once went through all the Gauss points in the same time incerement to build the nonlocal part of damage based on averaging over the neighboring Gauss points!
do you have any hint in this regard?
Best,
Hadi
Sanga Ha
Can u please upload any example code for it