Уменьшайте модель до фиксированных интерфейсных режимов в заданном частотном диапазоне и граничных интерфейсных степеней свободы.

Создайте переходную структурную модель для 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);

Сгенерируйте mesh.

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]