Re: [deal.II] Step-4 1D problem

[Praveen C]

In 1-D,

GridGenerator::hyper_cube(triangulation, -1, 1);

will assign different boundary indicators to left and right side, even though 
colorize=false is default, see

https://dealii.org/developer/doxygen/deal.II/namespaceGridGenerator.html#acea0cbcd68e52ce8113d1134b87de403
 


 I think left=0 and right=1. Due to this reason, bc may not be applied on both 
sides.

Do the apply bc twice, once with indicator=0 and once with indicator=1.

  std::map boundary_values;
  VectorTools::interpolate_boundary_values(dof_handler,
   0, // left boundary
   BoundaryValues(),
   boundary_values);
  MatrixTools::apply_boundary_values(boundary_values,
 system_matrix,
 solution,
 system_rhs);

  VectorTools::interpolate_boundary_values(dof_handler,
   1, // right boundary
   BoundaryValues(),
   boundary_values);
  MatrixTools::apply_boundary_values(boundary_values,
 system_matrix,
 solution,
 system_rhs);

Best
praveen

> On 07-Jun-2020, at 9:50 AM, 'Christoph Kammer' via deal.II User Group 
>  wrote:
> 
> Hi Deal.ii team,
> 
> I went through tutorial step-4 and I tried to use the solver for Poisson's 
> problem in 2D and rewrite it for a 1D problem. This is basically step-4 and 
> setting RHS f =1 and solving for homogeneous BCs. The program is compiling 
> just fine and executing, however, the solution I am getting from that doesn't 
> make any sense. The BCs aren't even satisfied at x = 1.
> 
> Could anyone let me know what's going on here? Setting f = 1 and changing 
> everything to solve for homogeneous BC in 1D seems like a very simple problem.
> 
> Thank you,
> 
> Christoph

-- 
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/C66D5394-FFCF-426E-A66F-9F3CA5A0F7B0%40gmail.com.

[deal.II] Step-4 1D problem

['Christoph Kammer' via deal.II User Group]

Hi Deal.ii team,

I went through tutorial step-4 and I tried to use the solver for Poisson's 
problem in 2D and rewrite it for a 1D problem. This is basically step-4 and 
setting RHS f =1 and solving for homogeneous BCs. The program is compiling 
just fine and executing, however, the solution I am getting from that 
doesn't make any sense. The BCs aren't even satisfied at x = 1.

Could anyone let me know what's going on here? Setting f = 1 and changing 
everything to solve for homogeneous BC in 1D seems like a very simple 
problem.

Thank you,

Christoph

-- 
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/5242ff6f-2add-4557-83b6-2905f9697dc4o%40googlegroups.com.

[deal.II] About metallic ductile damage

[itachi Ezio]

Hi,

I'm Abaqus user and fresh new about using deal.II. I'd like to ask a few 
questions base on deal.II functions.

I'm looking for the simulation for metallic ductile damage behavior with 
deal.II ( for example with triaxiality and lode angle), is it possible?

And how is it like to write user materials in deal.II, especially in 
plasticity behaviors? Is it also quit sophisticated as in UMAT in Abaqus?

Looking forward to the answers, thanks for that!  


-- 
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/8df55d20-2405-4e36-9be8-99d2b10d38a7o%40googlegroups.com.

[deal.II] Re: triangulation save not working for 1D domain

[Amaresh B]

Thank you peterrum for your reply.

I have tried the following commands. In both the cases, the triangulation 
was not saved. Please let me know if you have any suggestions.

1.GridGenerator::hyper_cube (triangulation, 0., 40.0);


2.   GridGenerator::hyper_rectangle(triangulation,
 Point(0.0),
 Point(40.0),
 true);


Thanks,
Amaresh


On Saturday, June 6, 2020 at 9:33:20 PM UTC+5:30, peterrum wrote:
>
> What type of triangulation are you using?
>
> Peter
>
> On Saturday, 6 June 2020 17:52:53 UTC+2, Amaresh B wrote:
>>
>>   Dear all,
>>
>> I am trying to save my triangulation and using the lines below. 
>> 'triangulation.save' commands works and saves the mesh if my domain is 
>> either 2D or 3D (i.e. dim=2 or 3). However, it does not save the 
>> triangulation if dim=1, nor it gives any error message. I will be thankful 
>> if someone can comment on it and give suggestions.
>>
>>
>>
>>std::vector 
>> sol_transfer_vectors;
>>sol_transfer_vectors.push_back();
>> 
>> 
>> parallel::distributed::SolutionTransfer 
>> sol_transfer(dof_handler);
>>  sol_transfer.prepare_for_serialization 
>> (sol_transfer_vectors);
>>
>>   triangulation.save("restart.mesh");
>>
>> Thank you,
>>
>> Regards
>> Amaresh
>>  
>>
>

-- 
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/d4f7eb57-9a58-489e-879a-d3a56906e568o%40googlegroups.com.

[deal.II] Re: triangulation save not working for 1D domain

['peterrum' via deal.II User Group]

What type of triangulation are you using?

Peter

On Saturday, 6 June 2020 17:52:53 UTC+2, Amaresh B wrote:
>
>   Dear all,
>
> I am trying to save my triangulation and using the lines below. 
> 'triangulation.save' commands works and saves the mesh if my domain is 
> either 2D or 3D (i.e. dim=2 or 3). However, it does not save the 
> triangulation if dim=1, nor it gives any error message. I will be thankful 
> if someone can comment on it and give suggestions.
>
>
>
>std::vector 
> sol_transfer_vectors;
>sol_transfer_vectors.push_back();
> 
> 
> parallel::distributed::SolutionTransfer 
> sol_transfer(dof_handler);
>  sol_transfer.prepare_for_serialization 
> (sol_transfer_vectors);
>
>   triangulation.save("restart.mesh");
>
> Thank you,
>
> Regards
> Amaresh
>  
>

-- 
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/00da9e61-551a-4a19-851a-3179fc4df131o%40googlegroups.com.

[deal.II] triangulation save not working for 1D domain

[Amaresh B]

  Dear all,

I am trying to save my triangulation and using the lines below. 
'triangulation.save' commands works and saves the mesh if my domain is 
either 2D or 3D (i.e. dim=2 or 3). However, it does not save the 
triangulation if dim=1, nor it gives any error message. I will be thankful 
if someone can comment on it and give suggestions.



   std::vector sol_transfer_vectors;
   sol_transfer_vectors.push_back();


parallel::distributed::SolutionTransfer 
sol_transfer(dof_handler);
 sol_transfer.prepare_for_serialization 
(sol_transfer_vectors);

  triangulation.save("restart.mesh");

Thank you,

Regards
Amaresh
 

-- 
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/1a081d39-c0b5-47c4-b58a-e227e62840c4o%40googlegroups.com.