40m QIL Cryo_Lab CTN SUS_Lab TCS_Lab OMC_Lab CRIME_Lab FEA ENG_Labs OptContFac Mariner WBEEShop
  Coating Ring-down Measurement Lab elog  Not logged in ELOG logo
Message ID: 368     Entry time: Fri Jul 14 16:43:24 2017
Author: Zach 
Type: Electronics 
Category: Modeling 
Subject: Force profile matlab script 

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
Attachment 1: forces.m  1 kB  | Hide | Hide all
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
... 41 more lines ...
Attachment 2: faster.m  24 kB  | Hide | Hide all
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');
... 403 more lines ...
ELOG V3.1.3-