example_helix_pipe Program

Uses

  • program~~example_helix_pipe~~UsesGraph program~example_helix_pipe example_helix_pipe module~forcad forcad program~example_helix_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 sweep a straight, pipe-like NURBS volume into a cylindrical helix with a prescribed radius, pitch, and number of turns.

The program: Creates a straight ring extruded along (a pipe segment), Refines it along the axial direction (knot insertion + degree elevation), Maps the control points onto a helix of radius rh, pitch pitch, with nturns turns, Exports the resulting NURBS volume to VTK, Displays the geometry.


Calls

program~~example_helix_pipe~~CallsGraph program~example_helix_pipe example_helix_pipe proc~build_helix build_helix program~example_helix_pipe->proc~build_helix proc~create~2 nurbs_volume%create program~example_helix_pipe->proc~create~2 proc~elevate_degree~2 nurbs_volume%elevate_degree program~example_helix_pipe->proc~elevate_degree~2 proc~export_xc~2 nurbs_volume%export_Xc program~example_helix_pipe->proc~export_xc~2 proc~export_xg~2 nurbs_volume%export_Xg program~example_helix_pipe->proc~export_xg~2 proc~export_xth_in_xg~2 nurbs_volume%export_Xth_in_Xg program~example_helix_pipe->proc~export_xth_in_xg~2 proc~insert_knots~2 nurbs_volume%insert_knots program~example_helix_pipe->proc~insert_knots~2 proc~set_ring~2 nurbs_volume%set_ring program~example_helix_pipe->proc~set_ring~2 proc~show~2 nurbs_volume%show program~example_helix_pipe->proc~show~2 none~get_knot~2 nurbs_volume%get_knot proc~build_helix->none~get_knot~2 none~get_nc~2 nurbs_volume%get_nc proc~build_helix->none~get_nc~2 none~get_wc~2 nurbs_volume%get_Wc proc~build_helix->none~get_wc~2 none~get_xc~2 nurbs_volume%get_Xc proc~build_helix->none~get_xc~2 none~set~2 nurbs_volume%set proc~build_helix->none~set~2 proc~is_rational~2 nurbs_volume%is_rational proc~build_helix->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
integer, parameter :: N = 50

interior knots to insert along .
real(kind=rk), parameter :: center(3) = [0.0_rk, 0.0_rk, 0.0_rk]

Pipe center .
integer :: i
real(kind=rk), parameter :: length = 1.0_rk

Length of the straight segment.
integer, parameter :: nturns = 3

Number of full turns.
real(kind=rk), parameter :: pitch = 1.5_rk

Helix pitch (axial rise per full turn).
real(kind=rk), parameter :: r1 = 0.1_rk

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

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

Helix radius.
type(nurbs_volume) :: ring

Straight pipe ring.
type(nurbs_volume) :: shape

Helical pipe shape.

Subroutines

pure subroutine build_helix(this, c, rh, p, n)

Map a straight pipe-like NURBS volume onto a cylindrical helix.

Read more…

Arguments

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

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

Center , helix radius , pitch .
real(kind=rk), intent(in) :: rh

Center , helix radius , pitch .
real(kind=rk), intent(in) :: p

Center , helix radius , pitch .
integer, intent(in) :: n

Number of turns.