ЦІЛЬ ЛЕКЦІЇ
· познайомити із задачами курсу "Числові методи та автоматизація проектування":
· викласти короткий зміст курсу;
· дати список рекомендованої літератури;
· ввести поняття «чисельні методи»;
· викласти основні положення апроксимації;
· навести приклад апроксимації табличних даних
Ключові слова: чисельний метод, апроксимація, сплайн, припасування кривої
ВСТУП
Мета
Метою викладання дисципліни “Чисельні методи та автоматизація проектування транспортних споруд” є формування знань та навиків, що необхідні для чисельного аналізу напружено-деформованого стану транспортних споруд, знайомство із сучасними дискретними моделями, знайомство із сучасними комп’ютерними технологіями проектування споруд. Значна увага приділяється формуванню навиків самостійної роботи з комп’ютером, як це передбачено моделлю інженера за спеціальністю 7.092106 “Мости і транспортні тунелі” .
Завдання вивчення дисципліни формулюються в наступному:
· знайомство з сучасними математичними дискретними моделями теорії споруд;
· знайомство із сучасними ідеями та принципами автоматизації проектування транспортних споруд;
· опанування технології роботи в програмних комплексах MathCAD, Excel, AutoCAD, SCAD.
Короткий зміст.
Задача предмету – познайомити з основними методами, використовуваними в теорії споруд, сучасними математичними моделями їх розрахунку, познайомити з ідеями і термінологією у області автоматизованого проектування споруд, надати можливість одержати навики використовування сучасних комп'ютерних технологій в проектуванні.
В курсі лекцій розглядаються наступні основні питання:
· апроксимація та інтерполяція;
· чисельна інтеграція;
· розв’язок систем лінійних рівнянь;
· розв’язок нелінійних рівнянь;
· матрична алгебра;
· математичні моделі статичного розрахунку споруд засновані на варіаційних принципах;
· метод скінчених елементів;
· принципи і програми автоматизованого проектування споруд.
Тут відводиться рівна кількість годинника на курс лекцій і на практичні заняття.
Практичні заняття ставлять свою за мету прищепити навики використовування сучасних інтегрованих програмних комплексів, вживаних в автоматизованому проектуванні. Йдеться про використовування, в першу чергу, таких програм: MathCAD, Excel, AutoCAD і SCAD.
2. Городецкий А.С., Заворицкий В.И., Лантух-Лященко А.И., Рассказов А.О. М.: Автоматизация расчетов транспортных сооружений. Транспорт, 1989. – 232 с.
3. Лантух-Лященко А.И. ЛИРА. Программный комплекс для расчета и проектирования конструкций. Учебное пособие М.-К.: Факт, 2001. – 359 с.
4. Варвак П.М. и др.. Метод конечних элементов. – К.: Вища школа, 1981. –176 с.
5. Плис И.А., Сливина Н.А. MathCAD 2000. Математический практикум. Учебное пособие. – М.: Финансы и статистика, 2000. – 656 с.
ЧИСЕЛЬНІ МЕТОДИ
Більшість задач розрахунку споруд не мають розв’язку у замкнутому вигляді за допомогою відомих функцій. В цих випадках вихід один – отримати розв’язок шляхом чисельної процедури, яка оперує не функціями а деякою множиною чисел. Ці процедури складають велику і дуже важливу гілку математики, яку називають чисельні методи.
¨
Визначення. Чисельними наз. методи наближеного або точного розв’язку задач прикладної математики, які ґрунтуються на побудові скінченної послідовності дій над скінченною множиною чисел.
Нерідко чисельні методи розрахунку споруд ототожнюють з дискретними методами. І це цілком вірно, тому що математична дискретизація континуального об’єкту за своєю суттю є чисельним розв’язком задачі розрахунку споруди.
¨
Визначення. Дискретними наз. методи наближеного розв’язку задач теорії споруд, які ґрунтуються на віртуальній (уявній) заміні суцільного середовища системою елементів скінчених розмірів які зв'язані між собою в скінченному числі точок.
Сьогодні застосування чисельних методів в задачах теорії споруд стало настільки поширеним, що аналітичні методи в складних задачах майже не застосовуються. Це пояснюється бурхливим розвитком обчислювальної техніки і тим фактом, що комп’ютерні програми для дискретних задач є більш простими
Апроксимація та інтерполяція
Апроксимація
Однією із важливих задач чисельних методів є математичний опис зв'язків в процесах і явищах. Універсальним способом такого опису є аналітична функціональна залежність
, (1.1.)
якою незалежним змінним , ставиться у відповідність одне значення функції . У кожному конкретному випадку формальне співвідношення (1.1.) має бути представлено певною функцією. Що ж стосується самою функції , як правило, вона невідома. Тоді і залучається апарат чисельних методів, який має назву апроксимація (наближення).
Наголосимо, що мова про апроксимацію йде в тому випадку, коли є потреба у визначенні функції за допомогою іншої функції .
Така заміна робиться, якщо нічого невідомо відносно функції або ж є недостатня кількість інформації для її точного визначення.
Наприклад, у проектуванні мостів, нам завжди потрібна функція, що описує зміну тиску на дану балку прогонової будови в залежності від положення одиничної сили, що переміщається в поперечному напрямку. Сьогодні в теорії споруд немає такої функції в загальному виді. Але її можна одержати апроксимацією дискретного чисельного розв’язку задачі або на основі таблиці експериментальних даних.
Або інший приклад: апроксимація часто потребується там де є підбір аналітичної залежності, що описує результати експерименту.
Не менш важливою областю апроксимації є апріорне наближене представлення шуканої функції в чисельних методах теорії споруд.
¨
Визначення. Апроксимація є наближене визначення невідомої функції F(x) функцією заданого виду f(x). Це наближення виконується на основі наявних дискретних значень в окремих точках заданого інтервалу х.
Як апроксимуючу функцію, часто застосовуються алгебраїчний поліном ступеня n. Така функція має вид:
, (1.2)
тобто парабола n– ого ступеня.
Полином (1.2) має n+1 коефіцієнт. Відповідним підбором коефіцієнтів можна задовольнити n+1 умов. Підбір виконується шляхом розв’язку системи n+1 лінійних рівнянь з вектором вільних членів з відомих значень функції у заданих точках (вузлах).
Звідси випливає, що якщо ми хочемо апроксимувати невідому функцію так, щоб значення збігалися з у трьох вузлах, досить параболи другого ступеня, у чотирьох – третього і т.д. (рис.1.1).
¡ ¡ ¡ F(x)
____ f(x)
Рис. 1.1. Апроксимація функції F(x) функцією f(x)
При цьому апроксимуюча крива проходить через кожний вузол. Ця загальна процедура припасування кривої і являє собою апроксимацію функції, яка задана таблицею.
Проте, в задачі припасування кривої, як правило, не ставиться вимога, щоби крива проходила через кожний вузол. Звичайно криву прагнуть провести так, щоб її відхилення від табличних даних були мінімальними. Найбільш поширеним способом мінімізації відхилення є метод найменших квадратів. За цим методом апроксимуюча функція визначається так, щоб звести до мінімуму суму квадратів різностей між табличними значеннями в вузлах та апроксимуючої кривої.
Іншим, дуже поширеним способом апроксимації є застосування сплайнів.
Сплайном називається функція виду:
, (1.3)
де - поліном ступеня не вище ;
- коефіцієнти які одержують з умов на кінцях інтервалу апроксимації.
У практиці апроксимації найчастіше застосовуються квадратні і кубічні сплайни. Так квадратний сплайн ( ) має вигляд:
(1.4)
Найбільш поширеним, і сьогодні навіть класичним, став кубічний сплайн (n=3).
Сплайн третього ступеня має яскраву механічну інтерпретацію: це тонка сталева лінійка, вигнута таким чином, щоб вона стикалася з заданими точками дискретних значень шуканої функції . Будучи деформована таким чином, лінійка здобуває форму, при якій запасена в ній пружна енергія мінімальна. Використовуючи теорію вигину бруса при малих деформаціях, можна строго показати, що сплайн – це група сполучених кубічних багаточленів, у місцях сполучення яких перша і друга похідні неперервні.
Теорія сплайнової апроксимації інтенсивно розвивалася в останні 30 років і тепер вважається універсальним методом. Апроксимація сплайнами є в усіх відомих математичних і графічних програмних комплексах.
В теорії споруд часто застосовується апроксимація тригонометричним поліномом.
Функція виду.:
, (1.5)
де - параметр поділу інтервалу апроксимації, ;
L– довжина відрізку апроксимації,
називається тригонометричним поліномом
Наближення тригонометричним поліномом часто застосовується для апроксимації функції в чисельних дискретних розв’язках задач теорії споруд.
Інші функції.
Тут тільки назвемо декілька найбільш відомих апроксимуючих функцій і прийомів:
· експоненціальні функції;
· гіперболічні функції;
· поліноми Лежандра;
· функції Бесселя;
· апроксимація поліномами Лагранжа.
Приклад апроксимації
Метою прикладу є підбір (припасування) аналітичної функції, яка відображає залежність отриману в результаті експерименту у вигляді табличних даних.
Дослідження впливу пластифікатора на міцність епоксидного клею дали такі результати:
Кількість пластифікатора, %
Зростання опору на розтяг, %
0,3
3,8
4,8
7,6
6,1
3,3
Потребується знайти апроксимуючу функцію у вигляді полінома, якою встановлюється приблизна аналітична залежність між кількістю пластифікатора в епоксидному клейові та зміною його опору на розтяг.
Виконаємо підбір функції апроксимації в програмному комплексі MathCAD.
За апроксимуючу функцію візьмемо поліном 3 ступеня:
P(x) = a0+a1x+a2x2+a3x3.
Визначимо апроксимуючу функцію P(x), аі, і=1,2, ..., 4, в програмному комплексі за методом найменших квадратів. Для цього скористаємося функціями MathCAD regress та interp. Обчислення та графіки заданої і апроксимуючої функцій показані на рис. 1.2. Як видно з графіка, крива P(x) не проходить через вузли апроксимації, проте знаходиться в достатній близькості від них.
В нашому прикладі не ставилася задача виписати функцію P(x) в явному вигляді, тобто з числовими коефіцієнтами. Зробити це нескладно, MathCAD має відповідну процедуру для виводу на екран значень коефіцієнтів аі.
U
Увага. Користуватись апроксимуючою функцією можна тільки в межах того інтервалу, на якому вона визначена.
Рис.1.2. Припасування кривої
Запитання
1 Дайте визначення поняття «чисельний метод»
2 Визначення понять «чисельні методи» та «дискретні методи»
3 Чим пояснюється широке застосування чисельних методів в теорії споруд?
4 Дайте визначення поняття «аппроксимація»
5 Сформулюйте математичну ідею апроксимації. Дайте синонім поняття «апроксимація»
6 Наведіть приклад апроксимації в розрахунках мостів
7 В чому полягає процедура апроксимації. Які є вимоги до кривої апроксимації відносно вузлів апроксимації. Дайте графік кривої апроксимації
8 Назвіть поширений спосіб припасування кривої апроксимації і дайте ідею способу. Дайте графік кривої апроксимації