Exponenta Banner

interpolateVelocity

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

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

Синтаксис

intrpVel = interpolateVelocity(structuralresults,xq,yq)
intrpVel = interpolateVelocity(structuralresults,xq,yq,zq)
intrpVel = interpolateVelocity(structuralresults,querypoints)

Описание

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.

Примеры

свернуть все

Открыть Live Script

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

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

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

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

Partial Differential Equation Toolbox производными

Поддержка