Научная статья на тему 'Твердотельная модель и конечно-элементный анализ в области суставного отростка нижней челюсти'

Твердотельная модель и конечно-элементный анализ в области суставного отростка нижней челюсти Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
178
38
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КОМПЬЮТЕРНА ТОМОГРАФИЯ / СУСТАВНОЙ ОТРОСТОК / РАСЧЕТНАЯ СХЕМА / ГЕОМЕТРИЧЕСКИЕ ПАРАМЕТРЫ / МЕХАНИЧЕСКИЕ ХАРАКТЕРИСТИКИ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Чуйко А. М., Погранична Х. Р.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Чуйко А. М., Погранична Х. Р.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Твердотельная модель и конечно-элементный анализ в области суставного отростка нижней челюсти»

© A.M. Чуйко, Х.Р. Погранична*

УДК 617.52

A.M. Чуйко, Х.Р. Погранична*

ТВЕРДОТІЛЬНА МОДЕЛЬ І СКІНЧЕННО-ЕЛЕМЕНТНИЙ АНАЛІЗ У ДІЛЯНЦІ СУГЛОБОВОГО ВІДРОСТКА НИЖНЬОЇ ЩЕЛЕПИ

М. Харків, Україна *Львівський національний медичний університет імені Данила Галицького (м. Львів)

Робота є фрагментом комплексної науково-дослідної роботи кафедри хірургічної стоматології та щелепно-лицевої хіруогії Львівського національного медичного університету імені Данила Галицького «Оптимізація діагностично-лікувального процесу хворих з кістковими і м’якотканинними дефектами та деформаціями різної етіології, травматичними і запальними ураженнями щелепно-лицевої ділянки», номер державної реєстрації 0110и008228.

Вступ. При сучасному рівні розвитку науки для вивчення таких складних систем, як нижня щелепа (НЩ) людини, найраціональніше використовувати твердотільне моделювання з подальшим звичайно-елементним аналізом [2].

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

Об’єкт і методи дослідження. Для зниження трудомісткості необхідних побудов проведемо

Рис.1. Віртуальна остеотомія щелепи.

віртуальну остеотомію, виділивши для аналізу тільки виростковий відросток правої гілки щелепи (рис.1).

Результати подальших побудов представлені на рис. 2,3,4. У кожній з моделей є свої особливості і переваги.

3D модель, представлена на рис.2,а,б якнайповніше відображає анатомію суглобового відростка, отриману на основі СТ після остеотомії.

Модель, представлена на рис.3,а,б, є стерео-литографічною 3D моделлю, побудованою програмою SolidWorks, на основі STL файлу, що згенеру-вав програмою MIMICS. При необхідності, файл STL може бути використаний при виготовленні реальної моделі STL за допомогою верстатів швидкого про-тотипирування. Тут відразу відзначаємо, що замість цієї трудомісткої і дорогої реальної моделі STL, в [2] рекомендовано використовувати STLV - віртуальну стереолитографічну модель, яка і точніша, і не вимагає спеціального устаткування, матеріалу і додаткового персоналу.

Модель на рис.4,а,б - це так звана, твердотільна модель (по термінології SolidWorks), побудована програмою SolidWorks, на основі STL файлу. Цим моделям можуть бути привласнені реальні механічні властивості кісткових і м’яких тканин і в рамках інтегрованої з SolidWorks програми COSMOSWorks, побудована скінченно-елементна модель, що дозволяє проводити статичні, динамічні, температурні і ін. дослідження.

Підкреслимо, що точність моделей, представлених на рис.2, визначається технічними характеристиками і налаштуваннями томографа. Точність моделей, представлених на рис.3, визначається налаштуваннями програми MIMICS. На рис.3 використана опція LOW, тобто найгрубіше налаштування.

Точність моделей, представлених на рис.4, визначається об’ємом «ручної роботи» по пере-строюванню STL моделі при створенні твердотільної моделі на екрані комп’ютера і залежить від досвіду і кваліфікації оператора. Тут можливі будь-які доопрацювання і виправлення моделей, представлених на рис.2,3, що і продемонстроване при моделюванні суглобової головки (рис.4).

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

Рис.2. 3D модель суглобового відростка справа (а) і зліва (б).

Рис.3. Віртуальна STL модель суглобового відростка справа (а) і зліва (б).

Рис.4. Твердотільна модель суглобового відростка справа (а) і зліва (б).

Рис.5. Внутрішні порожнини STL,

На рис.5 показано внутрішні порожнини суглобового відростка, в яких міститься губчаста кістка.

Моделі, представлені на рис.2-4, служать не тільки основою для твердотільного моделювання, але і, при переломі суглобового відростка, можуть бути використані для підгонки стандартного імплантата і проектування індивідуального.

в яких міститься губчаста кістка.

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

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

Рис.6. Скінченно-елементна сітка, граничні умови і один з варіантів навантаження.

Рис.7. Майданчики прикладання навантаження.

Рис.8. Циліндровий стрижень, як проста модель суглобового відростка.

остеотомія. На рис.6 відповідні граничні умови відмічені зеленими стрілками.

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

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

Розглянута вище методика визначення реакцій в суглобовій головці базується на плоскій розрахунковій схемі з урахуванням навантаження, що виникає при відкусуванні їжі. Природно, в житті спектр навантажень по величині і напряму різноманітніший і повноцінний напряженно-деформированное стан (ПДС) слід розглядати на базі 3D моделі. Але для глибшого і чіткішого розуміння подальших міркувань і результатів спочатку розглянемо простий циліндровий стрижень постійного перетину діаметром d і довжиною l з незмінними механічними властивостями (рис.8).

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

F

а' = a ■ <’>

де F - зовнішня сила (реакція в суглобовій головці);

A - площа поперечного перетину тіла.

Переміщення, в цьому випадку, визначаються по формулі

а

F

-A ■ <2)

У формулі (2) Е - модуль пружності матеріалу; А - площа поперечного перетину стрижня; І - довжина елементу. Крім того, утвір ЕА - називається жорсткістю стрижня на розтягування (стискування).

При дії поперечної сили елемент починає згинатися, як показано на рис.8,б. Нормальна напруга при вигині визначається по формулі

М

= -уУ> (3)

де Л - момент інерції поперечного перетину.

Напруга, що виникає при цьому, змінюється по лінійному закону, як по висоті поперечного перетину балки (залежно від координати у), так і по довжині залежно від величини моменту. Для консольної балки (рис.8, б) величина моменту М, що вигинає = Р(І-г), тобто залежить від плеча додатку сили.

Підкреслимо, що нормальна напруга від вигину на осі балки і на кінці, де прикладена сила, по всьому перетину, дорівнюють нулю (див. епюри на рис.8,б). Максимальна нормальна напруга від вигину виникає в жорсткому закріпленні балки і для круглого поперечного перетину можуть бути визначені по формулі

M

Fl

W it- d3 /32 ’ (4)

де враховано, що для круглого перетину момент опору

J ж- d3

W =

32

При сумісній дії вертикальної і горизонтальної складової, сили напруження складаються алгебраїчно (з урахуванням знаку), як показано на рис.8,в. При однаковій величині вертикальної і

горизонтальної складової сили, легко отримати співвідношення між максимальною напругою при вигині і стискуванні. В цьому випадку

4 = 8 1-. (5) а

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

Переміщення при поперечному вигині консольної балки (рис.8,б) можуть бути визначені по формулі

13

(6)

Визначимо співвідношення переміщень від вигину Аи і розтягування (стискування) Ас при дії однаковій по величині силі Р, скориставшись формулами (2), (4) і (6) (див. рис.8 а, б)

ґ і \2

• (7)

А = ■

А.

І

V у

16

Ас 7

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

При дії поперечної сили (бокової) виникає зрушення і поперечний вигин, що викликають дотичні напруження т, які при чистому зрушенні визначаються по формулі, аналогічній (1)

А

Т = —. (8)

А

При вигині, окрім нормальної напруги (3), виникає і дотична напруга, що розраховується по складнішій формулі Журавського, яку ми наводити не будемо.

При чистому крученні (рис.8,г) виникає дотична напруга тк що визначається по формулі

=

(9)

де Мк і Wk - величина моменту, що крутить, і момент опору крученню відповідно.

Особливий інтерес при крученні викликає напрям напруги, що обчислюється за формулою (9). При чистому крученні виникає так звана головна напруга, рівна по величині дотичній напрузі, але направлені під кутом 45° до осі елементу, як показано стрілками на рис.8,г. Досліди показують, що крихкі матеріали руйнуються по гвинтовій поверхні, нахиленій до осі під кутом 45°. Якщо ця напруга превалює в ПДС, то

перелом кістки буде косим. Це просте пояснення виникнення косих переломів в літературі не отримало достатнього освітлення.

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

аш = 1/о,5[(стх-ау)2 + (аг)2 + (ах-ах)2 + 6(^2 +т^2 +тА)] (10)

і характеризують загальний напружений стан у певній в точці. Сучасні програми, що реалізовують метод скінченних елементів (МСЕ), розраховують еквівалентну напругу в автоматичному режимі. У формулі (10) через ст і т позначені відповідно нормальна і дотична напруга, а індекси при них - напрями дії напруги уздовж осей х, у і и.

Всі позначення у формулі (10) - звичайні в україномовній літературі. У програмі COSMOSWorks, яка буде використана для проведення досліджень, ці ж величини позначаються таким чином: напруга по Мізесу уМ, як SM, вертикальна напруга уу, як SY і тому подібне.

Формула (1) служить основою для формулювання умови міцності, яка, після використання наведених формул, може бути розглянута конкретніше

СТс = А ~аТ , (11)

Властивості кістки конкретного пацієнта - це права частина нерівності (11). їхнє визначення залежить від статі, віку і особливо захворювання (пародонтит, остеопороз тощо) представляє не менш складне і актуальне завдання, яке до теперішнього часу також досліджене недостатньо [1].

Співвідношення (11) дозволяє вирішити наступні практично важливі завдання, що виникають при біомеханічному аналізі:

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

(12)

де стм - напруга, що діє або еквівалентні (напруга по Мізесу);

стТ - травмуюча напруга для кістки (межа текучості або межа втомної міцності для металів).

Відзначимо, що напруження по Мізесу (еквівалентні), визначувані по формулі (10), по своєму механічному сенсу припускають порівняння отриманої величини з межею текучості матеріалу. Тому тут і далі ми припускатимемо, що травмуюча напруга і межа текучості, ідентичні поняття.

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

к=?ї-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

. (13)

<7

M

Так само, як і напруга по Мізесу, програма COSMOSWorks обчислює це співвідношення в автоматичному режимі, позначаючи його як factor of safety, тобто

k = FOS. <13'а>

Легко побачити, що при k >1 умова міцності (11) задовільняється, при k <1, умова міцності не задовільняється і необхідно або зменшувати навантаження, або збільшувати площу поперечного перетину елементу, або підвищити міцність матеріалу.

2. Провести проектувальний розрахунок елемента системи, тобто при відомому максимальному навантаженні і травмуючому напруженні визначити розміри поперечного перетину

F (14)

A > max

(JT

3. Визначити навантаження, що допускається, при відомій площі поперечного перетину і міцності матеріалу

(15)

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

Ці три взаємозв’язані поняття, що виникають при біомеханічному аналізі, фактично відображають «умову міцності», математична і механічна суть якої визначена співвідношенням (11). Відзначимо, що

співвідношення (11) і (12) справедливі тільки для осьового навантаження за умови його рівномірного розподілу по перетину. У решті випадків виникатимуть зони концентрації напружень, які, при сучасному рівні розвитку науки можуть бути найуспішніше проаналізовані за допомогою МКЕ, що реалізовується на ЕОМ.

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

Після всіх попередніх міркувань приступимо до безпосереднього випробування і аналізу звичайно-елементної моделі виросткового відростка (рис.6).

Реакції в суглобовій головці визначалися залежно від величини жувального навантаження. Такий підхід природний при розгляді конкретного клінічного випадку. При загальному аналізі НДС в зоні суглобового відростка, для забезпечення порівнянності різних по напряму навантажень на суглобову головку, у всіх даних розрахункових випадках, прикладатимемо зовнішню силу Р = 100 Н. Такий підхід цілком виправданий, оскільки завдання вирішується в лінійній постановці і результат при будь-якій зовнішній силі може бути перерахований по відношенню до конкретного значення.

Як наголошувалося вище, найбільш повно напружений стан відображає напруження по Мізесу (10). Як випливає із змісту цієї формули, ця залежність завжди позитивна, що і відображається на діаграмі напруження. Тому для поглиблення аналізу, окрім напруження по Мізесу ^М), використовуватимемо і інші компоненти НДС, які складають основу напружень по Мізесу, наприклад, вертикальне напруження SY (або SZ залежно від орієнтації осей координат). Це напруження обчислюється з урахуванням знаку: «плюс» при розтягуванні і «мінус» при стискуванні, що може служити додатковим чинником при аналізі. Знак цього напруження відображається на діаграмі напруження.

Рис.9. Поле напружень по Мізесу (а) і полю вертикального напруження SZ (б).

Рис.10. Поля переміщень в недеформованому (а) і деформованому стані (б) суглобового відростка.

Спочатку розглянемо найпростіший розрахунковий випадок, коли стискуюча сила Р = 100 Н діє строго вертикально (рис.8,а, рис.9) на майданчику, практично співпадаючому з віссю суглобової головки. На рис.9,а представлено поле напружень по Мізесу, а на рис.9,б поля вертикального напруження SZ. Легко встановити, що конфігурація цих полів практично збігається як якісно, так і кількісно. Величина максимального значення SM = 5,755 Мпа, а SZ = -5,271 Мпа.

Не дивлячись на те, що результати розрахунків по формулах і скінченно-елементного аналізу можна зіставляти тільки якісно, спробуємо їх зіставити і кількісно. Припустимо, погодившись з рис. 2, 4 і 9, що з суглобової головки виділений стрижень з поперечним перетином у формі еліпса з розмірами ак = 7 мм і Ьк = 3,5 мм, площею Ак =19,24 мм. Тоді по формулі (1) отримаємо, що напруженняпри стисканні стс = - 5,19 МПа.

На рис.10 наиведені поля сумарних переміщень в недеформованому і деформованому стані суглобового відростка, значення яких також будуть використані при аналізі. Програма дозволяє виводити на друк всі компоненти переміщень по осях Х^ і Z, але для скорочення об’єму статті аналізуватимемо тільки сумарні переміщення УЯ. Оскільки у всіх випадках для аналізу використовується «нормативне» навантаження Р = 100 Н, те значення УЯ, що ділиться на 100, фактично виражає і величину податливості

системи - показник, який в [2] широко використовувався для аналізу рухливості зубів і імплантатів різного типу.

Як випливає з рис.10, точки суглобової головки переміщаються не тільки вертикально, але присутні і невеликі горизонтальні переміщення. Це вказує, що вертикальна сила не повністю збігається з умовною віссю головки.

Для зручності порівняльного аналізу всі числові значення основних параметрів ПДС зведені в табл. Базою для порівняння набуті значень напруги і переміщень при осьовому стискуванні.

У другій колонці табл. наводиться розташування майданчика (центральне або медіальне), до якого прикладається навантаження; у третій колонці -орієнтація навантажень (вертикальна, сагиттальная, бічна); у четвертій - значення максимального напруження по Мізесу SM (10); у п’ятій - значення максимальної вертикальної нормальної стискуючої напруги SZ (1) або дотичної напруги SXY (8,9); у шостій - відносна величина напруги (5); у сьомій

- коефіцієнт запасу міцності до (13); у восьмій -сумарні переміщення УЯ і в дев’ятій - відносна величина переміщень (7).

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

Таблиця

Основні параметри НДС, даних розрахункових випадків

№ Площадка Напрям Р БМ МПа(10) Б2( БХУ) МПа(1,8,9) Хст (5) к(13) иЯмм ХА (7)

1 Централ. Вертикал. 5,755 -5,271 1 37,24 2,98-10-3 1

2 Централ. Сагіттал. 31,03 19,24 17,11 5,39 6,9

3 Централ. Боковий 32,21 -30,64 5,59 6,65 1,23 -10-1 41,2

4 Медіальн. Вертикал. 22,31 9,06 -22,2 3,88 9,6 4,78-10-2 16

5 Медіальн. Сагіттал. 27,75 5,158 4,82 7,72 5,29-10-2 17,7

6 Медіальн. Боковий 31,5 30,6 5,47 6,8 1,18 -10-1 39,7

Рис.11. Напрям навантаження (а) і поле напруги по Мізесу (б).

жорсткості гілки, що відповідає дії жувального навантаження в нормі.

На рис.11 показаний напрям навантаження і поле напруги по Мізесу.

На рис.12 наведено основні компоненти напружень по Мізесу, рис.25,а представлене поле напружень по Мізесу, а на рис.25,6 поле вертикального напруження 87.: вертикальне напруження 87 в розтягнутій і стислій зонах.

Порівняння даних на рис.9 і 11 показує, що зона максимального напруження з головки перемістилася в шийку суглоьового відростка, що легко можна бути пояснити при порівняльному аналізі формул (1) і (3,4). На рис.12,б добре видно поле концентрації напруження в зоні різкої зміни геометрії. На цій особливості тут і далі ми зупинятися не будемо, оскільки основна мета роботи аналіз НДС виросткового відростка в нормі.

Наступним розрахунковим випадком є дія поперечного навантаження перпендикулярно до сагиттальной площини, тобто в напрямі, де вигинаюча жорсткість гілки НЩ мінімальна. На рис.13,а представлено поле напруженнь по Мізесу, а на рис.13,б поле вертикального напруження 87. Тут слід звернути увагу, що зона максимального напруження, якщо не враховувати концентрацію напруження, змістилася в ділянку нижче за шийку головки, туди, де вигинаючий момент, ближчий до максимуму, розташовуючись під невеликим кутом до вертикальної осі гілки. Абсолютна величина напруги збільшилася неістотно.

На рис.14 наведено поле сумарних переміщень ия.

По проведеному короткому аналізу трьох розрахункових випадків, цифрові дані по яких наведені у перших трьох стовпчиках табл., можна зробити додатковий порівняльний аналіз. Оскільки прийняте для аналізу навантаження Р = 100 Н співмірна з

Рис.12. Вертикальне напруження ЭЕ в розтягнутій (а) і стислій зонах (б).

Рис.13. Поле напружень по Мізесу БМ (а) і вертикальному напруженню Б2 (б).

Рис.14. Поле сумарних переміщень УР

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

Як відомо [1] «Нижня щелепа людини здійснює рухи в трьох напрямах: вертикальному (вгору і вниз), що відповідає закриванню і відкриттю рота; сагиттальном (вперед і назад); трансверсальному (управо і вліво). Кожен рух нижньої щелепи відбувається при одночасному ковзанні і обертанні суглобових головок. Різниця полягає лише в тому, що при одному русі в суглобах переважають шарнірні рухи, а при іншому - що ковзають». Навантаження при цих рухах можуть бути самі різні як по величині, так і по

Рис.15. Напрям навантаження (а) і поле напрження по Мізесу (б).

Рис. 16. Вертикальне напруження БЕ в розтягнутій (а) і стислій зонах (б).

напряму особливо при різноманітних побутових і спортивних ударах.

Для моделювання цієї особливості вище визначений майданчик, розташований медиальныше, ніж центральний, до якого прикладатимуться навантаження, величина яких як і раніше рівна Р = 100 Н. Послідовність проведення досліджень збережемо

- спочатку розглянемо вертикальну складову навантаження. На рис.15 показаний напрям навантаження і поле напруження по Мізесу.

На рис. 16 показано вертикальне напруження 87 в розтягнутій і стислій зонах.

На рис.17 наведено поле сумарних переміщень ия.

Отримані результати для цього розрахункового випадку (рис.15,а) по термінології опору матеріалів слід віднести до позацентрового стискування. Окрім осьової сили, тут виникає вигинаючий момент, постійний по довжині і рівний відтворенню сили на плече - відстані між центрами майданчиків. На рис.8 цьому виду навантаження відповідають

Рис.17. Поле сумарних переміщень УЯ.

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

Як і вище, наступним розглянемо розрахунковий випадок, коли навантаження на майданчику, розташованому медіальніше, направлене паралельно сагітальній площині (рис.18,а). При цьому виді вантаження на об’єкт одночасно діють вигин (рис.8,б), зрушення і кручення (рис.8,г). їх сумарну дію відображає напруження по Мізесу, представлені на рис.18,б. Дотичне напруження від кручення представлене на рис.19,а, а на рис.19.б - поле сумарних переміщень.

Порівнювати отримані результати слід, насамперед, з результатами другого рядка табл. і рис.11, 12. На рис.19,а, де представлено поле дотичного напруження від кручення, чітко видно їх «гвинтовий» характер, навіть не дивлячись на складність геометрії в цій зоні. Механіка цього явища проілюстрована на рис.8,г. Таким чином, якщо крутильная складова НДС буде істотною, в цьому випадку, слід чекати виникнення косого перелому.

Реверс навантаження (рис.20,а) загальну картину ПДВ, по відношенню до рис.18,б і 19,а не змінив.

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

Таким чином, наведений у статті скінченно-елементний аналіз основних параметрів НДС в ділянці суглобовового відростка дозволив виявити його основні особливості залежно від анатомічних особливостей конкретного пацієнта, а також величини, місця прикладання і напряму навантаження. Ці дані дозволяють прогнозувати можливість виникнення і характер травми і повинні обов’язково враховуватися при проектуванні і установці фіксуючих пристроїв, як стандартних, так і виготовлених індивідуально.

Рис.18. Напрям навантаження (а) і поле напружень по Мізесу (б).

Рис.19. Поле дотичного напруження БХУ (а) і сумарних переміщень ия (б).

Рис.20. Напрям навантаження в сагітальному (а) і бічному напрямі (б).

Висновки і практичні рекомендації.

1. Розроблена в [2] і адаптована в [3] для щелепно-лицьової хірургії, методика реконструкції травмованої щелепи без розтину операційного поля, як елементу технології СЇ/САО/САМ/САЕ, дозволила розглянути завдання міцності і жорсткості НЩ в зоні суглобового відростка конкретного пацієнта в нормі.

2. У основу біомеханічного аналізу покладена умова міцності або аналогічна йому умова жорсткості, які взаємно зв’язують КОНСТРУКЦІЮ (ЩЕЛЕПА В НОРМІ, ПРИ ПАТОЛОГИ АБО ПІСЛЯ ТРАВМИ) - ВЛАСТИВОСТІ КІСТКОВИХ ТКАНИН (КОНСТРУКЦІЙНИХ МАТЕРІАЛІВ) - НАВАНТАЖЕННЯ. У цьому триєдиному взаємозв’язку складових кожної з умов при розгляді конкретного пацієнта в нормі або після травми перші дві - КОНСТРУКЦІЯ і ВЛАСТИВОСТІ КІСТКОВИХ ТКАНИН аналізуються на основі даних СТ пацієнта. НАВАНТАЖЕННЯ, її

величина і напрям, визначається на основі статистичних даних.

3. Проведений скінченно-елементний аналіз основних параметрів НДС в зоні суглобового відростка, який дозволив виявити його основні особливості залежно від анатомії (конструкції) виросткового відростка конкретного пацієнта, величини, місця прикладання і напряму навантаження. Ці дані дозволяють прогнозувати можливість виникнення і вид травми і повинні обов’язково враховуватися при проектуванні і установці імплантатів, як стандартних, так і виготовлених індивідуально.

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

Список літератури

1. Гаврилов Е.И. Ортопедическая стоматология: Учебник 3-е изд., перераб. и доп. / Е.И. Гаврилов, А.С. Щербаков. - М. : Медицина, 1984. - 576 с., ил.

2. Чуйко А.Н. Биомеханика в стоматологии : Монография / А.Н. Чуйко, И.А. Шинчуковский. - X. : Изд-во «Форт». 2010. -516 с.

3. Чуйко А.Н. Компьютерная томография и биомеханическое сопровождение в челюстно-лицевой хирургии / А.Н. Чуйко, Д.К. Калиновский, К.Р. Пограничная // Ортопедия, травматология и протезирование. - 2011. - № 2, Ч. 1. - С. 23-35.

УДК 617.52

ТВЕРДОТІЛЬНА МОДЕЛЬ І СКІНЧЕННО-ЕЛЕМЕНТНИЙ АНАЛІЗ У ДІЛЯНЦІ СУГЛОБОВОГО ВІДРОСТКА НИЖНЬОЇ ЩЕЛЕПИ

Чуйко А.М., Погранична Х.Р.

Резюме. На основі аналізу уточненої розрахункової схеми, на базі даних комп’ютерної томографії пацієнта, яка дозволяє врахувати як геометричні параметри, так і механічні властивості кісткових і м’яких тканин, проведено аналіз міцності і жорсткості суглобового відростка нижньої щелепи в залежності від анатомії (конструкції), величини, місця прикладання і напряму навантаження.

Ключові слова: комп’ютерна томографія, суглобовий відросток, розрахункова схема, геометричні параметри, механічні характеристики.

УДК 617.52

ТВЕРДОТЕЛЬНАЯ МОДЕЛЬ И КОНЕЧНО-ЭЛЕМЕНТНЫЙ АНАЛИЗ В ОБЛАСТИ СУСТАВНОГО ОТРОСТКА НИЖНЕЙ ЧЕЛЮСТИ

Чуйко А.М., Пограничная К.Р.

Резюме. На основе анализа уточненной расчетной схемы, на базе данных компьютерной томографии пациента, которая позволяет учесть как геометрические параметры, так и механические свойства костных и мягких тканей, проведен анализ прочности и жесткости суставного отростка нижней челюсти в зависимости от анатомии (конструкции), величины, места приложения, и направления нагрузки.

Ключевые слова: компьютерна томография, суставной отросток, расчетная схема, геометрические параметры, механические характеристики.

UDC 617.52

Solid-State Model And Terminating Ultimate Analysis In The Field Of The Mandible Joint Shoot

Chuyko A.M., Pogranychna K.R.

Summary. On base of the analysis elaborated accounting scheme, on database computer tomography patient, which will allow to take into account both geometrics parameters, and mechanical characteristic bone and soft fabric, is organized analysis to toughness and acerbity mandibular condyle depending on anatomies (the designs), values, place of exhibit and directions of the load.

Key words: computer tomography, mandibular condyle, accounting scheme, geometrics parameters, mechanical characteristic.

Стаття надійшла 19.04.2012 р.

Рецензент - проф. Король М.Д.

i Надоели баннеры? Вы всегда можете отключить рекламу.