On 7/11/23 05:32, Mohammad Amir Kiani Fordoei wrote:
Screenshot from 2023-07-11 14-01-30.png
As I understood, I just need to specify a boundary_id on triangulation and
apply traction force on that with code like this:
{ - loop over cells
-loop over faces per cell && controlling being at_boundary
-loop over quadrature points of face
-loop over dof per cell && controlling intended boundary_id_NBC
-finally adding this term to rhs as
cell_rhs (i)+=
(fe_face_values.shape_value (i,q_point)
*traction_component
*fe_face_values.JxW(q_point));}
But, I didn't see noticble variation of displacement or norm of stress in my
solution, even for large traction force.
is it necessary to multiply traction_component to "* present_timestep *
velocity " (like Dirichlet BC displacement)? (I tried this too, but similar to
previous one no significant effect of load was observed.)
Is this a boundary where you are still *also* applying a prescribed
displacement? You can only do one or the other.
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/486aa979-8841-0075-e6ab-dd9bdf4e0106%40colostate.edu.
