Dear @Fanny.Callisaya,
I’m not aware that NREL has developed BModes models for the blades of the IEA Wind 15-MW or 22-MW reference wind turbines. I was not the one to develop these models, but I believe the mode shapes for these models were derived via Frame3DD rather than through BModes. That said, the inputs set in the beam section properties file of BModes are very similar to those of the ElastoDyn blade file, so, it should not be hard to create these files yourselves. For the columns in the beam section properties file of BModes not specified in ElastoDyn (e.g., flp_iner), see the guidance from the following forum topic: NREL 5MW Blade Mode Shape.
BModes does not have functionality to model TMDs. However, if you’ve already derived the mode shapes for ElastoDyn and are modeling the system in OpenFAST, you can see the effect of TMDs on the blade response through the OpenFAST results.
Best regards,