Re: [deal.II] Elasticity probelm: same matrices - different solution vector.
On 10/24/22 13:16, Мария Бронзова wrote: *** Caution: EXTERNAL Sender *** Dear all, I have been trying to implement orthotropic elasticity (similar as in step-18) and in order to verify my results I decided to set all E and G modulus values as in an isotropic case and to compare the outcome with the one for the same isotropic elastic problem but with another formulation (as in step-46 for the structure part). I created the simplest 1 cell models and compared the matrices in both cases (both LHS and RHS), they are absolutely identical. But weirdly enough I receive different solution vectors. Maria: That is of course a contradiction. If you have the same linear system, and you solve it the same way, you should get the same solution. We are generally happy to help with advice on this forum, but we can't debug other people's code. But I think you have everything you need to debug your codes yourself. Just output the elements of the matrix and rhs of the linear system for both codes. Make sure they are the exact same. If they are not, then you know that the problem lies with the assembly. If they are, then solve the linear system with the exact same solver and output the solution element by element. It should be the exact same. If it is not, you know that the problem lies with what must be different ways you are solving the linear system. If it is, then the problem lies in the part of the code that does something with the solution later on. 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/dad3cb12-9e95-33a3-1b1c-782ad6d68c29%40colostate.edu.
[deal.II] Elasticity probelm: same matrices - different solution vector.
Dear all, I have been trying to implement orthotropic elasticity (similar as in step-18) and in order to verify my results I decided to set all E and G modulus values as in an isotropic case and to compare the outcome with the one for the same isotropic elastic problem but with another formulation (as in step-46 for the structure part). I created the simplest 1 cell models and compared the matrices in both cases (both LHS and RHS), they are absolutely identical. But weirdly enough I receive different solution vectors. I have been struggling trying to identify the source of the problem. If you could please take a look at my small scripts - I would be very grateful..(you can find them attached). Thank you a lot in advance! Kind regards, Mariia -- 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/b3964286-359e-453b-a0ad-0b98aee9dc81n%40googlegroups.com. orth_1_cell_no_freq Description: Binary data elast_1_cell_no_freq Description: Binary data
[deal.II] Elasticity Preconditioner
Dear all, A student and me are trying to deal with (parallel, static) linearized elasticity similar to step-8. However, the problem is slightly different since the Lame parameters are functions instead of constants. Can anyone recommend good solvers & preconditioners for the resulting system? Our boundary conditions fix all degrees of freedom in the kernel. Best, Konrad -- 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/eba0cc0f-6028-4fde-9cee-307e09f8d2ban%40googlegroups.com.
[deal.II] Elasticity
Dear all, In step 18, the elasticity problem has been solved with this assumption that there was not initial strain in the equation. I want to add the initial strain with consideration of the fact that this tensor is a function of the solution. In other words, this tensor should be multiplied by a function that is related to the solution vectors. It would be very kind of you if you help what is the best way to introduce this additional term in the equation. I do really appreciate your kindness and your time you can answer to my question. Looking forward to hearing from you. Thanks, Benhour -- 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.
