On 11/20/22 07:20, Deepika Kushwah wrote:
I'm trying to solve elasticity problem refereeing step 8 and have got solution
for plane strain problem. Now I wish to solve the problem with plane stress
condition for which matrix D (constitutive matrix) is to be modified to get
stiffness (K = B^T*D*B).
Not clear where it has been defined in the code and how it can be modified.
Deepika:
assemble_system() is always the function where the system matrix is built. We
do not generally assemble systems via B^T*D*B, however, but rather based on
the weak formulation of the equation -- that is, the bilinear form that
corresponds to the equation you want to solve.
A good first step would be to understand what equation step-8 solves, and how
that corresponds to what that program's assemble_system() does. Then you'll
have to identify what the bilinear form is for the equation you care about.
The final step is to translate that in the same way as what the program
currently does.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c3a0d03d-93aa-d03c-8d88-87c0be10be1a%40colostate.edu.
