Поиск начального смещения и скорости, назначенных геометрической области
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')
Свойства GeometricStructuralIc | Свойства NodalStructuralICs | structuralIC | StructuralModel