[deal.II] Re: Assembly of material forces

The last piece of the puzzle is solved. The small problem I mentioned was due to the problem discussed here: https://groups.google.com/d/msg/dealii/5YXHbjPs3ls/bsGMAPwXBQAJ Thank you everyone. The material force computation is now working 100% correctly for any case. Kind regards, S. A. Mohsen

[deal.II] Re: Assembly of material forces

@Prof. Bangerth: The upper one where pseudocolor is written, is from VisIt and deal.II. Below as a reference I showed my results from the C++ code. Both simulations are done using Q4 quadrilateral elements. There are no triangular elements here and the meshes are identical. They are intentional

[deal.II] Re: Assembly of material forces

@Prof. Bangerth: The upper one where pseudocolor is written, is from VisIt and deal.II. Below as a reference I showed my results from the C++ code. Both simulations are done using Q4 quadrilateral elements. There are no triangular elements here and the meshes are identical. They are intentional

Re: [deal.II] Re: Assembly of material forces

On 01/28/2017 04:56 AM, 'Seyed Ali Mohseni' via deal.II User Group wrote: I somehow figured the assembly out by listening to Prof. Bangerth's advice that I should try to solve it the deal.II way. And after Jean-Paul told me to take a look at step-42 I realized how things are working in deal.II.

[deal.II] Re: Assembly of material forces

First of all, thank you for all your suggestions. They helped me a lot, especially from Prof. Bangerth and Jean-Paul. Many thanks ;) I solved the assembly by listening to Prof. Bangerth's advice that I should try to implement things in a deal.II way. And after Jean-Paul told me to take a look a

[deal.II] Re: Assembly of material forces

First of all, thank you for all your suggestions. They helped me a lot, especially from Prof. Bangerth and Jean-Paul. Many thanks ;) I somehow figured the assembly out by listening to Prof. Bangerth's advice that I should try to solve it the deal.II way. And after Jean-Paul told me to take a lo

[deal.II] Re: Assembly of material forces

Dear Seyed, On top of what Wolfgang has already said, I have a few more direct comments about what you posted a few days ago. > The computation of configurational forces is accomplished by the following > formula: > > > >

Re: [deal.II] Re: Assembly of material forces

Seyed, Is there a command that allows me to assign values to the DoF. Can you please explain this a bit more. The Gauss quadrature points never include the vertices [1] so saying that the values represent forces for each vertex and at the same time that the values are computed with

[deal.II] Re: Assembly of material forces

@Daniel: So you have a FESystem<2>(FE_Q<2>(1),2) element, i.e. linear continuous > elements with two components, and the vector represents corresponding local > DoF values? > Exactly. I defined first FESystem fe; and then template SolidMechanics::SolidMechanics(): ..., fe(FE_Q(1), dim), .

[deal.II] Re: Assembly of material forces

@Daniel: So you have a FESystem<2>(FE_Q<2>(1),2) element, i.e. linear continuous > elements with two components, and the vector represents corresponding local > DoF values? > Exactly. I defined first FESystem fe; and then template SolidMechanics::SolidMechanics(): ..., fe(FE_Q(1), dim), .

[deal.II] Re: Assembly of material forces

Dear Seyed, > My problem is not the theory. > Providing the theory is for our benefit (and therefore yours). No one here is necessarily an expert in what you're trying to accomplish. So, for all of your explanation, an equation or two might go a long way to help us understand exactly what y

[deal.II] Re: Assembly of material forces

Seyed, There is a vector called cell_cf for instance. It contains all integrated > material forces (summed up material forces from integration points) for > each cell. Each cell has 4 vertices and 8 DoFs in 2D. Hence, it has to be a > 8x1 vector. > [...] > Now imagine you have this vector which

[deal.II] Re: Assembly of material forces

My problem is not the theory. I am just not a deal.II expert and therefore it is still quite difficult for me to implement my approach. Going in to theoretical detail would take maybe much longer. Maybe it is better instead, if I reformulate my question like this: There is a vector called cell_

[deal.II] Re: Assembly of material forces

Seyed, What exactly are you trying to do? Which equations are you considering? How would what you want to do look in formulas? The information you provide are much too high-level to expect a useful answer. [...] > Conclusion: It would be great, if you can at least guide me to the correct > dir

[deal.II] Re: Assembly of material forces

Dear Prof. Bangerth, It fails, because I tried to use the same procedure for the standard assembly such as in step-18 or step-40 and I received some awkward contour plots which won't match with what I showed from the C++ code. Hence, fails is refered to my own failure in understanding the core

[deal.II] Re: Assembly of material forces

Dear Prof. Bangerth, It fails, because I tried to use the same procedure for the standard assembly such as in step-18 or step-40 and I received some awkward contour plots which won't match with what I showed from the C++ code. Hence, fails is refered to my own failure in understanding the core