Сокращение модели несущих конструкций
Rcb = reduce(structuralmodel,'FrequencyRange',[omega1,omega2])[omega1,omega2] и граничные границы степеней свободы.
Уменьшите модель до фиксированных режимов интерфейса в указанном диапазоне частот и степенях свободы граничного интерфейса.
Создайте переходную структурную модель для 3-D проблемы.
structuralmodel = createpde('structural','transient-solid');
Создайте геометрию и включите ее в модель. Постройте график геометрии.
gm = multicuboid(0.1,0.01,0.01); structuralmodel.Geometry = gm; pdegplot(structuralmodel,'FaceLabels','on','FaceAlpha',0.5)
Задайте модуль Юнга, коэффициент Пуассона и массовую плотность материала.
structuralProperties(structuralmodel,'YoungsModulus',70E9, ... 'PoissonsRatio',0.3, ... 'MassDensity',2700);
Создайте сетку.
generateMesh(structuralmodel);
Укажите концы балки в качестве интерфейсов несущих виражей. Метод модели с уменьшенным порядком сохраняет степени свободы на интерфейсах виража при конденсации степеней свободы на всех других границах. Для повышения производительности используйте набор кромок, ограничивающих каждую сторону балки, а не всю грань.
structuralSEInterface(structuralmodel,'Edge',[4,6,9,10]); structuralSEInterface(structuralmodel,'Edge',[2,8,11,12]);
Сокращение модели до фиксированных режимов интерфейса в диапазоне частот [-Inf,500000] и граничные границы степеней свободы.
R = reduce(structuralmodel,'FrequencyRange',[-Inf,500000])R =
ReducedStructuralModel with properties:
K: [166x166 double]
M: [166x166 double]
NumModes: 22
RetainedDoF: [144x1 double]
ReferenceLocations: []
Mesh: [1x1 FEMesh]
reconstructSolution | ReducedStructuralModel | solve | structuralBC | StructuralModel | structuralSEInterface