Стеснение деформаций в корневом сечении модели кессона крыла при кручении

УДК 539.319

СТЕСНЕНИЕ ДЕФОРМАЦИЙ В КОРНЕВОМ СЕЧЕНИИ МОДЕЛИ КЕССОНА

КРЫЛА ПРИ КРУЧЕНИИ

© 2012 В. И. Корольков

Воронежский государственный технический университет

В трёхмерной постановке представлено решение упругой задачи об определении напряжённо-деформированного состояния сплошного защемлённого стержня произвольного поперечного сечения, нагруженного сосредоточенным крутящим моментом на свободном конце.

Перемещения, деформации, напряжения, стеснённое сечение.

Ограничение перемещений в одном из поперечных сечений0 подверженного крутящему моменту стержня вызывает местный всплеск напряжений и носит название эффекта Сен-Венана. Так как эта зона с точки зрения прочности конструкций является наиболее опасной, то определение уровня напряжений необходимо для повышения несущей способности подобных деталей. Подобных зон в конструкции летательного аппарата значительное число, но наиболее опасной и высоконагруженной является зона крепления кессона крыла к центроплану. В связи с тем, что в этой зоне наблюдается значительное изменение жёсткости, то происходит сингулярный всплеск напряжений. Точное знание напряжений в этой зоне необходимо для правильного выбора конструктивного решения. Существующие методы расчёта дают лишь приближённую оценку. Численные методы дают в этой зоне неустойчивое решение. Для выбора оптимального решения необходимо обладать достаточно простым, достоверным решением, позволяющим также оценить и качественную картину напряжений и деформаций.

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

1. Постановка задачи

Пусть стержень имеет произвольное поперечное сечение. Боковая поверхность свободна от внешних усилий. Один торец стержня полностью неизменяем. К другому торцу приложена произвольная система сил, статически эквивалентная действию крутящего момента Мкр. Начало декартовой системы координат х0у расположено в центре тяжести неизменяемого сечения, а ось г прямолинейна и проходит через центры тяжести поперечных сечений.

При отсутствии массовых сил перемещения произвольной точки данного тела должны удовлетворять уравнениям Ламе:

л + (Л + л = 0, г, ] = 1,2,3, (1.1) где иг - перемещения произвольной точки стержня, вызванные его деформацией;

1 = Еу/(1 + у)(1 - 2У) , ц = Е/2(1 + У) -коэффициенты Ламе;

Е,  - модуль упругости и коэффициент Пуассона.

Также должны удовлетворяться граничные условия на боковой поверхности: о1]п] = 0, г, ] = 1,2,3

и на торцах стержня:

г = 0; иг = 0 , г, ] = 1,2,3

г = I; Цха3^П= Му , Ц

Л= Цст^О = = JJстззdQ = 0 ,

где ау - составляющие тензора напряжений,

П - направляющие косинусы нормали к поверхности;

- площадь поперечного сечения, I - длина стержня.

Требуется определить напряжённое и деформированное состояния стержня.

2. Метод решения

Известно, что вблизи стеснённого сечения возникает всплеск напряжений, имеющий затухающий экспоненциальный характер вдоль продольной координаты г. На удалении от стеснённого сечения существует однородное напряжённое состояние. Следовательно, используя концепцию пограничного слоя, можно получить асимптотическое решение, называемое внутренним, вблизи стеснённого сечения и срастить его по методу Ван-Дайка [1] с внешним решением для однородного напряжённого состояния. В качестве внешнего решения можно любое известное решение заданной геометрии.

использовать для стержня

внутреннего внутреннюю

характеризующий стрежня; п -

разложение обозначив

Для построения разложения введём

переменную

*=т (' -е -) •

где к - малый параметр, относительное сужение некоторая константа.

Запишем внутреннее перемещений в виде ряда, и = Щ, V = и2, = и3 :

и' = ки (X, у, *) + к 2и2 (х, у, *),

V' = ку;

w' = Щ (х, у, *) + к (х, у, *).

Затем, произведя замену переменных в уравнениях (1.1) и подставив в них разложения (2.2), получим, приравняв члены при одинаковых степенях к, систему шести уравнений.

Решение первых трёх уравнений, удовлетворяющее граничным условиям

Подставим (2.3) в оставшиеся уравнения системы и, удовлетворяя граничным условиям (1.3) на торце г=0, получим выражения для перемещений.

Неизвестные функции поперечных

можно

^ Ь2

И=иг,И-^И= И°, (2.4) где индекс "'" означает разложение по внутренней координате (2.1), а "0" - по внешней координате 2.

Записав условие сращивания (2.4) и приравняв члены одного порядка к, получим шесть уравнений относительно шести неизвестных функций, решением которых получим внутреннее разложение.

Далее, используя процедуру сращивания [1], можно получить составное равномерно пригодное для всего стержня решение. Однако составное решение может быть не единственным. Отличаться решения могут в зоне сращивания. Выбор определённым образом константы п, введённой во внутреннюю переменную (2.1), даёт возможность определить составное разложение единственным образом.

Логично использовать для этого принцип минимума потенциальной энергии деформации стержня [2]: дЭ/ дп = 0,

I I 3

где Э = ¿В' - потенциальная

0  2 ', 1=1

с,, с

12

координат а1, а2 определить из условий сращивания внешнего и внутреннего (2.2) разложений:

энергия деформации стержня; еь- -

составляющие тензора деформаций.

Уравнение (2.5) однозначно определяет константу п, которая является характеристикой эффекта стеснения.

3. Решение

Пространственная задача свободного кручения сплошных слабоконических стержней произвольного поперечного сечения решена Д.Ю.Пановым [3] методом малого параметра. Поэтому используем это решение в качестве внешнего разложения.

Так же, как и в работе [3], будем рассматривать стержень, боковая поверхность которого определяется уравнением:

/

/(£,ц) = 0 , если ввести вспомогательные переменные £ = х(1 - кг) и ц = у(1 - кг). Третья вспомогательная переменная С совпадает с осью г.

В качестве малого параметра используем величину относительного изменения поперечных размеров на единицу

длины стержня или относительное сужение стержня k. Будем считать в дальнейшем величину k настолько малой, что величиной k2 можно пренебрегать по сравнению с k, если это специально не оговорено.

Во внутреннее разложение перемещений вида (2.2), учитывая в нём k второго порядка малости, подставим компоненты перемещений и получим внутреннее разложение с шестью неизвестными функциями. Для определения этих неизвестных функций необходимо использовать условия сращивания (2.4). Для этого приравняем внутреннее разложение, выраженное через внешнюю переменную z, и внешнее разложение [3] по внутренней переменной % . Приравняв члены с одинаковыми степенями k, найдём неизвестные функции.

Составное решение,

удовлетворяющее поставленной задаче на всем отрезке [0, Г], будет определяться следующим образом [4]:

0 . г

и = и + и -

V = V + V -

и

V )0

0 , г (г V

ш = ш + ш - IШ I .

Подставив разложения перемещений, получим равномерно пригодное решение

задачи с точностью до к2:

и = -у +

т Г 2 1 +идФЛ

— у +--——

3п ,

3 и дх

(1 - е~т )е -

+ 2ку 2 + к тР

V = -тxz + -

3п

х -

- 2kтxz + ЫР2 (1 - е

2 1 + и дФ

3 и ду

(1 - е )е -

ш = 2 тф(1 - е ^)+ !тф(1 - е ^ )2 + ТР3 (х, у),

где т - относительное закручивание стержня на единицу длины;

1 / \ Г х ^

Р1

Р2

V

дФ дФ Р3 (х, у) = -2Ф(х, у) - х---у —;

дх ду

Ф

Условие (2.5) даёт возможность однозначно определить константу п.

Таким образом, выражениями перемещений (3.2) и напряжений полностью определяют решение задачи стеснённого кручения стержня заданной геометрии.

4. Призматический стержень Полученное решение показывают, что при стеснённом кручении стержней малой конусности все напряжения и перемещения отличны от нуля. Рассмотрим теперь предельный случай перехода

слабоконического стержня в

призматический, устремив в решениях параметр k к нулю.

Тогда получим для призматических стержней перемещения в виде:

1

и = -у +--т

3п ,

у

2 1 + и дФ

3 и дх

-т )е

1 Г 2 1 + и дФ

V = -х--т| х----

3п , 3 ¡л ду

ш = 2тФ(1 3

- е

')+^тф(1 - е 3

>

-е"")е,

)2.

Получаем, что и призматический стержень в зоне стеснения имеет сложное напряжённое состояние. На удалении от этой зоны е- ж стремится к нулю и выражения (3.3) дают решение Сен-Венана.

Если рассмотреть закреплённое сечение ^=0), то получим, что здесь нормальные напряжения могут даже превышать касательные, а касательное напряжение ст12, как и следовало ожидать, равно нулю. Наибольшим из нормальных напряжений будет продольное напряжение

Например, при коэффициенте Пуассона П=1/3 продольное нормальное напряжение П33 будет в два раза больше поперечных.

5. Стержни частных видов поперечных сечений

+

т

Рассмотрим подробнее решение задачи для некоторых конкретных областей поперечного сечения .

Пусть область  представляет собой эллипс с полуосями а и в в сечении 2=0 и определяется уравнением:

п2

IМ_^+п-1-

а Ь

где | = X (1 - К2), п _ Y (1 - К2).

Функция кручения для данной области, как известно, имеет вид ф(|,п)_ Сц$,

где С = -

Относительный угол закручивания стержня на единицу длины равен

_ Мг (а 2 + Ь 2 ) ж/а 3 Ь 3

Используя формулы (3.2), запишем выражения для перемещений:

_ -э2

+ КСу

3"

х1 -12 4 х2 - 3 у2

(1 - е-"2)-

-1 Кту

. _ «2(1 - 2К2)-тх[\ - 2^с\(1 - е-"2У"2 +

+ КСх

—--

1 V 2 12

у\\ у — х

2 А 3

(1 - е-"2 )+

+1 Ктх

со _ 2тСху(1 - е-"2)+1гСху{1 - е~"2) -

- 4KlCxyz. Для

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

Параметр " является важной характеристикой эффекта стеснения, обуславливающей не только величину зоны распространения краевого эффекта от стеснения, но и интенсивность и характер изменения напряжений в этой зоне. Кроме того, параметр " даёт возможность проводить качественный анализ

напряжённо-деформированного состояния материала в области влияния эффекта стеснения от геометрических размеров стержня, и характеристик материала. Параметр ", а следовательно и эффект стеснения, в случае призматического стержня не зависит от его длины, что подтверждает тем самым гипотезу Сен-Венана. Конусность же имеет некоторое, но не сильное, влияние на параметр ". Уменьшение же нормального напряжения а33 в стеснённом сечении 2=0 с увеличением конусности подтверждает вывод

Е.П.Гроссмана [5] о том, что эффект стеснения и эффект конусности создают при кручении нормальные напряжения разных знаков.

Зависимости коэффициента " от отношения полуосей эллипса при кручении призматического стержня дают возможность сделать некоторые качественные выводы: во-первых, коэффициент " имеет минимальное значение при Ь/а=0.85-0.9, т.е. стержни с таким соотношением полуосей имеют наибольшую зону распространения краевого эффекта от стеснения; во-вторых, с увеличением коэффициента Пуассона зона краевого эффекта, в связи с возрастанием ", уменьшается, а нормальные напряжения несколько возрастают; в-третьих, у круглого стержня (Ь/а^-1) и очень тонкой пластинки (Ь/а^-0), как и следовало ожидать, эффект стеснения исчезает, т.к. коэффициент "^<х>.

Сравнение графиков распределения нормального безразмерного напряжения

* 3

о33 = о33 а /М2 вдоль призматического стержня эллиптического сечения при следующих расчётных данных: у=0.32, Ь/а=0.5, а/Ь=0.25, х/а=0.75, у/а=-0.33 (рис. 1) показывает, что а33 имеет одинаковый характер поведения во всех трёх сравниваемых решениях. Расхождение наблюдается в численных значениях и достигает между представленным в работе решением и решением Н.В.Зволинского [6] в стесненном сечении 25%. С результатами А.Феппля [2] различие значительно меньше.

3 2,5 2

со

Г 1,5

Ь

1 0,5 0

данное решение

решение А.Фёппля решение Н.В.Зволинского

0 0,25 0,5 0,75 1 1,5 2 2,5 3 г/а

Рис. 1. Сравнительные графики распределения безразмерного продольного нормального напряжения а33 при кручении эллиптического стержня

Решения Н.В.Зволинского, А.Феппля и С.П.Тимошенко [7] не удовлетворяют всем уравнениям статической теории упругости. В первом решении не удовлетворяются точно уравнения равновесия, а в других двух решениях не выполняются условия совместности деформаций, т.е. они являются приближёнными решениями. Решение, представленное в данной работе, имеет более общий характер, т.к. оно справедливо для стержней произвольного поперечного

сечения, как призматических, так и имеющих отклонение от призматичности.

На рис.2 даны сравнительные графики распределения безразмерных нормальных и касательных напряжений ст*. = ст.  а3 / м г

вдоль оси г при кручении призматического эллиптического стержня. Расчётные данные прежние. Нормальное поперечное напряжение ст22 будет равно ст** , т.к. стержень имеет постоянное по длине сечение.

ст

2,25

1,5

0,75

\ * у-ст33 / К * ст13

N # \ -------- * СТ 23 ----12.3.__

ч * \ "ст12 .-к г ' — • !

0,75

1,5

2,25

3

Рис. 2. Сравнительные графики распределения нормальных и касательных напряжений вдоль оси г при кручении стержня эллиптического поперечного сечения

3

0

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

С целью проверки полученного аналитического решения было проведено экспериментальное изучение распределения напряжений, порождаемых эффектом стеснения, в сплошных стержнях. Использовался метод фотоупругости [8] в связи с тем, что с его помощью представляется возможность получения полной картины напряжений. Для получения картин полос использовался метод замораживания на моделях прямоугольного поперечного сечения из оптически чувствительного материала ЭД20-МТГФА.

Эксперименты полностью

подтвердили картину распределения и уровень напряжений вблизи стеснённого сечения.

Представленное решение позволяет оценить уровень всех напряжений и деформаций в зависимости от формы сечения и упругих характеристик материала.

Так как получено решение упругой задачи, то записав аналогично решения для других нагрузок и применив принцип суперпозиции, получим решение для реального случая нагружения.

CONSTRAINT STRAIN IN THE ROOT SECTION WING MODEL WITH TORSION

© 2012 V. I. Korolkov Voronezh State Technical University

In the three-dimensional formulation is presented the solution of the problem of determining the elastic stressstrain state of solid clamped rod of arbitrary cross section loaded by concentrated torque at the free end.

Displacement, strain, stress, constrained by section.

Информация об авторе

Корольков Владимир Иванович, доктор технических наук, профессор, заведующий кафедрой «Самолётостроение», Воронежский государственный технический университет. Email: korolkov_vi@bk.ru. Область научных интересов: математическое моделирование процессов локального циклического пластического формообразования тонкостенных оболочечных деталей, разработка методик расчёта на прочность конструкций с учётом зон сингулярного характера напряжений, моделирование браковочных признаков и изменений технологических параметров листовых материалов в процессе деформирования.

Библиографический список

1. Ван-Дайк, М. Методы возмущений в механике жидкости [Текст] / М. Ван-Дайк -М.: Мир, 1967. - 310с.

2. Феппль, А. Сила и деформация [Текст] / А. Феппль, Л. Феппль - М.; Л.:ОНТИ НКТП СССР, 1936. Т.2. - 408с.

3. Панов, Д.Ю. О кручении стержней, близких к призматическим [Текст] / Д.Ю. Панов // Докл. АН СССР.- 1939. Т.20,№4.-С.251-254.

4. Найфе, А.Х. Введение в методы возмущений [Текст] / А.Х. Найфе - М.:Мир, 1984.-536с.

5. Гроссман, Е.П. Расчёт на прочность лонжеронных крыльев [Текст] / Е.П. Гроссман // Труды ЦАГИ.- 1947.- Вып.628.-77с.

6. Зволинский, Н.В. Приближённое решение задачи кручения упругого цилиндрического бруса с одним неизменяемым сечением [Текст] / Н.В. Зволинский // Изв. АН СССР, ОТН. 1939 -№8. - С.91-100.

7. Тимошенко, С.П. Теория упругости [Текст] / С. П. Тимошенко - М.; Л.: Гостехиздат,1934. - 451.

8. Хесин, Г.Л. Метод фотоупругости [Текст] / Под ред. Г.Л. Хесина - М.: Стройиздат, 1975. Т.1. - 460 с.

Korolkov Vladimir Ivanovich, Dr. Sc. Sciences, head of the Department "Aircraft Manufacturing", Voronezh State Technical University, E-mail: korolkov_vi@bk.ru. Area of research: mathematical modeling of local cyclic Plastic Forming thin shell components, the development of methodologies for calculating the strength of structures subject areas of the singular nature of the stress, modeling for rejection and changes of process parameters in sheet metal deformation process.