Straathof R, Meijaard JP, van Vliet-Pérez SM, Kolkman-Deurloo IKK, Nout RA, Heijmen BJM, Wauben LSGL, Dankelman J, van de Berg NJ. Multibody dynamic modeling of the behavior of flexible instruments used in cervical cancer brachytherapy.
Med Phys 2024;
51:3698-3710. [PMID:
38226798 DOI:
10.1002/mp.16934]
[Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Revised: 10/24/2023] [Accepted: 12/09/2023] [Indexed: 01/17/2024] Open
Abstract
BACKGROUND
The steep radiation dose gradients in cervical cancer brachytherapy (BT) necessitate a thorough understanding of the behavior of afterloader source cables or needles in the curved channels of (patient-tailored) applicators.
PURPOSE
The purpose of this study is to develop and validate computer models to simulate: (1) BT source positions, and (2) insertion forces of needles in curved applicator channels. The methodology presented can be used to improve the knowledge of instrument behavior in current applicators and aid the development of novel (3D-printed) BT applicators.
METHODS
For the computer models, BT instruments were discretized in finite elements. Simulations were performed in SPACAR by formulating nodal contact force and motion input models and specifying the instruments' kinematic and dynamic properties. To evaluate the source cable model, simulated source paths in ring applicators were compared with manufacturer-measured source paths. The impact of discrepancies on the dosimetry was estimated for standard plans. To validate needle models, simulated needle insertion forces in curved channels with varying curvature, torsion, and clearance, were compared with force measurements in dedicated 3D-printed templates.
RESULTS
Comparison of simulated with manufacturer-measured source positions showed 0.5-1.2 mm median and <2.0 mm maximum differences, in all but one applicator geometry. The resulting maximum relative dose differences at the lateral surface and at 5 mm depth were 5.5% and 4.7%, respectively. Simulated insertion forces for BT needles in curved channels accurately resembled the forces experimentally obtained by including experimental uncertainties in the simulation.
CONCLUSION
The models developed can accurately predict source positions and insertion forces in BT applicators. Insights from these models can aid novel applicator design with improved motion and force transmission of BT instruments, and contribute to the estimation of overall treatment precision. The methodology presented can be extended to study other applicator geometries, flexible instruments, and afterloading systems.
Collapse