Трубопроводные системы. Точность расчета. Все ли так безоблачно?

В настоящее время, на рынке программных продуктов появилось достаточное количество расчетных программ, позволяющих проводить прочностной анализ трубопроводных систем. Перед пользователем возникает вопрос, какой продукт наиболее оптимально подходит под нужды предприятия, института, заказчика…
Стоит ли в этом вопросе опираться на всем известное соотношение цена/качество?
При всей схожести основных решаемых задач, которые ставятся перед программным комплексом, каждая из них имеет ряд отличительных особенностей. Не всегда покупатель может самостоятельно и достаточно быстро разобраться в них.
Причиной является то, что какие-то возможности программы лежат на поверхности, а о каких-то умалчивают.
Фразы о дружественном интерфейсе и простоте интерпретации результатов расчета не новы и знакомы многим. Однако мало кто задумывался о «механизме решения» программного комплекса, на каких алгоритмах он построен, насколько точны полученные результаты и какие граничные условия в нем имеются.
Однако именно это является основным показателем надежности программного продукта и гарантией адекватности полученных показателей.
Предыдущие статьи были посвящены необходимости использования модуля динамического анализа, как самого современного подхода для решения нестационарных режимов в трубопроводе. Не каждое программное обеспечение позволяет это сделать, что в свою очередь уже говорит о слабости алгоритмов, на основе которых происходит анализ.
Сегодня, мы приоткроем завесу и обратим внимание на некоторые задачи и возможности статического модуля Solver.
Трубопровод представляет собой неразрезную конструкцию, стержни которой обычно жестко соединены между собой. Для классификации таких систем существует два признака. Один из них, кинематический. По нему стержневые системы делятся на геометрически неизменяемые и геометрически изменяемые. В классическом понимании, геометрически неизменяемыми являются почти все строительные конструкции. Дополнительные вопросы появляются тогда, когда конструкция превращается в геометрически изменяемую, т.е. не предназначенную для восприятия нагрузок и передачу их на опоры. Достаточно часто отсутствие опыта расчетчика или простая человеческая ошибка, может привести к некорректности задания исходных данных в программе.
Наша задача выяснить, как алгоритм решения программы отреагирует на это. Подскажет ли об ошибке или попытается самостоятельно наложить дополнительные граничные условия, чтобы получить хоть какой-либо результат. Для проверки указанной задачи, проведем сравнение в двух программах по расчету трубопровода на прочность. Назовем их условно, программа 1 (от зарубежного разработчика) и программа 2 (отечественной разработки).
Смоделирован аппарат без мертвых опор. Первый вариант расчета: внешние нагрузки отсутствуют. Второй вариант расчета: наличие внешнего крутящего момента.


                                                   а)  программа №1                                                                                                  б) программа №2
Рисунок 1. Модель горизонтального аппарата
Модель, показанная на рисунке 1, кинематически неподвижна. При приложении к аппарату изгибающего момента, он может свободно вращаться вокруг продольной оси.

Таблица 1. Пример 1_ момент равен нулю. Перемещение узлов


Программа № 1

Программа № 2

Номер узла

Х,
мм

Y,
мм

Z,
мм

Номер узла

Х,
мм

Y,
мм

Z,
мм

RX, град

10

-5,74

0,00

0,00

10

-6,26

0,10

-0,10

15,07

20

-4,54

0,00

0,00

20

-4,96

0,04

0,04

15,07

30 опора

-3,59

0,00

0,00

30 опора

-3,91

-0,20

0,00

15,07

40 опора

4

0,00

0,00

40 опора

0,00

-0,20

0,20

15,07

50 точка приложения момента

0,84

0,00

0,00

50 точка приложения момента

0,91

-0,30

0,20

15,07

60

2,03

0,00

0,00

60

2,22

-0,20

0,25

15,07

Таблица 2. Пример 2_ момент равен 10 Н·м. Перемещение узлов


Программа № 1

Программа № 2

Номер узла

Х,
мм

Y,
мм

Z,
мм

RX, град

Номер узла

Х,
мм

Y,
мм

Z,
мм

RX, град

10

-5,74

0,00

0,00

5071

10

-6,26

0,02

-0,03

-4,53

20

-4,54

0,00

0,00

5071

20

-4,96

0,01

-0,01

-4,53

30 опора

-3,59

0,00

0,00

5071

30 опора

-3,91

0,00

-0,00

-4,53

40 опора

4

0,00

0,00

5071

40 опора

0,00

-0,05

0,05

-4,53

50 точка приложения момента

0,84

0,00

0,00

5071

50 точка приложения момента

0,91

-0,01

0,1

-4,53

60

2,03

0,00

0,00

5071

60

2,22

-0,01

0,1

-4,53

Таблица 3. Пример 3_ момент равен 20 Н·м. Перемещение узлов


Программа № 1

Программа № 2

Номер узла

Х,
мм

Y,
мм

Z,
мм

RX, град

Номер узла

Х,
мм

Y,
мм

Z,
мм

RX, град

10

-5,74

0,00

0,00

10142

10

-6,26

-0,05

0,04

6,01

20

-4,54

0,00

0,00

10142

20

-4,96

-0,02

0,02

6,01

30 опора

-3,59

0,00

0,00

10142

30 опора

-3,91

0,00

0,00

6,01

40 опора

4

0,00

0,00

10142

40 опора

0,00

0,10

-0,10

6,01

50 точка приложения момента

0,84

0,00

0,00

10142

50 точка приложения момента

0,91

0,10

-0,10

6,01

60

2,03

0,00

0,00

10142

60

2,22

0,10

-0,10

6,01

 

Модель, сделанная в программе №1, на рисунке 1а, и результаты ее расчета показывают, что в аппарате присутствуют только линейные расширения от температурного расширения, а при возникновении изгибающего момента, появляется угол поворота, который увеличивается совместно с приложенным моментом.
Другая модель, на рисунке 1б, показывает, что результаты получаются недостоверными еще на этапе первого примера (без приложенного изгибающего момента). Появление боковых линейных и угловых перемещений говорит о невозможности системы уравнений справиться с этой задачей. А при приложении изгибающего момента и его изменении, получаем непредсказуемые изменения перемещений.

Как показывают результаты расчета, во второй программе, при расчете таких систем, разрешающая система уравнений получается плохо обусловленной, в которой лишь малое изменение исходных данных приводит к непредсказуемым изменениям в результатах расчета. Это как нельзя лучше показывает ограничения и слабости «решателя».

Перейдем ко второму примеру. В нем мы так же рассмотрим проблемные/ошибочные участки: это установка осевых компенсаторов, без промежуточной неподвижной опоры, и близко расположенные тройники.

Рисунок 2. Модель трубопроводной системы в программе №1.

 

Рисунок 3. Модель трубопроводной системы в программе №2.
Рассматривается обычная система с рабочим давлением и температурой без дополнительных внешних воздействий. Результаты расчета приведены в таблицах ниже.

Таблица 4. Перемещение узлов_Рабочий режим


Программа № 1

Программа № 2

Номер узла

Х,
мм

Y,
мм

Z,
мм

Номер узла

Х,
мм

Y,
мм

Z,
мм

40 начало компенсатора

0,00

-0,27

4,61

40 точка компенсатора

0,00

-0,39

-1,76

50 окончание компенсатора

0,00

-0,27

-2,14

80 начало компенсатора

0,01

-0,32

1,56

80 точка компенсатора

0,01

-0,44

2,21

90 окончание компенсатора

0,01

-0,32

-4,31

130 тройник

0,20

0,90

0,61

130 тройник

0,17

0,88

0,82

140 тройник

0,13

1,33

0,76

140 тройник

0,05

0,97

2,31

Таблица 5. Перемещение узлов_Холодный режим


Программа № 1

Программа № 2

Номер узла

Х,
мм

Y,
мм

Z,
мм

Номер узла

Х,
мм

Y,
мм

Z,
мм

40 начало компенсатора

0,00

-0,26

-0,30

40 точка компенсатора

0,00

-0,39

0,73

50 окончание компенсатора

0,00

-0,26

0,13

80 начало компенсатора

0,01

-0,30

-0,11

80 точка компенсатора

0,01

-0,41

0,58

90 окончание компенсатора

0,01

-0,30

1,77

130 тройник

-0,28

0,60

1,45

130 тройник

-0,28

0,43

1,30

140 тройник

-0,49

0,89

1,44

140 тройник

0,49

0,66

1,30

Таблица 6. Напряжения


Программа № 1

Программа № 2

 

Холодное состояние

Рабочее состояние

 

Холодное состояние

Рабочее состояние

Номер узла

σдейст,
МПа

σдоп,
МПа

σдейст,
МПа

σдоп,
МПа

Номер узла

σдейст,
МПа

σдоп,
МПа

σдейст,
МПа

σдоп,
МПа

40 начало компенсатора

2,66

180,99

2,65

217,2

40 точка компенсатора

3,75

180,99

4,65

217,2

50 окончание компенсатора

2,66

180,99

2,66

217,2

80 начало компенсатора

2,82

180,99

2,88

217,2

80 точка компенсатора

3,73

180,99

4,91

217,2

90 окончание компенсатора

2,86

180,99

2,94

217,2

130 тройник

10,92

180,99

19,59

217,2

130 тройник

27,03

180,99

-

217,2

140 тройник

14,21

180,99

23,91

217,2

140 тройник

32,79

180,99

-

217,2

 

Результаты расчета так же показывают отсутствие полной сходимости результатов, особенно в плане перемещений. В программе №1 компенсатор моделируется двумя точками, что позволяет провести более точный расчет перемещений в начальном и конечном узле, для определения общей работы элемента.
Кроме этого, граничные условия решателя, накладываемые на другую расчетную модель, в плане мгновенно изменяемой системы (в пределах податливости), создают перемещения конструкции в противоположном направлении от перемещений, рассчитанных с помощью программы №1.

Не смотря на то, что, как уже говорилось ранее, модели являются тестовыми, а не рабочими, и они заведомо смоделированы с ошибками,  именно такие примеры позволяют как можно глубже узнать отличия и ограничения решателя программы, которые во второй программе присутствуют налицо.
У читателей может возникнуть логичный вопрос, а почему считается, что Solver программы №1 лучше и именно она проводит расчет грамотно, а не другое ПО. Ответ на этот вопрос лежит в логике полученных результатов расчета и верификационной проверке.

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

Не нашли нужную информацию?

Позвоните по телефону +7 (495) 744 00 11 или напишите нам на e-mail:  Адрес электронной почты защищен от спам-ботов. Для просмотра адреса в вашем браузере должен быть включен Javascript.
Наши специалисты ответят на ваши вопросы!


Яндекс.Метрика