Re: [deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-06-24 Thread Ramprasad R
Hello Bangerth,

The terms N_* are the resultant forces in the * direction and these forces
are calculated using the strains which in turn are calculated using the
displacements u. So the terms N_* are not constants, rather change with
each element. And these values directly depend on u.

Thank you.

Regards,
Ramprasad

On Mon, Jun 24, 2019 at 10:56 PM Wolfgang Bangerth 
wrote:

> On 6/24/19 9:24 AM, Ramprasad R wrote:
> >
> > The following image has the exact formulation of the geometric stiffness
> > matrix. Which is denoted by G in the image.
>
> Ramprasad,
> I don't think I quite understand the connection step-60. In the context
> of that problem, for which of the two variables (u and lambda) do you
> want to compute that matrix? Assuming you mean that it is supposed to be
> for u, it is not difficult to write down what you need to do based on
> what you already learn in step-3 or step-4, for example. In fact, these
> two early programs already implement the first two terms of G_ij if you
> assume that N_x=N_y=N_xy=1.
>
> I am not familiar with the terms used in the page you show. Are the N_*
> factors constants? Are they functions of space?
>
> 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/53a1ebd3-fdd4-485d-2486-b70913095c55%40colostate.edu
> .
> For more options, visit https://groups.google.com/d/optout.
>

-- 
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/CAKS2Wq58q8AU%2BJ%3DKR2Y7aExCd%2BFa1onWQMHR6EJv6-ADST8UKQ%40mail.gmail.com.
For more options, visit https://groups.google.com/d/optout.


Re: [deal.II] Fwd: Regarding Step 60 tutorial of deal.ii

2019-06-24 Thread Wolfgang Bangerth
On 6/24/19 9:24 AM, Ramprasad R wrote:
> 
> The following image has the exact formulation of the geometric stiffness 
> matrix. Which is denoted by G in the image.

Ramprasad,
I don't think I quite understand the connection step-60. In the context 
of that problem, for which of the two variables (u and lambda) do you 
want to compute that matrix? Assuming you mean that it is supposed to be 
for u, it is not difficult to write down what you need to do based on 
what you already learn in step-3 or step-4, for example. In fact, these 
two early programs already implement the first two terms of G_ij if you 
assume that N_x=N_y=N_xy=1.

I am not familiar with the terms used in the page you show. Are the N_* 
factors constants? Are they functions of space?

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/53a1ebd3-fdd4-485d-2486-b70913095c55%40colostate.edu.
For more options, visit https://groups.google.com/d/optout.


Re: [deal.II] running the "The 'Goal-oriented mesh adaptivity in elastoplasticity problems' code gallery program "

2019-06-24 Thread Muhammad Mashhood
Thank you Prof Wolfgang for your possible support and concern. I will be 
waiting for the further support from the author to successfully run it and 
to understand the workflow correctly. 

On Saturday, June 22, 2019 at 1:51:10 AM UTC+2, Wolfgang Bangerth wrote:
>
>
> Muhammad, 
> let me bring the original author of the program (Seyed Shahram Ghorashi) 
> into 
> the loop as well. 
>
> > I ran the case of "Cantiliver_beam_3d" and it returned the results after 
> > couple of minutes as mentioned in the attachement (where per time step, 
> > seemingly only one iteration was performed in Newton method and in all 
> > increamental load steps the result values are zero for both stress and 
> > displacement). 
>
> Good question. I can confirm that that is what the program generates 
> (though 
> at least the dual solution is non-zero), but as with all of the code 
> gallery 
> programs, we just make them available for others to look at -- we don't 
> generally know what specifically they are doing, or whether they are 
> correct. 
> You will have to ask Seyed Shahram Ghorashi this question. 
>
>
> > Secondly I tried to run the case "Thick_tube_internal_pressure" which 
> first 
> > popped up the dialog message of : 
>  > [...] 
> > Where as per my understanding of template based nature of deall.ii and 
> also 
> > this case is 2d so if I change the *dim = 2* from originally *dim = 3* 
> in main 
> > function of elastoplastic.cc 
>
> Yes, this is the correct solution. 
>
>
> > then the following error comes during compilation : 
> > [...] 
>
> Yes, I can confirm this as well. I will admit that I don't think this ever 
> worked. 
>
> In the end, this is really the limitation of the code gallery. We see it 
> as a 
> way for our users to share codes that may be of use to others in the 
> community, but we do not endorse them or verify the quality of these 
> codes. 
> It's also quite possible that they *used* to work correctly with older 
> deal.II 
> versions, but don't any more -- I could imagine that being the situation 
> with 
> the first problem above. 
>
> Let's hope that Seyed has something to add that help you! 
>
> Best 
>   WB 
>
> -- 
>  
> Wolfgang Bangerth  email: bang...@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/9cac3aa3-fa19-4dcd-acbe-29b251186028%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.