On 7/15/23 11:33, Mohammad Amir Kiani Fordoei wrote:
const Tensor<1,dim> n = fe_face_values.normal_vector(q_point);
cell_rhs (i)+=( (n*traction_component*fe_face_values.shape_value(i, q_point))*
fe_face_values.JxW(q_point));
}//for q_point
}//for dof_per_cell
But in computation of cell rhs i ge
Dear Wolfgang
Sorry for late response. Due to internet censorship in Iran, I didn't have
access to google group for few days.
I consider a specific region on triangulation for Neumann BC (a force
vector in 3D space).
There must be some mistakes in my program, but I couldn't figure out.
I tried t
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_boundar
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
