Static Kirchhoff Rods under the Action of External Forces: Integration via Runge-Kutta Method

<table>Frenet helices for circular cross section filaments and increasing values of <svg height="10.3875" id="M174" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 5.9250002 10.3875" width="5.9250002" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.125)"><path d="M324 430l-26 -36l-112 -4l-55 -265q-13 -66 7 -66q13 0 44.5 20t50.5 40l17 -24q-38 -40 -85.5 -73.5t-87.5 -33.5q-50 0 -21 138l55 262h-80l-2 8l25 34h66l25 99l78 63l10 -9l-37 -153h128z" id="x1D461"></path></g>
</svg>. As usual, <svg height="8.4250002" id="M175" style="vertical-align:-0.273pt" version="1.1" viewbox="0 0 8.5874996 8.4250002" width="8.5874996" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,8.013)"><path d="M475 445l-64 -69q-119 -128 -289 -280q27 6 81 -15q66 -26 102 -26q48 0 111 93l23 -19q-41 -84 -75.5 -116.5t-73.5 -32.5q-37 0 -147 53q-43 21 -63 7q-22 -15 -38 -41q-4 -7 -11 -7q-8 10 -8 26q0 38 47 83l292 273q-32 -5 -67 -1q-63 7 -105 7q-32 0 -53.5 -22
t-47.5 -67l-25 13q15 39 38 78q40 66 83 66q57 0 120 -12q52 -10 80 -10q29 0 67 42z" id="x1D467"></path></g>
</svg> is the vertical axis of the box.</table>

Journal of Computational Methods in Physics

fig3

Figure 3

Figure 3: Static Kirchhoff Rods under the Action of External Forces: Integration via Runge-Kutta Method 