CC BY
10
В1СНИК ПРИАЗОВСЬКОГО ДЕРЖАВНОГО ТЕХН1ЧНОГО УН1ВЕРСИТЕТУ 2008 р. Вип. № 18
УДК 539.3+620.178.15
Вовк Л.П.1, Лупаренко Е.В.2
ПРОСТРАНСТВЕННАЯ КОНТАКТНАЯ ЗАДАЧА ДИНАМИЧЕСКОГО ДЕФОРМИРОВАНИЯ НЕОДНОРОДНЫХ ЦИЛИНДРИЧЕСКИХ ДЕТАЛЕЙ МАШИН
Рассматриваются математические аспекты решения пространственной динамической задачи контактирования цилиндрических деталей машин. Приведена модификация метода суперпозиции, позволяющая провести аналитическое исследование напряженно-деформированного состояния и дать практические рекомендации по оценке интенсивности локальной концентрации напряжений в нерегулярных зонах контакта, что практически предопределяет прочностные параметры контактной пары.
Детали, воспринимающие контактные нагрузки, как правило, являются одними из наиболее ответственных, и от их прочности зависит работоспособность и надежность всего механизма в целом. Однако при рассмотрении систем контактирующих тел механическими явлениями в зоне контакта зачастую пренебрегают, проводя упрощение и схематизацию усилий, действующих на детали механизма. Такие упрощения могут приводить к неверным результатам, так как в большинстве реальных механизмов закон распределения контактных давлений оказывает существенное влияние на напряженно-деформированное состояние (НДС) механизма или конструкции в целом. В связи с этим возникает необходимость исследования НДС деталей механизма в местах их контактного взаимодействия при произвольной площадке контакта, размеры которой зависят от приложенных нагрузок.
Теоретические исследования в этой области сопряжены с решением краевых задач нестационарной динамики со смешанными граничными условиями. Одним из наиболее эффективных является аналитический подход, предложенный в [1] и развитый в работе [2]. Разработка аналитических подходов возможна при введении ряда гипотез [1, 2] или эмпирических коэффициентов [3]. В последние годы широкое распространение при решении упругопластических задач получили численные методы, среди которых наиболее универсальным является метод конечного элемента.
Цель работы - построение алгоритма метода суперпозиции, свободного от предположений раздельного рассмотрения волновых движений с различными типами симметрии относительно срединной плоскости цилиндра, распространения этого метода на трехмерные области с криволинейной границей и применения разработанной схемы для численно-аналитического расчета полей контактных напряжений в поршневом пальце (ПП) двигателя внутреннего сгорания (ДВС).
Постановка задачи. Пусть ПП, моделируемый конечным полым цилиндром, занимает в цилиндрической системе координат область D : — Н < z < 0 и b < R < аi.j 0 <(р < 2л .
Предполагаем, что материал цилиндра характеризуется упругими параметрами р -
плотность, Я, Ц - коэффициенты Ляме. Пусть вдоль границы области заданы следующие граничные условия в напряжениях { СО - круговая частота вибронагружения):
R = a: оя=Х(ф,2)е-ш, тRz = /2+ (фz)e~M, ткф = /3+ (фz)e**; R = b: <тЛ=/;(ф,г)е-ш, тRz = /" (ф, тщ = /; (ф ;
z = -h-gz=(p х{Я,ф)е-ш, tRz = у 2(Я,ф)е-ш, тфг = у 3(Д,0>Г,<В'; (1)
z = 0 :а=ср4(кф)е-ш, tRz = у 5(Я,ф)е-\ r^=<p 6Ш)е~ш.
1 ДонНТУ. д-р техн. наук, проф.
2ПГТУ, канд. техн. наук, доц.
Движение среды описывается уравнениями линейной теории упругости в смещениях, связь компонент векторов напряжений и деформаций определяется соотношениями закона Гука в цилиндрической системе координат [5].
Решение сформулированной краевой задачи строим с использованием принципа суперпозиции, используя тот факт, что данная область является пересечением подобластей
Вх\а<Я<Ь и В2\-И <2 < 0. Вектор смещения 0(11,(р,г) в области Д и В2 ищем в виде суммы векторов II (II,(р, г) = Ы1(К,(р,г) + и2{К,(р,г), определяющих решения вспомогательных задач для бесконечного цилиндрического слоя (- йх) и слоя с плоскопараллельными границами (/А —й2).
Решение вспомогательной задачи для бесконечного цилиндрического слоя. С этой целью рассмотрим установившиеся гармонические колебания области : а < Я< Ъ при следующих граничных условиях:
Я = а:<у я= Х;((р,г)е-ш,тЯ2 = Х+О^)^, тК(р = ХЦср^е^-
(2)
Функции напряжений
Хт, вообще говоря, неизвестны, но в рамках постановки
вспомогательной задачи предполагаем их заданными и принадлежащими к классу суммируемых вдоль всей границы рассматриваемой подобласти. На бесконечности задаются условия излучения энергии колебаний. Эти предположения являются естественными и необходимы для корректного использования при решении задачи методики интегральных преобразований Фурье [4]. В результате получено следующее представление для компонент вектора перемещений: щ{К,(р,г) = |мп = ит,и21 = и х,иъх = иЛ|
ик1 (3)
2 л /и р
+00 3
йк1(ЯЖа)= ^ В^ПЛ»«.® ЖЧ«»® )®фИ)
И=-СО / 1
Вид контура интегрирования Г в комплексной плоскости ос определяется использованием принципа предельного поглощения при удовлетворении условиям излучения на бесконечности [6] и имеет следующий вид: обходит положительные особенности подынтегральной функции снизу, отрицательные - сверху, а на остальной части совпадает с вещественной осью. Функции А-п(а,со), В-п(а,со) определяются при удовлетворении выражениями (3) граничным
условиям (2), ), (л* 2)(](,а ,0) ). Ан (а, 0) ) имеют довольно громоздкую
структуру, например,
\в)-а
, 7 = 1,2;
/ \ (4)
в1=Р^2а2 итп ¡л Л + 2/л
При подстановке представлений (3) в граничные условия (2) с использованием соотношений закона Гука в цилиндрической системе координат, в преобразованиях Фурье по координате 2 с разложением всех заданных и неизвестных функций в ряд Фурье по координате (р получаем для каждого значения параметра /7 систему шести функциональных
уравнений для определения неизвестных функций А ,п(а,со ), В,п(а,со ). Численные значения этих коэффициентов при заданных значениях параметров проще получить, реализуя численный
алгоритм их определения непосредственно на ЭВМ. Точность получаемых при этом результатов не ниже, чем при расчете по аналитическим представлениям. Итоговое представление решений этой задачи имеет вид (3), где
+СО 3
ик;(К,<р,а)=^ X[Ф'- '(Ка,ю Щп>(а) + Ф<>У (Л,а,ю Щ„>(а )]ехр(тср) (5)
и=-со у=1
71 +СО
^ /I тш
Здесь Х{г>{а) = — Г Í г/)) ехр [-/(//<р +az)\lz d<p Oír J J J
2л-
-71 -со
Решение вспомогательной задачи для бесконечного слоя с плоскопараллелъными границами. Рассмотрим задачу об установившихся гармонических колебаниях области /А : х, у се (— оо,+оо); — Н < г < 0 при следующих граничных условиях
г = -Я: а=¥1(х,у)е-ш, гХ2 = Г2 (х,у)е-ш, ту2=¥ъ{х,у)е^-
(6)
¿ = 0: о = У4 (х,у)е~ш, т„ =Г5 (х,уУш, туг=¥6(х,у)е^.
На бесконечности, как и ранее, задаем условия излучения энергии колебаний. Решение задачи строим методом интегральных преобразований с использованием принципа предельного поглощения [6]. В результате получаем следующее представление для амплитудной функции вектора смещения й2(х,у,г) = (м12,и22,иЪ2\.
и 2 (х, у, 2) = —^— [\ип (у,Р, г)} ехр[/(> X + ¡3 (7)
4 711 Ц р
Здесь Я]2(у,р,2)=^1]к(у,р,2)гк(у,р\ 7 = 1,2,3;
к=1
г) =—[(и2 +о22)22(ско ¿-ско }[Н-г)сИа 2Я)+%и2о хо ^Иа 2Н$коН-г) + ¡ли
п п п ( , ч \ (и2 +<т,)3л/гст,Я5/гсг 1Н-г)
+4и2(и2 +о22)(ско 2г-сИо 2(Н-г)ско }Н)+----!-^-¿1 ^
<У 1<У 2
2 2 п 2 2 2/л2 /л ®ат г2 Л+2 /Л 2 /Л
и2=у2+р2; а2=и2-в2; в }=—; V2 =-V2
р р
А1(м) = 8(м2 + о 22У(сИо хНска 2Н-\)-[{и2 +о 2)4м~2 + 16и2{а \а 2)](сг "'ст 2)$Ь<у ^ка 2Я; 4(;+з)(/ ./ = , получаем из г) меняя во всех формулах г и Я - г
местами,
%{/,Р)=\ JX (х, У )ехр(- ir х + Р y\fady-
Полученные соотношения позволяют проводить расчет амплитудных характеристик вектора смещения в слое при заданных на его гранях напряжениях Yk(x,y )ехр(— ico /).
Система интегральных уравнений пространственной динамической задачи для конечного цилиндра при однородных граничных условиях. Рассмотрим исходную краевую задачу о возбуждении колебаний в конечном цилиндрическом слое в пространственной постановке (граничные условия (1)). В соответствии с принципом суперпозиции ищем в виде суммы решений вспомогательных краевых задач для бесконечного цилиндрического слоя с граничными условиями (2) и плоскопараллельного слоя (граничные условия (6)):
й(х,У,2) = Щ(*,У,z) + и2(R, q>,z). (8)
Представления для компонент векторов й ■ определены соотношениями (3), (4), (5) - для
йх и (7) для й2. Функции напряжений, определенные граничными условиями вспомогательных краевых задач, вообще говоря, неизвестны. Для их определения используем граничные условия исходной краевой задачи (1) и соотношения закона Гука в цилиндрических и декартовых
координатах [5]. В результате получаем систему двенадцати интегральных уравнений относительно двенадцати неизвестных функций напряжений Х^р (ф, z), Yn(x,y), (j = 1,2,3; п = 1,2,..., 6), определяющих граничные условия вспомогательных задач.
Система интегральных уравнений пространственной динамической контактной задачи для поршневого пальца. При постановке контактной задачи предполагается, что на части поверхности цилиндра заданы компоненты вектора смещения, а на остальной части -компоненты вектора напряжения. Рассмотрим случай, когда на внешней поверхности цилиндра R = a в области Q: <р а [— <р р <р п ] (и z а [/?,, h2 ] заданы компоненты вектора смещения,
например uR (ср, z) — U(z, (р ) . При контакте без трения касательные напряжения в области Q полагаем равными нулю (/^ = ' = 0). Нормальные усилия в области контакта ( (R, (р )) неизвестны и подлежат определению. Получим систему 12 уравнений с 13 неизвестными. Для определения дополнительного неизвестного (нормального контактного давления используем условие равенства компоненты смещения поверхности цилиндра в области Q -uR (ср, z) смещению подошвы штампа (/(г, <р). Это дополнительное уравнение имеет следующую структуру
uR = ии(а,<р ) + yjuf2(acos(p,asm ср, z) + и 22(а eos ср, a sin cp,z) = U(z,cp) (9)
где ип определяется соотношениями (3) - (5); м12,м22 - соотношениями (7). В результате правая часть (9) является интегральным оператором от неизвестных функций напряжений X-,Y-, Получаемое интегральное уравнение замыкает СИУ. Алгоритм получения
представлений интегральных операторов достаточно прост при матричной записи. Однако при его реализации аналитические представления для подынтегральных функций операторов системы чрезвычайно громоздки и труднообозримы. Для исследования и построения решения систем предложено пользоваться численными методами, позволяющими производить расчет значений подынтегральных функций при фиксированных значениях параметров. Для заданных функций напряжений разработаны алгоритмы и программы расчета интегральных представлений, определяющих операторы системы.
Численные результаты. Рассмотрим приложение разработанной выше численно-аналитической методики к уточненному расчету полей контактных напряжений в деталях цилиндро-поршневой группы ДВС. Анализ проведем двумя способами: 1) при помощи метода конечных элементов; 2) решением системы интегральных уравнений с выделением математической особенности механических характеристик в нерегулярных точках области контактирования. В качестве основной расчетной модели был выбран 1111 406.1004020-01 для двигателей марок ЗМЗ-4061.10, ЗМЗ-4062.10 и 3M3-4063.10, шатун в сборе 24-1004045-02 с втулкой 406-1004052-10, поршень ЗМЗ 523-1004014.
При реализации первого способа воспользуемся методикой работы [8]. В связи с тем, что точное распределение контактных напряжений между 1111 и бобышкой и между головкой шатуна и 1111 неизвестны и подлежат определению, в зоне контакта моделировали сетку контактных элементов (общим количеством до 36000). Взаимодействие деталей должно удовлетворять условиям контактирования, в соответствии с которыми точки одного тела не должны проникать внутрь другого.
Для обнаружения условий проникновения задаются контактные узлы, а на поверхности поршневого пальца - контактные элементы. Вместе они образуют контактную пару. Контактный узел, попадая на поверхность контактного элемента, скользит вдоль нее с коэффициентом трения ц. При вычислении контактных сил использовался известный метод штрафных функций.
Результаты решения динамической задачи представлены в виде полей напряжений, деформаций и перемещений деталей поршневой группы. На рис. 1, 2 показаны нормальные напряжения оу,ох, соответствующие пятой моде колебаний и действующие в глобальных направлениях Y, X в принятой системе координат. Видно, что в зоне контакта бобышки
109Е-03 .450Е+09
450Е+09 3 8 0Е+09 310Е+09 2 40Е+09 170Е+09
.402Е+08 .110Е+09
поршня с поршневым пальцем возникает концентрация напряжений, обусловленная влиянием граничных эффектов. Поля напряжений имеют характерный пик в окрестности особой точки и
ниспадающий характер на удалении. Похожая картина возникает в зоне контакта поршневого пальца с головкой шатуна. На рис. 1, 2 видно, что кроме сжимающих напряжений, обусловленных контактным взаимодействием, возникают значительные растягивающие напряжения, обусловленные изгибом, и приводящие к овализации пальца. Попеременное воздействие этих напряжений приводит к появлению усталостных трещин и разрушению. В настоящее время в инженерной практике принято использовать приближенные зависимости, Рис. 1 - Поле нормальных напряжений . предложенные в работе [7].
Предполагается, что нагрузка по длине поршневого пальца распределяется по кусочно-линейному закону. Ясно, что такая схема нагружения будет давать весьма приближенные результаты расчета НДС пальца, так как не учитывает концентрации напряжений в особых точках, характерные черты которой присутствуют на рис. 1,2.
При расчете по предлагаемой в данной работе численно-асимптотической методике предполагаем, что поле напряжений имеет корневую особенность вдоль границы зоны контакта. Это позволяет выделить сингулярные члены в разложениях механических характеристик и улучшить сходимость функциональных рядов для контактных напряжений.
Для примера детально рассмотрим интенсивность локальной концентрации напряжений,
возникающей в зоне контакта поршневого пальца с головкой шатуна (см. рис. 1, 2). Подсчитаем относительное значение внутренней энергии Е [5], накопленной в нерегулярной зоне контакта Г2 возле кромки головки шатуна. Значение Е представляет собой процентное отношение энергии, накопленной в области Г2 к энергии, накопленной во всей области контакта головки шатуна с поршневым пальцем. Если взять начало координат в плоскости кромки головки шатуна, то область Г2 задается неравенствами (координата ^ измеряется в миллиметрах): П: (ра [— 7т,7г]^ г с [—1Д].
ВМХ =.109Е-03 ЭМН =-.450Е+09 ЗМХ =.3 61Е+08 ^^ -. 450Е+09 СИ - . 39 6Е+09 ИН -. 342Е+09
|-1 -.28 8Е+09
|-1 -.23 4Е+09
-. 180Е+09
|-1 -.12 6Е+09
|-1 -.719Е+08
^^ -. 179Е+08 . 3 61Е+08
Рис. 2 - Поле нормальных напряжений а
В таблице 1 подсчитаны значения внутренней энергии на первых шести резонансных частотах, полученные при помощи метода конечных элементов (Ех) и при помощи решения
определяющей системы интегральных уравнений
Таблица 1 - Относительная внутренняя энергия, накопленная в области концентрации
№ моды Значение собственной частоты, Гц Ег е2
1 30.7 51.3 56.4
2 113.6 46.7 50.5
3 326.2 46.1 47.7
4 648.7 40.8 41.8
5 864.5 38.8 38.5
6 1023.6 38.0 34.9
можно принять за концентрации
(Е2). Величины Е1 и Е показатели интенсивности напряжений в области П.
Из данных таблицы 1 следует: 1. Для любой рассмотренной моды колебаний наблюдается существенная концентрация напряжений в нерегулярной области контакта Г2, которая составляет лишь примерно 4 % от общей площади контакта поршневого пальца с головкой шатуна. Это еще раз свидетельствует о невозможности применять в уточненных прочностных расчетах инженерные гипотезы [7].
2. С увеличением порядкового номера моды колебаний интенсивность концентрации напряжений несколько уменьшается, приближаясь к постоянной величине.
3. Результаты расчетов, проведенных по обеим методикам, практически совпадают, что свидетельствует о достоверности полученных данных.
4. Для достаточно больших значений частоты нагружения со целесообразно применять конечно-элементные методики расчета, поскольку весь алгоритм решения определяющей системы интегральных уравнений контактной задачи базируется на том, что длина волны много больше размеров поршневого пальца.
Большой интерес представляет оценка роли граничных условий в формировании спектра резонансных частот и форм колебаний полых цилиндрических деталей. С этой точки зрения в работе проведено рассмотрение динамической смешанной контактной задачи, которая имеет выраженное прикладное значение. Получено интегральное уравнение, определяющее неизвестные напряжения в области контакта, которое замыкает определяющую систему интегральных уравнений относительно неизвестных функций вспомогательной задачи. Таким образом, получены возможности численного исследования всех особенностей контактного деформирования полых цилиндрических деталей при общем случае нагружения их внешней поверхности.
В качестве перспектив дальнейших исследований следует рекомендовать распространение предложенной численно-аналитической методики для анализа полей контактных напряжений для существенно анизотропных тел и тел с трещинами.
Выводы
Предложенная в работе методика обобщения метода суперпозиции для построения решения пространственной динамической задачи вибрационного деформирования полых цилиндрических тел обладает достаточным элементом научной новизны. Суть ее состоит в возможности представления общего решения в виде суммы решений вспомогательных задач для бесконечного цилиндрического слоя и бесконечного слоя с плоскопараллельными границами. Применение интегральных преобразований для решения вспомогательных задач позволяет получить их решение в замкнутой форме.
Перечень ссылок
1. Александров В.М. Контактные задачи в машиностроении / В.М. Александров, Б. Л. Ромалис. -М.: Машиностроение, 1986. - 176 с.
2. Кадомцев И. Г. Осесимметричное упругопластическое соударение двух тел, одно из которых коническое / И.Г. Кадомцев II Изв. СКНЦ ВШ. - 1990. - № 4. - С. 50 - 54.
3. Матлин М.М. Определение параметров первоначально точечного упругопластического контакта по физико-механическим свойствам контактирующих тел / М.М. Матлин II Проблемы машиностроения и автоматизации. - 1993. - № 5. - С. 11 - 20.
4. Уфлянд Я. С. Интегральные преобразования в задачах теории упругости I Я. С. Уфлянд. - М.: Наука, 1968. - 402 с.
5. Гринченко В.Т. Гармонические колебания и волны в упругих телах / В.Т. Гринченко, В.В. Мелешко. - Киев: Наук, думка. 1981. - 284 с.
6. Ворович И.И. Динамические смешанные задачи теории упругости для неклассических областей IИ.И. Ворович, В.А. Бабешко. - М.: Наука, 1979. - 319 с.
7. Кинасошвили P.C. Расчет поршневого пальца авиационного двигателя / P.C. Кинасошвили. -М.: Оборонгаз, 1977. - 148 с.
8. Вовк Л. П. Применение акустического метода свободных колебаний к диагностике технологически неоднородных деталей / Л.П. Вовк II Известия вузов. Северо-Кавказский регион. Технические науки. - 2004. - № 3. - С. 55 - 58.
Рецензент: Ю.Е. Коляда
д-р физ.-мат. наук, проф., ПГТУ
Статья поступила 27.03.2008
