Abstract
We propose an effective and flexible way to assemble tangent stiffness matrices in MATLAB. Our technique is applied to elastoplastic problems formulated in terms of displacements and discretized by the finite element method. The tangent stiffness matrix is repeatedly assembled in each time step and in each iteration of the semismooth Newton method. We consider von Mises and Drucker-Prager yield criteria, linear and quadratic finite elements in two and three space dimensions. Our codes are vectorized and available for download. Comparisons with other available MATLAB codes show, that our technique is also efficient for purely elastic problems. In elastoplasticity, the assembly times are linearly proportional to the number of integration points.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
