Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Алгебра нечётких чисел и её применение в инженерной практике. Препринт № 1-2013

.pdf
Скачиваний:
0
Добавлен:
11.08.2026
Размер:
643 Кб
Скачать

 

Å

Å

Å

Å

 

Свойство 3.

A( A(xɶ1 , xɶ2 ), xɶ3 )

A(xɶ1 , A(xɶ2 , xɶ3 ))

сочетательный

(ассоциативный) закон сложения.

 

 

 

 

 

 

 

 

Свойство 4.

M (M (xɶ1 , xɶ2 ), xɶ3 ) M (xɶ1 , M (xɶ2 , xɶ3 ))

– сочетатель-

ный (ассоциативный) закон умножения.

 

 

Å

Å

 

Свойство 5.

M (xɶ1 , A(xɶ2 , xɶ3 )) A(M (xɶ1 , xɶ2 ), M (xɶ1 , xɶ3 ))

Å

 

 

Å

 

A(M (xɶ2 , xɶ1 ), M (xɶ3 , xɶ1 ))

M ( A(xɶ2 , xɶ3 ), xɶ1 )

– распределительный

(дистрибутивный) закон умножения относительно сложения.

 

 

Å

Å

 

Свойство 6.

M (xɶ1 , A(xɶ2 , xɶ3 ))

A(M (xɶ1 , xɶ2 ), M (xɶ1 , xɶ3 ))

Å

 

Å

 

 

A(M (xɶ2 , xɶ1 ), M (xɶ3 , xɶ1 )) M ( A(xɶ2 , xɶ3 ), xɶ1 )

– распределитель-

ный (дистрибутивный) закон умножения относительно вычитания.

 

 

Å

 

 

 

 

), xɶ-1 )

Å

 

 

 

 

 

, xɶ-1 ))

Свойство 7.

M ( A(xɶ

 

, xɶ

 

A(M (xɶ

 

, xɶ

-1 ), M (xɶ

 

 

 

1

2

3

 

1

3

2

3

Å

 

 

 

 

 

 

 

 

Å

 

 

 

 

 

A(M (xɶ-1

, xɶ

), M (xɶ

-1

, xɶ

)) M (xɶ

-1

, A(xɶ

, xɶ

))

– распределитель-

3

1

3

 

2

 

3

1

 

2

 

 

 

ный (дистрибутивный) закон умножения относительно сложения.

 

 

Å

 

 

xɶ

 

Å

 

 

 

 

, xɶ-1 ))

Свойство 8.

M ( A(xɶ

 

,

), xɶ

-1 ) A(M (xɶ

, xɶ-1 ), M (xɶ

 

 

2

 

3

1

2

1

 

3

1

Å

 

 

 

 

 

 

 

Å

, xɶ

 

 

 

A(M (xɶ-1

, xɶ ), M (xɶ

-1

, xɶ

)) M (xɶ-1

, A(xɶ

))

– распредели-

1

3

1

 

2

 

1

2

3

 

 

 

тельный (дистрибутивный) закон умножения относительно вычитания.

51

Автоматическая проверка справедливости этих свойств осуществлена экспериментально с использованием разработанных и реализованных численно-символьных алгоритмов.

Целесообразно отобразить следующие записи математических выражений, которые возможно доказать методом индукции:

x1

; x1

 

Å

, x1 ) = 2x1;

x1 +¼+ x1

= nx1 (n Î ).

+ x1 ¬ A(x1

ɶ

ɶ

ɶ

ɶ

ɶ

ɶ

ɶ

 

ɶ

 

ɶ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

ɶ2

 

 

 

ɶn

(n Î )

 

ɶ

ɶ ɶ

 

ɶ ɶ

ɶ

ɶ

 

 

x ; x ·x

¬ M (x ,x ) = x

; x ·¼·x

= x

 

1

1 1

 

1 1

1

1

1

 

1

 

n

Имеют место и справедливы также следующие операции с нечёткими числами:

Операция 1.

 

ɶα

 

β

 

 

ɶα +β

.

 

 

 

 

 

 

ɶ

 

 

 

 

 

 

 

M (x

, x

) @ x

 

 

 

 

 

 

Операция 2.

ɶα

)

β

 

ɶα·β

.

 

 

 

 

 

 

 

 

(x

 

= x

 

 

 

 

 

 

 

 

 

Операция 3.

 

 

ɶ

ɶ

 

p

 

 

 

ɶp

ɶp

).

 

M (x, y)

= M (x

, y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Операция 4.

 

 

ɶ

ɶ

-1

p

 

ɶp

ɶ

- p

).

M (x, y

 

 

)

= M (x

, y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Операция 5.

 

ɶp

 

-q

 

ɶp-q

.

 

 

 

 

 

 

ɶ

 

 

 

 

 

 

M (x

, x

 

 

) = x

 

 

 

 

 

 

 

4 Решение задач строительства в нечёткой постановке

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

52

железобетонного образца ζɶ и на оценке момента трещинообразо-

вания Mɶcrc как нечёткой величины.

В работе [4] предложен новый метод построения нечётких аппроксимаций диаграмм деформирования материалов, в частности бетона. Нечёткие аппроксимации применяются к решению задачи о возникновении, стабилизации и росте трещин в нечёткой постановке. Методами кластеризации определены интерполяционные узлы аппроксимирующих функций как нечёткие величины, затем построены «функция-мажоранта» и «функция-миноранта». Разработан итерационный метод определения узлов интерполяции как нечётких величин. В качестве исходного материала использованы результаты экспериментальных исследований деформирования бетонных и железобетонных образцов – « восьмёрок» при центральном растяжении балок и при четырёхточечном изгибе.

Работа [5] посвящена построению и реализации алгоритма оп-

ределения изгибающего момента Mɶcrc , соответствующего возник-

новению трещины в железобетонном элементе в нечёткой постановке. Исходные данные для решения задачи представляются в виде нечётких величин, моды и носители которых определяются по результатам экспериментов.

Формулы для вычисления Mɶcrc и ζɶ (нечёткой величины) имеют вид:

 

 

 

 

 

ζɶ =

 

μɶ

αɶ

s

+

 

 

μɶ

2αɶ2

+

 

μɶ αɶ

 

 

+

 

Eɶ

 

Wɶ

Eɶ

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

s

s

 

 

 

 

s

 

·

 

rt

 

 

s

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɶ

 

ɶ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

4

 

 

 

 

 

 

s

 

s

 

 

 

 

 

 

2

εɶ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

E

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

bt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɶ

 

 

 

εɶbt

 

Bɶ

 

 

 

Aɶεɶbt

 

 

Eɶbεɶbt

 

 

ζɶ3

 

 

 

 

 

 

 

1ɶ

ɶ

ɶ

ɶ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M crc

=

 

 

 

 

ɶ

+

 

 

ɶ

 

+

 

 

 

 

 

 

 

 

 

 

 

 

ɶ

+

ɶ

 

 

ɶ

;

ɶ

ɶ

2

 

 

 

ɶ

 

 

 

ɶ

 

 

2

 

 

 

 

μs Esεbt

 

 

 

 

 

 

3

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

1

 

 

 

 

 

 

 

 

 

(1

)

 

 

 

 

 

 

 

 

 

 

 

(1

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɶ

 

 

 

ɶ

 

 

 

·σɶ

 

 

 

ɶ

ɶ2

.

 

 

 

 

 

 

 

 

 

 

 

(4.1)

 

 

 

 

 

 

 

 

 

 

 

M

crc

= M

crc

02

·b ·h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

53

В расчётах использованы нечёткие величины и их дефаззифицированные значения методом центра тяжести, полученные в работе [4]:

Rɶbt = [0, 7197; 1, 0092]0 [0,82; 0,82]1;

Wɶ

 

=

 

 

 

 

 

 

 

−5

 

 

 

−5

 

 

 

 

 

 

−5

 

 

−5

;

 

4, 3986·10

; 5, 721·10

5, 4797·10

; 5, 4797·10

 

rt

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

1

 

 

εɶ

 

 

 

 

 

 

 

−5

 

 

−5

 

 

 

 

 

−5

 

−5

;

 

 

 

 

 

 

 

= 7, 58·10

; 8,15·10

7, 95·10

; 7, 95·10

 

 

 

 

 

 

 

bt

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

1

 

 

 

 

 

Eɶ

=

 

 

 

 

 

 

4

 

 

 

4

 

 

 

 

 

 

 

4

 

4

;

 

2, 4865·10

; 2, 7692·10

 

2, 7692·10

; 2, 7692·10

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

1

 

Aɶ =

 

 

 

 

 

 

 

 

8

 

 

 

 

8

 

 

 

 

 

 

8

 

 

 

 

 

8

 

−4, 3833·10 ; − 2, 4814·10

− 2, 9757·10

;− 2, 9757·10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

1

Bɶ =

 

 

 

 

 

 

4

 

 

 

4

 

 

 

 

 

 

4

 

4

 

 

;

2, 6726·10

; 4, 2064·10

3,1242·10 ; 3,1242·10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

R

 

= Rɶd

= (0,85, 0, 33);

 

 

Eb = Eɶbd

= (2, 675·104 , 0, 33);

 

 

 

 

bt

 

rt

bt

(

 

 

 

 

)

 

 

 

 

(

 

 

 

)

 

 

 

 

 

 

rt

 

 

 

 

 

−5

 

 

 

ɶd

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɶ d

 

 

 

 

, 0, 33

 

;

A =

=

 

 

 

8

 

 

 

;

 

 

W =W

 

= 5, 2·10

 

A

(

−3, 28·10 , 0, 33

 

 

 

 

bt

 

 

bt

 

(

 

 

 

−5

 

)

 

ɶd

 

 

 

 

)

 

 

 

 

 

ε

 

= εɶ

d

 

 

 

 

 

, 0, 33 ;

 

=

 

 

 

4

 

 

 

 

 

 

 

 

 

= 7,893·10

 

B = B

 

3, 334·10 , 0, 33 .

 

В таблицах 4 и 5 и на рисунках 3 и 4 приводятся результаты определения параметров ζɶ и Mɶcrc и их дефаззифицированных значений ζ и M crc (методом центра тяжести [2]) для балки с размерами: b = 0,1 м ; h0 = 0,16 м ; l =1, 2 м .

В таблице 4 приведены значения нечётких величин ζɶ и Mɶcrc ,

полученных с применением алгебры нечётких чисел по формулам

(4.1), где исходными данными являлись нечёткие значения Rɶbt ,

Wɶrt , εɶbt , Eɶb , Aɶ , Bɶ . Расчёт выполнен для характерных коэффици-

ентов армирования.

54

Таблица 4

– Результаты расчётов для различных

значений

μs

с применением алгебры нечётких чисел

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μs

 

ζɶ

 

 

 

 

 

 

Mɶcrc

 

 

 

 

 

 

 

ζɶ1

= [0, 489; 0, 633]0

 

 

 

 

 

Mɶcrc,1

=

 

 

 

 

 

0, 000

 

 

 

−3

; 2, 380 ×10

−3

 

 

 

[0, 560; 0, 560]1

= 1, 061×10

 

 

 

0

 

 

 

 

×10

−3

; 1, 502

×10

−3

 

 

 

 

 

 

 

 

 

1, 502

 

 

1

 

ζɶ2

= [0, 602; 0, 749]0

 

 

 

 

 

Mɶcrc,2

=

 

 

 

 

 

0, 005

 

 

 

 

−3

; 2, 312 ×10

−3

 

 

 

[0, 668; 0, 668]1

= 1, 089 ×10

 

 

0

 

 

 

 

×10

−3

; 1, 518

×10

−3

 

 

 

 

 

 

 

 

 

 

1, 518

 

 

1

 

ζɶ3

= [0, 721; 0,873]0

 

 

 

 

 

Mɶcrc,3 =

 

 

 

 

 

0, 010

 

 

 

−3

; 2, 292 ×10

−3

 

 

 

[0, 782; 0, 782]1

= 1,132 ×10

 

 

 

0

 

 

 

 

×10

−3

; 1, 558

×10

−3

 

 

 

 

 

 

 

 

 

1, 558

 

 

1

В первой колонке таблицы 5 приводятся дефаззифицирован-

ные значения (методом центра тяжести [4]) нечётких величин ζɶ и

Mɶcrc , найденных с использованием алгебры нечётких чисел. Де-

фаззификация проводилась на последнем этапе после определения нечётких величин ζɶ и Mɶcrc при разных коэффициентах армирова-

ния. Во второй колонке таблицы 5 приведены результаты решения задачи в чёткой постановке (определены величины ζ и M при разных коэффициентах армирования), где в качестве исходных

55

данных приняты дефаззифицированные значения нечётких величин

Rɶ

, Wɶ

, εɶ

, Eɶ

, Aɶ , Bɶ .

 

 

 

 

 

 

 

bt

rt

bt

 

b

 

 

 

 

 

 

 

 

 

Таблица 5 –

Сравнение величин ζ и

M после дефаззификации

на различных этапах вычислений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Результаты решения задачи

 

Дефаззификация нечётких

 

в чёткой постановке

величин ζɶ

и Mɶ , найденных

при исходных данных −

 

с использованием алгебры

дефаззифицированных

 

значениях нечётких величин

 

 

нечётких чисел

 

 

Rɶ

 

, Eɶ

, εɶ , Wɶ

, Aɶ , Bɶ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

bt

b

 

bt

rt

 

ζ1 = 0, 560

 

M1 = 1, 648 ×10−3

ζ1 = 0, 559

 

M1 = 1, 592 ×10−3

ζ2 = 0, 673

 

M 2

= 1, 640 ×10−3

ζ2 = 0, 610

 

M 2

= 1, 654 ×10−3

ζ3 = 0, 792

 

M 3

= 1, 661×10−3

ζ3 = 0, 660

 

M 3

= 1, 721×10−3

Примечание:

 

 

ζ1 = ζ (μ = 0),

 

 

ζ2 =ζ (μ = 0, 005),

ζ3 = ζ (μ = 0, 01);

M1 = M = 0),

 

M 2 = M = 0, 005),

M3 = M = 0, 01).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчётные функции принадлежности дают вполне приемлемые для практического использования в инженерной практике оценки момента трещинообразования. Интересным для практического применения является отношение моды нечёткой величины Mɶcrc

(значение момента, для которого функция принадлежности равна 1) к значению момента, соответствующего левой границе носите-

ля Mɶcrc (см. рисунок 4). Это отношение для всех значений коэф-

фициента армирования приближенно можно считать равным 1, 4 ,

соответственно обратное отношение ≈ 0,72. Полученные резуль-

56

таты доказывают возможность установить единый коэффициент запаса по моменту образования трещины.

Рисунок 3 – Кортеж функций принадлежности нечёткой величины ζɶ и их дефаззифицированные значения

57

Рисунок 4 – Кортеж функций принадлежности нечёткой величины Mɶ и их дефаззифицированные значения

58

ЗАКЛЮЧЕНИЕ

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

Показан результат работы алгоритма определения изгибающего момента, соответствующего возникновению трещины в железобетонном элементе. Исходные данные для решения задачи представляются в виде нечётких величин, моды и носители которых определяются по результатам экспериментов.

Расчётные функции принадлежности дают вполне приемлемые для практического применения оценки момента трещинообразования.

59

СПИСОК ЛИТЕРАТУРЫ

1.Stefanini L. Fuzzy Numbers and Fuzzy Arithmetic / L. Stefanini, L. Sorini, M. L. Guerra // Handbook of granular computing / eds. W. Pedrycz, A. Skowron, V. Kreinovich. – A John Wil ey & Sons, Ltd, Publication, 2008. – Р. 249–285.

2.Штовба С. Д. Проектирование нечётких систем средствами MATLAB / С. Д. Штовба. – М. : Горячая линия– Телеком, 2007. – 288 с.

3.RILEM TС QFS ‘Quasibrittle fracture scaling and size effect’ – Final report // Materials and Structures. – October 2 004. – Vol. 37. – Р. 547–568.

4.Адищев В. В. Применение кластерного анализа для построения нечётких аппроксимаций диаграмм деформирования бето-

на / В. В. Адищев, Д. С. Шмаков // Известия вузов. Строитель-

ство. – 2012. – № 4. – C. 60–71.

5.Адищев В. В. Определение момента образования трещины как нечёткой величины по экспериментальным данным / В. В. Адищев, Д. С. Шмаков // Известия вузов. Строительство. – 2012. –

№ 5. – C. 3–14.

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]