Exponenta Banner
exponenta event banner

interpolateAcceleration

Интерполяция ускорения в произвольных пространственных местоположениях для всех временных или частотных шагов для динамической структурной модели

свернуть все на странице

Синтаксис

intrp Accel = interpolateAcceleration (структурный результат, xq, yq)
intrp Accel = interpolateAcceleration (структурный результат, xq, yq, zq)
intrp Accel = interpolateAcceleration (структурные результаты, точки запроса)

Описание

intrpAccel = interpolateAcceleration(structuralresults,xq,yq) возвращает интерполированные значения ускорения в точках 2-D, указанных в xq и yq для всех временных или частотных шагов.

пример

intrpAccel = interpolateAcceleration(structuralresults,xq,yq,zq) использует точки 3-D, указанные в xq, yq, и zq.

intrpAccel = interpolateAcceleration(structuralresults,querypoints) использует точки, указанные в querypoints.

Примеры

свернуть все

Открыть сценарий в реальном времени

Интерполировать ускорение в геометрическом центре пучка при гармоническом возбуждении

Создайте переходную динамическую модель для 3-D проблемы.

structuralmodel = createpde('structural','transient-solid');

Создайте геометрию и включите ее в модель. Постройте график геометрии.

gm = multicuboid(0.06,0.005,0.01);
structuralmodel.Geometry = gm;
pdegplot(structuralmodel,'FaceLabels','on','FaceAlpha',0.5)
view(50,20)

Figure contains an axes. The axes contains 3 objects of type quiver, patch, line.

Задайте модуль Юнга, коэффициент Пуассона и массовую плотность материала.

structuralProperties(structuralmodel,'YoungsModulus',210E9, ...
                                     'PoissonsRatio',0.3, ...
                                     'MassDensity',7800);

Закрепите один конец балки.

structuralBC(structuralmodel,'Face',5,'Constraint','fixed');

Применить синусоидальное смещение вдоль y- направление на конце, противоположном фиксированному концу балки.

structuralBC(structuralmodel,'Face',3,'YDisplacement',1E-4,'Frequency',50);

Создайте сетку.

generateMesh(structuralmodel,'Hmax',0.01);

Задайте нулевое начальное смещение и скорость.

structuralIC(structuralmodel,'Displacement',[0;0;0],'Velocity',[0;0;0]);

Решите модель.

tlist = 0:0.002:0.2;
structuralresults = solve(structuralmodel,tlist);

Интерполяция ускорения в геометрическом центре балки.

coordsMidSpan = [0;0;0.005];
intrpAccel = interpolateAcceleration(structuralresults,coordsMidSpan);

Постройте график y- компонент ускорения геометрического центра балки.

figure
plot(structuralresults.SolutionTimes,intrpAccel.ay)
title('Y-Acceleration of the Geometric Center of the Beam')

Figure contains an axes. The axes with title Y-Acceleration of the Geometric Center of the Beam contains an object of type line.

Входные аргументы

свернуть все

Решение задачи динамического структурного анализа, определяемого как TransientStructuralResults или FrequencyStructuralResults объект. Создать structuralresults с помощью solve функция.

Пример: structuralresults = solve(structuralmodel,tlist)

точки запроса координат x, заданные как вещественный массив. interpolateAcceleration вычисляет ускорения в 2-D координатных точках [xq(i),yq(i)] или в точках координат 3-D [xq(i),yq(i),zq(i)]. Поэтому xq, yq, и (при наличии) zq должно иметь одинаковое количество записей.

interpolateAcceleration преобразует точки запроса в векторы столбцов xq(:), yq(:), и (при наличии) zq(:). Функция возвращает ускорения как FEStruct объект со свойствами, содержащими векторы того же размера, что и эти векторы столбцов. Чтобы убедиться, что измерения возвращенного решения согласуются с измерениями исходных точек запроса, используйте reshape функция. Например, использовать intrpAccel = reshape(intrpAccel.ux,size(xq)).

Типы данных: double

точки запроса координат y, заданные как вещественный массив. interpolateAcceleration вычисляет ускорения в 2-D координатных точках [xq(i),yq(i)] или в точках координат 3-D [xq(i),yq(i),zq(i)]. Поэтому xq, yq, и (при наличии) zq должно иметь одинаковое количество записей. Внутри, interpolateAcceleration преобразует точки запроса в вектор столбца yq(:).

Типы данных: double

точки запроса координат z, заданные как вещественный массив. interpolateAcceleration вычисляет ускорения в 3-D координатных точках [xq(i),yq(i),zq(i)]. Поэтому xq, yq, и zq должно иметь одинаковое количество записей. Внутри, interpolateAcceleration преобразует точки запроса в вектор столбца zq(:).

Типы данных: double

Точки запроса, заданные как вещественная матрица с двумя строками для 2-D геометрии или тремя строками для 3-D геометрии. interpolateAcceleration вычисляет ускорения в координатных точках querypoints(:,i), таким образом, каждый столбец querypoints содержит только одну 2-D или 3-D точку запроса.

Пример: Для 2-D геометрии querypoints = [0.5,0.5,0.75,0.75; 1,2,0,0.5]

Типы данных: double

Выходные аргументы

свернуть все

Ускорения в точках запроса, возвращаемые как FEStruct объект со свойствами, представляющими пространственные компоненты ускорения в точках запроса. Для точек запроса, которые находятся вне геометрии, intrpAccel прибыль NaN. Свойства FEStruct объекты доступны только для чтения.

См. также

| | | | | | | | | | | |

Представлен в R2018a

Документация по набору инструментов для дифференциальных уравнений в частных производных

Поддержка