Hi daeho,
Based on the description of your approach used to perform finite element analysis from the ANSYS data seems logical. Nonetheless, the issue with numerical instability and unrealistic results can be factored with some workaround and improvements:
1. You can ensure that the initial condition on which the ANSYS data was captured remains consistent with “ode45”. Along with condition, unit consistency can also be looked at since those increase the numeric values by factors of 10 as in this case.
2. A physical system with large sized matrices (2800 X 2800) has the possibility to become a stiff system for which “ode45” is not the most suited. MATLAB recommends “ode15s” for such cases. Additionally, there was a fixed step of 0.02 seconds while “ode45” has an adaptive size of step that might lead to discrepancies. Please refer to the below documentation to know more about the different types of solver and their applications:
3. The reverse calculation of forces with the equation mentioned f=mx’’ +kx should be applied appropriately especially, x’’ for acceleration and x for deformation, since any mistake in these can change the factors/scale of the results.
Hope this helps!
