Dear Yuqi,
I would not expect a perfect match to the mode shapes in the “NRELOffshrBsline5MW_Onshore_ElastoDyn_Tower.dat” file because I derived those mode shapes using ADAMS (with a flexible rotor) in place of BModes, as discussed in the following forum topics:
*NREL 5MW_Onshore: Using Bmodes & ModeShapePolyfitting.xls - #8 by Jason.Jonkman
*BModes : Input parameters about tower support subsystem
*Tower Eigenfrequencies of NREL 5MW Turbine
I see that you are using BModes v1.03 rather than the alpha version (BModesJJ). You may want to see if upgrading to the alpha version provides better results: nwtc.nrel.gov/BModes.
Best regards,