Dear @Ran.Tu
As Jason mentioned, we have recently looked into the topic of extracting 6x6 matrices from the linearization outputs of OpenFAST. I’m enclosing an alternative document that derives the relationships between the Jacobians at a source node and destination node: MappingLoadJacobians
There are two ways to obtain the 6x6 matrix from MAP at the ElastoDyn reference point:
-
using MAP++ (note that to transfer from the origin to the reference point of elastodyn, one needs the value of the force at the operating point. This was added very recently to MAP++. I therefore recommend using the dev branch of pyMAP. Look for instance in example 1 or 3 here
-
by extracting the outputs from a linearization analysis with OpenFAST (requires LinInputs = 2, LinOutputs = 2, LinOutJac = True), and applying the relationship between load jacobians. A example of implementation is provided here
A (non-running) example of both approaches is available here
I hope that helps,
Emmanuel