【机械】板的线性和非线性弯曲分析附MATLAB 代码

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，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.

⛄ 完整代码

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料