example_twist_taper Program

Uses

  • program~~example_twist_taper~~UsesGraph program~example_twist_taper example_twist_taper module~forcad forcad program~example_twist_taper->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 apply a progressive twist and linear taper to a hexahedral NURBS volume along its axial (z) direction.

The program: Creates a straight hexahedral block (control points nc = [7,7,9] over a box of size L), Copies it to shape, Applies a z-dependent twist up to twist_deg and a linear taper to factor taper, Exports the resulting NURBS volume to VTK, Displays the geometry.


Calls

program~~example_twist_taper~~CallsGraph program~example_twist_taper example_twist_taper proc~build_twist_taper build_twist_taper program~example_twist_taper->proc~build_twist_taper proc~create~2 nurbs_volume%create program~example_twist_taper->proc~create~2 proc~export_xc~2 nurbs_volume%export_Xc program~example_twist_taper->proc~export_xc~2 proc~export_xg~2 nurbs_volume%export_Xg program~example_twist_taper->proc~export_xg~2 proc~export_xth_in_xg~2 nurbs_volume%export_Xth_in_Xg program~example_twist_taper->proc~export_xth_in_xg~2 proc~set_hexahedron nurbs_volume%set_hexahedron program~example_twist_taper->proc~set_hexahedron proc~show~2 nurbs_volume%show program~example_twist_taper->proc~show~2 none~get_knot~2 nurbs_volume%get_knot proc~build_twist_taper->none~get_knot~2 none~get_wc~2 nurbs_volume%get_Wc proc~build_twist_taper->none~get_wc~2 none~get_xc~2 nurbs_volume%get_Xc proc~build_twist_taper->none~get_xc~2 none~set~2 nurbs_volume%set proc~build_twist_taper->none~set~2 proc~is_rational~2 nurbs_volume%is_rational proc~build_twist_taper->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~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 proc~set_hexahedron->none~set~2 proc~hexahedron_xc hexahedron_Xc proc~set_hexahedron->proc~hexahedron_xc 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_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~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 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~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->set interface~compute_multiplicity compute_multiplicity proc~cmp_nc~2->interface~compute_multiplicity proc~cmp_tgc_3d->interface~kron proc~cmp_tgc_3d->proc~basis_bspline proc~repelem repelem proc~compute_knot_vector->proc~repelem proc~compute_multiplicity1 compute_multiplicity1 interface~compute_multiplicity->proc~compute_multiplicity1 proc~compute_multiplicity2 compute_multiplicity2 interface~compute_multiplicity->proc~compute_multiplicity2 proc~get_multiplicity~2->interface~compute_multiplicity

Variables

Type Attributes Name Initial
real(kind=rk), parameter :: L(3) = [1.0_rk, 1.0_rk, 3.0_rk]

Domain extents in .
type(nurbs_volume) :: hexa
integer, parameter :: nc(3) = [7, 7, 9]

Control point counts per direction.
type(nurbs_volume) :: shape
real(kind=rk), parameter :: taper = 0.1_rk

Target in-plane scale at the top face (0<taper≤1).
real(kind=rk), parameter :: twist_deg = 360.0_rk

Total twist angle (degrees) at the top face.

Subroutines

pure subroutine build_twist_taper(this, Length, nc, twist_deg, taper)

Apply a z-progressive twist and linear taper to a NURBS hexahedron.

Read more…

Arguments

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

Volume to be transformed.
real(kind=rk), intent(in) :: Length(3)

Box lengths .
integer, intent(in) :: nc(3)

Control points sizes.
real(kind=rk), intent(in) :: twist_deg

Total twist at top face (degrees).
real(kind=rk), intent(in) :: taper

In-plane scale at top face (0<taper≤1).