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

2017-01-25 Thread Wolfgang Bangerth
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: question about a lage deformation model based on the updated Lagrangian foemolation

2017-01-25 Thread Bruno Turcksin
H.B., On Wednesday, January 25, 2017 at 3:02:01 PM UTC-5, ha.badn...@gmail.com wrote: > > > The image of the mapping applied to cell with center [0.295666 0.850622] > is distorted. The cell geometry or the mapping are invalid, giving a > non-positive volume fraction of -0.000119973 in quadratur

[deal.II] Re: question about a lage deformation model based on the updated Lagrangian foemolation

2017-01-25 Thread Jean-Paul Pelteret
Dear H.B. For this mesh and deformation field, you can expect to run into this problem regardless of whether you use the total or updated Lagrangian formulation. As you can see, with your updated formulation you run into an issue when the mapping Jacobian becomes singular (due to element skewn

[deal.II] Re: Assembly of material forces

2017-01-25 Thread 'Seyed Ali Mohseni' via deal.II User Group
@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

2017-01-25 Thread 'Seyed Ali Mohseni' via deal.II User Group
@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), .

Re: [deal.II] Re: Add heterogeneous parameters controlled by function on cells

2017-01-25 Thread Wolfgang Bangerth
On 01/25/2017 05:37 AM, Jean-Paul Pelteret wrote: Dear Lam, There are a number of ways that you can do this. One approach would be to define a class that derives from the function class ; this is alluded to in step-8

[deal.II] Re: Assembly of material forces

2017-01-25 Thread Jean-Paul Pelteret
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

2017-01-25 Thread Daniel Arndt
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: Add heterogeneous parameters controlled by function on cells

2017-01-25 Thread Jean-Paul Pelteret
Dear Lam, There are a number of ways that you can do this. One approach would be to define a class that derives from the function class ; this is alluded to in step-8

[deal.II] Re: Assembly of material forces

2017-01-25 Thread 'Seyed Ali Mohseni' via deal.II User Group
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] Add heterogeneous parameters controlled by function on cells

2017-01-25 Thread Lam DANG
Dear All I want to add some heterogeneous parameters (for example heterogeneous permeability which is varied on each cell by a function, etc...) on all cells. Could you let me know how to deal it? PS: As I know in the step 28, cell->material_id() only can help to get the parameters which were

[deal.II] Re: Assembly of material forces

2017-01-25 Thread Daniel Arndt
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

2017-01-25 Thread 'Seyed Ali Mohseni' via deal.II User Group
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

2017-01-25 Thread 'Seyed Ali Mohseni' via deal.II User Group
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