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

книги / Несущая способность конструкций в условиях теплосмен

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
4.03 Mб
Скачать

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

0

 

1

 

2

 

3

 

4

 

5

 

6

 

7

 

8

 

9

 

10

 

20

 

30

0,1 0,9

Рис. 1.3. Числовая ось и распределение вещественных чисел, представимых на условном одноразрядном компьютере

Поскольку значения всех вещественных чисел ограничены приведенными выше значениями ω и Ω, можно сделать вывод, что на числовой оси вблизи нулевого значения имеется область чисел, которые невозможно представить в ЭВМ в нормализованном виде

(рис. 1.4).

0 ω

Рис. 1.4. Область [0, ω] чисел, которые невозможно хранить

вЭВМ в нормализованном виде

1.3.4.Погрешности округления чисел в ЭВМ

Оценка погрешности результата вычисления арифметического выражения при наличии ошибки в представлении данных может быть выполнена следующим образом. Пусть x – точное (как правило, неизвестное) значение, а x – приближенное значение некоторой величины. Абсолютная погрешность δ представления x приближенным значением x определяется разностью

δ = x x.

Поскольку точное значение x неизвестно, вводится «верхняя» оценка для погрешности δ ≤ . Определяется величина относи-

тельной погрешности

21

ε= δ x = (x x)x = 1 x x.

Вприведенном выражении абсолютная погрешность делится на приближенное значение x , поскольку точное значение величи-

ны x неизвестно.

Поскольку не все вещественные числа можно представить с использованием для записи мантиссы конечного числа двоичных разрядов, хранящиеся в ЭВМ данные являются «округленными». В рассматриваемом случае округлением называется операция замены заданного числа x другим числом x , первые S значащих цифр которого совпадают с соответствующими цифрами исходного числа x,

аначиная с S + 1 разряда содержат нули.

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

Например, если x = 0,123 456 789, тогда при сохранении S = 7 значащих (десятичных) цифр округленное число примет значение x = 0,123 456 700. В этом случае модуль погрешности операции ок-

ругления δ = x x = 0, 000 000 089 < 0,000 000100 = 107 , то есть не

превышает единицы (с соответствующим порядком) в младшем сохраненном разряде округленного числа. В рассмотренном при-

мере «верхняя» оценка погрешности округления = 107 = 10S . Для определенности, упрощения и однозначности дальнейше-

го изложения способа представления чисел в ЭВМ принимается p = 10 , 0,1 ≤ a < 1,0 . При таком предположении в представлении

чисел на персональных компьютерах IBM достоверными оказываются 6–7 значащих цифр10 (для хранения вещественного числа используются 4 байта оперативной памяти).

10 При использовании вещественных переменных длиной 8 байт достоверными считаются 15 цифр; 19 цифр – при использовании вещественных переменных длиной 10 байт.

22

При S = 7 округленное число представляется в виде

±0, XXXXXXX 10b , где под символом X понимается любая цифра от 0 до 9. Модуль абсолютной погрешности определяется значением

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

=

 

x x

 

≤ = 107 10b = 10bS ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

модуль относительной погрешности

 

ε

 

=

 

x x

 

 

 

x

 

10bS 0, XXXXXXX 10b 10bS 0,1 10b = 101S .

 

 

 

 

 

 

Для некоторых частных случаев относительная погрешность ε представления вещественных чисел оценивается значениями:

S = 7, ε ≤ 106 = 0, 0001 %,

S = 15,

 

 

 

ε

 

1014 = 0, 000 000 000 001 %,

 

 

 

 

S = 19,

 

ε

 

10

18 = 0,000 000 000 000 000 1 %.

 

 

Необходимо подчеркнуть, что для всех чисел, представимых в электронно-вычислительной машине, относительная погреш-

ность одна и та же, что очень существенно для получения оценок погрешностей математических моделей.

1.3.5. «Потеря порядка» и «переполнение» при выполнении арифметических операций

Используются ранее принятые обозначения: ω – наименьшее положительное и Ω – наибольшее положительное числа, представимые в оперативной памяти вычислительного устройства11.

Действительные числа из интервала −ω < x < ω нельзя представить в ЭВМ в нормализованном виде («потеря порядка»). В этом случае приближенное значение числа x в компьютере – x = 0. Но это означает, что имеет место катастрофическая потеря точности арифметических вычислений:

 

ε = 1x x → ± ∞, x 0.

 

 

 

11 ω ≈ 2,9 1039

при использовании вещественных переменных длиной

4 байта, и ω ≈ 5,56 10309

при8-байтовом представлениивещественнойпеременной.

 

23

Все действительные числа x > Ω также нельзя представить

в нормализованной форме («переполнение»). Приближение числа x определяется как x = Ω . Очевидно, что и в этом случае имеет место потеря точности:

ε = 1x x → ± ∞,

x = Ω, x → ± ∞.

1.3.6. Оценка погрешности вычислительного эксперимента

Рассматривается задача об оценке погрешности вычисления

1000

произведения ai . В соответствии с определением относительной

i=1

погрешности приближенное значение переменной можно представить в виде x = x(1 − ε ) . Приближенный (округленный) результат

произведения двух первых чисел в указанном выражении может быть оценен формулой (a1a2 )(1 − ε ) . Для получения значения всего

выражения следует многократно повторить эту операцию:

a1a2

 

 

 

 

a3

 

a4

 

a1000

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1− ε

1− ε

 

− ε

 

1− ε

 

Для ε = 106

значение выражения

1000

=ai (1− ε)999 .

i=1

(1 − ε )999 равно 1,000 999 5,

что соответствует относительной погрешности вычисления всего результата, то есть перемножения 1000 сомножителей, равной 0,099 85 %. Несложно убедиться, что при перемножении 1 000 000 сомножителей относительная погрешность, рассчитанная аналогичным образом, превышает 63 %. На практике такого катастрофического нарастания погрешности, как правило, не наблюдается в силу того, что погрешности округления имеют разные знаки и, следовательно, компенсируют друг друга.

1.4. Верификация и проверка адекватности модели

Вычислительная модель исследуемого объекта, процесса или явления должна проходить верификацию, а именно:

– оценку погрешности аппроксимации разностной схемой дифференциальных уравнений;

24

устойчивость получаемого численного решения относительно возмущения входных данных (в том числе шагов дифференцирования и интегрирования);

сходимость последовательности численных решений, полученных с различными шагами интегрирования и дифференцирования;

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

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

Пример. На материальную точку с массой m, находящуюся на гладкой горизонтальной поверхности, действует горизонтальная гармоническая сила F (t ) = Asin t. Найти закон движения точки, если

ее начальная скорость v0.

Уравнение, описывающее движение точки, имеет вид

md 22x = Asin t. dt

Вкачестве начальных условий примем

x(0) = 0, v(0) = dxdt(0) = v0 .

Интегрирование исходного уравнения с заданными начальными условиями дает решение

v(t ) = v0 + mA (1cost );

x(t ) = v0 + mA t mA sin t,

25

непрерывно зависящее от начальных условий

x (0) = 0 ,

v (0) = v0 и

коэффициентов A, m.

 

 

 

 

 

 

Предполагается, что v

= − A

. В этом частном случае решение

 

0

m

 

 

 

 

исходной задачи имеет вид

 

 

 

 

 

 

 

 

 

 

x(t ) = −

A sin t,

 

 

 

 

 

 

m

 

 

 

соответствующий случаю гармонического колебания точки вблизи

ее начального положения. Иными словами, t > 0

x(t)

A .

 

 

 

 

 

 

m

Очевидно, что при малом отклонении величины v0 от выбран-

ного значения, например за счет погрешности округления данных

в ЭВМ, при численном решении задачи получаем (рис. 1.5)

 

x (t ) → ∞.

 

 

 

 

 

t → ∞

 

 

 

Таким образом, даже при корректно поставленной задаче ее

численное решение может быть получено со значительной погреш-

ностью.

 

 

 

 

 

 

x(t)

 

 

 

 

 

 

10

 

 

 

 

 

 

8

 

 

1

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

4

 

 

 

 

 

 

2

 

 

 

2

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-2

 

 

 

 

 

 

0

10

 

20

30

 

t

Рис. 1.5. Действительное решение (нижняя кривая 2) и возмущенное

решение (верхняя кривая 1) задачи из рассмотренного примера

26

2. ВВЕДЕНИЕ В МЕХАНИКУ СПЛОШНОЙ СРЕДЫ

Механика сплошной среды – обширная часть механики, посвященная движению газообразных, жидких и твердых деформируемых тел. В отличие от теоретической механики, изучающей движения материальной точки, систем материальных точек и абсолютно твердых тел, в механике сплошной среды рассматриваются движения таких материальных тел, которые заполняют пространство непрерывно, сплошным образом и расстояния между точками которых во время движения меняются. Иными словами, механика сплошной среды имеет дело с деформируемыми объектами, то есть с жидкостями, газами, твердыми телами и конструкциями.

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

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

2.1. Актуальные проблемы механики сплошной среды

Существует множество прикладных вопросов, на которые можно дать ответы только на основании специальных знаний. Например:

какова скорость вытекания газа из отверстия в баллоне, в котором газ находится в сжатом состоянии;

как будет двигаться в атмосфере воздушный циклон;

как можно снизить воздушное сопротивление автомобиля

исамолета или водяное сопротивление корабля;

27

как построить здание или телевизионную башню высотой

всотни метров, мост с пролетом между двумя ближайшими опорами более километра;

что произойдет с увеличением или уменьшением диаметра воздушного винта на самолете;

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

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

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

Примерами актуальных в настоящее время проблем являются следующие:

снижение сопротивления тел при движении в воде с большими скоростями;

создание и удержание плазмы с температурой в миллионы градусов;

выяснение особенностей поведения материалов при больших нагрузках и высоких температурах;

определение сил, действующих на здания и сооружения при взрывах;

объяснение общей циркуляции воздуха в атмосфере, построение прогнозов погоды;

изучение процессов в растениях и живых организмах;

эволюции звезд, явлений, происходящих на Солнце,

и многие другие.

Некоторые наиболее существенные разработанные в настоящее время проблемы механики сплошной среды [23]:

1. Воздействие жидкости и газа на движущиеся в них тела; силы, действующие со стороны газа или жидкости на тело, определяются скоростью движения газа или жидкости, их плотностью и вязкостью, формами и размерами взамиодействующих тел. Ввиду этого изучение движения тел непосредственно связано с изучением

28

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

2. Движение жидкости и газа по трубам внутри различных машин. В этих задачах основное значение отводится законам взаимодействия жидкости с границами потока, в частности сопротивлению подвижных и неподвижных твердых стенок; явлениям неравномерности в распределении скоростей и др. Эти задачи имеют важное значение для проектирования газопроводов, нефтепроводов, насосов, турбин и других гидравлических машин.

3. Фильтрация, то есть движение жидкости сквозь почву и другие пористые среды. Движение воды в почве необходимо учитывать при постройке фундаментов различных сооружений (плотин, опор мостов, гидростанций), при создании подземных туннелей. Большое значение фильтрация имеет в нефтяном деле.

4.Гидростатика, или равновесие жидкостей и тел, плавающих внутри и на поверхности жидкости; равновесие вращающихся масс жидкости под действием сил тяготения.

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

6.Неустановившиеся движения газов с химическими превращениями при взрывах, детонации и горении, например в потоке воздуха в цилиндрах поршневых машин или камерах реактивных двигателей.

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

29

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

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

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

прочность частей машин.

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

тельного решения.

 

 

 

 

Можно

упомянуть

о

механических

проблемах,

связанных

с движением

гомогенных

и

гетерогенных

смесей, песков, снега

и различных

грунтов, сплавов, жидких

растворов,

суспензий

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

2.2. Основные гипотезы

Механика материалов, по определению С.П. Тимошенко12 [28], представляет собой раздел механики сплошной среды, в котором изучается поведение твердых деформируемых тел при различных

12 Тимошенко Степан Прокофьевич (11.12.1878–29.05.1972) – российский и американский ученый-механик, внес значительный вклад в развитие теории упругости. Иностранный член АН СССР, профессор Мичиганского и Стэнфордского университетов, один из основателей Национальной академии наук Украины.

30

Соседние файлы в папке книги