Dear Deal.II group
I am trying to add boundary force (b) to step-18 as below
[image: 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.)
instead of procedure described in step-18, can I use functions "
create_right_hand_side" and "create_boundary_right_hand_side" of
VectorTools for applying body and traction force, respectively or they just
belongs to Laplace equation?
Thanks for considering this!
Amir
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c5c74fbf-7bf5-47e9-906c-6e7db0d23a76n%40googlegroups.com.
