Кубоиды, цилиндры и сферы

В этом примере показано, как создать 3-D геометрии, образованные одной или несколькими кубическими, цилиндрическими и сферическими камерами, с помощью multicuboid, multicylinder, и multisphere функций, соответственно. С помощью этих функций можно создать сложенные или вложенные геометрии. Можно также создать геометрии, в которых некоторые камеры пусты; для примера, полых цилиндров, кубов или сфер.

Все камеры в геометрии должны быть одного типа: или кубоиды, или цилиндры, или сферы. Эти функции не объединяют камер разных типов в одной геометрии.

Единая сфера

Создайте геометрию, которая состоит из одной сферы и включите эту геометрию в модель PDE.

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

gm = multisphere(5)
gm = 
  DiscreteGeometry with properties:

       NumCells: 1
       NumFaces: 1
       NumEdges: 0
    NumVertices: 0
       Vertices: []

Создайте модель УЧП.

model = createpde
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: []
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Включите геометрию в модель.

model.Geometry = gm
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: [1x1 DiscreteGeometry]
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Постройте график геометрии.

pdegplot(model,'CellLabels','on')

Figure contains an axes. The axes contains 2 objects of type quiver, patch.

Вложенные кубоиды той же высоты

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

Создайте геометрию при помощи multicuboid функция. Получившаяся геометрия состоит из трёх камер.

gm = multicuboid([2 3 5],[4 6 10],3)
gm = 
  DiscreteGeometry with properties:

       NumCells: 3
       NumFaces: 18
       NumEdges: 36
    NumVertices: 24
       Vertices: [24x3 double]

Создайте модель УЧП.

model = createpde
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: []
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Включите геометрию в модель.

model.Geometry = gm
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: [1x1 DiscreteGeometry]
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Постройте график геометрии.

pdegplot(model,'CellLabels','on','FaceAlpha',0.5)

Figure contains an axes. The axes contains 3 objects of type quiver, patch, line.

Сложенные цилиндры

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

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

gm = multicylinder(10,[1 2 3 4],'ZOffset',[0 1 3 6])
gm = 
  DiscreteGeometry with properties:

       NumCells: 4
       NumFaces: 9
       NumEdges: 5
    NumVertices: 5
       Vertices: [5x3 double]

Создайте модель УЧП.

model = createpde
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: []
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Включите геометрию в модель.

model.Geometry = gm
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: [1x1 DiscreteGeometry]
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Постройте график геометрии.

pdegplot(model,'CellLabels','on','FaceAlpha',0.5)

Figure contains an axes. The axes contains 3 objects of type quiver, patch, line.

Полый цилиндр

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

Создайте полый цилиндр при помощи multicylinder функция со Void аргумент. Получившаяся геометрия состоит из одной камеры.

gm = multicylinder([9 10],10,'Void',[true,false])
gm = 
  DiscreteGeometry with properties:

       NumCells: 1
       NumFaces: 4
       NumEdges: 4
    NumVertices: 4
       Vertices: [4x3 double]

Создайте модель УЧП.

model = createpde
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: []
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Включите геометрию в модель.

model.Geometry = gm
model = 
  PDEModel with properties:

           PDESystemSize: 1
         IsTimeDependent: 0
                Geometry: [1x1 DiscreteGeometry]
    EquationCoefficients: []
      BoundaryConditions: []
       InitialConditions: []
                    Mesh: []
           SolverOptions: [1x1 pde.PDESolverOptions]

Постройте график геометрии.

pdegplot(model,'CellLabels','on','FaceAlpha',0.5)

Figure contains an axes. The axes contains 3 objects of type quiver, patch, line.