Программирование МКЭ в МATLAB. Даутов P.З. - 21 стр.

§ 2. Определение геометрии области, построение P1 сеток                                                        21


        Имеется две возможности задания геометрии:
   1) при помощи матрицы (если все сегменты являются отрезками
      прямых или сегментами окружностей или эллипсов);
   2) при помощи m.файла.
1 способ. Cледующая матрица, формируемая функцией gtm, опре-
деляет геометрию области для нашего примера. Матрица содержит
15 столбцов (по числу сегментов) и каждый столбец определяет 1
сегмент. Первый элемент в столбце определяет тип сегмента (2 —
отрезок прямой, 4 — эллипса); 2, 3 (4, 5) строки определяют x (y) ко-
ординаты начальной и конечной точки соответственно; 6, 7 строки —
номер области слева и справа по направлению обхода. Для отрезков
прямых следующие строки не нужны и задаются нулями; для эллип-
са 8, 9 строки определяют x и y координаты центра эллипса, а 10, 11
строки — его x и y полуоси (для окружности они совпадают и равны
радиусу окружности); 12 строка определяет угол поворота эллипса
вокруг центра против часовой стрелки (в радианах).1)
f u n c t i o n g=gtm
% geometry matrix
g =[2          2    2 2         2         2 2       2       2       2   2        4      4        4      4
     2 1 .25 1 .25 2            2   0 .75 0         0       0   0 .75   0 0 .75         1 1 .25         1
     0 1 .25 0 .75 2            2   1 .25 0         0       2   0 .75   2        1 1 .25         1 0 .75
   −1 0 . 7 5 0 . 2 5 −1        0   0 . 7 5 −1      0       1   0 .25   0 −0 . 5 −0 . 7 5 −0 . 5 −0 . 2 5
   −1 0 . 2 5 0 . 2 5 0         1   0 .75 0         1       1   0 .75   0 −0 . 7 5 −0 . 5 −0 . 2 5 −0 . 5
     0         3    3 1         3         3 0       0       0       3   3        0      0        0      0
     1         2    2 0         0         2 1       3       3       2   1        1      1        1      1
     0         0    0 0         0         0 0       0       0       0   0        1      1        1      1
     0         0    0 0         0         0 0       0       0       0   0 −0 . 5 −0 . 5 −0 . 5 −0 . 5
     0         0    0 0         0         0 0       0       0       0   0 0 .25 0 .25 0 .25 0 .25
     0         0    0 0         0         0 0       0       0       0   0 0 .25 0 .25 0 .25 0 .25
     0         0    0 0         0         0 0       0       0       0   0        0      0        0      0 ];

Следующие команды строит сетку, представленную на рис. 3:
g=gtm ;
[ p , e , t ]= i n i t m e s h ( g , 'Hmax' , 0 . 1 5 ) ;
pdemesh ( p , e , t ) , a x i s e q u a l

2 способ. К рис. 3 приводят также команды
g=' g t f ' ;
[ p , e , t ]= i n i t m e s h ( g , 'Hmax' , 0 . 1 5 ) ;
pdemesh ( p , e , t ) , a x i s e q u a l
  1)
       если все сегменты прямые, то матрица может иметь лишь 7 строк.