Essentially, there is no need in changing anything in assemble_system function from step-8, at least for stiffness matrix. You should change boundary conditions and triangulation to suit for your task, that is all. However you should not mix up Young's modulus with Lame parameters which are used in the program. if lambda and mu stay there, they have to stay there :). There is a great table with identities between different modulus and parameters of elasticity at https://en.wikipedia.org/wiki/Lam%C3%A9_parameters.
-- 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. For more options, visit https://groups.google.com/d/optout.