Российский журнал биомеханики

В XXI веке J. Wolff отметил, что кость здорового человека или животного адаптируется к тем нагрузкам, которым подвергается. На уровне ткани в кости просматриваются слабопористые участки компактного вещества и сильнопористые - трабекулярного (губчатого). Известно, что в трабекулярной костной ткани механизм адаптации реализуется посредством выстраивания трабекул (костных балочек, образующих твердый матрикс) вдоль линии действия преобладающей нагрузки. При достижении трабекулярной тканью оптимальной структуры для существующего в локальной области нагружения кость переходит в состояние равновесия (гомеостаза). В своих работах S. Cowin предложил описывать положение трабекул в каждый дискретный момент времени главными направлениями тензора структуры, отыскиваемыми из решения системы кинетических уравнений. На основе предложенных соотношений авторами решены многие задачи динамики трабекулярной костной ткани, результатом которых являются функции некоторой величины, зависящей от времени. Например, рассмотрено изменение пористости, компоненты девиатора тензора структуры. В данной работе кинетические уравнения используются с целью определения структуры, напряженного состояния или упругих свойств костного образца, находящегося в состоянии равновесия. Математическая постановка представлена задачей теории упругости анизотропного тела. Все численные расчеты выполнены с использованием программного продукта ANSYS Mechanical APDL на примере растяжения прямоугольной пластины. Проведено сравнение результатов задачи, полученных аналитическим методом и методом конечных элементов. Предполагается, что полученные соотношения могут быть полезны для определения реализуемого в кости известной структуры напряженно-деформированного состояния или, наоборот, для предсказания оптимальной структуры кости при заданных условиях нагружения.

Статья посвящена исследованию нестационарного течения и оседания запыленного воздуха (многофазной смеси газов и твердых частиц) в крупных нижних дыхательных путях в рамках разработки математической модели дыхательной системы человека. Исследуется содержание во вдыхаемом воздухе взвешенных частиц PM 10 (диаметр частиц менее 10 мкм) оксида алюминия на уровне норматива предельно допустимых среднесуточных концентраций. Несущая фаза описывается моделью вязкой жидкости; твердые частицы являются отдельными включениями различного размера, для которых записываются отдельные уравнения; полагается, что твердые частицы имеют сферическую форму. Расчеты параметров течения были выполнены с использованием программного продукта ANSYS CFX . Получены параметры течения газовой смеси с твердыми частицами, приведены поля скоростей и траектории движения твердых частиц различных размеров во время вдоха в различные моменты времени. Результаты, полученные с помощью подмодели движения воздуха в крупных нижних воздухоносных путях человека, являются входной информацией для моделирования движения газовзвеси в легких человека. Кроме того, полученные результаты могут быть использованы для решения задач по оценке риска здоровью человека от действия взвешенных частиц, осевших в нижних воздухоносных путях человека. Дальнейшее развитие модели предполагает совместное решение задач газодинамики в дыхательных путях человека и распространения воздуха в упругодеформируемой насыщенной пористой среде легких человека.

Целью исследования является создание методики построения расчетной конечно-элементной модели по данным томографии. Для оценки модели проводились расчеты бедренной кости. Актуальность данного исследования продиктована необходимостью считывать в расчетах распределение механических свойств костного органа по объему, что позволяет индивидуализировать подход к моделированию. Численные исследования выполнены с помощью метода конечных элементов в пакете Ansys , обработка данных компьютерной томографии проводилась в пакете Avizo . В работе рассматривалась задача линейного негомогенного упругого тела. Для определения модуля Юнга и предельного напряжения использовались степенные функции от оптической плотности, в свою очередь оптическая плотность определялась по линейным соотношениям в зависимости от чисел Хаунсфилда. Для конечно-элементной дискретизации производилось распределение механических свойств материала по каждому элементу согласно данным томографии. После решения задачи напряженно-деформированного состояния в каждом узле определялся коэффициент запаса с учетом свойств материала по данным томографии. В работе для рассматриваемого образца были проведены расчеты трех модельных задач в приведенной постановке и для осредненных значений механических параметров. Численные результаты в данных задачах наглядно иллюстрируют значимые различия в результатах напряженно-деформированного состояния органа и позволяют судить о локальной прочности костной ткани.

Одним из ключевых факторов успеха дентальной имплантации является остеоинтеграция - срастание имплантата с костью. Чрезмерная микроподвижность на поверхности имплантат - кость при нагружении нарушает остеоинтеграцию. Возникает естественный вопрос: нельзя ли минимизировать микроподвижность за счет управления конструкционными характеристиками имплантата, в частности резьбы? В работе изучается влияние характеристик резьбы (профиля, глубины и шага) на стабильность дентальных имплантатов при немедленном нагружении (первичную стабильность), когда остеоинтеграция еще не произошла и на интерфейсе имплантат - кость нет полного сцепления. Изменение профиля резьбы моделировалось вариацией угла наклона сторон профиля с последовательным переходом от треугольной через трапециевидную и квадратную резьбу к шиповой. Минимумы максимальных перемещений на интерфейсе имплантат - кость получены для квадратного и шипового профилей. Их величины составляют микроны, что соответствует данным других работ. влияние глубины и шага резьбы изучалось на квадратном профиле. Установлено, что увеличение глубины резьбы с 0,1 до 0,4 мм (как и уменьшение шага резьбы от 2,0 до 0,4 мм) приводит к уменьшению осадки имплантата (глобальной подвижности) на 3-6% и к значимому (кратному) монотонному падению перемещений на интерфейсе. Таким образом, профиль резьбы оказывает существенное влияние на микроподвижность на поверхности имплантат - кость при нагружении: при одинаковых значениях окклюзионной нагрузки и наблюдаемой макроподвижности зуба он может менять микроподвижность в разы. Минимумы микроподвижности получены для квадратного и шипового профилей. Максимумы локальных перемещений возникают на апексном витке резьбы.

Разрывы передней крестообразной связки являются самым распространенным повреждением среди изолированных травм связочных структур коленного сустава, которое требует раннего оперативного вмешательства, особенно у молодых трудоспособных пациентов. Существует большое количество способов восстановления передней крестообразной связки, при которых формируют каналы в бедренной и большеберцовой костях, а трансплантат фиксируют при помощи канюлированных интерферентных винтов. Особый интерес представляют операции, когда в качестве трансплантата используются синтетические протезы на основе полиэтилентетрафталата. Задачей данного исследования явилось изучение жёсткости фиксации имплантата передней крестообразной связки при различных видах его фиксации интерферентными винтами. На примерах экспериментальной модели «имплантат - кость» было показано, что применение искусственных протезов передней крестообразной связки на основе полиэтилентетрафталата позволяет достичь изначальной жёсткости фиксации, которая даст возможность ранней полной нагрузки на коленный сустав. Наибольшей прочностью обладает система, когда имплантат фиксирован четырьмя интерферентными винтами. При фиксации имплантата двумя интерферентными винтами происходит его проскальзывание в большеберцовом канале, что несёт потенциальную опасность для больного в случае начала ранней реабилитации.