|
a |
|
b/Application/Examples/BikeSpring/Model/objective.m |
|
|
1 |
function [F] = objective(x) |
|
|
2 |
|
|
|
3 |
F = 0; |
|
|
4 |
|
|
|
5 |
% AnyVar L0 = 0.24; //Spring slack length |
|
|
6 |
% AnyVar F1 = 100; //Spring force at 50% strain |
|
|
7 |
% AnyVar SpringRad = 0.15; //Radius of spring attachment on the crank |
|
|
8 |
% AnyVar SpringAngle = 45; //Angle of spring attachment points on the crank |
|
|
9 |
% AnyVar SpringX = -0.10; //X coordinate of spring attachment on frame |
|
|
10 |
% AnyVar SpringY = 0.30; //Y coordinate of spring attachment on frame |
|
|
11 |
|
|
|
12 |
% Create include file with variable values |
|
|
13 |
fid=fopen('AnyInc.any','wt'); |
|
|
14 |
fprintf(fid,'AnyVar L0 = %12.4f;\n',x(1)); |
|
|
15 |
fprintf(fid,'AnyVar F1 = %12.4f;\n',x(2)); |
|
|
16 |
fprintf(fid,'AnyVar SpringRad = %12.4f;\n',x(3)); |
|
|
17 |
fprintf(fid,'AnyVar SpringAngle = %12.4f;\n',x(4)); |
|
|
18 |
fprintf(fid,'AnyVar SpringX = %12.4f;\n',0.0); % Taken out of the problem |
|
|
19 |
fprintf(fid,'AnyVar SpringY = %12.4f;\n',x(5)); |
|
|
20 |
fclose(fid); |
|
|
21 |
|
|
|
22 |
fprintf(1,'\n Running simulation...\n'); |
|
|
23 |
dos('anybodycon springeval.anymcr'); |
|
|
24 |
%[s,w]=dos('c:\program files\anybody\anybodycon springeval.anymcr'); |
|
|
25 |
%while ans~=0%length(w)<5747 |
|
|
26 |
% pause(1) |
|
|
27 |
%end |
|
|
28 |
|
|
|
29 |
Data = csvread('MaxActivity.csv',13,0); |
|
|
30 |
time = Data(:,1); |
|
|
31 |
b = Data(:,2); |
|
|
32 |
F=max(b); |