В этом примере показано, как создать 3-D геометрии, образованные одной или несколькими кубическими, цилиндрическими и сферическими ячейками с помощью multicuboid, multicylinder, и multisphere соответственно. С помощью этих функций можно создавать сложенные или вложенные геометрии. Можно также создать геометрию, в которой некоторые ячейки пусты; например, полые цилиндры, кубы или сферы.
Все ячейки в геометрии должны быть одного типа: кубоиды, цилиндры или сферы. Эти функции не объединяют ячейки различных типов в одной геометрии.
Создайте геометрию, состоящую из одной сферы, и включите ее в модель PDE.
Используйте multisphere для создания единой сферы. Результирующая геометрия состоит из одной ячейки.
gm = multisphere(5)
gm =
DiscreteGeometry with properties:
NumCells: 1
NumFaces: 1
NumEdges: 0
NumVertices: 0
Vertices: []
Создайте модель 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
Vertices: [24x3 double]
Создайте модель 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
Vertices: [5x3 double]
Создайте модель 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
Vertices: [4x3 double]
Создайте модель 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)
