针对目前商业有限元软件中缺少挠曲电模块且有限元自主编程数值模拟挠曲电问题的复杂性问题，提出在COMSOL软件中的偏微分方程模块(PDE模块),对考虑挠曲电效应的力电耦合问题直接进行数值模拟的方法.首先，通过数值计算固定边界位移的挠曲电圆筒，验证有限元数值计算方法的正确性.接着，研究挠曲电系数和不同材料梯度指数对圆筒力电耦合的影响.数值结果表明：在微米尺度下，增大挠曲电系数对径向位移的影响很小，但是对圆筒内部的电势和电位移分布影响很大.增大材料梯度指数对挠曲电圆筒内部力电耦合影响非常明显，圆筒内部径向位移相对边界处变化范围逐渐减小，圆筒内部电势呈现负值逐渐减小后转为正数的趋势，而圆筒内部电位移值负值逐渐增大.本工作可以为挠曲电效应数值计算及挠曲电结构设计提供帮助.