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

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

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

Одна сфера

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

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

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

       NumCells: 1
       NumFaces: 1
       NumEdges: 0
    NumVertices: 0

Создайте модель PDE.

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')

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

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

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

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

       NumCells: 3
       NumFaces: 18
       NumEdges: 36
    NumVertices: 24

Создайте модель PDE.

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)

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

Создайте геометрию, которая состоит из трех сложенных цилиндров, и включайте эту геометрию в модель 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

Создайте модель PDE.

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)

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

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

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

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

       NumCells: 1
       NumFaces: 4
       NumEdges: 4
    NumVertices: 4

Создайте модель PDE.

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)