example_bend_pipe Program

Uses

  • program~~example_bend_pipe~~UsesGraph program~example_bend_pipe example_bend_pipe module~forcad forcad program~example_bend_pipe->module~forcad module~forcad_kinds forcad_kinds module~forcad->module~forcad_kinds module~forcad_nurbs_curve forcad_nurbs_curve module~forcad->module~forcad_nurbs_curve module~forcad_nurbs_surface forcad_nurbs_surface module~forcad->module~forcad_nurbs_surface module~forcad_nurbs_volume forcad_nurbs_volume module~forcad->module~forcad_nurbs_volume module~forcad_nurbs_curve->module~forcad_kinds fordebug fordebug module~forcad_nurbs_curve->fordebug module~forcad_utils forcad_utils module~forcad_nurbs_curve->module~forcad_utils module~forcad_nurbs_surface->module~forcad_kinds module~forcad_nurbs_surface->fordebug module~forcad_nurbs_surface->module~forcad_utils module~forcad_nurbs_volume->module~forcad_kinds module~forcad_nurbs_volume->fordebug module~forcad_nurbs_volume->module~forcad_utils module~forcad_utils->module~forcad_kinds

Example program demonstrating how to bend a straight pipe-like NURBS volume into a circular arc.

The program: Creates a straight pipe segment (as a ring extruded in ), Refines the shape, Applies bending with different bend angles (90°, 270°, 360°), Exports the resulting NURBS volumes to VTK, Displays the geometry.


Calls

program~~example_bend_pipe~~CallsGraph program~example_bend_pipe example_bend_pipe proc~bend_pipe bend_pipe program~example_bend_pipe->proc~bend_pipe proc~create~2 nurbs_volume%create program~example_bend_pipe->proc~create~2 proc~elevate_degree~2 nurbs_volume%elevate_degree program~example_bend_pipe->proc~elevate_degree~2 proc~export_xc~2 nurbs_volume%export_Xc program~example_bend_pipe->proc~export_xc~2 proc~export_xg~2 nurbs_volume%export_Xg program~example_bend_pipe->proc~export_xg~2 proc~export_xth_in_xg~2 nurbs_volume%export_Xth_in_Xg program~example_bend_pipe->proc~export_xth_in_xg~2 proc~insert_knots~2 nurbs_volume%insert_knots program~example_bend_pipe->proc~insert_knots~2 proc~set_ring~2 nurbs_volume%set_ring program~example_bend_pipe->proc~set_ring~2 proc~show~2 nurbs_volume%show program~example_bend_pipe->proc~show~2 none~get_knot~2 nurbs_volume%get_knot proc~bend_pipe->none~get_knot~2 none~get_nc~2 nurbs_volume%get_nc proc~bend_pipe->none~get_nc~2 none~get_wc~2 nurbs_volume%get_Wc proc~bend_pipe->none~get_wc~2 none~get_xc~2 nurbs_volume%get_Xc proc~bend_pipe->none~get_xc~2 none~set~2 nurbs_volume%set proc~bend_pipe->none~set~2 proc~is_rational~2 nurbs_volume%is_rational proc~bend_pipe->proc~is_rational~2 interface~compute_xg~2 compute_Xg proc~create~2->interface~compute_xg~2 interface~ndgrid ndgrid proc~create~2->interface~ndgrid proc~create~2->proc~is_rational~2 set set proc~create~2->set proc~elevate_degree~2->none~get_knot~2 proc~elevate_degree~2->none~set~2 proc~elevate_degree_a_5_9 elevate_degree_A_5_9 proc~elevate_degree~2->proc~elevate_degree_a_5_9 proc~elevate_degree~2->proc~is_rational~2 proc~cmp_elem_xc_vis~2 nurbs_volume%cmp_elem_Xc_vis proc~export_xc~2->proc~cmp_elem_xc_vis~2 proc~export_vtk_legacy export_vtk_legacy proc~export_xc~2->proc~export_vtk_legacy proc~export_xc~2->set proc~cmp_elem_xg_vis~2 nurbs_volume%cmp_elem_Xg_vis proc~export_xg~2->proc~cmp_elem_xg_vis~2 proc~export_xg~2->proc~export_vtk_legacy proc~export_xg~2->set proc~export_xth_in_xg~2->interface~compute_xg~2 interface~unique unique proc~export_xth_in_xg~2->interface~unique proc~export_xth_in_xg~2->proc~export_vtk_legacy proc~export_xth_in_xg~2->proc~is_rational~2 proc~export_xth_in_xg~2->set interface~compute_multiplicity compute_multiplicity proc~insert_knots~2->interface~compute_multiplicity proc~insert_knots~2->none~get_knot~2 proc~findspan findspan proc~insert_knots~2->proc~findspan proc~insert_knot_a_5_1 insert_knot_A_5_1 proc~insert_knots~2->proc~insert_knot_a_5_1 proc~insert_knots~2->proc~is_rational~2 s_loc s_loc proc~insert_knots~2->s_loc proc~insert_knots~2->set proc~set_ring~2->none~set~2 proc~compute_multiplicity1 compute_multiplicity1 interface~compute_multiplicity->proc~compute_multiplicity1 proc~compute_multiplicity2 compute_multiplicity2 interface~compute_multiplicity->proc~compute_multiplicity2 proc~compute_xg_bspline_3d compute_Xg_bspline_3d interface~compute_xg~2->proc~compute_xg_bspline_3d proc~compute_xg_bspline_3d_1point compute_Xg_bspline_3d_1point interface~compute_xg~2->proc~compute_xg_bspline_3d_1point proc~compute_xg_nurbs_3d compute_Xg_nurbs_3d interface~compute_xg~2->proc~compute_xg_nurbs_3d proc~compute_xg_nurbs_3d_1point compute_Xg_nurbs_3d_1point interface~compute_xg~2->proc~compute_xg_nurbs_3d_1point proc~ndgrid2 ndgrid2 interface~ndgrid->proc~ndgrid2 proc~ndgrid3 ndgrid3 interface~ndgrid->proc~ndgrid3 proc~unique_integer unique_integer interface~unique->proc~unique_integer proc~unique_real unique_real interface~unique->proc~unique_real proc~get_knot_all~2 nurbs_volume%get_knot_all none~get_knot~2->proc~get_knot_all~2 proc~get_knoti~2 nurbs_volume%get_knoti none~get_knot~2->proc~get_knoti~2 proc~get_nc_all~2 nurbs_volume%get_nc_all none~get_nc~2->proc~get_nc_all~2 proc~get_nc_dir~2 nurbs_volume%get_nc_dir none~get_nc~2->proc~get_nc_dir~2 proc~get_wc_all~2 nurbs_volume%get_Wc_all none~get_wc~2->proc~get_wc_all~2 proc~get_wci~2 nurbs_volume%get_Wci none~get_wc~2->proc~get_wci~2 proc~get_xc_all~2 nurbs_volume%get_Xc_all none~get_xc~2->proc~get_xc_all~2 proc~get_xcid~2 nurbs_volume%get_Xcid none~get_xc~2->proc~get_xcid~2 proc~get_xci~2 nurbs_volume%get_Xci none~get_xc~2->proc~get_xci~2 proc~set1~2 nurbs_volume%set1 none~set~2->proc~set1~2 proc~set2~2 nurbs_volume%set2 none~set~2->proc~set2~2 proc~set3~2 nurbs_volume%set3 none~set~2->proc~set3~2 proc~set4~2 nurbs_volume%set4 none~set~2->proc~set4~2 interface~elemconn_c0 elemConn_C0 proc~cmp_elem_xc_vis~2->interface~elemconn_c0 proc~cmp_elem_xg_vis~2->interface~elemconn_c0 proc~elevate_degree_a_5_9->interface~compute_multiplicity proc~bincoeff bincoeff proc~elevate_degree_a_5_9->proc~bincoeff proc~cmp_elemconn_c0_l cmp_elemConn_C0_L interface~elemconn_c0->proc~cmp_elemconn_c0_l proc~cmp_elemconn_c0_s cmp_elemConn_C0_S interface~elemconn_c0->proc~cmp_elemconn_c0_s proc~cmp_elemconn_c0_v cmp_elemConn_C0_V interface~elemconn_c0->proc~cmp_elemconn_c0_v proc~factln factln proc~bincoeff->proc~factln interface~kron kron proc~compute_xg_bspline_3d->interface~kron proc~basis_bspline basis_bspline proc~compute_xg_bspline_3d->proc~basis_bspline proc~compute_xg_bspline_3d_1point->interface~kron proc~compute_xg_bspline_3d_1point->proc~basis_bspline proc~cmp_tgc_3d cmp_Tgc_3d proc~compute_xg_nurbs_3d->proc~cmp_tgc_3d proc~compute_xg_nurbs_3d_1point->interface~kron proc~compute_xg_nurbs_3d_1point->proc~basis_bspline proc~get_nc_dir~2->interface~compute_multiplicity proc~set1~2->set proc~cmp_degree~2 nurbs_volume%cmp_degree proc~set1~2->proc~cmp_degree~2 proc~cmp_nc~2 nurbs_volume%cmp_nc proc~set1~2->proc~cmp_nc~2 proc~set2~2->set proc~set2~2->proc~cmp_nc~2 proc~compute_knot_vector compute_knot_vector proc~set2~2->proc~compute_knot_vector proc~set3~2->set proc~set3~2->proc~cmp_degree~2 proc~set4~2->set proc~kron3 kron3 interface~kron->proc~kron3 proc~kron_t1_t1 kron_t1_t1 interface~kron->proc~kron_t1_t1 proc~kron_t1_t2 kron_t1_t2 interface~kron->proc~kron_t1_t2 proc~cmp_degree~2->set proc~get_multiplicity~2 nurbs_volume%get_multiplicity proc~cmp_degree~2->proc~get_multiplicity~2 proc~cmp_nc~2->interface~compute_multiplicity proc~cmp_nc~2->set proc~cmp_tgc_3d->interface~kron proc~cmp_tgc_3d->proc~basis_bspline proc~repelem repelem proc~compute_knot_vector->proc~repelem proc~get_multiplicity~2->interface~compute_multiplicity

Variables

Type Attributes Name Initial
real(kind=rk), parameter :: a(3) = [90.0_rk, 270.0_rk, 360.0_rk]

Bend angles (degrees).
real(kind=rk), parameter :: c(3) = [0.0_rk, 0.0_rk, 0.0_rk]

Center of the ring.
real(kind=rk), parameter :: l = 1.0_rk

Length of the straight segment.
real(kind=rk), parameter :: r1 = 0.3_rk

Inner radius of the pipe.
real(kind=rk), parameter :: r2 = 0.5_rk

Outer radius of the pipe.
real(kind=rk), parameter :: rb = 1.5_rk

Bend radius.
type(nurbs_volume) :: ring

Straight pipe ring
type(nurbs_volume) :: shape

Bent pipe shape

Subroutines

pure subroutine bend_pipe(this, center, radius1, radius2, length, rbend, angle_deg)

Bend a straight pipe-like NURBS volume into a circular arc. Each control point is transformed relative to the center . Let the total bend angle in degrees be (input), and define the total bend angle in radians For each control point: Mapping: Knots and weights are preserved; only the control lattice is updated.

Arguments

Type IntentOptional Attributes Name
type(nurbs_volume), intent(inout) :: this

NURBS volume to be bent.
real(kind=rk), intent(in) :: center(3)

Pipe center coordinates .
real(kind=rk), intent(in) :: radius1

Inner radius of the pipe (not used in bending).
real(kind=rk), intent(in) :: radius2

Outer radius of the pipe (not used in bending).
real(kind=rk), intent(in) :: length

Length of the straight pipe segment before bending.
real(kind=rk), intent(in) :: rbend

Bend radius , i.e. distance from the bend centerline.
real(kind=rk), intent(in) :: angle_deg

Bend angle in degrees .