I believe Ansys, COMSOL, and Abaqus are all good software packages to use for this application, although COMSOL may be the easiest to setup a piezoelectric Multiphysics.
The paper below is an example of a simulation (performed in Ansys) that provides a description of the model setup and key parameters of interest.
Article Kinetic Energy Harvesting for Wearable Medical Sensors
Since in energy harvesting applications you are dealing with Multiphysics problems (conversion from one form of energy to another form of energy), the most popular software in this area is COMSOL. I agree that you can use Ansys or Abaqus. But COMSOL is more straightforward to use.
For the commercial software, I recommend the COMSOL. It is more straightforward, and easy to learn. As far as the ANSYS is concerned, you should be careful with the IEEE and ANSYS forms about the piezoelectric material properties.
The Matlab is suitable for secondary development or structural optimization (eg, topology optimization research, see my published paper).@Abbas Homayouni Amlashi
[1] He M , Zhang X , Fernandez L , et al. Multi-material topology optimization of piezoelectric composite structures for energy harvesting[J]. Composite Structures, 2021, 265(1):113783.
[2] Homayouni-Amlashi A , Schlinquer T , · Abdenbi Mohand-Ousaid, et al. 2D Topology Optimization MATLAB Codes for Piezoelectric Actuators and Energy Harvesters[J]. Structural and Multidisciplinary Optimization, 2020.
For COMSOL
[1] Wang L , Zhao L , Jiang Z , et al. High accuracy comsol simulation method of bimorph cantilever for piezoelectric vibration energy harvesting[J]. AIP Advances, 2019, 9(9):095067.
COMSOL is very good, but you have to start with the theoretical basis of this this topic, means starting from the Fourth Order Partial Differential Equation.