example_toroidal_pipe Program

Uses

  • program~~example_toroidal_pipe~~UsesGraph program~example_toroidal_pipe example_toroidal_pipe module~forcad forcad program~example_toroidal_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 onto a toroidal pipe (donut) with an optional progressive cross-section twist and a sinusoidal wobble in the global z-direction as it travels around the torus.

The program: Creates a straight pipe segment (ring extruded along ), Refines it along the axial direction (knot insertion + degree elevation), Maps the control points onto a torus of major radius R, Applies a cross-section twist of twist_turns full rotations over one loop, Superimposes a vertical wobble , Exports the resulting NURBS volume to VTK and displays it, Prints the approximate enclosed volume (Gaussian integration).


Calls

program~~example_toroidal_pipe~~CallsGraph program~example_toroidal_pipe example_toroidal_pipe proc~create nurbs_volume%create program~example_toroidal_pipe->proc~create proc~elevate_degree nurbs_volume%elevate_degree program~example_toroidal_pipe->proc~elevate_degree proc~export_xc nurbs_volume%export_Xc program~example_toroidal_pipe->proc~export_xc proc~export_xg nurbs_volume%export_Xg program~example_toroidal_pipe->proc~export_xg proc~export_xth_in_xg nurbs_volume%export_Xth_in_Xg program~example_toroidal_pipe->proc~export_xth_in_xg proc~insert_knots nurbs_volume%insert_knots program~example_toroidal_pipe->proc~insert_knots proc~map_to_torus_sinez map_to_torus_sineZ program~example_toroidal_pipe->proc~map_to_torus_sinez proc~set_ring nurbs_volume%set_ring program~example_toroidal_pipe->proc~set_ring proc~show nurbs_volume%show program~example_toroidal_pipe->proc~show interface~compute_xg compute_Xg proc~create->interface~compute_xg interface~ndgrid ndgrid proc~create->interface~ndgrid proc~is_rational nurbs_volume%is_rational proc~create->proc~is_rational set set proc~create->set none~get_knot nurbs_volume%get_knot proc~elevate_degree->none~get_knot none~set nurbs_volume%set proc~elevate_degree->none~set proc~elevate_degree_a_5_9 elevate_degree_A_5_9 proc~elevate_degree->proc~elevate_degree_a_5_9 proc~elevate_degree->proc~is_rational proc~cmp_elem_xc_vis nurbs_volume%cmp_elem_Xc_vis proc~export_xc->proc~cmp_elem_xc_vis proc~export_vtk_legacy export_vtk_legacy proc~export_xc->proc~export_vtk_legacy proc~export_xc->set proc~cmp_elem_xg_vis nurbs_volume%cmp_elem_Xg_vis proc~export_xg->proc~cmp_elem_xg_vis proc~export_xg->proc~export_vtk_legacy proc~export_xg->set proc~export_xth_in_xg->interface~compute_xg interface~unique unique proc~export_xth_in_xg->interface~unique proc~export_xth_in_xg->proc~export_vtk_legacy proc~export_xth_in_xg->proc~is_rational proc~export_xth_in_xg->set interface~compute_multiplicity compute_multiplicity proc~insert_knots->interface~compute_multiplicity proc~insert_knots->none~get_knot proc~findspan findspan proc~insert_knots->proc~findspan proc~insert_knot_a_5_1 insert_knot_A_5_1 proc~insert_knots->proc~insert_knot_a_5_1 proc~insert_knots->proc~is_rational s_loc s_loc proc~insert_knots->s_loc proc~insert_knots->set proc~map_to_torus_sinez->none~get_knot none~get_nc nurbs_volume%get_nc proc~map_to_torus_sinez->none~get_nc none~get_wc nurbs_volume%get_Wc proc~map_to_torus_sinez->none~get_wc none~get_xc nurbs_volume%get_Xc proc~map_to_torus_sinez->none~get_xc proc~map_to_torus_sinez->none~set proc~map_to_torus_sinez->proc~is_rational proc~set_ring->none~set 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->proc~compute_xg_bspline_3d proc~compute_xg_bspline_3d_1point compute_Xg_bspline_3d_1point interface~compute_xg->proc~compute_xg_bspline_3d_1point proc~compute_xg_nurbs_3d compute_Xg_nurbs_3d interface~compute_xg->proc~compute_xg_nurbs_3d proc~compute_xg_nurbs_3d_1point compute_Xg_nurbs_3d_1point interface~compute_xg->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 nurbs_volume%get_knot_all none~get_knot->proc~get_knot_all proc~get_knoti nurbs_volume%get_knoti none~get_knot->proc~get_knoti proc~get_nc_all nurbs_volume%get_nc_all none~get_nc->proc~get_nc_all proc~get_nc_dir nurbs_volume%get_nc_dir none~get_nc->proc~get_nc_dir proc~get_wc_all nurbs_volume%get_Wc_all none~get_wc->proc~get_wc_all proc~get_wci nurbs_volume%get_Wci none~get_wc->proc~get_wci proc~get_xc_all nurbs_volume%get_Xc_all none~get_xc->proc~get_xc_all proc~get_xci nurbs_volume%get_Xci none~get_xc->proc~get_xci proc~get_xcid nurbs_volume%get_Xcid none~get_xc->proc~get_xcid proc~set1 nurbs_volume%set1 none~set->proc~set1 proc~set2 nurbs_volume%set2 none~set->proc~set2 proc~set3 nurbs_volume%set3 none~set->proc~set3 proc~set4 nurbs_volume%set4 none~set->proc~set4 interface~elemconn_c0 elemConn_C0 proc~cmp_elem_xc_vis->interface~elemconn_c0 proc~cmp_elem_xg_vis->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->interface~compute_multiplicity local local proc~ndgrid2->local proc~ndgrid3->local proc~set1->set proc~cmp_degree nurbs_volume%cmp_degree proc~set1->proc~cmp_degree proc~cmp_nc nurbs_volume%cmp_nc proc~set1->proc~cmp_nc proc~set2->set proc~set2->proc~cmp_nc proc~compute_knot_vector compute_knot_vector proc~set2->proc~compute_knot_vector proc~set3->set proc~set3->proc~cmp_degree proc~set4->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->set proc~get_multiplicity nurbs_volume%get_multiplicity proc~cmp_degree->proc~get_multiplicity proc~cmp_nc->interface~compute_multiplicity proc~cmp_nc->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->interface~compute_multiplicity

Variables

Type Attributes Name Initial
real(kind=rk), parameter :: Az = 0.1_rk

amplitude of vertical wobble
integer, parameter :: Nref = 30
real(kind=rk), parameter :: R = 1.00_rk

major radius
real(kind=rk), parameter :: center(3) = [0.0_rk, 0.0_rk, 0.0_rk]

center of the ring
integer :: i
real(kind=rk), parameter :: length = 1.0_rk

length of the pipe
integer, parameter :: nwaves_z = 6

number of sine waves per full loop
real(kind=rk), parameter :: phase_z = 0.1_rk

phase shift (radians)
real(kind=rk), parameter :: r1 = 0.12_rk

inner radius
real(kind=rk), parameter :: r2 = 0.20_rk

outer radius
type(nurbs_volume) :: ring

straight pipe segment
type(nurbs_volume) :: sh

toroidal pipe shape
real(kind=rk), parameter :: twist_turns = 1.0_rk

cross-section twist (turns over full loop)

Subroutines

pure subroutine map_to_torus_sineZ(this, c, R, twist_turns, Az, nwaves_z, phase_z)

Map a straight pipe (ring extruded in ) onto a toroidal pipe with: • cross-section twist: , • vertical wobble: .

Read more…

Arguments

Type IntentOptional Attributes Name
type(nurbs_volume), intent(inout) :: this
real(kind=rk), intent(in) :: c(3)
real(kind=rk), intent(in) :: R
real(kind=rk), intent(in) :: twist_turns
real(kind=rk), intent(in) :: Az
integer, intent(in) :: nwaves_z
real(kind=rk), intent(in) :: phase_z