68 ДОБЫЧА
УДК 621.77.07
Научные основы моделирования процесса упрочнения впадины резьбы бурильных труб обкатыванием роликом
м.В. Песин
к.т.н., доцент1, первый заместитель директора по нефтепромысловому оборудованию2 M.Pesin@mail.ru
1ПНИПУ, Пермь, Россия
2ЗАО «Торговый дом ПКНМ», Пермь, Россия
Выполнено компьютерное моделирование процесса упрочнения впадины резьбы обкатыванием роликом. Приведена концептуальная и математическая постановка задачи. В результате серии вычислительных экспериментов с использованием пакета ABAQUS получены оценки распределений остаточных напряжений в зависимости от степени вдавливания ролика.
материалы и методы
Использованы методы математического моделирования.
Ключевые слова
упрочнение резьбы обкатыванием, остаточные напряжения, метод конечных элементов
Одним из методов поверхностного упрочнение впадин конической замковой резьбы бурильных труб является обкатывание этих впадины роликом. В этом процессе приповерхностные слои материала впадины резьбы трубы испытывают пластические деформации, а возникающие при этом остаточные напряжения и поверхностное упрочнение материала (наклеп) повышают усталостную прочность труб, увеличивая срок их эксплуатации. Экспериментальное изучение остаточных напряжений вблизи впадины резьбы затруднено в силу сложного профиля поверхности резьбы и сильной неоднородности остаточных напряжений в малом по толщине слое материала резьбы. Поэтому целью данной работы являлась оценка с помощью методов математического моделирования распределений остаточных напряжений, возникающих после обкатывания. На рис. 1 приведена резьбовая часть ниппеля бурильной трубы, а также профиль и параметры резьбы. Геометрия ролика для обкатывания (рис. 2).
Ролик устанавливается во впадине резьбы без наклона к оси обрабатываемой детали и прижимается усилием P. Затем труба начинает вращаться, заставляя прижатый ролик катиться по впадине резьбы.
С точки зрения механики деформируемого твердого тела задача моделирования процесса обкатывания резьбы относится к трехмерным нестационарным контактным задачам упругопластического деформирования. Сложная геометрия моделируемых тел исключает возможность использования аналитических методов для решения подобной задачи. Далее представлена система гипотез для построения математической модели.
В качестве объекта моделирования выбираем область материала ниппеля трубы с резьбой. Введем некоторые упрощения.
Так как нагрузки при обкатывании вызывают только локальные пластические деформации во впадине резьбы под роликом, то в качестве моделируемого объема можно рассмотреть впадину с двумя витками резьбы. При этом будем рассматривать резьбу не на конической, а на цилиндрической поверхности. Учитывая достаточно большой (по сравнению с глубиной впадины резьбы) диаметр трубы, локальность деформирования и оценочный характер расчетов можно моделируемый объем развернуть с цилиндрической поверхности на плоскость. Все эти упрощения позволяют представить моделируемый объем в виде фрагмента одного витка длиной 20*10-3 м, так как показано на рис. 3. В начальном положении ролик располагается над канавкой резьбы в 5*10-3 м от левого торца витка (рис. зб). Ролик считаем абсолютно твердым телом. Деформации материала резьбы считаем малыми, а физико-механические свойства материала — изотропными. В начальном состоянии материал резьбы принимаем находящимся в естественном состоянии (ненапряженным и ненагруженным). Процесс нагружения принимаем изотермическим. Влияние массовых сил пренебрегаем. Принимаем гипотезу об аддитивности скоростей упругих и пластических деформаций. Задача рассматривается как контактная, причем контактные условия постоянно уточняются в процессе решения. Нагружение выполняется в 2 шага. На первом шаге ролик вдавливается во впадину резьбы на величину На втором шаге весь фрагмент резьбы смещается влево (рис. зб) на 10*10-3 м. При этом вдавленный ролик катится по впадине, вращаясь относительно своей неподвижной оси.
Перейдем к математической постановке данной задачи. Согласно сформулиро-
Рис. 1 — Нипель резьбовой части бурильной трубы (а), профиль и размеры резьбы (б)
Рис. 2 — Геометрия ролика
ванным выше гипотезам на всех этапах рассматриваемого технологического процесса с точки зрения математического моделирования задачу можно представить как трехмерную термовязкоупругопла-стическую нестационарную квазистатическую контактную задачу. Её общая постановка приведена в [1, 2].
Пусть в некоторый момент времени X £ (О, ^ деформируемое тело занимает область г £ '0х, '0х £ к3 (трехмерное пространство) с границей З^З^З^О^З^Х), $у(х) — часть границы, подвергающаяся кинематическому воздействию, З (О — свободная граница, Зкон(х) — граница контакта (смешанные граничные условия), '0х= '0х ^З'.
Необходимо определить поля скоростей перемещений V, напряжений о и деформаций е для любого X £ (О, го), удовлетворяющие следующей системе уравнений: уравнениям равновесия в скоростях
, (1)
+ Л7а> = 0 геП'
определяющим
соотношениям (ОС)
(2)
статическим и кинематическими граничным
условиям
VX=0, 0ХУ =Oxz = 0 f е 5"! ; (4а)
V=0, 0ху = Oyz = О Г 6 s;, ; (46)
°x = 0, = 0z = 0 ; (4в)
v =0, о =0 =0 г«А""'1 (4г)
z xz yz ■ vt v
Vz=-0,1 м/сек, 0xz = 0yz = о Г e.V.';1'1; (4Д)
условия для оси ролика v = =0, v =-0,1 м/сек, для шага 1,
х VZ ' у ' ' 7 П 7
vx=vy=vz=o, для шага 2 при начальных условиях о(г) = е(г) = о, v(r) = о, геГ5" ; (5) Для первого этапа техпроцесса соотношения (5.1о)-(5.1з) выполнятся автоматически согласно принятым гипотезам, для последующих — исходя из непрерывности процесса (начальный момент нового этапа соответствует концу предыдущего).
Перемещения u и напряжения о определяются интегрированием:
(6)
(7)
в лрочивноч с
кинематическим соотношениям
I) ^(Cv.vC), reti
(з)
я|гу)- <г„(г} ■ |о-[Г,Г к/Г
По полю перемещений (6) можно построить произвольную меру деформации.
В уравнениях (1)-(7): С — набла-опе-ратор (векторный дифференциальный оператор Гамильтона в текущем лагран-жевом базисе), верхний индекс к означает
Рис. з — Модель обкатки: (а — часть витка; б — вид с боку; в — вид с вершины оси г)
коротационную производную, й, — тензоры полной и пластической деформации скорости, п — вектор единичной нормали к границе, оЗ — упрочнение материала, ои — интенсивность напряжений.
Соотношениям (1-7) должны удовлетворять поля (V, о, е) при моделировании на каждом этапе рассматриваемого технологического процесса. Отличия касаются начального состояния материала на каждом этапе, начальной для этапа конфигурации области и задаваемых граничных условий. Отметим, что определяющие соотношения (2) задают поведение материала резьбы согласно ма-крофеноменологической теории пластического течения.
Данная задача является существенно нелинейной, в частности — вследствие наличия контактных граничных условий. Рассмотрим более подробно контактные граничные условия.
Следует отметить, что контактная граница Зкон(') априори неизвестна. Для её определения используется следующий подход.
Обозначая через З ('1) и З
г возм.кон.рез4 у возм.
('1) поверхность исследуемой области
кон.рол4 ' Г
резьбы и ролика, где контакт между этими телами возможен, а через Зс('1) — контактную поверхность резьбы для момента времени Область контакта можно определить как область пересечения поверхностей Зв ('1) следующим образом:
(t) и S
рол v г возм.кон.рез
возм.кон.
ЗД)={г е _ „ ц .„(Г,) №<д} (8)
Заметим, что Б (Г) представляет собой область «геометрического контакта», в которой возможны подобласти отсутствия «физического контакта» (т.е. отсутствия взаимодействия геометрически касающихся поверхностей друг с другом). В дальнейшем будем устанавливать граничные условия для некоторого фиксированного момента времени, поэтому индекс ^ будем опускать. Выделим в области ЗС две зоны:
={ге Л'й^Дг.О =0)} - зона «отлипания» (9.1)
х'. = (гея(|(гг,(г,н)с0)| — зона контакта (9.2)
где оп=п о п — нормальное напряжение. При скоростной постановке задачи подразделяется еще на две подобласти, в которых скорость нормальных напряжений либо оп < о, либо оп =0, первая из которых в дальнейшем относится к З0с, а на второй задаются граничные условия свободной поверхности.
Зона Зс также подразделяется на две части: зону прилипания Зса и зону проскальзывания Зс5. Не уменьшая общности, будем считать справедливым закон трения Зибе-ля-Кулона. Тогда
где Г — коэффициент трения, принимаемый в данном случае 0,2,
Рис. 4 — Исследуемая область
- касательная составляющая (к поверхности ролика) вектора напряжения,
— величина модуля вектора действующего касательного напряжения,
= <т — предел текучести
материала на сдвиг.
Для получения приближенного решения поставленной задачи (1)-(ю) был использован программный комплекс ЛБАЦиБ, использующий традиционный для механики деформируемого твердого тела метод конечных элементов [3, 4]. В частности, в силу существенной нелинейности задачи был
использован модуль Abaqus/Explicit. Для описания исследуемой области конечными элементами была построена конечно-элементная сетка (рис. 5), состоящая из 44232 элементов С3081^ и 39593 узлов.
Для исследования распределения остаточных напряжений после обкатки резьбы роликом были проведены вычислительные эксперименты при различных величинах вдавливания ролика Результаты исследований были внедрены в производство [5].
Итоги
С точки зрения механики деформируемого твердого тела сделана постановка трехмерной нестационарной контактной задачи упругопластического деформирования.
Рис. 5 — Конечно-элементная сетка
ENGLISH
Выводы
Подводя итог можно заключить, что в процессе обкатывания резьбы роликом в приповерхностных слоях материала межвитковой впадины резьбы формируются области сжимающих напряжений. Данное обстоятельство, без сомнения, приводит к упрочнению приповерхностного слоя и препятствует возникновению микротрещин. С другой стороны, в результате данного процесса в подповерхностных слоях материала резьбы могут формироваться области значительных растягивающих напряжений, что может приводить к возникновению внутренних микро- и макротрещин, приводящих к разрушению конструкции. Подбирая оптимальные параметры процесса (геометрия ролика и сила его прижатия), можно получать благоприятные для данной конструкции распределения остаточных напряжений.
Список используемой литературы
1. Поздеев А.А., Няшин Ю.И., Трусов П.В. Остаточные напряжения: Теория и приложения. М.: Наука, 1982. 112 с.
2. Поздеев А.А., Трусов П.В., Няшин Ю.И. Большие упругопластические деформации: теория, алгоритмы, приложения. М.: Наука, 1986. 232 с.
3. Зенкевич О. Метод конечных элементов в технике. М.: Мир, 1975. 542 с.
4. Сегерлинд Л. Применение метода конечных элементов. М.: Мир, 1979. 392 с.
5. Песин М.В. Повышения надежности бурильных труб на стадии проектирования путем использования математического моделирования процесса упрочнения резьбовой поверхности // Экспозиция Нефть Газ. 2013. №2 (27). С. 56-57.
OIL PRODUCTION
Scientific the bases of the simulation of the process of strengthening the bottom of thread of drill pipes by deep roll
Authors:
Mikhail V. Pesin — dr. sci. tech., senior lecturer1; deputy director2; M.Pesin@mail.ru
•State National Research Polytechnical University of Perm, Perm, Russian Federation 2CJSC Trading house PKNM (Perm Oil Engineering Company) , Perm, Russian Federation
Abstract
UDC 621.77.07
Performed computer simulations of the hardening of the hollow carving deep roll. Provides conceptual and mathematical formulation of the problem. As a result of a series of numerical experiments using ABAQUS package evaluation of residual stresses distributions are obtained depending on the degree mounting of the roll.
Materials and methods
Methods of mathematical modelling are used.
Results
From the point of view of the mechanics solid of body being deformed is made setting the three-dimensional nonstationary contact problem of elastic-plastic deformation.
Conclusions
It is possible to conclude that in the process of deep roll thread by roller in the near-surface layers of the material of interturn bottom of thread the regions of the compressive stresses are formed. This circumstance undoubtedly leads to strengthening of near-surface layer and it prevents the appearance of microscopic
cracks. From other side, as a result this process in the sub-surface layers of the material of thread can be formed the regions of significant tensile stresses, which can lead to the appearance of the internal micro- and of the macrofissures, which lead to the structural failure. Selecting the optimum parameters of process (geometry of roller and the force of its pressing), it is possible to obtain the favorable for the given construction distributions of residual stresses.
Keywords
hardening of deep roll thread, residual stresses, finite element method
References
1. Pozdeev A.A., Nyashin Yu.I., Trusov P.V. Ostatochnye napryazheniya: Teoriya i prilozheniya [Residual stresses: Theory and application]. Moscow: Nauka, 1982. 112 p.
2. Pozdeev A.A., Trusov P.V., Nyashin Yu.I. Bol'shie uprugoplasticheskie deformatsii: teoriya, algoritmy, prilozheniya [Large elastic-plastic deformations: theory,
algorithms, application]. Moscow: Nauka, 1986, 232 p.
3. Zenkevich O. Metod konechnykh elementov vtekhnike [Method of final elements in the technology]. Moscow: Mir, 1975, 542 p.
4. Segerlind l. Primenenie metoda konechnykh elementov [Application of a method of final elements]. Moscow: Mir, 1979, 392 p.
5. Pesin M.V. Povysheniya nadezhnosti
buril'nykh trub na stadiiproektirovaniya putem ispol'zovaniya matematicheskogo modelirovaniya protsessa uprochneniya rez'bovoy poverkhnosti [Improving reliability of drill pipe at the design stage through the use of mathematical modelling of the process of hardening the threaded surface] Exposition Oil Gas, 2013, issue 2, pp. 56-57.