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