MATLABでの超音波ホーンの最適設計法について

超音波を増幅させる鉄製のホーンの最適設計で、現状MATLABの中心複合計画で計画した因子と水準で、AbaqusでFEM解析を実施して、

ホーン先端の変位（振幅）を最大化するように、応答曲面でホーンの寸法を最適化しようとしています。３６回も解析モデルを作成し、計算させるのは非効率と思えてきました。MATLABでFEM解析ができればもっと効率があがると思います。MATLABでのFEM解析のToolboxなどあればご教授ください、また応答曲面法以外の最適化法もご教授頂けますと幸いです。よろしくお願いいたします。