search

GazzolaLab / PyElastica

Python implementation of Elastica, an open-source software for the simulation of assemblies of slender, one-dimensional structures using Cosserat Rod theory.
https://www.cosseratrods.org
MIT License
228 stars 110 forks source link

The values of Q_i, e_i , D_i and \kappa_i for the last spatial node i=n+1. #418

Closed Srishti-Siddharth-iitb closed 2 months ago

Srishti-Siddharth-iitb commented 2 months ago

I am trying to simulate a single Cosserat rod in MATLAB using your codes and the numerical scheme discussed in the following 2018 paper: https://doi.org/10.1098/rsos.171628. To stimulate the discrete equations (3.8) and (3.9) of the paper, the values of Q_i, e_i, D_i and \kappa_i are needed for the last spatial node i.e. i=n+1. However, according to equations (3.3), (3.4) and (3.5) of the paper, the above quantities are not described for i=n+1. I wanted to understand how these quantities have been computed in your Python code. It would be helpful if you could provide the line number in the Python code which can be used as a reference to look at the computation of these quantities for the last spatial node. Also, do ghost nodes play a role here?

armantekinalp commented 2 months ago

You can take a look to this file to see how rod equations are computed. I am also moving this issue to discussion.