Thanks a lot that helped, I applied symmetric<2,dim> tensor because it is a 
heat problem. I also implemented the convection boundary condition. I have 
a problem in the results, it seems the contact between surfaces is not 
detected automatically, I want to model such contact.Please I will be 
grateful for your guidance.[image: contact problem.png]
On Monday, August 16, 2021 at 9:43:25 PM UTC+2 Wolfgang Bangerth wrote:

> On 8/16/21 10:37 AM, ibrahim zawra wrote:
> > 
> >  real problem for graphical board . I used GMSH and identified different 
> > material_ids on different parts and made sure it is correct using the 
> function 
> > print mesh info.
> > I wrote the class Material data to substitute the coefficient function. 
> but I 
> > tried several things but it didn't work.
>
> You don't say what specifically you have tried and how or why it didn't 
> work. 
> But the usual way to do this would be in the assemble_system() function, 
> where 
> you would do
>
> for (auto &cell : dof_handler.active_cell_iterators())
> {
> fe_values.reinit (cell);
>
> SymmetricTensor<4,dim> strain_stress_tensor;
> if (cell->material_id() == ...)
> strain_stress_tensor = ...whatever is appropriate for this material...;
> else
> strain_stress_tensor = ...whatever is appropriate for other materials..;
>
> for (unsigned int q=...)
> for (unsigned int i=...)
> for (unsigned int j=...)
> cell_matrix(i,j) += ...something that uses the strain stress tensor;
>
> ...
> }
>
> Best
> W.
>
> -- 
> ------------------------------------------------------------------------
> 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/c3c85b48-2c12-41c8-bb9f-86a2eccacd38n%40googlegroups.com.

Reply via email to