findStructuralIC

Найдите начальное смещение и скорость присвоенными геометрической области

Описание

пример

sica = findStructuralIC(structuralmodel.InitialConditions,RegionType,RegionID) возвращает начальное смещение и скорость, присвоенную заданной области.

Примеры

свернуть все

Найдите начальное смещение и скорость присвоенными ячейкам 3-D геометрии.

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

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

Создайте геометрию, состоящую из трех вложенных цилиндров, и включайте ее в модель. Постройте геометрию.

gm = multicylinder([5 10 15],2);
structuralmodel = createpde('structural','transient-solid');
structuralmodel.Geometry = gm;
pdegplot(structuralmodel,'CellLabels','on','FaceAlpha',0.5)

Установите начальные условия для каждой ячейки. Когда вы задаете только начальную скорость или начальное смещение, structuralIC принимает, что не использованный параметр является нулем.

structuralIC(structuralmodel,'Displacement',[0;0;0],'Velocity',[0;0;0],'Cell',1);
structuralIC(structuralmodel,'Displacement',[0;0.1;0],'Cell',2);
structuralIC(structuralmodel,'Velocity',[0;0.2;0],'Cell',3);

Проверяйте начальную спецификацию условия на ячейку 1.

SICACell1 = findStructuralIC(structuralmodel.InitialConditions,'Cell',1)
SICACell1 = 
  GeometricStructuralICs with properties:

             RegionType: 'Cell'
               RegionID: 1
    InitialDisplacement: [3x1 double]
        InitialVelocity: [3x1 double]

SICACell1.InitialDisplacement
ans = 3×1

     0
     0
     0

SICACell1.InitialVelocity
ans = 3×1

     0
     0
     0

Используйте статическое решение в качестве начального условия для динамической структурной модели. Проверяйте и постройте начальное смещение.

Создайте статическую модель.

staticmodel = createpde('structural','static-solid');

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

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

Задайте модуль Молодежи, отношение Пуассона и массовую плотность.

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

Примените граничное условие и статическую загрузку.

structuralBC(staticmodel,'Face',5,'Constraint','fixed');
structuralBoundaryLoad(staticmodel,'Face',3,'SurfaceTraction',[0;1E6;0]);
generateMesh(staticmodel,'Hmax',0.02);
Rstatic = solve(staticmodel);

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

dynamicmodel = createpde('structural','transient-solid');
gm = multicuboid(0.06,0.005,0.01);
dynamicmodel.Geometry = gm;

Примените граничное условие.

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

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

generateMesh(dynamicmodel,'Hmax',0.02);
structuralIC(dynamicmodel,Rstatic)
ans = 
  NodalStructuralICs with properties:

    InitialDisplacement: [113x3 double]
        InitialVelocity: [113x3 double]

Проверяйте начальную спецификацию условия на dynamicmodel.

sica = findStructuralIC(dynamicmodel.InitialConditions,'Cell',1)
sica = 
  NodalStructuralICs with properties:

    InitialDisplacement: [113x3 double]
        InitialVelocity: [113x3 double]

Постройте z- компонент начального смещения.

pdeplot3D(dynamicmodel,'ColorMapData',sica.InitialDisplacement(:,3))
title('Initial Displacement in the Z-direction')

Входные параметры

свернуть все

Начальные условия переходной структурной модели, заданной как InitialConditions свойство StructuralModel объект.

Геометрический тип области, заданный как 'Face', 'Edge', или 'Vertex' для 2D или 3-D модели модели или 'Cell' для 3-D модели.

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

Геометрический ID области, заданный как вектор положительных целых чисел. Найдите идентификаторы области при помощи pdegplot.

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

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

свернуть все

Структурное начальное условие для конкретной области, возвращенной как GeometricStructuralICs или NodalStructuralICs объект. Для получения дополнительной информации смотрите GeometricStructuralICs Properties и NodalStructuralICs Properties.

Введенный в R2018a

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