4.4 Пример использования команды

В качестве примера рассмотрим расчет собственных форм и частот пружины жесткостью k=1e+5 [кг/с^2] с грузом весом в m=1 [кг]. На груз действует гармоническая сила амплитудой 100 [Н]. Сила задана через абстрактную функцию f.

Рассчитывается частотная характеристика датчика длины пружины от величины силы. Анализ проводится для списка частот от 10 [Гц] до 70 [Гц] с шагом 1 [Гц].

Текст проекта

body body1=body( color = RGB( 102, 255, 204 ) );

body body2=body( color = RGB( 204, 204, 204 ) );

point p=point( 0 [m], 0 [m], 0 [m] );

solid solid1=sphere( p, 1 [m], mass = 1 [kg] );

body body2 < ( solid1 );

force force1=spring( body1, p, body2, p, 1e5 [kg/s2] );

function f=abstract( "[]", list( "[s]" ) );

function func(t [s])=100 [N] * f(#t);

sensor input=func(time);

force force2=force( body2, p, projectX, func, list( time ) );

sensor output=bodyDisplacement( body1, p, projectX, body2, p );

command Vibration=vibrationResponse( list( output ), range( 10 [Hz], 70 [Hz], 1 [Hz] ), freqDamping = yesDamping: );

command VibrationNL=vibrationResponseNL( f, input, list( output ), range( 10 [Hz], 70 [Hz], 1 [Hz] ), 10, comDynamics );

/\///////////////////////////////////////////////////////////////////////////////////

set ground = body1;

/\///////////////////////////////////////////////////////////////////////////////////

/\ Единицы измерения;

set units = SI;

Результатом выполнения команды являются частотные характеристики системы: амплитудно-частотная характеристика (рис. 19), фазо-частотная характеристика (рис. 20) и амплитудно-фазовая частотная характеристика (рис. 21).

Рис. 19. Амплитудно-частотная характеристика системы

Рис. 20. Фазо-частотная характеристика системы

Рис. 21. Амплитудно-фазовая частотная характеристика системы

Пик АЧХ приходится на частоту 50 [Гц], что примерно совпадает с теоретической частотой колебаний пружины:

[Гц].