алгоритма, следует, что в процессе эксплуатации происходит смена ведущегоэлемента. Результаты решения представлены в табл. 2. Решение 1 получено приигнорировании смены ведущего элемента, решение 2 – с учётом его смены.Таблица 2 - Значения долговечности и вид активного ограниченияМетод решенияАктивное ограничениеR 240 МПа f 1400LЧисленное решение 1. t* 8,507 лет(9,33%)*t 7,314 лет(0,04%)t0,1года n 85 ( n 73 ) 2. t* 7,816 лет(0,45%)Эталонное решение **t 7,781 летt 7,311 лет t 0,005 лет n 1557( n 1460)Следует отметить, что при активном ограничении по жёсткостирасхождение результатов находится в рамках приемлемой погрешности. Этообусловлено тем, что предельное состояние в данном случае определяется ненапряжением, а перемещением, и смена ведущего элемента не оказывает влиянияна результат решения задачи.Таким образом, игнорирование смены ведущего элемента, очевидно, можетпривести к серьезному искажению полученных результатов. Гарантироватьдостоверность решения задачи долговечности при использовании методов типаРунге-Кутта без усложнения логики алгоритма или увеличения размерностейрабочих массивов невозможно.Устойчивость алгоритма относительно входных данных. При решенииоптимизационной задачи заранее неизвестны конкретные значения варьируемыхпараметров (ВП), которые будут реализованы на данном шаге поискаоптимального проекта. Как правило, величина шага интегрирования СДУявляется параметром, который не изменяется во время работы программы.Следовательно, число итераций при постоянном значении заданнойдолговечности также постоянно. В этом случае погрешность численного решениябудет существенно зависеть от конкретных значений ВП на данном шаге задачиНЛП. Таким образом, алгоритм, основанный на использовании методов типаРунге-Кутта, при всех вышеперечисленных недостатках может оказатьсянеустойчивым относительно входных данных.Для численной иллюстрации нарушения устойчивости известныхалгоритмов относительно входных данных, предлагается решение задачидолговечности стержня растянутого силой Q . На рис. 5 приведены кривые ростанапряжений для стержней с различными площадями поперечных сечений иосевой нагрузкой. Рассматривались размеры поперечного сечения и величинарасчётной нагрузки, для которых долговечность стержня, определённая по*аналитической формуле, была равна t 8, 00 годам. Затем задача решаласьчисленно. Оценивались значения напряжений в момент времени t* . Прииспользовании методов типа Рунге-Кутта для вычисления долговечности привеличине шага по времени t 0, 25 года погрешность численного решениясоставила 0,34% при R 0 2 (кривая 1), 2,44% при R 0 6 (кривая 2) и 25,06%ISSN 2079.3944. Вісник НТУ “ХПІ». 2012. №44(950)8
при R 0 16 (кривая 3). При этом дляA0 52,89 см 2 ; 0, 49 см-1 (кривая 1)значения напряжений оказались равными 239, 18 МПа; для A0 9,36 см 2 ; 1, 16 см -1(кривая 2) – 234, 14 МПа; для A0 4,83см 2 ; 1,61см -1 (кривая 3) – 179, 85 МПа. Такоерасхождение численных решений связано,очевидно, с использованием равномернойвременной сетки.Таким образом, традиционныеРис. 5. Зависимость напряжений отвремениалгоритмы, основанные на одношаговых методах, могут оказатьсянеустойчивыми относительно входных данных. Этой ситуации можно избежатьпри увеличении числа итераций (уменьшении шага интегрирования t ) илииспользования неравномерного шага t .Эффективность. Если при решении задачи расчёта НДС поведениеконструкции в агрессивной среде исследуется путём численного решения задачиКоши для СДУ (3) с помощью одношаговых методов, а для обеспеченияприемлемой точности решения задачи долговечности, шаг по времени tприходится принимать достаточно малым, причём в каждом узле временнойсетки решается задача с помощью МКЭ, то общее число обращений к процедуреМКЭ может достигать нескольких десятков, а в некоторых случаях и сотен. Вэтом случае возможно накопление погрешности вычислений при решении системлинейных алгебраических уравнений итерационными методами [5]. Для того,чтобы получить для всех случаев, рассматриваемых в предыдущем примере,результат с погрешностью не более 0,35%, необходимо выполнить не менее 200итераций методом Эйлера. Следовательно, возникает проблема уменьшениячисла итераций при численном решении задачи Коши и одновременно сниженияразмерности задачи МКЭ.Задача долговечности (блоки {А} и {В}) решается на каждом шаге поискаоптимального проекта при проверке функций ограничений. Число же шаговитерационного процесса при решении задачи НЛП исчисляется десятками тысяч.Проблема повышения эффективности вычислительных алгоритмов можетоказаться не менее значимой, чем вышеперечисленные проблемы.Совершенно очевидно, что вычислительные алгоритмы, основанные нанепосредственном использовании методов типа Рунге-Кутта в той их реализации,которая использовалась до настоящего времени, зачастую, ненадёжны,неэффективны, неустойчивы относительно входных данных, а в ряде случаев негарантируют получение приемлемого результата.Основные результаты и выводыНа основании проведенного анализа проблем, вызванных влияниемнапряжений на скорость коррозионного процесса, можно сделать следующийвывод. Непосредственное использование существующих численных методоврешения СДУ применительно к поставленной задаче не только не являетсяISSN 2079.3944. Вісник НТУ “ХПІ». 2012. №44(950)9
- Page 1 and 2: МІНІСТЕРСТВО ОСВІТ
- Page 3 and 4: УДК 624.04:539.4:519.853ТЕХН
- Page 5 and 6: Расчёт НДС и долгов
- Page 7: Долговечность конс
- Page 11 and 12: УДК 629.7.036.001Б. Щ. МАМ
- Page 13 and 14: полном отсутствии
- Page 15: газодинамического
- Page 19 and 20: Согласно выведенно
- Page 21 and 22: этапа: образование
- Page 23 and 24: окисных плёнок, раз
- Page 25 and 26: УДК 681.518.3ИНФОРМАЦИ
- Page 27 and 28: програмне забезпеч
- Page 29 and 30: характеризують ста
- Page 31 and 32: ситуацій {K S }, що ві
- Page 33 and 34: Процеси сценарного
- Page 35 and 36: Рис. 3. Узагальнений
- Page 37: ВступДля моделюван
- Page 40 and 41: ТомуGopttH H (3 V V ) tQ Q V V
- Page 42 and 43: является исчерпыва
- Page 44 and 45: выбор частных коэф
- Page 46 and 47: внутри кластеров. В
- Page 48 and 49: 3. Суммарная внутри
- Page 50 and 51: 6x3 0,540780 10,45922 33,84673 0,00
- Page 52 and 53: 23Продолжение табли
- Page 54 and 55: наибольшее примене
- Page 56 and 57: Отримання точних р
- Page 58 and 59:
маршрутних мереж м
- Page 60 and 61:
моделей маршрутних
- Page 62 and 63:
у вільних каналах,
- Page 64 and 65:
R 1R1N 2kn1jNXP xk Re. n 1nn210n20N
- Page 66 and 67:
Разработка методов
- Page 68 and 69:
проблем в протипож
- Page 70 and 71:
Таблиця 2- Стан поже
- Page 72 and 73:
УДК 579.61ТЕХНОЛОГИИ
- Page 74 and 75:
Перелік найпопуляр
- Page 76 and 77:
штамів, які є їх осн
- Page 78 and 79:
Микробиол. журн. 2002.
- Page 80 and 81:
превышением предел
- Page 82 and 83:
водопользования, п
- Page 84 and 85:
концентрацию либо
- Page 86 and 87:
серпня 2009 р. за № 767/
- Page 88 and 89:
Анализ последних и
- Page 90 and 91:
Подобная инертност
- Page 92 and 93:
В то же время изопр
- Page 94 and 95:
N-Ацетокси-N-метокси
- Page 96 and 97:
Ключові слова: N-хло
- Page 98 and 99:
MeO 2COAcNOMeMeOHt-BuOHMeO 2COMeNOM
- Page 100 and 101:
алкоксимочевин тре
- Page 102 and 103:
H(2b)…O(1) 2,21 Å, которы
- Page 104 and 105:
диалкокси-N'-арилмо
- Page 106 and 107:
упарили в вакууме 25
- Page 108 and 109:
(6H, NOCHMe 2 , 3 J = 6,4 Гц), 3,
- Page 110 and 111:
9,48 (с 1H, NH), 9,93 (с, 1H, NHO
- Page 112 and 113:
Выводы: Таким образ
- Page 114 and 115:
Ключевые слова: ник
- Page 116 and 117:
аккумуляторов, %; ω[N
- Page 118 and 119:
2. С повышение темпе
- Page 120 and 121:
К- коэффициент гран
- Page 122 and 123:
Данные погрешности
- Page 124 and 125:
очистки от величин
- Page 126 and 127:
S 2 = 10,55рН 3 - 226,5рН 2 + 1
- Page 128 and 129:
природних вод, є по
- Page 130 and 131:
Алколін-МХ, Алколін
- Page 132 and 133:
Як видно із таблиці
- Page 134 and 135:
МПА-А-1МПА-А-2МПА-ДМ-1
- Page 136 and 137:
Розроблено спосіб
- Page 138 and 139:
производства полик
- Page 140 and 141:
кремния от 50 до 300 м
- Page 142 and 143:
В процессе водород
- Page 144 and 145:
10 -12 10 - 11 c [5]. Это об
- Page 146 and 147:
продиффундировавш
- Page 148 and 149:
Emets B., Almazova E. // Bulletin o
- Page 150 and 151:
nw n KБыли сделаны вы
- Page 152 and 153:
y methods of pressure casting for p
- Page 154 and 155:
составляет 32000 л/(м 2
- Page 156 and 157:
В дальнейшем необх
- Page 158 and 159:
Следует отметить в
- Page 160 and 161:
Развитие гигиениче
- Page 162 and 163:
ЗМІСТЛяшенко О. А.,
- Page 164:
НАУКОВЕ ВИДАННЯВІС