CC BY
35
УДК 531/534: [57+61]
Российский Журнал
www.biomech.ru
ОБ ОПТИМИЗАЦИИ УПРУГОГО ЭЛЕМЕНТА ПРОТЕЗА СТОПЫ
М.А. Осипенко, Ю.И. Няшин
Кафедра теоретической механики Пермского государственного технического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: [email protected]
Аннотация. Приведен обзор новых результатов по оптимизации конструкции упругого элемента (УЭ) протеза стопы. Протез моделируется в виде консольно закрепленной многолистовой рессоры. Между листами имеет место контакт с возможным отставанием; задача изгиба поэтому относится к контактным задачам теории балок. Найдена переменная толщина равнонапряженного однолистового УЭ при сильном изгибе. Показано, что такой УЭ имеет эллиптический профиль. Найдены переменные толщины листов равнонапряженного двухлистового УЭ при слабом изгибе для некоторых значений отношения длин листов. Показано, что при неравных длинах листов более длинный лист такого УЭ имеет внутренний шарнир. Доказано, что при дополнительных слабых ограничениях равнонапряженный УЭ может не существовать. Рассмотрен коэффициент использования материала УЭ. Сопоставлены задачи максимизации этого коэффициента и минимизации массы УЭ и выяснено, что связь этих задач не является простой. Обе задачи могут не иметь решений; первая задача может иметь решения, а вторая - не иметь; первая задача может иметь много решений, одно из которых является единственным решением второй.
Ключевые слова: протез стопы, упругий элемент, листовая рессора, линейный и нелинейный изгиб, контактная задача, равнонапряженный упругий элемент, коэффициент использования материала, минимизация массы.
Введение
Проектирование и оптимизация протеза стопы являются одним из актуальных направлений протезостроения. Протез стопы должен воспроизводить ряд функций естественной стопы - опорную, функцию амортизации переднего толчка, функцию накопления потенциальной энергии, рессорную (распластывание под действием вертикальных сил), балансировочную (регуляция позы), адаптационную (способность ходьбы по наклонным поверхностям). Протез стопы должен также обладать достаточной прочностью.
Перечисленным требованиям наиболее полно удовлетворяют конструкции, содержащие так называемый упругий элемент (УЭ), который представляет собой консольно закрепленную многолистовую рессору. УЭ может быть стальным. На рис. 1, а показан серийный образец, изготовленный на экспериментальном заводе средств протезирования, г. Реутово, Московская область. Такие протезы сейчас наиболее распространены. Однако более перспективными являются УЭ, изготовленные из легких композиционных материалов. На рис. 1, б показан опытный образец УЭ, изготовленный в Пермском НИИ композиционных материалов из углепластика КМУ-4Л. На рис. 1, в представлен УЭ протеза Carbon Copy, изготовленный фирмой
© Осипенко М.А., Няшин Ю.И., 2011
Осипенко Михаил Анатольевич, к.ф.-м.н., доцент кафедры теоретической механики, Пермь Няшин Юрий Иванович, д.т.н., профессор, завкафедрой теоретической механики, Пермь
09806267
м
а
б
в
Рис. 1. Разновидности упругого элемента протеза стопы
Ohio Willow Wood Company, США [6]. На рис. 1, г показан УЭ протеза стопы Quantum Foot, изготовленный фирмой Hanger Inter Med, Великобритания.
Такие конструкционно сложные элементы протеза, будучи более эффективными при эксплуатации, требуют применения более новых методов проектирования, в частности, метода математического моделирования, который позволяет также оптимизировать конструкцию протеза стопы. В настоящей статье изложен ряд результатов по оптимизации УЭ протеза стопы, полученных авторами в последние годы. Эти результаты представляют интерес для специалистов в области расчета рессор и общей теории листовых рессор.
Простейшая модель УЭ [2, 3, 7] представляет собой пачку N консольно закрепленных балок (листов) одинаковой ширины м>, но с различными длинами ¡х,.., ¡м
и различными переменными толщинами //(х),..., hN(х) (х - координата, отсчитываемая
вдоль балки) (рис. 2).
В естественном состоянии (без нагрузки) листы являются плоскими и плотно прилегают друг к другу; трение между листами отсутствует. К краю самого длинного листа, перпендикулярно листу, приложена сосредоточенная следящая сила F. Форма изогнутого листа (в равновесии) задается функцией ф(х), где ф - угол между профилями изогнутого и неизогнутого листов (рис. 3). Если лист изогнут, то в точке с координатой х со стороны правой (по отношению к данной точке) части листа на левую действует сосредоточенный момент
где Е - модуль Юнга. При этом напряжение (единственная существенная компонента а хх ) на нижней поверхности листа
Оптимальными считаются [1, 5] равнонапряженные УЭ, т.е. удовлетворяющие условию
Математическая модель УЭ
(1)
о(x) = Eh(х)ф'(х)/ 2.
(2)
\on (х) = а = const, 1 < n < N.
(3)
Рис. 2. Модель упругого элемента
Рис. 3. Схема изгиба одного листа
Равнонапряженность может быть достигнута подбором функций hn (х). Рассмотрим вопрос об отыскании этих функций в некоторых частных случаях.
Оптимальный профиль однолистового УЭ при нелинейном изгибе
Пусть УЭ однолистовой. Уравнение равновесия листа следует из (1) с учетом сделанного предположения о нагрузке (сосредоточенная сила):
Ewh\х )ф'( х)/12 = F J1 cos (ф(/) -ф(s)) ds. (4)
Из (2)-(4) и условия защемления листа (при х = 0) находим
|ф"/ф'3 = А cos (ф(/) -ф),
{ ф(0) = 0, ф'(/) = +»,
(5)
где А = 3Е2F^4ст3w. Нелинейность задачи (5) означает, что рассматривается сильный изгиб. Решение (5) существует (и единственно) при I/А < 2 и имеет вид
ф(x) = arccos (1 -I/А) - arccos (1 - (l - x)/А).
Тогда из (2), (3) получаем
h(x) = ( 2стА/E)^ 1 -(1 - (l - x)/А)2 ,
т.е. однолистовой равнонапряженный УЭ имеет эллиптический профиль. Практически при нахождении профиля задается не ст, а прогиб УЭ
г i
А = J sin ф(x)dx.
Тогда из (6) находим, что l/ А = 1 - cos а, где 0 < а < п - корень уравнения
(sin а -а cos а)/(1 - cos а) = 2 А/1.
(6)
(7)
(8)
Легко установить, что при Д/1 < я/4 это уравнение имеет единственное решение.
Если Д/1 << 1 (слабый изгиб; эквивалентное условие: |ф(х)| << 1), то из (7), (8)
получаем известный результат [5] - И(х) = 21 (^Ем ’Д)13 ^1 - ~х[1 - однолистовой равнонапряженный УЭ при слабом изгибе имеет параболический профиль.
Оптимальные профили двухлистового УЭ при линейном изгибе
Теперь рассмотрим двухлистовой УЭ при слабом изгибе. Предположим, что а) листы УЭ контактируют (кроме точки защемления) только в точке, расположенной на краю более короткого листа; б) стп (х) = ст (а не модуль, что требуется в (3)). Тогда уравнения равновесия листов УЭ имеют вид
Г h3( ) ' ( )/12 j-P(/2 - x) + F(/\ -x) (0 - x - 2X (9)
Ых{ ф\ (х)\1 = •! (9)
[ F (/\ - x) (/2 - x - AX
Ewh2(x^2(x)/\2 = P(/2 -x), (\0)
ф\ (0) = Ф2(0) = 0, (\\)
где P - сила взаимодействия листов, которая находится из условия контакта листов в точке x = /2 :
y\(/2) = y2 (/2 ) ; (\2)
<• x
Уп (x) = J0 Фп (s)ds . (\3)
Из (9)-(\3), (2), (3) (с учетом предположения (б)) можно получить уравнение для P в виде
5 + 2ß = (i; + ß)7Vb4r,
где ß = yj \ - /2//\ , 5 = yj \ - P/2 /F/\ . Легко установить, что при 0 - ß < \ это уравнение имеет единственное решение 5 . Однако находя из (9)-(\\), (2), (3) ф\(x), ф2(x) и затем из (\3) y\(x), y2(x) , можно показать, что при /2 < /\ limy2(x)/y\(x) < \ (верхний лист
x^0
«проникает» под нижний), а при /2 = /\ y\(x) = y2(x) . Таким образом, при /2 = /\ = / получаем одинаковые параболические профили листов
h\ (x) = h2 (x) = 2/ (F/2EwA)13 ^/\ - x// , а при /2 < /\ равнонапряженный УЭ, для которого выполняются принятые выше предположения (а) и (б), не существует.
Докажем, что при сохранении предположения (б) при /2 < /\ равнонапряженный УЭ не существует, даже если ослабить предположение (а) и допустить, что листы могут контактировать в произвольном конечном числе точек 0 < xM < ... < x\ < x0 = /2 < /\.
Учитывая, что при 0 - x - xM
Ewh\ (x^ (x)/\2 = F(/\ - x) - £ M=0 Pm (xm - x) ,
Ewh2 (X)Ф2 (x)/12 = Z M=0 Pm (xm - x) , ф\(0) = Ф2(0) = 0,
где Pm - силы взаимодействия листов, получим с учетом (2), (3), (\3) из условий
ф\ (xM ) = ф2(xM ) и y\ (xM ) = у2 (xM ) соответственно:
А\ + А2 = А3 + А4 , (A\ + 2А2 )К A\ + A2 ) =(А3 + 2А4 )/(A3 + А4 ) , (\4)
где
А= P(h - XM )-Z Mo Pm X “ XM ) , ¿2 = M Mo P^m , (15)
^3 = ^ Мт:1 Рт (хт - хм ), А = ^ !=0 Рт • (16)
Из (14)-(16) вздует, что ^Мт= 0рт =F/2 и 2М= 0Ршхш/11 =F/2. но эти два равенства
противоречивы, так как Рт > 0 и хт < /1. Поэтому высказанное выше утверждение доказано.
Откажемся тогда от предположения (б), сохранив (а). Из (9), (2) следует, что
если
Р/2 > F/l, (17)
то а1(х) меняет знак в точке х* = (Р/2 - F/1 )/(Р - F) . Тогда из (9)-(13), (2), (3) получаем уравнение для Р в виде
(1+л)2 =(2 -Тл(л +3) )>/л + У/ (у-1)
где у = /1//2 , л = (Р^ — у')/(у-1) . Легко установить, что при у> 1 это уравнение имеет единственное решение ^, которому соответствует Р, удовлетворяющее условию (17). Тогда профиль нижнего листа
/гі( х) =
\/б(Р-Р)(х~Х)7 (0 < х < х*),
у[в{Р—Е)(х—х*^о^ (х* < х < /2), ■у/бЕ^/^х^а^ (/2 < х < /1).
Заметим, что \(х*) = 0 ; это означает, что в точке х = х* нижний лист имеет внутренний шарнир. Профиль верхнего листа параболический:
1^2(х) = д/6Р(/2 - х)/ аw .
Для проверки непротиворечивости предположения (а) и требования равнонапряженности УЭ нужно исследовать знак разности у2 (х) - у1 (х). Это
исследование проводилось численно. Установлено, что У2(х) > У1(х) при С<к//1 <^ где С = 2(4 - 42)17, а при /2 //1 < С существуют 0 < х < /2, для которых у2 (х) < у1 (х) . Тем самым равнонапряженный УЭ с найденными выше профилями существует только при /2//1 > С. Вопрос о существовании равнонапряженных двухлистовых УЭ при слабом изгибе для /2//1 < С остается открытым и требует дальнейшего исследования.
О КОЭФФИЦИЕНТЕ ИСПОЛЬЗОВАНИЯ МАТЕРИАЛА УЭ
К вопросу об оптимизации УЭ протеза стопы можно подойти и с другой стороны. В литературе [4], [5] качество УЭ характеризуется коэффициентом
использования материала (КИМ). Пусть УЭ изготовлен из материала с модулем упругости Е и находится в равновесии под некоторой нагрузкой; и - потенциальная энергия упругой деформации УЭ, V - объем УЭ, а тах - максимальный модуль главного напряжения в УЭ. Тогда по определению КИМ - это величина
2Еи
Т =
V о2
^тах
Для вычисления КИМ однолистового УЭ обычно используется несколько более общая, чем рассмотренные выше, модель листа. Эта модель показана на рис. 4. Поперечное сечение - зависящая, вообще говоря, от х область D(х). Главные центральные оси у, z сечения имеют постоянную (не зависящую от х) ориентацию. Изгиб листа происходит в плоскости х z под действием поперечной нагрузки с заданной эпюрой моментов М(х). Для многолистового УЭ при вычислении КИМ обычно используется модель рис. 2 с листами постоянной толщины.
Рис. 4. Модель однолистового упругого элемента, используемая при вычислении
В рамках этих моделей можно найти [4], что для однолистового УЭ с произвольным переменным сечением
max
С (X) M (x)
0 < X < l J (x )
0 J(x)
где S (x) - площадь сечения,
J (x) = if z2 dydz , С (x) = max Izl.
JJ (y, z)e D (x)
D(x)
Для двухлистового УЭ T =--------------------
ц (1 + ц3) [А 2(3 + А) + 4ц3 ]
9(1 + Ц-А) j где А = 1 -l2 /11, ц = hxjh2 .
max
ц2(2 + А), 2ц3 -А ,2А(1 + ц3)'
По отношению к КИМ можно ставить две различные задачи оптимизации. Задача 1. Найти параметры УЭ, для которых T = max.
Задача 2. При заданных E, l, M(х) (для двухлистовго УЭ - E, l1, P ) и при условиях
f = fo, CTmax (18)
где f - прогиб УЭ (f0, ст0 заданы), найти параметры УЭ, для которых V = min
(плотность материала считаем постоянной, поэтому минимизация объема эквивалентна минимизации массы).
Можно показать [4], что для решений этих задач справедливы следующие утверждения.
Утверждение 1. Для однолистового УЭ с произвольным переменным сечением задача 1 не имеет решения: T < 1 и существует последовательность УЭ, для которой T ^ 1. В этом случае и задача 2 не имеет решения: существует число V1, такое, что
V > V1,, и существует последовательность УЭ, для которой (18) выполнены и V ^ V1, .
Утверждение 2. Для однолистового УЭ с произвольным постоянным сечением задача 1 не имеет решения: существует число T, такое, что T < T,, и существует последовательность УЭ, для которой T ^ T,. В этом случае и задача 2 не имеет решения: существует число V2, такое, что V > V2,, и существует последовательность УЭ, для которой (18) выполнены и V ^ V2,.
2
2
I
Утверждение 3. Для однолистового УЭ с сечением в виде прямоугольника постоянной ширины и произвольной переменной высоты задача 1 имеет
неединственное решение: ^х) = const^JMr(x), причем Т = 1/3. В этом случае
существование решения задачи 2 зависит от значения величины
В =
хйх
л/М(Х) ■
Если В < да, то задача 2 имеет единственное решение:
2ст0 В^М (х )
/0 Е
3/№_
2а0 В 2
причем
V =
3/с Е
ст0В о
йх.
Если В = да, то задача 2 не имеет решений: существует последовательность УЭ, для которой (18) выполнены и V ^ 0.
Утверждение 4. Для двухлистового УЭ задача 1 имеет неединственное решение: Х = 9/13, ц = 3/2 (неединственность состоит в том, что определяются только отношения /2/11 и h Х1 h2 ), причем Т = 71/423. В этом случае задача 2 имеет единственное решение:
¡2 = 4/1/13, \ = (142/169)•( /2ст^Е/о), h2 = 2^/3,
щ = (59319/10082 ) • (Р/02 Е2 / ст0 ¡1),
причем
V = (423/71)-(Р/0Е/ст2) .
Аналогично результату утверждения 4, можно найти оптимальное значение КИМ для Глистового УЭ с листами постоянного прямоугольного сечения. Соответствующий результат приведен на рис. 5. При больших N величина Т близка к значению 1/3, соответствующему однолистовому УЭ параболического профиля.
Рис. 5. Зависимость оптимального КИМ от числа листов упругого элемента
Заключение
Ряд новых результатов оптимизации конструкции УЭ протеза стопы позволяет перейти к проектированию новых нестандартных конструкций. При сильном изгибе следует использовать не параболический, а эллиптический профиль листа. Также можно применить конструкцию с внутренним шарниром. Установлена связь задач максимизации КИМ и минимизации массы УЭ. Эта связь не является простой, что говорит о необходимости с некоторой осторожностью использовать данный коэффициент для оценки оптимальности УЭ.
Список литературы
1. Брынских С.И., Осипенко М.А., Няшин Ю.И. Оптимизация двухлистового упругого элемента протеза стопы с использованием линейной и нелинейной теорий изгиба // Российский журнал биомеханики. - 2003. - Т. 7, № 2. - С. 9-19.
2. Няшин Ю.И., Осипенко М.А., Рудаков Р.Н. К теории изгиба листовой рессоры // Известия РАН. Механика твердого тела. - 2002. - № 6. - С. 134-143.
3. Осипенко М.А., Брынских С.И. Отыскание толщин листов равнонапряженных листовых рессор // Вычислительная механика. - Пермь, 2004. - № 2. - С. 51-54.
4. Осипенко М.А., Таланцев Н.Ф. О свойствах коэффициента использования материала листовой рессоры // Изв. вузов. Поволжский регион. Технические науки. - 2009. - № 2(10). - С. 134-144.
5. Пархиловский И.Г. Автомобильные листовые рессоры. - М.: Машиностроение, 1978. - 227 с.
6. Geil M.D., Pamianpour M., Berme N. Significance of nonsagittal power terms in analysis of a dynamic elastic response prosthetic foot // Journal of Biomechanical Engineering. - 1999. - Vol. 121. - P. 521-524.
7. Osipenko M.A., Nyashin Yu.I., Rudakov R.N. A contact problem in the theory of leaf spring bending // International Journal of Solids and Structures. - 2003. - Vol. 40. - P. 3129-3136.
ON OPTIMIZATION OF THE FOOT PROSTHESIS ELASTIC ELEMENT
М^. Osipenko, Y.I. Nyashin (Perm, Russia)
The survey of the new results on the optimization of the foot prosthesis elastic element (EE) is given. The prosthesis is modelled as the cantilever multileaf spring. There is the unbonded contact between the leaves, therefore the bending problem belongs to the beam contact problems. The variable thickness of the uniform-strength of the one-leaf EE is found in the case of the strong bending. It appears that such EE has the elliptic profile. The variable thicknesses of the uniform-strength two-leaf EE leaves are found in the case of the weak bending for some values of the leaves lengths ratio. It appears that in the case of different leaves lengths the longer leaf has the internal joint. It is proved that the uniform-leaf EE may not exist under the additional weak constraints. The material utilization coefficient of the EE is considered. The problems of this coefficient maximization and EE mass minimization are compared. It turns out that the relations of these problems are not simple. The both problems can have no solutions; the first problem can have the solution, while the second problem has no solutions; the first problem can have many solutions, one of them being the unique solution of the second problem.
Key words: foot prosthesis, elastic element, leaf spring, linear and nonlinear bending, contact problem, uniform-strength elastic element, material utilization coefficient, mass minimization.
Получено 11 мая 2011
