1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
| unit Main;
interface
uses
Formulas, Variables, Input_Function;
{…}
var
Dsn:mas;
matr_A: mas2; {основная матрица СЛАУ}
stolb_b:mas; {столбец правой части СЛАУ}
implementation
procedure border; {разбивка на граничные элементы}
var
i:integer;
begin
Setlength(x,N+1);
Setlength(y,N+1);
Setlength(a,N+1);
Setlength(beta,N+1);
for i:=1 to M do
begin
a[i]:=AO/M/2;
beta[i]:=0;
x[i]:=-AO+a[i]+2*a[i]*(i-1);
y[i]:=0;
end;
for i:=M+1 to M+n_OB do
begin
a[i]:=L/n_OB/2;
beta[i]:=0;
x[i]:=a[i]+2*a[i]*(i-M-1);
y[i]:=0;
end;
for i:=M+n_OB+1 to N do
begin
a[i]:=CD/n_CD/2;
beta[i]:=alfa;
x[i]:=L+BH+(-CD/2+a[i]+2*a[i]*(i-M-n_OB-1))*cos(beta[i]);
y[i]:=(-CD/2+a[i]+2*a[i]*(i-M-n_OB-1))*sin(beta[i]);
end;
end;
//___________
procedure base; {формирование и решение СЛАУ}
var
i,j:integer;
begin
{формирование СЛАУ в зависимости от заданных граничных условий}
{задаем размеры массивов, нумерация элементов начинается с 0, нулевые элементы не исполь-зуем}
Setlength(matr_A,2*N+1,2*N+1); {размер основной матрицы 2N*2N}
Setlength(stolb_b,2*N+1); {размер правой части 2N}
Setlength(Dsn,2*N+1); {размер вектора неизвестных 2N}
{формирование основной матрицы СЛАУ}
for i:=1 to N do
begin
for j:=1 to N do
begin
{
на участке АО
}
if check1 then
if i<=M then
begin
matr_A[i,j]:=Bss(i,j); {первое уравнение СЛАУ}
matr_A[i,j+N]:=Bsn(i,j)
end;
if check2 then
if i<=M then
begin
matr_A[i,j]:=Ass(i,j); {второе уравнение СЛАУ}
matr_A[i,j+N]:=Asn(i,j);
end;
if (check3)or(check10) then
if i<=M then
begin
matr_A[i+N,j]:=Bns(i,j); {третье уравнение СЛАУ}
matr_A[i+N,j+N]:=Bnn(i,j);
end;
if check4 then
if i<=M then
begin
matr_A[i+N,j]:=Ans(i,j); {четвертое уравнение СЛАУ}
matr_A[i+N,j+N]:=Ann(i,j);
end;
{
на участке ОВ
}
if check5 then
if (i>M)and(i<=M+n_OB) then
begin
matr_A[i,j]:=Bss(i,j); {первое уравнение СЛАУ}
matr_A[i,j+N]:=Bsn(i,j)
end;
if check6 then
if (i>M)and(i<=M+n_OB) then
begin
matr_A[i,j]:=Ass(i,j); {второе уравнение СЛАУ}
matr_A[i,j+N]:=Asn(i,j);
end;
if check7 then
if (i>M)and(i<=M+n_OB) then
begin
matr_A[i+N,j]:=Bns(i,j); {третье уравнение СЛАУ}
matr_A[i+N,j+N]:=Bnn(i,j);
end;
if check8 then
if (i>M)and(i<=M+n_OB) then
begin
matr_A[i+N,j]:=Ans(i,j); {четвертое уравнение СЛАУ}
matr_A[i+N,j+N]:=Ann(i,j);
end;
//________ |
(
