✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。
🍎个人主页：Matlab科研工作室
🍊个人信条：格物致知。
⛄ 内容介绍
当今世界,无论是国防尖端武器,航空航天飞行器,高端建筑设备,机器人,还是在日常生活中家电的控制器及摄影设备中,都记载着智能材料的点滴发展.因此,智能材料已经在明显提高着人民大众的生活质量.由于具有独特的压电效应,压电材料现已在智能结构的研究与工程应用中占有重要的地位. 由于位移场和电场的耦合作用,当涉及到非线性压电效应时,反映物理问题的数学模型的求解变得非常复杂.
⛄ 部分代码
function [pb]=PlateInplaneNonlinear(nnel,dshapedx,dshapedy,u_p,v_p,w_p)
%--------------------------------------------------------------------------
% Purpose:
% determine the kinematic matrix expression relating nonlinear inplane action
% Synopsis:
% [pb]=PlateInplaneNonlinear(nnel,dshapedx,dshapedy,u_p,v_p,w_p)
%
% Variable Description:
% nnel - number of nodes per element
% dshapedx - derivatives of shape functions with respect to x
% dshapedy - derivatives of shape functions with respect to y
%--------------------------------------------------------------------------
dudx=dshapedx(1)*u_p(1)+dshapedx(2)*u_p(2)+dshapedx(3)*u_p(3)+dshapedx(4)*u_p(4);
dudy=dshapedy(1)*u_p(1)+dshapedy(2)*u_p(2)+dshapedy(3)*u_p(3)+dshapedy(4)*u_p(4);
dvdx=dshapedx(1)*v_p(1)+dshapedx(2)*v_p(2)+dshapedx(3)*v_p(3)+dshapedx(4)*v_p(4);
dvdy=dshapedy(1)*v_p(1)+dshapedy(2)*v_p(2)+dshapedy(3)*v_p(3)+dshapedy(4)*v_p(4);
dwdx=dshapedx(1)*w_p(1)+dshapedx(2)*w_p(2)+dshapedx(3)*w_p(3)+dshapedx(4)*w_p(4);
dwdy=dshapedy(1)*w_p(1)+dshapedy(2)*w_p(2)+dshapedy(3)*w_p(3)+dshapedy(4)*w_p(4);
for i=1:nnel
i1=(i-1)*5+1;
i2=i1+1;
i3=i2+1;
i4=i3+1;
i5=i4+1;
pb(1,i3)=dwdx*dshapedx(i);
pb(2,i3)=dwdy*dshapedy(i);
pb(3,i3)=dwdy*dshapedx(i)+dwdx*dshapedy(i);
pb(4,i4)=0;
pb(5,i5)=0;
pb(6,i5)=0;
end
⛄ 运行结果
⛄ 参考文献
[1]顾春霞. 强电场作用下压电层合板的非线性弯曲分析[D]. 广西大学, 2013.
[2]石艺娜. 压电层合圆板传感器件的非线性弯曲与振动分析[D]. 兰州大学, 2001.
