Аналитическое решение контактной задачи об упругих деформациях тонкой пластины, помещенной в цилиндрическую полость Текст научной статьи по специальности «Физика»

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408

электронный научно-технический журнал

Аналитическое решение контактной задачи об упругих

деформациях тонкой пластины, помещенной в цилиндрическую

полость

# 01, январь 2013

DOI: 10.7463/0113.0517977

Дмитриев С. Н.

УДК 539.3:62-27

Россия, МГТУ им. Н.Э. Баумана dim. sm2@yandex.ru

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

Рассмотрим тонкую прямоугольную пластину (рис. 1) с размерами а х Ь в плане и толщиной И << min(a, Ь).

к-а-А

Рис. 1

Пусть каким-либо образом, например, путем приложения распределенного изгибающего момента постоянной интенсивности вдоль двух противоположных кромок длиной Ь, пластина свернута в цилиндр радиуса Я и помещена в цилиндрическую полость того же радиуса (рис. 2).

А' А

Рис. 2

При таком сворачивании перемещения точек пластины оказываются сопоставимыми с габаритными размерами пластины. В работе [1] В.Д. Черненко получил аналитическое решение данной задачи в предположении, что пластина плотно прилегает к полости по всей поверхности. Техническое приложение решения дано в [2], где определены величины контактных усилий, возникающих при изгибе и сворачивании в цилиндр резинокордной конвейерной ленты. Предположение о плотном прилегании применимо в тех случаях, когда изгибная жесткость пренебрежимо мала по сравнению с жесткостью на растяжение-сжатие, что имеет место для резинокордной ленты. В случае если пластина изготовлена из изотропного материала будет наблюдаться неплотное прилегание к внутренней поверхности цилиндра, и возникнет зона отслоения. Такая задача относится к классу контактных задач с неизвестной заранее областью контакта [3]. В статье [4] В.Л. Бидерман и С.В. Бояршинов получили аналитическое решение контактной задачи о посадке упругого кольца на вал. При этом упругое кольцо контактирует с жесткой цилиндрической поверхностью, зона контакта заранее неизвестна. Ими были определены контактные усилия и найдено значение угла отслоения. Перемещения кольца при посадке на вал малы и поэтому задача, рассмотренная в [4] является геометрически линейной. В докладах [5, 6] был предложен и развит численный

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

Условия равновесия на отслоенном и прилегающем участках

Следуя [2], будем считать, что срединная поверхность изогнутой пластины мало отличается от поверхности кругового цилиндра радиуса R, поэтому целесообразно деформацию пластины рассматривать в системе координат, связанной с поверхностью цилиндра. Считаем, что пластина находится в условиях цилиндрического изгиба, так что можно ограничиться рассмотрением деформаций в одной плоскости. При этом в отличие от [2], предполагается наличие зоны отслоения, размер которой определяется некоторым углом а . Обычно [3, 4], распределение контактных усилий в подобных задачах принимается в соответствии с рис. 3: здесь д - распределенное по площади контактное давление постоянной интенсивности, Q1 и Q2 - усилия, распределенные вдоль границ зон контакта.

Рис. 3

Допустимость принятой силовой схемы должна быть подтверждена в ходе решения. При этом, как отмечалось в [7], принимаемое в подобных задачах распределение сил, http://technomag.bmstu.ru/doc/517977.html 3

противоречит предположению о плотном прилегании пластины к поверхности цилиндра на участке BC. Здесь пластина должна находиться в условиях чистого изгиба и при использовании теории пластин Кирхгофа-Клебша поперечная нагрузка q должна быть равна нулю, а нагрузка должна быть представлена изгибающими моментами, приложенными к наружному контуру пластины. При принятой схеме нагрузок, должны каким-то образом, проявится противоречия при решении уравнений в перемещениях и удовлетворении граничным условиям. Как будет показано ниже, предположение о плотном прилегании на участке BC выполняется с некоторой погрешностью. Здесь будут перемещения, но они оказываются существенно меньше, чем в зоне отслоения AB.

Так как в направлении образующей цилиндра, в силу принятых допущений, пластина не изгибается, внутренние силовые факторы могут быть найдены из условий равновесия методом сечений. Силовые факторы действующие на отслоенном участке АВ показаны на рисунке 4.

А

Рис. 4

Записывая уравнения равновесия, получим:

Nx = Q\sin ф

Qx = Qicos^ (1)

Mx = Q1R sin ф

Здесь Nx - погонная продольная сила, Qx - погонная поперечная сила, Mx - погонный изгибающий момент, ф - некоторый текущий угол. Принятые положительные направления силовых факторов, показаны на рис. 4.

Для прилегающего участка силовые факторы показаны на рис. 5:

А

Рис. 5

Уравнения равновесия для участка ВС:

Nx = -Q1 sin р- Q2 sin(p -а) + qR[cos(p - а) -1] Qx = Q1 cos р + Q2 cos(p - а) + qR sin(p - а) Mx = Q1R sin р + Q1R sin(p - а) + qR2[1 - cos(p - а)] а - угол, характеризующий размер зоны отслоения.

Определение контактных сил и давления

На участке BC пластина в силу плотного прилегания к цилиндру находится в условиях чистого изгиба и момент на этом участке постоянный, чтобы на участке ВС выполнялось условие M x = M = const, из третьего уравнения системы (2) получим

а = -f-

sin а

q2 = -qR ■ ct§a (з)

M = qR2

При чистом изгибе изгибающий момент связан с кривизной соотношением

1 M qR2

R D D

(4)

E ■ h3

Где D = (-2) цилиндрическая жесткость, E - модуль упругости материала,

/и - коэффициент Пуассона. Из (4) следует

Б Я

Таким образом, оказываются определенными значения всех контактных усилий.

q = (5)

Определение перемещений на отслоенном участке и размера зоны отслоения

Для определения неизвестного угла а потребуется рассмотрение деформаций отслоенного участка АВ пластины. Как и при записи уравнений равновесия, будем считать, что срединная поверхность изогнутой пластины мало отличатся от поверхности кругового цилиндра радиуса R и используем геометрические и физические уравнения цилиндрической оболочки [8]. В случае изгиба в одном направлении продольное усилие Nx и изгибающий момент Mx связаны с компонентами перемещения u и w соотношениями:

Eh 1 ,du

Nx =--т~(~г + w)

1 - ^ R dp

„ D„ 1 d2 w 1 du. (6)

Mx = — (1---^ +--)

x R R dp R dp

Первое слагаемое D в формуле для момента Mx учитывает момент, возникающий при R

сворачивании пластины в цилиндр. Положительные направления перемещений u и w показаны на рис. 3.

Используя уравнения равновесия для зоны отслоения (1) и найденные значения силовых факторов (3) перепишем (6) в виде:

Eh 1 ,du . _ . D

--(— + w) = -Q1 sinp = --—-sinp

1 -v R dp R sina

D „ 1 d2 w 1 du ч ^ . D (7)

— (1---- +---) = Q1R ■ sinp =-sinp

R R dp R dp Rsina

После преобразований приходим к системе из двух дифференциальных уравнений для определения неизвестных перемещений u и w в зоне отслоения:

du sin p

--+ w =-a -

7 u

dp sina

,2 , • (8)

dwdu sinp

-2---= R ■ (1--)

dp dp sina

h2

где au =-. Из системы (8) получим дифференциальное уравнение для определения w :

12R

d2w ^ ,, sin y. sin®

— + w = R ■ (1 ^ - (9)

dy sin a sin a

Общее решение уравнения (9):

R + a

w = C1 ■ siny + C2 ■ cosy + R +---cosy (10)

2sina

Неизвестные константы C1 и C2, а так же угол отслоения a определяются из граничных условий. Потребуем выполнения на концах участка АВ равенства нулю поперечных перемещений, а в зоне сопряжения с прилегающим участком BC равенства нулю угла поворота:

w(o) = w{a) = 0

^ (a)=o (11)

dy

Удовлетворяя граничным условиям, получим значения констант:

R + a

C1 = -R sin а ' cos2 а 2sin а

C2 =- R

и трансцендентное уравнение для определения угла а :

1 1 + ««

(12)

R R а —cosа+----= 1

cosa+ —--;-=1 (13)

2 2 sin a

Ввиду малости толщины пластины h, коэффициент au << R, поэтому (13) можно приближенно переписать в виде:

a ,

cosa+ —-= 1 (14)

sin a

Уравнение (14) совпадает с полученным в [4] аналогичным уравнением в контактной задаче о посадке упругого кольца на вал. Значение угла a, удовлетворяющее (14) будет таким же, как и в (4): a = 122035'. При малых значениях au решения (13) и (14) отличаются

незначительно, так при R = 1000 мм и h = 1 мм, значения угла a = 122034'. Значения внутренних силовых факторов могут быть найдены по уравнению (6), а поперечного перемещения w по уравнению (10) с учетом значений констант (12). При вычислениях можно пренебречь au по сравнению с радиусом цилиндрической полости R . Максимальное значение поперечного смещения в относительных единицах составляет приблизительно:

_ w — R

wmax = max « 0,1, максимальное значение изгибающего момента Mxmax = Mxmax— « 1,2.

max x max x max p 1

Хотя значения перемещений велики, они, так же как и размер зоны отслоения получаются близкими к найденным в работе [5], где применялся численный метод для нелинейных задач. Поэтому данное решение может быть использовано в качестве первого приближения.

Перемещения на прилегающем участке

Оценим так же величину продольного смещения u на свободных кромках свернутой в цилиндр пластины. Используя первое уравнение из системы (8) и решение (10) уравнения (9) запишем:

du sinp ■ п v R + au sinp пч\ — = -w - au-= -C1 • sinp- C2 • cosp- R---p- cosp- au--(15)

dp sin a 2sina sin a

Где константы C1 и C2 определяются формулами (12). При интегрировании уравнения (15) появляется еще одна константа C3, которая должна быть определена из граничных условий

на отслоенном участке. Учитывая, что наше решение является приближенным, пренебрежем продольным перемещением на прилегающем участке BC и примем в качестве граничного условия u(a) = 0. После интегрирования и удовлетворения граничному условию получим u(p):

u =

г a \

C1 + . u I • (cos (p - cos a)- C 2 • (sin (p- sin a) - R •(p-a)-sina )

(16)

R + au ( ■ ■ \ R + au ( \

—:-(p • sin pp - a • sin a)--:-(cos p- cos a)

2sina 2sina

Подставим в (16) значения констант C1, C 2 и упростим конечное выражение при au << R, полагая a. « 0:

u = - R •

^sin2 a + 2 1 ^

1 ' (cos (p - cos a) - (sin pp - sin a) + (p-a) + —:-(p • sin p - a • sin a) (17)

-:--(cos p- cos a) - (sin p- sin a) + ( - a) + —:-(p • sin p - a • sin a

2 • sin a 2sina )

На свободной кромке пластины для значения угла отслоения а = 122035' получаем

и (о) = ОД (18)

При таких перемещениях свободные кромки пластины будут наползать друг на друга (что фактически имеет место при сворачивании конвейерной ленты), если рассмотреть не пластину, а тонкую проволоку, сворачиваемую в кольцо, то будет наблюдаться перехлест концов проволоки. Условия, при которых реализуется решение данной задачи, можно выполнить, если взять ширину пластины Ь < 2 -ж - R, учитывая (18), Ь « 1,8 -ж - R.

Оценим теперь величины перемещений на прилегающем участке ВС. Здесь продольная сила и изгибающий момент (2) после подстановки контактных сил (3) будут:

N = -дК

М г = дЯ> <19>

Подставляя в уравнения (6), связывающие усилия и перемещения при изгибе цилиндрической оболочки в одном направлении получим систему из двух дифференциальных уравнений для перемещений в зоне прилегания, аналогичную системе (8) для зоны отслоения:

du

— + w = -аи

1 и

dy

Л (20)

dwdu

--= 0

dy dy

Система содержит очевидное противоречие с исходным предпосылками, принятыми при решении задачи. В силу прилегания пластины к поверхности цилиндра на участке BC должно быть w = 0 и, как отмечалось выше, в рамках классической теории изгиба при чистом изгибе q = 0. Ясно, что для правомерности приведенного выше решения, перемещения, полученные из (20) должны быть малыми. Решая систему уравнений в той же последовательности, что и (8), при граничных условиях

w(n) = 0, — (п) = 0, п{ж) = 0 (21)

dy

получим значения перемещений в зоне прилегания

u = a„ ■ sin y

U (1 ) (22) w = -au ■ (1 + cos y)

Таким образом, в зоне прилегания перемещения оказываются малыми порядка 0,1h, в то время как, на отслоенном участке большими - порядка 0,1R . Видно, что предположение о плотном прилегании на участке BC выполняется приближенно.

Выводы

1.В статье получено решение контактной задачи о пластине, свернутой в цилиндр и помещенной в цилиндрическую полость, причем в отличие от известных решений [1,2] размер зоны прилегания заранее неизвестен и определяется в ходе решения.

2.Известно [7], что обычно принимаемое в подобных задачах распределение сил [3], противоречит (при использовании классической теории изгиба) предположению о чистом изгибе части пластины, плотно прилегающей к поверхности цилиндра. В данной статье

показано, что это противоречие дает перемещения в зоне прилегания порядка 0,1 от толщины пластины.

З.Распространение полученного решения, учитывающего наличие зоны отслоения, на ортотропные пластины [1,2] не представляет сложности и требует только коррекции параметров в уравнениях (6), связывающих перемещения с внутренними силовыми факторами.

Список литературы

1. Черненко В.Д. Изгиб ортотропной цилиндрической пластины в цилиндрическую оболочку // Прикладная механика. 1972. Т. 15, № 4. С. 49-53.

2. Черненко В.Д. Расчет средств непрерывного транспорта. СПб.: Политехника, 2008. 386 с.

3. Моссаковский В.И., Гудрамович В.С., Макеев Е.М. Контактные задачи теории пластин и оболочек. М.: Машиностроение, 1978. 248 с.

4. Бидерман В.Л., Бояршинов С.В. Расчет храпового механизма пружинного типа // Труды кафедры сопротивления материалов МВТУ им. Н.Э. Баумана, 1947. С. 127-142.

5. Дмитриев С.Н., Солодовникова С.Ю. Нелинейная контактная задача об изгибе тонкой пластины, помещенной в цилиндрическую полость // Ракетно-космическая техника: фундаментальные и прикладные проблемы : труды 2-й Междунар. науч. конф. (Москва, 18-21 ноября 2003 г.) : в 4 ч. 2005. Ч. 2. С. 56.

6. Дмитриев С.Н., Солодовникова С.Ю. Контактная задача о равновесии пластины, свернутой в цилиндр // Аэрокосмические технологии, 2004-2007 : сб. тр. М.: Изд. МГТУ им. Н.Э. Баумана, 2008. С. 348-351.

7. Феодосьев В.И. Избранные вопросы и задачи по сопротивлению материалов. М.: Наука, 1967. 376 с.

8. Погорелов В.И. Строительная механика тонкостенных конструкций. СПб.: БХВ-Петербург, 2007. 528 с.

SCIENTIFIC PERIODICAL OF THE RAIJMAN MS TU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Analytical solution of the Hertzian contact problem of elastic

deformations of a thin plate positioned in a cylindrical cavity

# 01, January 2013

DOI: 10.7463/0113.0517977

Dmitriev S.N.

Russia, Bauman Moscow State Technical University

dim. sm2@yandex.ru

The author proposes an analytical solution of the contact problem of elastic deformations of a thin plate rolled up into a cylinder and placed into a cylindrical cage of the same radius. In this case, the plate is only partially adjacent to the surface of the cylinder, on the other part there is delamination. For linearization of basic correlations it is assumed that the plate deviates from the surface of the cylindrical cavity to a very little degree; therefore, it is possible to use the equilibrium equation, as well as physical and geometric equations of cylindrical shell. It is considered that the shell bending is cylindrical. The rolling of the plate is taken into account as a pre-stressed state. It is supposed that it is a state of pure bending in plane. As a result of solution, the angle of delamination, displacement and internal power factors are determined.

Publications with keywords: bending, slice, shell, deformations Publications with words: bending, slice, shell, deformations

References

1. Chernenko V.D. Izgib ortotropnoi tsilindricheskoi plastiny v tsilindricheskuiu obolochku [Bending of orthotropic cylindrical plate into the cylindrical shell]. Prikladnaia mekhanika [Applied mechanics], 1972, vol. 15, no. 4, pp. 49-53.

2. Chernenko V.D. Raschet sredstv nepreryvnogo transporta [Calculation of means of continuous transport]. St. Petersburg, Politekhnika, 2008. 386 p.

3. Mossakovskii V.I., Gudramovich V.S., Makeev E.M. Kontaktnye zadachi teoriiplastin i obolochek [Contact problems of the theory of plates and shells]. Moscow, Mashinostroenie, 1978. 248 p.

4. Biderman V.L., Boiarshinov S.V. Raschet khrapovogo mekhanizma pruzhinnogo tipa [Calculation of the ratchet mechanism of spring type]. Trudy kafedry soprotivleniia materialov MVTU im. N.E. Baumana [Works of the Department of Strength of Materials of the Bauman MHTU], 1947, pp. 127-142.

5. Dmitriev S.N., Solodovnikova S.Iu. Nelineinaia kontaktnaia zadacha ob izgibe tonkoi plastiny, pomeshchennoi v tsilindricheskuiu polost' [Nonlinear contact problem of the bending of

a thin plate, placed in a cylindrical cavity]. Raketno-kosmicheskaia tekhnika: fundamental'nye i prikladnye problemy : trudy 2-i Mezhdunar. nauch. konf. [Rocket-space technology: fundamental and applied problems : proc. of the 2th International scientific conference], Moscow, 18-21 November 2003, 2005, pt. 2, pp. 56.

6. Dmitriev S.N., Solodovnikova S.Iu. Kontaktnaia zadacha o ravnovesii plastiny, svernutoi v tsilindr [Contact problem of the equilibrium of the plate, rolled into a cylinder]. Aerokosmicheskie tekhnologii, 2004-2007 : sb. tr. [Aerospace technologies, 2004-2007 : the collection of works]. Moscow, Bauman MSTU Publ., 2008, pp. 348-351.

7. Feodos'ev V.I. Izbrannye voprosy i zadachi po soprotivleniiu materialov [Selected questions and tasks in Strength of Materials]. Moscow, Nauka, 1967. 376 p.

8. Pogorelov V.I. Stroitel'naia mekhanika tonkostennykh konstruktsii [Structural mechanics of thin-walled structures]. St. Petersburg, BKhV-Peterburg, 2007. 528 p.