Уменьшайте структурную модель
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);
Сгенерируйте 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]
reconstructSolution | StructuralModel | solve | structuralBC | structuralSEInterface | ReducedStructuralModel