2017-07-14

• I have completed a rough, but functioning script that calculates the modal force profiles. The force values are still coming out incorrect (on the order of 10^14) but the script can take in my model as a .m file and return an array with a force value per mode. I am attaching both the .m file and the matlab script
• I have done very little work with the numerical integration itself, based on the 2D numerical integration code I received I just appended a z component and left it at that so when I return from Livingston I will  fix that component

par.a = 75e-3/2;    % radius [m]
par.h = 1.004e-3;   % thickness [m]
par.E = 73.2e9;     % Young's modulus [Pa]
par.nu = 0.155;     % Poisson's ratio
par.rho = 2202;     % density [kg/m^3]

%Calculate fundamental modes of the disk
[freqs, modes, shapes, x, y] = disk_frequencies(par, 10000, 1, 'shapes', 0.5e-3);

%Now we extract the force profile from the COMSOL model

function out = model
%
% faster.m
%
% Model exported on Jul 14 2017, 14:47 by COMSOL 5.2.1.262.

import com.comsol.model.*
import com.comsol.model.util.*

model = ModelUtil.create('Model');