Dear deal.II experts, I have a question regarding the assembly of material forces also known as configurational forces.
First of all, I have constructed my own solid mechanics code within deal.II which I already validated with our own C++ code. Hence, the deal.II code works until here. <https://lh3.googleusercontent.com/-b6_TF8L_Fo0/WIdw2UhIGgI/AAAAAAAAADc/w5h60oEjoxokH2YN5kP2tMLyPeaiIjSAQCLcB/s1600/Bild1.png> <https://lh3.googleusercontent.com/-j5HV0dZ71M8/WIdxR5AXVQI/AAAAAAAAADg/slPwCpDBMbc0cECXZaoIAw7e3PZud_18ACLcB/s1600/Bild2.png> <https://lh3.googleusercontent.com/-gULxuXw7QTk/WIdxluJIZQI/AAAAAAAAADk/cwhKZpBz53QuTzn0iF3XNzyh1TZPiEYCgCLcB/s1600/Bild3.png> Attached you find my results on material forces for one single element. It also works perfectly. <https://lh3.googleusercontent.com/-X5WzN-koW9c/WIdx6PKxJgI/AAAAAAAAADo/k5KEGY9Dycc18RbmbZM9wmLcdQgnDZ7YACLcB/s1600/Bild4.png> <https://lh3.googleusercontent.com/-vsT6hLQk1BI/WIdx_obfu9I/AAAAAAAAADs/J_eThAgNPKweXvMiu2Kgm1ct-uG4mgTZwCLcB/s1600/Bild5.png> Now I try to compute the material forces for the whole structure. Hence, I tried using the same procedure of assembly for our LHS or RHS matrices in deal.II, but it fails. Do I have to take a special approach in order to make it work? I assume my local material forces have to be assembled to the global nodes, but this assembly is it just a simple mapping to the nodes? Since there are quadtree constraints, it will not be that simple, I think. <https://lh3.googleusercontent.com/-XONr8pFcP-M/WIdyv4aRCEI/AAAAAAAAADw/0nEOzI78WEQ4_5IAgEY4e58I6ekyhnJlQCLcB/s1600/Bild6.png> <https://lh3.googleusercontent.com/-FR50q3yJRgU/WIdyym0YwyI/AAAAAAAAAD0/S648FlQjWi07u9KSfUags45PMKMW1FZkgCLcB/s1600/Bild7.png> The successful results are computed within our simple c++ code and I just show them for comparison reasons. Extra question: Would my work be of interest for your code gallery maybe in the near future? Kind regards, S. A. Mohseni -- 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.
