Найдите начальное смещение и скорость присвоенными геометрической области
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')
structuralIC | Свойства GeometricStructuralICs | Свойства NodalStructuralICs | StructuralModel