【有限元动态水波模拟】基于有限元的三维动态水波运动模拟matlab仿真

1.软件版本

MATLAB2021a

2.核心代码

% Finite element solution of wave equation
 
close all;
load meshdata;
% c(x,y)===================================================================
c=@(X)1;
% =========================================================================
numberofnodes=size(p,1);
numberofelements=size(t,1);
% Initialisation of K, M ==================================================
K=zeros(numberofnodes,numberofnodes);
M=zeros(numberofnodes,numberofnodes);
% Computation of Ke, Me & assemby of K, M==================================
for e=1:numberofelements
   nodes=t(e,:); % row with nodenumbers of element e
   P=[ones(1,3);p(nodes,:)'];
   Coeffs=inv(P);
   areaelemente=abs(det(P))/2; % area of element e
   diffPhi=Coeffs(:,2:end);
   Ke=diffPhi*diffPhi'*areaelemente;
   Me=1/3^2*ones(3)*1/(c(mean(p(nodes,:))))^2*areaelemente;
   K(nodes,nodes)=K(nodes,nodes)+Ke;
   M(nodes,nodes)=M(nodes,nodes)+Me;
end
Mlumped=diag(M(1:size(M,1)+1:end));
M=Mlumped;
% Dirichlet boundaries=====================================================
% =unique(boundedges(p,t));
Dirichlet=[

...............................................

3.操作步骤与仿真结论

 5.参考文献

[1]杨永丰, 王玉玮, 樊瑜波. 三维有限元动态分析方法的建立[C]// 中华口腔医学会第三次全国口腔修复学术会议. 0.

D215

6.完整源码获得方式

方式1：微信或者QQ联系博主

方式2：订阅MATLAB/FPGA教程，免费获得教程案例以及任意2份完整源码