lisa.pizzo/LisaPizzoExercises
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Delete Project/AddMass_mult_rot.m

Browse source
This commit is contained in:
Lisa Pizzo 2026-01-26 20:38:07 +01:00
commit 358bc227a3
1 changed files with 0 additions and 35 deletions
Download patch file Download diff file Expand all files Collapse all files

35
Project/AddMass_mult_rot.m
View file

@ -1,35 +0,0 @@
function M = AddMass_mult_rot(model, M, c_wall, c_fluid, c_air)
mesh = model.Mesh;
nodes = mesh.Nodes;
elements = mesh.Elements;
Nelems = size(elements,2);
regions = zeros(Nelems,1);
regions(findElements(mesh,'region','Face',1)) = 1; % wall
regions(findElements(mesh,'region','Face',2)) = 2; % fluid
regions(findElements(mesh,'region','Face',3)) = 3; % air
for e = 1:Nelems
vert = elements(:,e);
x = nodes(1,vert); % r-coordinates
y = nodes(2,vert); % z-coordinates
Ae = polyarea(x,y); % element area
rbar = mean(x); % mean radius (axisymmetric weight)
switch regions(e) % Select heat capacity
case 1
c = c_wall;
case 2
c = c_fluid;
case 3
c = c_air;
end
% Axisymmetric element mass matrix
Me = rbar * c * Ae / 12 * [2 1 1; 1 2 1; 1 1 2];
% Assemble
M(vert,vert) = M(vert,vert) + Me;
end
end