2D cantilever beam: reference to step-7,8,11,22

[priyank...@iitgn.ac.in]

I am solving for a 2D cantiliever beam  with a rectangular grid and referring to step-8. The problem is because boundary id is set to 0 by default, though I am appliying tensile force, could not obsrve any elongation towards right end. This is seen in code vector_pbm1.cc

So tried to change boundary id on which I want to apply dirichlet boundary condition(zerofunction). I want dirichlet to be applied on to left side of the rectangular grid and with force on to right side, beam should elongate. How do I achieve this??

Any suggestion is highly appreciated.
Thanks in advance

[Wolfgang Bangerth]

On 10/05/2017 12:19 AM, priyank...@iitgn.ac.in wrote:
> I am solving for a 2D cantiliever beam with a rectangular grid and referring
> to step-8. The problem is because boundary id is set to 0 by default, though I
> am appliying tensile force, could not obsrve any elongation towards right end.
> This is seen in code vector_pbm1.cc

There is really not enough information here for us to help you. If I run the
program, I do see an elongation that results from the solution vector.

> So tried to change boundary id on which I want to apply dirichlet boundary
> condition(zerofunction). I want dirichlet to be applied on to left side of the
> rectangular grid and with force on to right side, beam should elongate. How do
> I achieve this??

What have you tried so far? What do you observe? How have you simplified the
problem to isolate where the problem may be? Have you applied a case where you
prescribe Dirichlet boundary conditions all around, and does that work? If it
does, how do you apply the tensile force? Via non-zero Dirichlet boundary
conditions, or via a contribution to the right hand side of the weak form?

Best
W.

--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

[Jie Cheng]

Hi there

Like Wolfgang said, you need to elaborate your question. I ran your code and I did see positive x-displacement in a certain area. I guess your question was that you expected the grid to be actually deformed? If so, that is because your results were plotted on the original mesh. Look at the move_mesh function in step-18 or the output_results function in step-44.

Best

Jie

[Wolfgang Bangerth]

On 10/07/2017 12:12 PM, Jie Cheng wrote:
>
> Like Wolfgang said, you need to elaborate your question. I ran your code and I
> did see positive x-displacement in a certain area. I guess your question was
> that you expected the grid to be actually deformed? If so, that is because
> your results were plotted on the original mesh. Look at the move_mesh function
> in step-18 or the output_results function in step-44.

Or you can do this in Paraview using the "warp by vector" plot.

[Priyanka Nemani(15310051)]

Thanks all for the suggestions. I am working on this.

Best,

Priyanka Nemani

Research Scholar

Discipline of Chemical Engineering

IIT Gandhinagar | +91 9913578357

--
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 a topic in the Google Groups "deal.II User Group" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/dealii/AnhMEjEYMcQ/unsubscribe.
To unsubscribe from this group and all its topics, send an email to dealii+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.