CC BY
62
I
Орипнальы досл1дження
Original Researches
Травма
УДК 616.717.51-089.28
МАКАРОВ В.Б.1, ЛЕВАДНЫЙ Е.В.2, СТРАФУН А.С.3
1ГУ «Специализированная медико-санитарная часть № 6», г. Днепропетровск 2Днепропетровский национальный университет
3Отдел микрохирургии и реконструктивно-восстановительной хирургии верхней конечности НИИ травматологии и ортопедии НАМН Украины, г. Киев
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КОНТАКТНЫХ НАПРЯЖЕНИЙ И ОБЪЕМА ДВИЖЕНИЙ В ЛОКТЕВОМ СУСТАВЕ ПРИ ПЕРЕЛОМЕ БЛОКА ПЛЕЧЕВОЙ КОСТИ
СО СМЕЩЕНИЕМ
Резюме. Авторами построена необходимая геометрическая модель локтевого сустава, соответствующая требованиям к проведению математического моделирования. В программной среде ABAQUS на основе построенной геометрической модели был смоделирован перелом 13В2 (AO/ASIF) дистального отдела плечевой кости с различными величинами смещения суставной поверхности блока плечевой кости — 2, 3, 4 и 5 мм. На основе геометрической модели создана расчетная конечно-элементная модель и выполнен численный эксперимент. Установлено, что с ростом величины смещения отломка плечевой кости растет величина контактного напряжения суставных поверхностей в плечелучевом суставе с одновременным уменьшением площади контакта. Контактные напряжения суставных поверхностей резко возрастают при смещении фрагментов на 4-5 мм. Установлена зависимость угла сгибания и величины смещения отломка в зависимости от роста контактных напряжений, которая показывает, что смещение отломка на 2-3 мм в случае перелома 13В2 (AO/ASIF) существенно не уменьшает угол сгибания в локтевом суставе, несмотря на повышение контактных напряжений до 20-30 % по сравнению с физиологическими напряжениями. Данный факт может свидетельствовать о возможности консервативного лечения и последующей активной реабилитации. Таким образом, смещение внутрисуставных отломков на 2-3 мм можно рассматривать как максимально допустимое смещение, которое не будет в значительной степени ограничивать движения в локтевом суставе.
Ключевые слова: перелом дистального отдела плеча, объем движений, площадь контактных напряжений, напряженно-деформированное состояние, конечно-элементная сетка.
Введение
Лечение внутрисуставных переломов дистального отдела плечевой кости — сложная и многогранная проблема. Переломы данной локализации являются сравнительно редкой травмой в структуре повреждений опорно-двигательной системы и составляют до 2 % от всех переломов [1]. При этих повреждениях чрезвычайно трудно достичь максимально возможного функционального результата. Спорными остаются вопросы показаний и противопоказаний к консервативному лечению, особенностей хирургического доступа, методов фиксации при открытой репозиции [2, 3]. Особую проблему данный вид повреждения представляет у пациентов старше 60 лет. Длительная иммобилизация у таких пациентов (более 14 суток) вызывает резкое ухудшение функции локтевого сустава и снижение потенциала для восстановления движений в минимально необходимом физиологическом объеме по Моггеу 100°
[4]. В связи с этим остается нерешенным вопрос, при каком минимальном смещении отломков, составляющих суставную поверхность, можно получить достаточную функцию локтевого сустава при активной реабилитации. Для ответа на данный вопрос необходимо изучить контактные взаимодействия между суставными поверхностями при моделировании внутрисуставного перелома дистального метаэпифиза плечевой кости. Ряд авторов исследовали контактные напряжения на трупном материале [5—7]. Несмотря на максимальную объективность физического эксперимента, весьма трудно учесть разницу между минеральной плотностью разных препаратов, а также невозможность повторения опыта на одном и том же материале из-за стрессо-
© Макаров В.Б., Левадный Е.В., Страфун А.С., 2015 © «Травма», 2015 © Заславский А.Ю., 2015
Рисунок 1. Геометрическая модель: а — локтевого сустава; б — внутрисуставного перелома дистального метаэпифиза плечевой кости, тип B2 (AO/ASIF)
вых перегрузок и его полного или частичного разрушения. Математическое моделирование, или численный эксперимент, лишено этих недостатков [8].
Цель исследования — с помощью математического моделирования определить зависимость между контактными напряжениями, а также площадью контакта в плечелучевом и плечелоктевом суставах и возможной амплитудой движений в локтевом суставе при внутрисуставном переломе дистального отдела плечевой кости с различной величиной смещения — 2, 3, 4 и 5 мм.
Материалы и методы
На основе данных, предоставленных компанией Zygote Media Group, Inc., U.S. (http://www.3dscience. com), построена необходимая геометрическая модель (рис. 1), отвечающая требованиям к проведению математического моделирования.
Для последующего анализа был выбран тип перелома 13B2 (AO/ASIF) дистального отдела плечевой кости. Переломы типа B дистального метаэпифиза плечевой
кости возникают при падении на кисть с разогнутым локтевым суставом, в сочетании с абдукцией предплечья. При этом нередко ломается и локтевой отросток. Перелом типа В2 может быть результатом прямого удара в область локтевого сустава, а именно медиального надмыщелка.
Этот тип перелома является внутрисуставным, что, безусловно, влияет на функциональный исход, и достаточно простым для моделирования смещения относительно суставной поверхности блока плечевой кости величиной 2, 3, 4 и 5 мм (рис. 2).
В программной среде ABAQUS на основе построенной геометрической модели для проведения математического моделирования (численного анализа) создана расчетная модель в виде конечно-элементной сетки, отражающей геометрию и физические характеристики тел; контактные элементы, позволяющие средствами нелинейного расчета отразить динамическое взаимодействие различных тел между собой; нагрузки и закрепления, уточняющие физическую картину работы расчетной модели и позволяющие определить объем движения в локтевом суставе при различных величинах смещения блока плечевой кости. Также конечно-элементная модель использовалась для прогноза площади контакта суставных поверхностей при различных смещениях отломков. Конечно-элементная сетка (рис. 3, 4) расчетной модели состоит из 396 607 элементов, 390 353 из которых являются C3D4 (4-узло-вой элемент, линейный тетраэдр).
Так как анализ данной задачи подразумевает расчет модели при различных углах поворота костей в локтевом суставе, то, чтобы учесть взаимное влияние разных тел, входящих в модель, задача ставилась как контактная, с учетом возможности возникновения контакта между отдельными участками контактирующих тел. Влияние связок и капсулы локтевого сустава в данном расчете не учитывалось.
Нагрузки и закрепления. Мы допустили, что положение локтевой и плечевой костей при пронации-супинации остается неизменным. Чтобы сохранить
Рисунок 2. Геометрическая модель локтевого сустава с переломом дистального метаэпифиза плечевой кости типа В2 (AO/ASIF): а — со смещением 2 мм; б — 3 мм; в — 4 мм; г — 5 мм
обратное влияние на лучевую и локтевую кость со стороны контактирующих участков зафиксированных костей, фиксация была выполнена в узлах конечно-элементной сетки, наиболее удаленных от мест контакта. Так как предполагается анализ результатов расчетов статических задач уравновешенных моделей при различных углах сгибания в локтевом суставе, то каждый расчет подразумевает наличие уравновешивающей системы сил, действующих на модель. Нагружение лучевой и локтевой костей предполагает получение уравновешенного состояния за счет осевого давления на диафизарную часть с одной стороны и реакции опоры со стороны головки и блока плечевой кости.
Таким образом, задавая осевую нагрузку на лучевую и локтевую кость в направлении локтевого сустава, зафиксировав при этом положение плечевой кости, мы получим установившееся статическое положение системы.
Для оценки осевой нагрузки на лучевую и локтевую кость примем, что эта нагрузка появляется вследствие отжимания на руках. Отметим, что лучевая кость принимает 70—80 % нагрузки [9]. Предположим, что вес человека составляет 100 кг:
О = 100 (кг), т.е. сила притяжения будет составлять:
F веса = О• 9,81 = 100• 9,81 = 981 (Н).
Рисунок 3. Конечно-элементная модель локтевого сустава
С учетом того, что при отжимании человек перераспределяет нагрузку от силы тяжести на 4 конечности, сила реакции опоры на одну руку будет составлять:
К = 981: 4 = 245,25 (Н).
Так как предплечье состоит из 2 костей, то сила, действующая вдоль одной лучевой и локтевой кости, будет составлять:
К лучевой = 245,25 • 0,7= 171,67 (Н); К локтевой = 245,25 • 0,3 = 73,57 (Н).
Следует заметить, что если величина осевой нагрузки на лучевую и локтевую кость приведет к достаточно большим напряжениям, то на практике это вызовет ее разрушение. Результаты расчета не смогут показать разрушение как таковое. Но причины, приводящие к дальнейшему разрушению, могут быть обнаружены. Поэтому после получения результатов расчетов и до начала анализа следует оценить уровень напряжений по отношению к максимально допустимым значениям, которые приведены в соответствующих свойствах материалов.
В расчетную модель входит три суставных поверхности плечевой, лучевой и локтевой костей. Кроме того, контактные участки моделировались с учетом снижения уровня модуля Юнга, отражая особенности наличия хрящевой ткани на контактных поверхностях костей. Расчет включал в себя анализ напряженно-деформированного состояния системы при углах сгибания в локтевом суставе 30, 60, 90 и 120° и в среднем положении между пронацией и супинацией.
Свойства материалов, а также характер нагружения модели для всех расчетов были приняты постоянными и одинаковыми (табл. 1, 2) [10, 11]. Величина коэффициента Пуассона для всех материалов была принята равной 0,3.
Результаты и их обсуждение
В результате расчета были получены контактные напряжения на суставных поверхностях локтевого сустава в норме, а также различные площади контакта суставных поверхностей и величины контактных напряжений в зависимости от угла возможного
Рисунок 4. Конечно-элементная модель локтевого сустава с переломом дистального метаэпифиза плечевой кости типа В2 (AO/ASIF): а — со смещением 2 мм; б — 3 мм; в — 4 мм; г — 5 мм
Таблица 1. Прочностные свойства костей
Локализация Модуль Юнга, ГПа Прочность на растяжение, МПа Прочность на сжатие, МПа
Плечевая кость 17,2 30 132
Лучевая кость 18,6 149 114
Локтевая кость 18,0 148 117
сгибания локтевого сустава при внутрисуставном переломе дистального метаэпифиза плеча типа 13В2 (рис. 5).
Анализ контактных площадей. Установлено, что площадь контакта между локтевой и плечевой костями в случае смещения отломков на 2 мм по линии перелома больше на 30 %, чем при смещении на 5 мм и угле сгибания 90° (рис. 5).
Анализируя данные табл. 3, следует отметить, что при переломе блока плечевой кости со смещением от 2 до 5 мм площадь контакта головки лучевой кости с головкой мыщелка плечевой кости изменяется незначительно в зависимости от угла сгибания (рис. 6), при этом чем больше угол сгибания, тем больше площадь контакта. При рассмотрении площади контакта блока плечевой кости с локтевым отростком увеличение смещения отломков приводит к уменьшению площади контакта, а увеличение угла сгибания уменьшает величину площади контакта (рис. 7).
Данные о величинах контактных напряжений на головке лучевой кости, локтевом отростке и плечевой кости при различных углах сгибания в локтевом суставе и различной величине смещения представлены на графиках (рис. 8—10).
Анализ данных контактных напряжений. Приняв среднефизиологическое контактное напряжение на суставных поверхностях в норме 5,5—8,5 МПа [12, 13], установили, что до угла сгибания 90° напряжение медленно растет. В среднем при смещении отломков на 2—3 мм и до угла сгибания до 90° рост контактных напряжений составляет 20—30 %. При угле сгибания 90° в локтевом суставе напряжения на суставной поверхности лучевой кости возрастают с 23 МПа при смещении отломков на 2 мм до 26 МПа, а при смещении на 5 мм
Таблица 2. Модуль Юнга хрящей суставных поверхностей
Локализация суставного хряща Модуль Юнга, ГПа
На головке плечевой кости 9,6
На головке лучевой кости 9,3
Со стороны локтевой кости 9,0
1
Таблица 3. Данные контактных площадей на головке лучевой кости и суставной поверхности плечевой кости в зависимости от угла сгибания и величины смещения отломков
Угол сгибания в градусах Площадь контакта на головке лучевой кости, мм2 Площадь контакта на суставной поверхности блока плечевой кости, мм2
Смещение отломков на 2 мм
0 52 117
30 55 108,225
60 57 96,525
90 60 72,54
120 63 65,52
Смещение отломков на 3 мм
0 53,56 104
30 56,65 96,2
60 58,71 85,8
90 61,8 64,48
120 64,89 58,24
Смещение отломков на 4 мм
0 54,08 91
30 57,2 84,175
60 59,28 75,075
90 62,4 56,42
120 65,52 50,96
Смещение отломков на 5 мм
0 54,6 78
30 57,75 72,15
60 59,85 64,35
90 63 48,36
120 66,15 43,68
Рисунок 5. Площади контакта при переломе дистального метаэпифиза плечевой кости: а — со смещением 2 мм и угле сгибания 90°; б — со смещением 5 мм и угле сгибания 90°; в — со смещением 5 мм и угле сгибания 90°; г — со смещением 5 мм и угле сгибания 90°
на локтевом отростке с 21 до 32 МПа, на плече — с 23 до 32 МПа (рис. 11). При угле сгибания 120° контактные напряжения возрастают с 33—37 до 45—50 МПа.
Таким образом, исходя из вышепредставленных данных, а именно из уменьшения площади контакта суставных поверхностей локтевой и плечевой костей при увеличении смещения отломков и угла сгибания в локтевом суставе, а также увеличения контактных
напряжений на суставных поверхностях, можно предположить, что ограничение угла сгибания происходит как из-за нарушения конгруэнтности суставных поверхностей, так и из-за нарастания болевого синдрома, вызванного повышением давления на суставных поверхностях.
Приняв контактное напряжение на суставных поверхностях, при котором начинается резкий рост
140
я 120
о >
я 100 а
» 80 J 60
о
40 20 0
50 100
Площадь (мм2)
150
— Смещение 2 мм
— Смещение 4 мм
■Смещение 3 мм Смещение 5 мм
0
Рисунок 6. Зависимость площади контакта головки лучевой кости от угла сгибания и величины смещения отломков
Рисунок 7. Зависимость площади контакта блока плечевой кости с суставной поверхностью локтевой кости от угла сгибания и величины смещения отломков
150
х
а
о >
д
а а
н
а
\о
л о
100
50
И I
20 40
Напряжение (МПа)
60
-ф— Смещение 2 мм — Смещение 4 мм
Смещение 3 мм Смещение 5 мм
0
140
я 120
о >
if 100 р
» 80 и
ан 60
л о
40 20 0
0
20
40
60
Напряжение (МПа) -ф—Смещение 2 мм —■—Смещение 3 мм —Смещение 4 мм —•—Смещение 5 мм
Рисунок 8. Зависимость величины контактного напряжения на головке лучевой кости от угла сгибания в локтевом суставе при различных величинах смещения блока плечевой кости
Рисунок 9. Зависимость величины контактного напряжения на локтевом отростке от угла сгибания в локтевом суставе при различных величинах смещения блока плечевой кости
140
я 120
с §
я 100 р
в я и н
а
\о
л о
80 60 40 20
И I
20 40
Напряжение (МПа)
60
Смещение 2 мм
Смещение 3 мм
—±—Смещение 4 мм —•—Смещение 5 мм
Рисунок 10. Зависимость величины контактного напряжения на суставной поверхности плечевой кости от угла сгибания в локтевом суставе при различных величинах смещения блока плечевой кости
контактных напряжений (в среднем 24 МПа), за начальное условие для построения зависимости угла сгибания и величины смещения отломка в зависимости от роста контактных напряжений, можно построить график (рис. 12). Из графика на рис. 13 видно, что смещение отломка в случае перелома В2 (AO/ASIF) на 2—3 мм значительно не уменьшает угол сгибания в локтевом суставе, что может косвенно свидетельствовать о возможности консервативного лечения (иммобилизации не более 14 суток) и последующей активной реабилитации.
Выводы
1. С помощью разработанной трехмерной конечно-элементной модели локтевого сустава возможно выполнять исследование напряженно-деформированного состояния как в норме, так и при моделировании перелома.
2. При рассмотрении площади контакта блока плечевой кости с локтевым отростком увеличение смещения отломков приводит к уменьшению площади контакта, а увеличение угла сгибания уменьшает площадь контакта.
3. Установлено, что площадь контакта между локтевой и плечевой костями в случае смещения отломков на 2 мм по линии перелома больше на 30 %, чем при смещении на 5 мм и угле сгибания 90°.
4. Из-за уменьшения площади контакта суставных поверхностей локтевой и плечевой костей при увеличении смещения отломков и угла сгибания в локтевом суставе, а также увеличении контактных напряжений на суставных поверхностях можно предположить, что ограничение угла сгибания происходит как из-за нарушения конгруэнтности суставных поверхностей, так и из-за нарастания болевого синдрома, вызванного повышением давления на суставных поверхностях.
5. Установлено, что с ростом величины смещения отломка плечевой кости увеличивается контактное напряжение суставных поверхностей в плечелучевом суставе с одновременным уменьшением площади контакта. Контактные напряжения суставных поверхностей резко возрастают при смещении фрагментов на 4-5 мм (до 45-50 МПа).
6. Установлена зависимость угла сгибания и величины смещения отломка от роста контактных напряжений, которая показывает, что смещение отломка на 2-3 мм в случае перелома 13В2 (AO/ASIF) значительно не уменьшает угол сгибания в локтевом суставе, несмотря на повышение контактных напряжений до 20-30 % в сравнении с физиологическими напряжени-
Рисунок 11. Контактные напряжения при угле сгибания в локтевом суставе 90°: а — при смещении 2 мм; б — при смещении 5 мм
0
Рисунок 12. Зависимость величины контактного напряжения на суставной поверхности плечевой кости от угла сгибания в локтевом суставе при различных величинах смещения блока плечевой кости
ями. Данный факт может косвенно свидетельствовать о возможности консервативного лечения и последующей активной реабилитации. Таким образом, смещение внутрисуставных отломков на 2—3 мм можно рассматривать как максимально допустимое смещение, которое не будет в значительной мере ограничивать движения в локтевом суставе.
Список литературы
1. Garcia J.A., Mykula R., Stanley D. Complex fractures of the distal humerus in the elderly // J. Bone Joint Surg [Br]. — 2002. — 84-B. — 812-6.
2. Wang Y., Zhuo Q, Tang P, Yang W. Surgical interventions for treating distal humeral fractures in adults // Cochrane Database Syst. Rev. — 2013 Jan 31.
3. Korner J., Lill H., Müller L.P., Hessmann M., Kopf K., Goldhahn J., Gonschorek O., Josten C., Rommens P.M. Distal humerus fractures in elderly patients: results after open reduction and internal fixation // Osteoporos Int. — 2005Mar. — 16, Suppl. 2. — S73-9. Epub 2004 Oct 29.
4. Morrey B. The elbow and its disorders / B. Morrey. — 3ed. — Philadelphia; London; New York: W.B. Saunders company, 2000. — 934p.
5. Shepard M.F., MarkolfK.L., Dunbar A.M. Effects of radial head excision and distal radial shortening on load-sharing in cadaver forearms // J. Bone Joint Surg. Am. — 2001 Jan. — 83-A(1). — 92-100.
6. Rao Z.T., Yuan F, Li B., Ma N. Effect of elbow flexion angles on stress distribution of the proximal ulnar and radius bones under a vertical load: measurement using resistance strain gauges// J. Orthop. Surg Res. — 2014 Jul. — 31, 9. — 60.
7. Eckstein F, Lohe F, Hillebrand S, Bergmann M., Schulte E, Milz S, Putz R. Morphomechanics of the humero-ulnar joint: I. Joint space width and contact areas as a function of load and flexion angle // Anat. Rec. — 1995 Nov. — 243 (3). — 318-26.
8. Edmonds E.W., Santago A.C. 2nd, Saul K.R. Functional Loss With Displacement of Medial Epicondyle Humerus Fractures: A Computer Simulation Study // J. Pediatr. Orthop. — 2014 Dec 9.
9. Pfaeffle H.J., Fischer K.J., Manson T.T., Tomaino M.M., Woo S.L., Herndon J.H. Role of the forearm interosseous ligament: is it more than just longitudinal load transfer? // J. Hand Surg. Am. — 2000. — 25. — 683-688.
10. Карлов А.В. Системы внешней фиксации и регулятор-ные механизмы оптимальной биомеханики / А.В. Карлов, В.П. Шахов. — Томск: STT, 2001. — 480 с.
11. Трушко А.В. АСМ-анализ хрящевых тканей / А.В. Трушко, С.А. Чижик: Мат-лы VIIМеждунар. семинара (1—3 ноября 2006 г., Минск). — Минск: ИТМО НАН Беларуси, 2006. — С. 162-168.
12. Бойко И.В. Биомеханическое изучение особенностей распределения нагрузок в головке и шейке лучевой кости / И.В. Бойко, В.Б. Макаров, Д.Е. Щербаков и др. // Всник ортопеда, травматологи та протезу-вання. — 2013. — № 1. — С. 22-27.
13. Kim. Sunghwan Contact Stress Analysis of the Native Radial Head and Radial Head Implants: Doctoral Dissertation. — University of Pittsburgh, 2014. —119p.
Получено 05.01.15 ■
Макаров В.Б.1, Левадний £.В.2, Страфун О.С.3
1ДУ «Спецалзована медико-сан'парна частина № 6»», м. Днпропетровськ
Днпропетровський нацюнальний унверситет
3В^л мкро^рургП'та реконструктивно-в'дновлювально!xipypni верхньоiюн^вки HAI травматологи та ортопедаНАМН Украни, м. ¡Мв
МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ КОНТАКТНИХ НАПРУЖЕНЬ ТА ОБСЯГУ PYXiB В ЛКТЬОВОМУ СУГЛОБi ПРИ ПЕРЕЛОМi БЛОКУ ПЛЕЧОВО1 ЮСТКИ i3 ЗМЩЕННЯМ В^ЛАМЮВ
Резюме. Авторами побудована необхщна геометрич-на модель лжгьового суглоба, що вщповщае вимогам до проведения математичного моделювання. У програмному середовищi АВАрИ8 на основi побудовано! геометрично! моделi був змодельований перелом 13В2 (АО/А8Ш) дис-тального вщдщу плечово! истки з рiзними величинами змщення суглобово! поверхш блоку плечово! истки — 2, 3, 4 i 5 мм. На основi геометрично! моделi створена розра-хункова инцево-елементна модель та виконано чисельний
експеримент. Встановлено, що 3i збшьшенням величини змщення уламка плечово! истки збиьшуеться величина контактного напруження суглобових поверхонь в плече-променевому суглобi з одночасним зменшенням площi контакту. Контактш напруження суглобових поверхонь рiзко зростають при змщенш фрагменив на 4—5 мм. Встановлено залежшсть кута згинання i величини змщення уламка залежно вщ зростання контактних напружень, яка показуе, що змщення уламка на 2—3 мм в разi перелому
13В2 (AO/ASIF) icTOTHO не зменшуе кут згинання в лжтьо-вому суглоб^ незважаючи на пщвищення контактних на-пружень до 20—30 % у порiвняннi з фiзiологiчними напру-женнями. Даний факт може свщчити про можливкть консервативного лiкування i подальшо! активно! реабииацп. Отже, змщення внутрiшньоcуглобових уламкiв на 2—3 мм
можна розглядати як максимально допустиме змщення, що не буде значною мiрою обмежувати рухи в лжгьовому cуглобi.
Ключовi слова: перелом дистального вщдшу плеча, обсяг ру-хш, площа контактних напружень, напружено-деформований стан, юнцево-елементна ciтка.
Makarov V.B.1, Levadnyi Ye.V.2, Strafun A.S.3
1State Institution «Specialized Industrial Hospital № 6»», Dnipropetrovsk
2Dnipropetrovsk National University, Dnipropetrovsk
3Department of Microsurgery and Reconstructive Surgery of the Upper Extremity of Research Institute of Traumatology and Orthopedics of National Academy of Medical Sciences of Ukraine, Kyiv, Ukraine
MATHEMATICAL MODELING OF CONTACT STRESS AND THE RANGE OF MOTION IN THE ELBOW JOINT IN DISPLACED
FRACTURE OF HUMERAL TROCHLEA
Summary. The authors constructed a necessary geometrical model of the elbow joint that meets sufficient requirements for mathematical modeling. In ABAQUS software, based on created geometric model, there has been simulated 13B2 (AO/ASIF) fracture of the distal humerus with different values of displacement of the articular surface of the humeral trochlea — 2, 3, 4 and 5 mm. On the basis of geometrical model, finite element model has been created and numerical experiment has been carried out. It was found that with increase in the displacement of the fragment of the humerus, the value of contact stress in the articular surfaces in the humeroradial joint increases with simultaneous reduction of the contact area. Contact stress on joint surfaces rise sharply at 4—5 mm displacement of frag-
ments. There has been shown a dependence of the angle of flexion and the value of displacement of fragments on the growth of contact stresses, which shows that the amount of displacement of fragment by 2—3 mm in case of13B2 fracture (AO/ASIF) did not significantly reduce the elbow angle, despite an increase in contact stresses up to 20—30 % compared with the physiological stresses. This fact may indirectly indicate a possibility of conservative treatment and further active rehabilitation. Thus, displacement of intraarticular fragments by 2—3 mm can be considered as maximal allowable displacement that will not greatly restrict the movements of the elbow joint.
Key words: distal humerus fracture, range of motion, contact stress area, stress-strain state, finite-element mesh.
