Exponenta Banner
exponenta event banner

interpolateVelocity

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

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

Синтаксис

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

Описание

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

пример

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

intrpVel = interpolateVelocity(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];
intrpVel = interpolateVelocity(structuralresults,coordsMidSpan);

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

figure
plot(structuralresults.SolutionTimes,intrpVel.vy)
title('Y-Velocity of the Geometric Center of the Beam')

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

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

свернуть все

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

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

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

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

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

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

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

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

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

Точки запроса, заданные как вещественная матрица с двумя строками для 2-D геометрии или тремя строками для 3-D геометрии. interpolateVelocity оценивает скорости в координатных точках 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 объект со свойствами, представляющими пространственные компоненты скорости в точках запроса. Для точек запроса, которые находятся вне геометрии, intrpVel прибыль NaN. Свойства FEStruct объекты доступны только для чтения.

См. также

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

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

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

Поддержка