Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Автоматизація нелінійного розрахунку складних просторових систем споруд методом скінченних елементів

Автореферат

дисертації на здобуття наукового ступеня доктора технічних наук

Загальна характеристика роботи

Актуальність теми. Робота базується на новому методі дослідження складних технічних об'єктів, який знаходиться на стику теоретичних і експериментальних підходів, з використанням сучасних комп'ютерних технологій автоматизованого проектування. Новий підхід включає ряд послідовних взаємозалежних етапів дослідження: створення математичного забезпечення аналізу і синтезу проектних рішень та програмного забезпечення у системі автоматизованого проектування (САПР) об'єктів будівництва на основі адекватної і ефективної математичної моделі та чисельного методу нелінійного аналізу; проведення комплексу чисельних досліджень; узагальнення результатів досліджень; отримання оптимальних проектних рішень. У результаті розвитку комп'ютерних технологій цей стиль теоретичного чисельного аналізу трансформувався в нову інформаційну технологію дослідження у вигляді чисельного експерименту. Ідеологія автоматизованого проектування і чисельного експерименту є методологічним фундаментом даної роботи.

Одним із основних етапів роботи є розробка адекватної математичної моделі об'єкта дослідження. Математична модель є традиційним інструментом теоретичного дослідження. Проте застосування нової запропонованої в роботі комп'ютерної технології дослідження суттєво підвищує можливості цього інструмента і вимоги до нього. Адекватність моделі визначається тим, наскільки задовільно вона відображає якісні і кількісні характеристики об'єкта дослідження.

Просторові пластинчаті і стержневі системи та їхні комбінації є широко розповсюдженими несучими конструктивними системами сучасних складних будівельних об'єктів різного призначення в звичайних та складних (70% території України) грунтових умовах. Це пояснюється їхньою високою міцністю, значною жорсткістю і порівняно малою вагою. Розвиток сучасних будівельних об'єктів характеризується збільшенням інтенсивності силових та деформаційних навантажень. У сполученні з прагненням до зниження матеріалоємності це приводить до розгляду проектних рішень, які в процесі експлуатації допускають деформації і напруження у системах споруда-основа, що перевищують їх умовно пружні значення. У таких системах протікають процеси, характерні для нелінійних механічних систем, у яких в залежності від параметрів навантажень реалізуються принципово різні режими деформування. Процес деформування може супроводжуватися порушенням контакту споруди і основи, що приводить до різкої зміни просторової конфігурації статичного стану, що суттєво змінює традиційні уявлення про розподіл напружень і деформацій у системі.

Розглянутий клас об'єктів має суттєво просторовий характер роботи, особливо, при «екзотичних» режимах деформаційних навантажень у вигляді кривизни опуклості (увігнутості), ступінчатих деформацій земної поверхні. Тому дослідження виділеного класу об'єктів необхідно проводити в просторовій постановці.

Існуючі численні математичні моделі і методи аналізу споруд в складних ґрунтових умовах, як правило, реалізовані в рамках плоскої постановки. Ряд найважливіших для теорії і практики задач у просторовій постановці ще не вирішені, або отримані результати за спрощеними чисельними моделями, а також із застосуванням «домашинних» методів досліджень не можуть задовольняти сучасним вимогам точності результатів.

Створення сучасного математичного і програмного забезпечення автоматизованого проектування на основі універсальної, адекватної і ефективної математичної моделі, чисельного методу нелінійного аналізу складних просторових систем споруд на деформованій основі, що забезпечують суттєве уточнення існуючих моделей і методів аналізу, прийняття оптимальних проектних рішень, представляє серйозну і актуальну наукову проблему САПР об'єктів будівництва в складних ґрунтових умовах, що має важливе народногосподарське значення.

Зв'язок роботи з науковими програмами. Дослідження виконані відповідно до плану науково-дослідних робіт Київського національного університету будівництва і архітектури Міністерства освіти і науки України на 2002-2004 р., зокрема, по темі «Теорія і методи дослідження динамічних хвильових процесів деформування, механізмів руйнування та втрати стійкості просторових конструкцій» (номер державної реєстрації 0102U000928).

Мета і задачі досліджень. Мета роботи - створення, адекватної і ефективної математичної моделі, універсального методу нелінійного аналізу у складі сучасного програмного забезпечення САПР об'єктів будівництва в звичайних і складних ґрунтових умовах на основі методу скінченних елементів (МСЕ), суперелементів, шагового методу навантаження та методу оптимізації з єдиним законом деформування матеріалів конструкцій та ґрунтів; встановлення на основі чисельного аналізу у САПР загальних закономірностей деформування у всьому діапазоні силових і деформаційних навантажень у процесі зведення, експлуатації і реконструкції об'єктів для отримання оптимальних проектних рішень.

Досягнення мети здійснюється розв'язанням наступних основних задач:

- створення універсальної, ефективної і адекватної математичної моделі у вигляді математичного забезпечення САПР просторової системи споруда-основа з єдиним законом деформування матеріалів конструкцій і ґрунтової основи у всьому діапазоні силових і деформаційних навантажень;

- розробка загального методу автоматизації нелінійного аналізу просторових систем споруда-основа у всьому діапазоні силових і деформаційних навантажень на основі теорії МСЕ, суперелементів, методу оптимізації і шагового методу послідовного навантаження;

- реалізація у вигляді алгоритмів та програм загального методу нелінійного аналізу просторової системи споруда-основа у складі сучасного програмного забезпечення САПР для проведення чисельних досліджень і вирішення практичних задач оптимального проектування, безпеки зведення, експлуатації і реконструкції споруд у звичайних і складних ґрунтових умовах;

- обґрунтування достовірності запропонованої математичної моделі та загального методу автоматизації нелінійного розрахунку на основі розв'язання ряду тестових задач і проведення комплексу чисельних досліджень експериментальної споруди при нерівномірних деформаціях земної поверхні;

- встановлення на основі чисельних досліджень загальних закономірностей розподілу напружень і деформацій у системі споруда-основа при впливі нерівномірних деформацій земної поверхні в складних ґрунтових умовах у процесі зведення, експлуатації і реконструкції для отримання оптимальних проектних рішень;

- реалізація результатів досліджень у нормативних документах, навчальному процесі, при автоматизованому проектуванні об'єктів будівництва у звичайних і складних ґрунтових умовах.

Об'єктом дослідження є автоматизація проектування споруд складної просторової структури при взаємодії з нерівномірно-деформованою основою в звичайних та складних ґрунтових умовах у всьому діапазоні силових та деформаційних навантажень при їх зведенні, експлуатації та реконструкції.

Предмет дослідження: математичне моделювання, чисельний аналіз і синтез проектних рішень на основі математичного та програмного забезпечення САПР складних об'єктів будівництва у звичайних і складних ґрунтових умовах.

Методи дослідження. Створений і реалізований у складі програмного забезпечення САПР загальний метод автоматизації нелінійного розрахунку складної просторової системи об'єкта в складних режимах силових і деформаційних навантаженнях базується на застосуванні теорії МСЕ для визначення компонентів напружено-деформованого стану (НДС) системи на кожнім етапі навантажень, методу послідовного навантаження для розв'язання нелінійної задачі і методу оптимізації для отримання оптимальних проектних рішень.

Наукова новизна одержаних результатів:

- вперше запропонована і створена у рамках вирішення проблеми САПР нова універсальна, адекватна і ефективна математична модель нелінійного деформування складних просторових систем споруд на основі єдиного експонентного закону деформування матеріалів конструкцій і грунта;

- створено у складі сучасного програмного забезпечення САПР загальний метод автоматизації нелінійного розрахунку на основі методу скінченних елементів, суперелементів і методу послідовного навантаження складних просторових систем споруд із застосуванням критеріальних параметрів, що повною мірою характеризують особливий стан і еволюцію статичного процесу у всьому діапазоні силових і деформаційних навантажень;

- на основі отриманих нових результатів комплексу чисельних досліджень виявлено суттєві ефекти при врахуванні фактора нелінійності взаємодії споруди і основи на якісні і кількісні оцінки параметрів, що дозволяє багаторазово знижувати напруження в системі і давати вичерпний прогноз змін структури статичних станів у всьому діапазоні навантажень;

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

Практичне значення одержаних результатів. Практична цінність роботи полягає у використанні результатів виконаних досліджень:

- у проектних та науково-дослідних інститутах для автоматизованого проектування, чисельного аналізу та обґрунтування безпеки зведення, експлуатації і реконструкції споруд у звичайних і складних ґрунтових умовах України та зарубіжжя. Зокрема, при вирівнюванні автоматизованими домкратними системами житлових споруд: 14-ти поверхової у Москві, 9-ти поверхових у Тбілісі та Катовіце;

- у навчальному процесі для автоматизованого проектування будівельних об'єктів у складних ґрунтових умовах України та зарубіжжя із застосуванням створеного програмного забезпечення САПР на основі програмного комплексу (ПК) ЛІРА;

- при розробці програмного забезпечення САПР, зокрема, при реалізації загального методу нелінійного розрахунку у вигляді алгоритмів та програм ПК ЛІРА, що застосовується численними вітчизняними і зарубіжними організаціями;

- при розробці нормативних документів, зокрема, «СНиП 2.01.09-91. Здания и сооружения на подрабатываемых территориях и просадочных грунтах», «Методические указания по расчету и проектированию комплекса защитных мероприятий каркасных зданий, строящихся на просадочных грунтах».

Особистий внесок здобувача полягає в створенні у складі математичного та програмного забезпечення САПР на основі ПК ЛІРА об'єктів будівництва у звичайних і складних ґрунтових умовах універсальних, адекватних і ефективних скінченноелементних та суперелементних моделей споруда-основа, як єдиних просторових систем, загального методу нелінійного аналізу і синтезу проектних рішень складних просторових систем споруд з єдиним експонентним законом деформування матеріалів конструкцій і ґрунтів на основі теорії МСЕ, суперелементів, методу оптимізації та шагового методу навантаження; проведенні комплексу чисельних досліджень і встановленні загальних закономірностей деформування систем у всьому діапазоні силових та деформаційних навантажень у процесі зведення, експлуатації та реконструкції.

Апробація результатів дисертації. Основні положення дисертаційної роботи доповідалися і схвалені на науково-технічному семінарі в Челябінському політехнічному інституті (Челябінськ, 1988); на міжнародних конференціях в університетах Польщі (Криніца, 1994; Щецін, 1996); на міжнародних семінарах в Науково-дослідному інституті механіки і прикладної математики Ростовського державного університету (Ростов-на-Дону, 1989, 1998); на міжнародному симпозіумі в Грузинському технічному університеті (Тбілісі, 2001); на міжнародній науково-технічній конференції Національного авіаційного університету (Київ, 2002); на науково-практичних конференціях Київського національного університету будівництва і архітектури (Київ, 2002,2003, 2004).

Публікації. Основні положення дисертації опубліковані в 28 статтях.

Структура та обсяг дисертації. Дисертація складається з вступу, п'яти розділів, основних висновків, списку використаних джерел та додатка. Обсяг дисертації - 323 стор., у тому числі ілюстрації - 82 стор., таблиці - 2 стор., список використаних джерел з 295 найменувань - 29 стор.

Основний зміст роботи

автоматизація деформаційний математичний

У вступі обґрунтована актуальність, наукова новизна та практичне значення роботи, дана її загальна характеристика.

РОЗДІЛ 1 присвячений огляду робіт по темі дисертаційної роботи і аналізу стану проблеми математичного моделювання, чисельного аналізу у САПР об'єктів будівництва в звичайних і складних грунтових умовах.

В теорії математичного моделювання основ споруд існують три основні напрямки.

Перший напрямок характеризується застосуванням методів аналізу несучих систем споруд, що базуються на моделі Фусса-Вінклера. Спосіб коефіцієнта постелі був детально розроблений у роботах С.П. Тимошенко, П.Л. Пастернака, Н.Л. Пузиревського, А.Н. Крилова, А.А. Уманського, В.А. Кисильова та ін. Завдяки відносній простоті розв'язань задач, основаних на гіпотезі Вінклера, вона знайшла особливо широке застосування в світовій практиці, а також у нашій країні в розрахунках конструкцій споруд на деформованій основі.

Другий напрямок характеризується застосуванням методів аналізу несучих систем споруд, що базуються на моделях пружного півпростору і комбінованих моделей поверхні основи. Спостереження за деформаціями ґрунтової основи показали, що осідання відбувається не тільки в навантаженому місці, але і поруч з ним. З цією метою була запропонована Г.Е. Проктором гіпотеза пружної напівплощини і пружного півпростору, яка одержала широке поширення тому, що математичний апарат для неї був добре розроблений раніше в теорії пружності. Н.М. Герсеванов, В.А. Флорин, Б.Н. Жемочкін, В.І. Кузнецов, М.І. Горбунов-Посадов і ін. внесли суттєвий вклад у розробку нової теорії розрахунку конструкцій на пружній основі. М.І. Горбуновим-Посадовим, О.Я Шехтером, К.Е. Єгоровим та ін. розвинута модель пружного шару. Ця модель, як правило, дає кращий збіг чисельних і експериментальних даних у порівнянні з моделлю півпростору, в тому випадку, коли найбільш деформованим є тонкий поверхневий шар ґрунту. Однак, через те, що модель пружного півпростору суттєво перебільшує розподільну здатність ґрунту і експериментами не підтверджується П.Л. Пастернаком, М.М.Філоненко-Бородичем та ін. запропоновані комбіновані моделі ґрунтової основи, що дають проміжний результат між розрахунками за гіпотезами Вінклера і пружного півпростору.

Третій напрямок характеризується застосуванням методів аналізу несучих систем споруд, що базуються на уточнених моделях поверхні основи з урахуванням нелінійних деформацій грунтів, які виникають при значних вертикальних і горизонтальних переміщення земної поверхні, характерних для складних ґрунтових умов. В результаті теоретичних та унікальних експериментальних досліджень споруд у складних грунтових умовах С.М. Клєпіковим суттєво розвинута модель Вінклера, на основі якої запропонована модель перемінного коефіцієнта жорсткості з застосуванням гіперболічного закону деформування поверхні основи, методу скінченних різниць і методу послідовних наближень для розв'язання контактних задач взаємодії конструкцій споруд з основою.

Математичні моделі деформування матеріалів конструкцій, зокрема, бетону і залізобетону почали активно розроблятися стосовно до скінченноелементного підходу для аналізу сучасних об'єктів різного призначення. Способи моделювання утворення тріщин при скінченноелементному аналізі, запропоновані різними авторами, можна розділити на дві групи. До першої групи відносяться дискретні моделі, в яких утворення тріщин фіксується на границях між скінченними елементами. Цей спосіб зв'язаний з необхідністю перенумерації вузлів після утворення тріщин і пошуку на кожнім кроці розрахунку відповідної схеми розчленовування конструкцій на скінченні елементи. До другої групи відносяться більш ефективні інтегральні моделі, у яких тріщина розглядається в тілі скінченного елемента. Аналіз розвитку теорій бетону і залізобетону показує, що завдяки ефективності, найбільше застосування мають інтегральні моделі, в тому числі моделі Карпенко Н.І. для просторових стержневих, пластинчастих і масивних конструкцій і моделі Генієва Г.А., що описують поводження матеріалів конструкцій з позицій теорії пластичності для загального випадку трьохосьового напруженого стану. В роботах Баженова В.А., Гуляра А.И., Козака А.Л., Сахарова А.С. на основі узагальнення теорії Н.І. Карпенко запропоновано моментну теорію деформування бетону і залізобетону з тріщинами з хорошим узгодженням отриманих результатів з експериментальними даними у всьому діапазоні трьохосьового напруженого стану.

Аналіз розвитку чисельних методів розрахунку просторових несучих систем споруд на деформованій основі дозволяє виділити два основних напрямки: чисельне моделювання несучих конструктивних систем споруд без врахування та з врахуванням деформованої основи з використанням відповідно роздільного і спільного методів розрахунку. У випадках, коли величини нерівномірності деформацій основи не перевищують граничних нормативних значень, у практиці проектування застосовується роздільний метод розрахунку. Однак, встановлені нормами граничні величини деформацій основи є не досить обґрунтованими і тому мають обмежену область застосування. Для виконання роздільних і спільних розрахунків споруд на деформованій основі широке застосування одержав метод скінченних елементів і суперелементів. МСЕ сполучає у собі переваги варіаційних і різницевих методів із гнучкістю одержання дискретних моделей інваріантних стосовно будь-якої системи координат і ефективно реалізується на сучасних персональних комп'ютерах. Ідея апроксимації континуальної області дискретними елементами належить Л. Ейлеру і М. Пуасону, яка одержала розвиток в основних роботах Р. Куранта, В. Прагера, Д. Аргіріса, О. Зенкевича, Р. Клафа та ін. У розвиток МСЕ вагомий вклад внесли роботи: А.В. Александрова, В.А. Баженова, П.М. Варвака, О.С. Городецького, Є.О. Гоцуляка, О.І. Гуляра, Є.С Дєхтярюка, В.С. Дейнека, І.Д. Євзєрова, В.В. Киричевського, В.Н. Кислоокого, В.Г. Корнєєва, А.І. Лантух-Лященко, Ю.І. Немчинова, В.Г. Піскунова, В.А. Постнова, О.О. Расказова, Л.А. Розіна, О.С. Сахарова, І.В. Сергієнко, О.Л. Синявського, І.Я. Хархуріма, Н.Н. Шапошнікова та інших.

Становленню та розвитку САПР об'єктів будівництва присвячені наукові праці В.А. Баженова, Д.В. Вайнберга, О.С. Городецького, Є.О. Гоцуляка, В.М. Глушкова, Л.Г. Дмітрієва, В.С. Дейнека, В.В. Киричевського, В.Н. Кислоокого, В.Г. Кобієва, А.І. Козачевського, А.І. Лантух-Лященко, П.П. Лізунова, Ю.І. Немчинова, В.І. Скуріхіна, І.В. Сергієнко, П.М. Сосіса, О.Л. Синявського, В.С. Шмуклера та інших.

Узагальнюючи огляд фундаментальних робіт, присвячених створенню комп'ютерних технологій розрахунку і проектування об'єктів будівництва, слід зазначити суттєвий вклад відомих Київських шкіл:

- Д.В. Вайнберга (В.А. Баженов, Є.О. Гоцуляк, О.І. Гуляр, Є.С Дєхтярюк, О.Л. Синявський, О.С. Сахаров, В.Н. Кислоокий, В.В. Киричевський, В.Г. Кобієв та інші), де створена автоматизована система ПРОЧНОСТЬ для дослідження і проектування будівельних та машинобудівельних конструкцій;

- П.М. Сосіса (О.С. Городецький, Г.Б. Гільман, І.Д. Глікін, Л.Г. Дмітрієв та інші). Сучасна школа САПР заснована О.С. Городецьким, у якій створені інтелектуальні програмні комплекси МОНОМАХ, ЛІРА для розрахунку та проектування будівельних об'єктів широкого призначення;

- І.В. Сергієнко, В.С. Дейнека (Інститут кібернетики ім. В.М. Глушкова), де за участю їх учнів створені математичне забезпечення та програмні комплекси NDS, НАДРА, ПЕОМ-Кластер для дослідження НДС ґрунтового масиву, просторового руху рідини в ґрунтовому середовищі на суперкоп'ютері із тридцяти двома 32-х розрядними процесорами для реалізації інформаційних технологій (СКІТ). Створено комплекс програм для розв'язання систем лінійних алгебраїчних рівнянь великих порядків (більше 0,5 мільйона невідомих) із симетричними матрицями у послідовно-паралельному режимі.

Розгляд загального стану зазначеної проблеми створення сучасних інформаційних технологій розрахунку і проектування об'єктів будівництва показує, що у переважній більшості опублікованих робіт комп'ютерний аналіз просторових систем споруд у складних грунтових умовах здійснюється, як правило, в спрощеній пружній постановці, або в плоскій нелінійно-пружній постановці. Практично відсутні публікації, зв'язані із створенням комп'ютерних технологій нелінійного аналізу та проектування просторових систем споруд у складних ґрунтових умовах з урахуванням всього комплексу найважливіших факторів: фізичної нелінійності матеріалів конструкцій та ґрунтового масиву, порушення контакту споруди і основи у процесі зведення, експлуатації та реконструкції аварійних споруд.

На основі аналізу стану проблеми визначені мета і задачі досліджень.

РОЗДІЛ 2 присвячений створенню універсальної математичної моделі системи споруда-основа на основі теорії МСЕ та суперелементів для чисельного аналізу у САПР статичних процесів фізично-нелінійного деформування складних просторових систем споруд при силових та деформаційних навантаженнях із застосуванням запропонованого і обґрунтованого єдиного експонентного закону деформування конструкційних матеріалів і ґрунту, який використовується при лінеаризації фізично нелінійної задачі.

Дискретизація суперелементами (ССЕ) системи споруда-основа дозволяє представити схему побудови вирішального рівняння аналогічно СЕ за виключенням побудови матриці жорсткості і вектора вузлових навантажень ССЕ, які визначаються не в явному вигляді як для СЕ, а в результаті чисельного розрахунку і побудови апроксимуючих функцій ССЕ .

Розглядаючи ССЕ як СЕ, для яких матриця жорсткості будується не аналітично, а за допомогою чисельного розрахунку, можна застосувати дві схеми побудови матриці жорсткості.

Схема 1 побудови матриці жорсткості ССЕ полягає в тім, що на супервузли розглянутої підсистеми накладаються в'язі, а елементи матриці жорсткості визначаються як реакції в цих в'язях від послідовних одиничних переміщень. Для ССЕ побудова матриці жорсткості і вектора вузлових навантажень виконується як для звичайного СЕ з підстановкою замість .

Апроксимація багатошарового ґрунтового масиву моделлю поверхні основи споруди в умовах просторового навантаження. У загальному випадку при дії вертикальних і горизонтальних навантажень , реакції основи в кожній точці (x, y) контакту основи з спорудою представимо у вигляді контактних нормальних і дотичних напружень , .

Для отримання законів розподілу зазначених контактних напружень по підошві споруди запропонована просторова скінченноелементна модель поверхні основи, що дозволяє замінити багатошаровий ґрунтовий масив. Для визначення вищезазначених вертикальних (w) і горизонтальних (u, v) переміщень існують наступні методи.

1. Теоретичний метод, заснований на розрахунку для кожної вертикалі геологічного розрізу осадок одним зі способів механіки ґрунтів, методом пошарового підсумовування, результати яких збігаються з даними натурних спостережень за осадками споруд.

2. Експериментальний метод випробовувань ґрунту штампами в умовах природного залягання застосовується для одержання фізико-механічних властивостей ґрунтів і визначення осадки штампа при дії рівномірно розподіленого навантаження.

3. Метод натурних випробовувань ґрунтів, заснований на експериментальному вдавленні фундаментів споруд у ґрунт в умовах природного залягання, площі і розміри яких практично збігаються з розмірами запроектованих фундаментів споруд.

В рамках дослідження для побудови фізично-нелінійної просторової моделі поверхні основи вперше запропоновані і обгрунтовані дотичні коефіцієнти жорсткості основи ,, при дії вертикальних і горизонтальних , навантажень на основі експонентних діаграм деформування поверхні грунту.

Фізично нелінійний СЕ поверхні основи споруди в умовах просторового навантаження із запропонованим експонентним законом деформування грунту.

Для побудови вирішального рівняння фізично нелінійного СЕ поверхні основи споруди в умовах просторового навантаження (рис. 2.4) при дії вертикальних та горизонтальних навантажень ,, (силових), ,, (деформаційних) з урахуванням запропонованих експонентних діаграм деформування ґрунту при стиску в напрямку Z і зсуві в напрямках X, Y (рис. 2.5) розглянемо статичну, геометричну і фізичну сторони задачі та їхній синтез:

Статичні рівняння: ; ; .

Геометричні рівняння: ; ; .

Фізичні рівняння: ; ; .

Лінеаризація нелінійних рівнянь виконується шаговим методом навантаження. У випадку приростів навантажень на нескінченно малі величини , , , розкладаючи в ряд Тейлора функції ,, , виключаючи члени вищих порядків малості.

З огляду на те, що , , , дотичні коефіцієнти жорсткості поверхні основи дорівнюють:

при стиску в діапазоні , :

- при зсуві в діапазоні , :

- при зсуві в діапазоні ,:

Дотичні жорсткістні характеристики , , обчислюються на кожнім m кроці послідовного навантаження за значеннями дотичних коефіцієнтів жорсткості, що визначаються по отриманих з попереднього m-1 кроку розрахунку сумарних деформаціях , , .

Напружено-деформований стан СЕ поверхні основи споруди на m етапі навантаження визначається МСЕ у переміщеннях.

Фізично нелінійний СЕ стержневих систем в умовах просторового навантаження із запропонованим експонентним законом деформування матеріалів конструкцій. Вирішальне рівняння фізично нелінійного СЕ стержня в умовах просторового навантаження з урахуванням запропонованих і обгрунтованих експонентних діаграм для матеріалів стержневих конструкцій отримаємо на основі наступної схеми, вперше використовуваної в дослідженнях О.С. Городецького.

Дотичні жорсткістні характеристики , , , , , обчислюються на кожному m кроці навантаження для перерізів, розташованих у точках інтегрування по довжині стержня, за значеннями дотичних модулів деформації , у дискретних точках поперечного перерізу відповідно до заданої сітки дискретизації перерізів. Поточні значення дотичних модулів деформації , у точці визначаються на основі гіпотези плоских перерізів по отриманій з попереднього m-1 кроку розрахунку сумарній деформації в точці по експонентних діаграмах для бетону й арматури (рис. 2.7).

Напружено-деформований стан СЕ стержня на m етапі навантаження визначається МСЕ у переміщеннях.

Апроксимуючі функції переміщень , , визначаються в результаті розв'язанням системи лінеаризованих рівнянь (2.20-2.22) відносно .

Фізично нелінійний СЕ пластинчатих систем в умовах плоского напруженого стану і згину із запропонованим експонентним законом деформування матеріалів конструкцій. Вирішальне рівняння фізично нелінійного СЕ тонкої пластини (рис. 2.8) в умовах плоского напруженого стану і згину з урахуванням запропонованих експонентних діаграм для матеріалів пластинчатих конструкцій (рис. 2.7) отримаємо на основі схеми, використовуваної у дослідженнях О.С. Городецького, В.С. Здоренко.

Апроксимуючі функції переміщень по області прямокутного 4-х вузлового СЕ

- для плоского напруженого стану - у вигляді полілінійних законів, що еквівалентно гіпотезі про лінійний розподіл переміщень уздовж координатних осей

- для згину - у вигляді неповного полінома 4-й степені від двох змінних (аналогічно СЕ Клафа)

Вирішальне рівняння у приростах r-го СЕ пластини в умовах плоского напруженого стану і згину має вигляд:

РОЗДІЛ 3 присвячений створенню загального методу нелінійного аналізу у САПР просторових систем споруд в складних грунтових умовах (рис. 3.1).

Запропонований загальний метод нелінійного аналізу дозволяє простежити дійсне поводження системи у всьому діапазоні навантаження, починаючи з ненавантаженого стану, включаючи стан нелінійного деформування матеріалів конструкцій і ґрунту і завершуючи процес навантаження руйнуванням системи. При цьому величина руйнівного навантаження визначається в процесі нелінійного розрахунку на основі прийнятих критеріїв руйнування по досягненню граничних деформацій або граничних напружень в елементах системи.

Універсальний алгоритм розв'язання фізично нелінійної задачі. Розв'язання нелінійної задачі виконується шаговим методом навантаження. Для визначення компонентів НДС системи на кожному ступіні навантаження використовується МСЕ у переміщеннях.

Схема методу послідовного навантаження. Призначається число видів навантажень і число ступенів навантаження для кожного -го виду навантаження.

Установлюються закони навантаження і визначаються прирости вузлових навантажень для кожного -го виду навантаження . При рівномірному навантаженні коефіцієнт навантаження .

Обчислюються дотичні жорсткістні характеристики і на їхній основі коефіцієнти матриці жорсткості r-го СЕ , що відповідають .

Підсумовуючи відповідні коефіцієнти, будується .

Оптимізація просторових систем споруд при багатьох навантаженнях.

Задача оптимального армування перетинів стержневих систем споруд від дії багатьох навантажень у загальному випадку розрахунку по міцності у термінах математичного програмування формулюється в такий спосіб.

Формування виразів дозволяє сформулювати задачу лінійного програмування, як кроку оптимізації армування від дії декількох розрахункових комбінацій зусиль. Задача оптимального армування на кожнім кроці оптимізації зводиться до розв'язання наступної задачі лінійного програмування.

РОЗДІЛ 4 присвячений обґрунтуванню адекватності математичної моделі та загального методу автоматизації нелінійного розрахунку систем споруда-основа на основі зіставлення даних чисельних досліджень з результатами відомих вітчизняних і зарубіжних експериментальних досліджень нелінійного деформування стержнів, пластин, ґрунтів і натурних випробовувань споруд.

Стержневі системи. Розрахункові моделі експериментальних зразків побудовані на основі фізично нелінійних 2-х вузлових СЕ стержня із застосуванням експонентних діаграм деформування бетону й арматури при стиску і розтягуванні і рівномірних законів навантаження. В якості параметрів діаграм прийняті нижчеподані експериментальні дані.

Балки. Серія 2. Зразки Б2, Б2P (експерименти Зайцева). Вихідні експериментальні дані (рис. 4.1.): L = 3 м; b = 0.15 м; h = 0.4 м Кубікова міцність бетону R=44.5 МПа. E_b = 26600 МПа. A_s = 7.9см²; A_s = 3.28см². R_s = 750 МПа; R_sc=402 МПа; E_s = 210000 МПа. Навантаження зосереджене. Граничні умови: опори защемлені (для зразка Б2 на одній опорі дозволене лінійне переміщення уздовж осі Х).

Пластинчаті системи. Розрахункові моделі експериментальних зразків побудовані на основі фізично нелінійних прямокутних 4-х вузлових СЕ тонкої пластини із застосуванням експонентних діаграм деформування бетону й арматури при стиску і розтягуванні. За параметри діаграм прийняті нижченаведені експериментальні дані. Закони навантаження: рівномірні.

Балка-стінка. Зразок 103 (експерименти Нільсона). Вихідні експериментальні дані (рис. 4.2): L = 2 м; H = 1 м; h = 0.08 м. Кубікова міцність бетону R=44.5 МПа. Крок арматурних стержнів s_x = s_z = 10 см, f_ax = f_az = 0.0282 см. R_s = 680 МПа, R_sc = 410 МПа. E_s = 210000 МПа, навантаження зосереджене. Граничні умови: шарнірні опори.

Плити. Зразки 825-827 (експерименти Баха і Графа). Вихідні експериментальні дані (рис. 4.3): a = b = 2 м; h = 0.122 м. Кубікова міцність бетону R=26.5 МПа, крок арматурних стержнів на нижній поверхні плити s_x = s_y = 10 см, f_ax = f_ay = 0.0407 см, R_s = 592.9 МПа, R_sc = 408 МПа. Навантаження у вигляді 16 вантажів рівномірно розподілених по полю плити. Граничні умови: шарнірні опори на контурі плити.

РОЗДІЛ 5 Можливості запропонованого загального методу чисельного аналізу у САПР розглянутих класів об'єктів ілюструються шляхом розв'язання ряду складних практично важливих задач автоматизованого проектування. До їхнього числа відносяться чисельні дослідження 9-ти поверхової експериментальної споруди на підроблюваних територіях Донбасу та 14-ти поверхового аварійного жилого будинку, зведеного на слабих ґрунтах у Москві.

Загальна характеристика експериментальної 9-ти поверхової споруди серії 179, запроектованої КиївЗНДІЕП за участю Донецького ПромбудНДІпроекта. Несуча просторова система споруди (рис. 5.1, 5.2) довжиною 26.4 м утворена із системи подовжніх і поперечних несучих стін, об'єднаних дисками перекриттів. Висота поверху - 2.8 м, проліт перекриттів в осях - 5.7 м. По периметрі всіх стін споруди запроектовані поверхові залізобетонні пояси.

Для проведення комплексу чисельних досліджень побудовані наступні три комбінації чисельних просторових моделей експериментальної споруди при ступінчатих деформаціях земної поверхні (рис. 5.1, 5.2): модель 1 лінійна споруда - лінійна основа; модель 2 лінійна споруда - нелінійна основа; модель 3 нелінійна споруда - нелінійна основа.

Лінійна просторова модель 1 побудована на основі наступних типів суперелементів з використанням пластинчатих СЕ в умовах плоского напруженого стану і вигину: суперелемент 1 підвалу; суперелемент 2 1-го поверху; суперелемент 3 типових 2, 4, 6, 8 - го поверхів; суперелемент 4 типових 3, 5, 7, 9 - го поверхів; суперелемент 5 технічного поверху; суперелемент 6 перекриттів поверхів. Основна схема включає модель поверхні основи та модель балкового фундаменту з використанням стержневих СЕ в умовах просторового навантаження.

Нелінійна просторова модель 2 побудована на основі одного типу суперелемента, що включає лінійні пластинчаті і стержневі СЕ несучої системи споруди (балкові фундаменти, зовнішні і внутрішні стіни, перекриття і покриття) і фізично-нелінійної скінченноелементної моделі поверхні основи в умовах просторового навантаження.

Нелінійна просторова модель 3 побудована на основі фізично нелінійних пластинчатих і стержневих СЕ несучої системи споруди і фізично нелінійних стержневих СЕ поверхні основи в умовах просторового навантаження. Схема навантажень моделей 1, 2, 3 прийнята у відповідності зі схемою натурних випробовувань споруди, яка складається із 11 етапів навантажень: 1-го силового та 10-ти деформаційних навантажень. Траєкторія навантажень нелінійних моделей 2, 3 прийнята рівномірною з величиною ступіні навантаження, що дорівнює 5% від величини навантаження, тобто 20 ступенів кожного виду навантаження:

Результати чисельного аналізу отримані у вигляді детальної картини розподілу максимальних N1 і мінімальних N3 головних напружень в просторових моделях 1, 2, 3 на всіх 11 етапах силових та деформаційних навантажень і ілюструються для найбільш навантаженої конструкції системи - внутрішньої подовжньої стіни при найбільш невигідній комбінації навантажень: уступ довжиною L_u=8.50 м, висотою H_u=0.132 м.

Аналіз отриманих результатів чисельних досліджень свідчить про суттєву якісну і кількісну відмінність у розподілі напружень у просторових моделях 1, 2, 3. Суттєве вирівнювання по висоті стіни і зниження максимальних по абсолютній величині головних напружень отримано в просторовій моделі 2 за рахунок нелінійного моделювання поверхні основи. При використанні моделі 3, нелінійна споруда-нелінійна основа, величини головних напружень N1 і N3 практично вирівнюються по висоті споруди, а також суттєво знижуються їхні максимальні значення по абсолютній величині в порівнянні з даними лінійної моделі 1.

Отримані багаторазові зниження значень головних напружень, у 2-3 рази в порівнянні з даними лінійно-пружного розрахунку, дозволили забезпечити допустимий за умовами експлуатації споруди рівень розкриття тріщин при заданих параметрах деформацій земної поверхні. Наведена на рис. 5.4 розрахункова схема розкриття тріщин і зон їхнього розподілу по подовжній стіні практично збігається з аналогічними даними натурного випробовування споруди.

Чисельний аналіз у САПР на основі ПК ЛІРА безпеки реконструкції складних будівельних об'єктів. Аварійний об'єкт у Москві. Основною метою чисельного аналізу у САПР було встановлення граничної величини нерівномірності вертикальних переміщень домкратів автоматизованої домкратної системи для безпеки вирівнювання й експлуатації аварійного об'єкта. Особливістю об'єкта, зведеного на слабих ґрунтах Москви, було те, що 14-ти поверховий жилий будинок з поперечними несучими стінами, і розмірами в плані 15х56 м був запроектований на пальових фундаментах (довжина паль 6 і 7 м) без врахування можливих нерівномірних осадок основи. Після зведення споруди в результаті навантажень від нерівномірних деформацій основи максимальна величина відхилення споруди від вертикалі в поперечному напрямку досягла 75 см (16.7 мм/м), яка більш ніж на порядок перевищувала граничну нормативну величину.

Для чисельного аналізу у САПР побудована просторова скінченноелементна модель 14-ти поверхової споруди відповідно до проектних геометричних і фізичних характеристик несучих елементів конструкцій, їхніх з'єднань і конструкцій посилення фундаментно-підвальної частини. Чисельні дослідження НДС просторової моделі споруди проводилися при найбільш невигідних комбінаціях навантажень: силових - від власної ваги споруди та деформаційних - від нерівномірних переміщень домкратів у вигляді прогинів та згинів у поперечному та подовжньому напрямках споруди. На основі аналізу компонентів НДС просторової системи споруди встановлена гранична величина нерівномірності переміщень домкратів - 1.2 мм/м, що дозволило вперше вирішити унікальну науково-практичну задачу реконструкції автоматизованою домкратною системою аварійного об'єкта для його безпечної експлуатації.

Чисельний аналіз у САПР безпеки складних технічних об'єктів атомної енергії. Мета - побудова адекватної скінченноелементної моделі у САПР та чисельний аналіз безпеки експлуатації типової захисної оболонки АЕС України в екстремальних режимах навантажень.

Дискретизація захисної оболонки виконана на основі стержневих, пластинчастих і просторових скінченних елементів. Циліндрична і купольна частини оболонки мають ортогональну сітку СЕ, що повторює траєкторії армоканатів і перехідну зону, що забезпечує стикування з анкерним карнизом. Моделювання системи попередньо напружених армоканатів виконано на основі стержневих скінченних елементів, що дозволяє адекватно відображати складні умови взаємодії армоканатів з поверхнею каналів.

Таким чином, запропоновані та створені адекватні математичні моделі, загальний метод аналізу у САПР на основі ПК ЛІРА дозволяють отримати оптимальні проектні рішення, що охоплюють життєвий цикл розглянутих об'єктів будівництва у звичайних та складних грунтових умовах в процесі зведення, експлуатації та реконструкції.

Висновки

У дисертаційній роботі вирішена актуальна наукова проблема розроблення науково-обґрунтованих основ математичного моделювання, створення та використання автоматизованих систем проектування, прийняття оптимальних проектних рішень в ситуаціях невизначеності, обумовлених реальним станом об'єктів будівництва в звичайних та складних ґрунтових умовах України. На основі застосування системного підходу, ідей і методів математичного моделювання, сучасних інформаційних технологій у роботі створені математичні моделі, алгоритмічні та програмні засоби САПР, застосування яких дозволило суттєво поліпшити якісні та кількісні показники автоматизованого проектування об'єктів будівництва в Україні.

У дисертації одержані такі основні теоретичні та практичні результати.

1. Вперше створено математичне та програмне забезпечення САПР на основі ПК ЛІРА із застосуванням запропонованих універсальних, ефективних і адекватних математичних моделей та загального методу нелінійного аналізу складних просторових систем споруд у всьому діапазоні складних режимів силових і деформаційних навантажень. У рамках вирішеної проблеми САПР здійснено вибір і обґрунтування єдиного закону деформування матеріалів конструкцій і ґрунту для розв'язання нелінійної задачі, що дозволяє отримати оптимальні проектні рішення складних просторових систем споруд.

2. На базі методу скінченних елементів і суперелементів, методу оптимізації, шагового методу послідовного навантаження розроблені і реалізовані ефективні алгоритми та програмне забезпечення САПР для аналізу станів складних просторових систем споруда-основа при силових і деформаційних навантаженнях. Інваріантість схеми чисельної дискретизації щодо конкретних особливостей просторового обрису об'єктів дозволяє розглянути з єдиних позицій широке коло задач автоматизації нелінійного аналізу і проектування складних просторових систем споруда-основа в екстремальних умовах статичних режимів навантажень.

3. Запропоновано та реалізовано на основі єдиного алгоритму загальний метод автоматизації нелінійного розрахунку для дослідження процесу нелінійного деформування просторової системи споруда-основа для можливих і принципово різних комбінацій лінійних та нелінійних станів споруди та основи у всьому діапазоні силових та деформаційних навантажень, який дозволяє оптимізувати обчислювальний процес у залежності від конкретних особливостей станів об'єктів дослідження.

4. Встановлено, що застосування запропонованої суперелементної моделі системи споруда-основа підвищує більш ніж на порядок ефективність загального методу автоматизованого розрахунку, що суттєво знижує трудомісткість і об'єм обчислювальних операцій і підвищує якість і економічність проектних рішень.

5. Запропонований загальний метод автоматизації нелінійного розрахунку на основі шагового методу послідовного навантаження є ефективною алгоритмічною основою комплексного аналізу нелінійного деформування систем споруда-основа при силових і деформаційних навантаженнях, у тому числі, у випадках виникнення зон порушення контакту в процесі еволюції режимів деформаційних навантажень.

6. Створений універсальний математичний метод реалізований у вигляді алгоритмів і програм у сучасному промисловому програмному комплексі ЛІРА, що дозволяє ефективно і достовірно розв'язувати нелінійні задачі статики складних просторових систем споруд при силових і деформаційних навантаженнях.

7. Адекватність запропонованих математичних моделей і загального методу чисельного аналізу складних просторових систем споруда-основа підтверджена високою збіжністю результатів тестових прикладів з даними відомих експериментальних і натурних досліджень нелінійного деформування пластинчатих і стержневих систем та основних видів ґрунтів.

8. Виконано комплекс чисельних досліджень, зокрема, експериментальної 9-ти поверхової споруди, на основі яких виявлено суттєвий вплив при врахуванні найважливішого фактора нелінійності взаємодії споруди і основи на якісні і кількісні оцінки параметрів, що характеризують протікання статичних процесів при різних варіантах силових і деформаційних навантажень. Виявлено суттєві ефекти при врахуванні просторової схеми моделювання і процесу зведення споруд при деформаційних навантаженнях.

9. Встановлено загальні закономірності розподілу напружень і деформацій у просторовій системі споруда-основа при різних комбінаціях моделей і режимів навантажень. Показано, що застосування фізично нелінійних моделей споруди і основи дозволяє багаторазово знижувати напруження в системі і давати вичерпний прогноз змін структури статичних станів при варіюванні параметрів зовнішніх навантажень. Цей прогноз передбачає еволюційний розвиток напружено-деформованого стану системи споруда-основа в процесі силових і деформаційних навантажень та перетворення його у вихідний стан, що відповідає проектному положенню об'єкта дослідження.

10. Результати виконаних досліджень використані при розробці ПК ЛІРА у САПР, що використовується численними вітчизняними і зарубіжними організаціями, у навчальному процесі, при розробці нормативних документів, в процесі проектування, а також для проведення чисельних досліджень напружено-деформованого стану і автоматизованого розрахункового обґрунтування безпеки зведення, експлуатації і реконструкції особливо відповідальних об'єктів АЕС України в умовах екстремальних навантажень, в тому числі аварійних споруд у звичайних і складних ґрунтових умовах країни і зарубіжжя. Зокрема, при вирівнюванні автоматизованими домкратними системами житлових споруд: 14-ти поверхової у Москві, 9-ти поверхових у Тбілісі та Катовіце.

Публікації за темою дисертації

1. Клепиков С.Н., Слободян Я.Е. Взаимодействие балки с нелинейно-неупругим неоднородным основанием во всем диапазоне нагружения // Основания, фундаменты и механика грунтов. - 1989. - №4. - С. 22-24.

2. Клепиков С.Н., Слободян Я.Е. Расчет зданий при многократной подработке // Строительная механика и расчет сооружений. - 1990. - №2. - С. 62-65.

3. Программный комплекс ЛИРА-Windows. Теоретические основи / Городецкий А.С., Слободян Я.Е. и др. - Киев: НИИАСС. - 1997. -90 с.

4. Болотов Ю.К., Слободян Я.О., Мілявський В.Г., Садовий В.С. Нові технології при відновлюванні експлуатаційної придатності аврійного житлового будинку // Будівництво України. -1997. - №4. - С. 17-18.

5. Козачевский А.И., Слободян Я.Е. Оптимизация площади рабочей арматуры элементов железобетонных конструкций // Ресурсосберегающие технологии автоматизированного проектирования и управления в строительстве. - Киев: НИИСП. - 1987. - С. 45-50.

6. Клепиков С.Н., Слободян Я.Е., Козачевский А.И. Напряженно-деформированное состояние рамных каркасов с балочными фундаментами при воздействиях неравномерных оседаний основания // Строительные конструкции. - Киев. - 1988. - Вып. 41. - С. 31-35.

7. Клепиков С.Н., Слободян Я.Е. Расчет перекрестной системы на воздействия просадки основания // Основания и фундаменты. - Киев. - 1989. - Вып. 22. - С. 34-38.

8. Слободян Я.Е., Кисиль А.И. Распределение усилий в стенах высотного здания при подработке // Методы расчета и исследования зданий и оснований на подрабатываемых территориях и просадочных грунтах. - Киев: НИИСК. - 1990. - С. 12-16.

9. Слободян Я.Е. Расчет конструкций фундаментов на многократные воздействия горизонтальных деформаций основания // Строительные конструкции. - Киев. - 1990. - Вып. 43. - С. 21-27.

10. Слободян Я.Е. Расчет жилых зданий на нелинейно-деформируемом основании при карстовых проявлениях // Строительные конструкции. - Киев. - 1993. - Вып. 45-46. - С. 107-110.

11. Слободян Я.О. Експериментальне обґрунтування чисельного розв'язання нелінійно-непружних контактних задач взаємодії споруд з нерівномірно-деформованими основами // Опір матеріалів і теорія споруд. - Київ: КНУБА. - 2000. - Вип. 68. - С. 147-155.

12. Козачевский А.И., Слободян Я.Е. К вопросу оптимизации несущих пространственных стержневых систем сооружений на деформируемом основании // Конструкции гражданских зданий. - Киев: КиевЗНИИЭП. - 2001. - С. 61-67.

13. Слободян Я.О. Оцінка точності чисельного розв'язання нелінійних контактних задач взаємодії споруд з основами // Опір матеріалів і теорія споруд. - Київ: КНУБА. - 2002. - Вип. 70. - С. 136-144.

14. Козачевский А.И., Слободян Я.Е. К вопросу построения расчетных моделей систем «сооружение-основание» // Конструкции гражданских зданий. - Киев: КиевЗНИИЭП. - 2003. - С. 79-88.

15. Баженов В.А., Слободян Я.О. Моделювання нелінійного деформування складних просторових систем споруд при силових та деформаційних впливах // Опір матеріалів і теорія споруд. - Київ: КНУБА. - 2003. - Вип.72. - С. 41-73.

16. Баженов В.А., Слободян Я.О. Чисельні дослідження напружено деформованого стану споруд при ступінчатих деформаціях земної поверхні // Опір матеріалів і теорія споруд. - Київ: КНУБА. - 2003. - Вип..73. - С. 91-99.

17. Баженов В.А., Слободян Я.О. Застосування сучасних комп'ютерних технологій розрахунку і проектування для обгрунтування безпеки об'єктів атомної енергії // Опір матеріалів і теорія споруд. - Київ: КНУБА. -2004. - Вип..74. - С. 35-43.

18. Баженов В.А., Слободян Я.О. Автоматизація оптимального проектування просторових систем споруд // Опір матеріалів і теорія споруд. - Київ: КНУБА. - 2004. - Вип..75. - С. 96-100.

19. Слободян Я.О., Білеуш А.І, Євтодій О.О, Шаравара Н.П. Застосування інформаційних технологій автоматизації проектування для дослідження конструкцій на нелінійно деформованій основі // Вісник НАУ. - Київ. - 2004. - №4. - С. 131-133.

20. Слободян Я.О., Гузь А.В., Кудря Ю.С. Нові інформаційні технології автоматизації проектування і дослідження нелінійного деформування пластинчатих конструкцій // Вісник НАУ. - Київ. - 2004. - №4. - С. 128-130.