Уменьшите структурную модель
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 | ReducedStructuralModel | solve | structuralBC | StructuralModel | structuralSEInterface