|
Технически университет гр. Габрово
|
СТУДЕНТСКА НАУЧНА СЕСИЯ
20-21 октомври 2011 година
|
МАТЕМАТИЧЕН МОДЕЛ нА ОСНОВНИТЕ ХАРАКТЕРИСТИКИ НА високоточЕН цангов патронник за револверни стругове и многовретенни хоризонтални прътови автомати
инж. Радомир Атанасов
Технически университет - Габрово
Катедра МТТ
Научен ръководител: доц. д-р инж. Пенка Неделчева
Резюме: Разработен е обемен параметричен изчислителен модел на високоточен цангов патронник, предназначен за револверни стругове и многовретенни хоризонтални прътови автомати, като са симулирани основните му характеристики с прилагане метода на крайните елементи. Проведен е планиран многофакторен числен експеримент с цел създаване на математичен модел.
Ключови думи: CAD/CAE, цангов патронник, метод на крайните елементи, статичен анализ, математичен модел, стъпкова регресия
ВЪВЕДЕНИЕ
При цанговите затягащи механизми предаването и затварянето на силовия поток се свежда до преобразуване на осовата сила от превода на машината в радиална сила на затягане на заготовката от затягащите елементи.
На фиг. 1 е показана конструкция на високоточен цангов патронник (ВЦП) [1], който има стандартна затегателна цанга ЗЦ и междинна цанга МЦ с напречни пера, разположени в равнина, перпендикулярна на оста на патронника. Челюстта на ЗЦ преобразуваща осовото движение под въздействие на теглещото усилие SΣ от превода в радиално преместване на челюстта на МЦ, осъществяваща радиална сила на затягане ТΣ на заготовката (TΣ = Т z, където Т е радиална сила на затягане, създавана от една челюст, а z – брой челюсти). Равнините на симетрия на челюстите на двете цанги – затягаща и междинна се сливат.
Напречните пера на МЦ са прави, формирани чрез изрези във фланеца и осъществяват еластичната връзка с челюстите. Те имат затворен контур и в съществена степен влияят на основните характеристики на тези патронници.
Основни задачи по усъвършенстване на ВЦП са повишаване стабилността на затягащата радиална сила Т, осовото изместване на заготовката в процес на затягане и радиалната стабилност на системата патронник-детайл в рамките на допуските на диаметралния размер на затяганата заготовка при минимални напрежения в перата на цангата, които водят до загуба на еластичните им свойства [1].
а б
Фигура 1. Конструктивна схема на ВЦП с МЦ: Вр. – вретено; ЗЦ – затегателна цанга; МЦ – междинна цанга; ЕФ – еластичен фланец; Заг. – затягана заготовка
Цел на настоящата работа е създаване на математичен модел, описващ основните характеристики на високоточен цангов патронник, предназначен за револверни стругове и многовретенни хоризинтални прътови автомати.
ИЗЛОЖЕНИЕ
Изследваната конструкция ВЦП (фиг. 1) е със ЗЦ изп. II 93 x d 65 mm БДС 6100-74, работеща на опън. Разработването на конструкцията на МЦ, като нестандартизиран елемент от ВЦП, се подчинява на следните ограничения (вж. фиг. 1): еластичната част на фланеца на МЦ се разполага в пръстеновидна област с размери D и d (фиг. 1б), съответстващи на стандартните размери D и d на ЗЦ изп. ІІ от БДС 6100-74 (фиг. 1а); прорезите във фланеца, с цел по-голяма еластичност на перата, трябва да бъдат с максимална дължина, т.е. радиусът им на закръгление тангира с външният граничен диаметър D на пръстеновидната област (фиг. 1б).
В среда на хибридната параметрична система SolidWorks Simulation [5] е създаден обемен параметричен крайно-елементен изчислителен модел на патронника [1, 2], като са наложени разгледаните по-горе геометрични ограничения. По този начин се намалява броят на независимите параметри, управляващи модела.
Изследвани са следните основни характеристики на ВЦП с МЦ (вж. фиг. 2) за диаметър на затяганата заготовка 65 mm:
- максималното еквивалентно напрежение (по IV якостна теория) в перата на междинната цанга екв. мах ;
- радиалната сила на затягане Т, създавана от една челюст, се отчита като контактна сила, възникваща по повърхнината на затяганата заготовка. Въз основа на нея се определя разсейването на коефициента на усилване на патронника – Ку = T / S, където S = 15000 N е осовата теглеща сила действаща на една челюст;
- осовото изместване x на заготовката при затягане (в следствие увличането ѝ) се отчита спрямо ос (например ос х), която е успоредна на оста на патронника;
- отчитаните радиални премествания yр на затегнатата заготовка, възникващи вследствие приложена радиална сила Fр, позволяват да бъде определена радиалната стабилност на системата патронник-детайл jp = Fр / yр и нейното разсейване в зависимост от параметрите на патронника.
Фигура 2. Изчислителен модел на основни характеристики на ВЦП с МЦ:
Т - радиална сила на затягане; Fp – радиална сила приложена върху заготовката; екв. мах - максимално еквивалентно напрежение в перата на МЦ; x - осово изместване на пръта при затягане; Kу - коефициент на усилване на патронника; yp - радиалнo преместване на затегнатата заготовка
Еластичните премествания на заготовката се отчитат на определено разстояние от патронника, което в случая е 120 mm. Радиалната сила Fр= 5000 N е приложена на същото разстояние.
ПЛАНИРАНЕ НА МНОГОФАКТОРЕН ЧИСЛЕН ЕКСПЕРИМЕНТ
За реализирането на математичния модел на ЗМ е избран ротатабилен централно-композиционен план на експеримента (РЦКП) от типа 23 с ядро на пълен факторен експеримент (ПФЕ) [3]. РЦКП осигурява еднаква дисперсия на равни разстояния от центъра на плана (табл. 1). Чрез подходящ избор на повторните наблюдения в центъра на плана може да се получи почти равномерно разпределение на дисперсията в цялата област, ограничена от сфера (или окръжност) с единичен радиус [4]. Този тип планиране се прилага за получаването на регресионни математични модели от втори ред, а също и за описание на локални участъци от факторното пространство най-често при брой на факторите n 4.
Таблица 1. Данни необходими за построяване на плана на експеримента
n
|
N1
|
Nα
|
N0
|
N
|
звездно рамо α
|
ядро
|
брой
фактори
|
брой опити за ядрото
|
брой на звездните точки
|
брой на
централните точки
|
общ брой
на опитите
|
3
|
8
|
6
|
6
|
20
|
1,682
|
ПФЕ
|
Нелинейните математични модели се използват в случаите когато не е възможно да се получи адекватен линеен модел. С моделите от по-висок ред се описва по-точно повърхнината на реакциите и областта на оптимума е по-подробна. Броят на наблюденията в центъра на плана N0 се избира така, че да се получи униформ-планиране. Общият вид на нелинеен полиномен модел от втори ред, е [3, 4]:
. (1)
В разглеждания случай, управляващи параметри са геометричните размери, формиращи еластичната част на фланеца на МЦ - b, l, t (фиг. 3а), които са основни фактори на взаимодействие в кибернетичния модел, а изходни параметри на изследването (целеви функции) са посочените по-горе основни характеристики (екв. max, x, Ky, jp) (фиг. 3б).
a
|
б
|
Фигура 3. Управляващи параметри (а) и кибернетичен модел (б) на многофакторния експеримент
Изборът на локална област за провеждане на многофакторния експеримент, се състои в определянето на основните нива и интервалите на вариране на управляващите параметри.
Таблица 2. Локална област на факторното пространство
Характеристики
|
Фактори
|
наименование
|
кодирана стойност
|
b
|
l
|
t
|
mm
|
mm
|
mm
|
основно ниво
|
0
|
5.75
|
14
|
19
|
интервал на вариране
|
|
0.25
|
0.5
|
1
|
горно ниво
|
+1
|
6
|
14.5
|
20
|
долно ниво
|
-1
|
5.5
|
13.5
|
18
|
кодово означение
|
xi
|
X1
|
X2
|
X3
|
Проведени са експерименти за уточняване на диапазоните на вариране на входните фактори съгласно избрания план на експеримента (табл. 2). Следва съставяне на матрицата на експеримента (табл. 3) за РКЦП при брой на факторите n = 3.
Таблица 3. Матрица в кодиран вид на РКЦП за n = 3 и опитни стойности
N
|
x1
|
x2
|
x3
|
x1x2
|
x1x3
|
x2x3
|
x12
|
x22
|
x32
|
y1
|
y2
|
y3
|
y4 (*106)
|
1
|
+1
|
+1
|
+1
|
+1
|
+1
|
+1
|
+1
|
+1
|
+1
|
189,47
|
0,0364
|
1,698
|
2,7543657
|
2
|
-1
|
+1
|
+1
|
-1
|
-1
|
+1
|
+1
|
+1
|
+1
|
181,76
|
0,0421
|
1,711
|
1,4575135
|
3
|
+1
|
-1
|
+1
|
-1
|
+1
|
-1
|
+1
|
+1
|
+1
|
204,74
|
0,0402
|
1,700
|
1,3716292
|
4
|
-1
|
-1
|
+1
|
+1
|
-1
|
-1
|
+1
|
+1
|
+1
|
183,98
|
0,0398
|
1,692
|
1,4988459
|
5
|
+1
|
+1
|
-1
|
+1
|
-1
|
-1
|
+1
|
+1
|
+1
|
214,36
|
0,0468
|
1,680
|
1,4391181
|
6
|
-1
|
+1
|
-1
|
-1
|
+1
|
-1
|
+1
|
+1
|
+1
|
200,4
|
0,0534
|
1,714
|
0,8914245
|
7
|
+1
|
-1
|
-1
|
-1
|
-1
|
+1
|
+1
|
+1
|
+1
|
220,18
|
0,0473
|
1,686
|
1,4665982
|
8
|
-1
|
-1
|
-1
|
+1
|
+1
|
+1
|
+1
|
+1
|
+1
|
204,72
|
0,0518
|
1,701
|
1,9504584
|
9
|
+1,682
|
0
|
0
|
0
|
0
|
0
|
+2,829
|
0
|
0
|
203,25
|
0,0434
|
1,694
|
2,1443122
|
10
|
-1,682
|
0
|
0
|
0
|
0
|
0
|
-2,829
|
0
|
0
|
190,01
|
0,04899
|
1,718
|
1,8307306
|
11
|
0
|
+1,682
|
0
|
0
|
0
|
0
|
0
|
+2,829
|
0
|
195,97
|
0,04319
|
1,680
|
1,3495641
|
12
|
0
|
-1,682
|
0
|
0
|
0
|
0
|
0
|
-2,829
|
0
|
194,86
|
0,0455
|
1,717
|
1,5360275
|
13
|
0
|
0
|
+1,682
|
0
|
0
|
0
|
0
|
0
|
+2,829
|
180,19
|
0,0358
|
1,698
|
1,0725241
|
14
|
0
|
0
|
-1,682
|
0
|
0
|
0
|
0
|
0
|
-2,829
|
214,64
|
0,0527
|
1,713
|
2,7739251
|
15
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
204,68
|
0,0421
|
1,709
|
2,7634997
|
16
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
198,64
|
0,04259
|
1,711
|
1,6307893
|
17
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
205,86
|
0,04219
|
1,707
|
1,9969646
|
18
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
199,7
|
0,04342
|
1,685
|
1,7831669
|
19
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
197,62
|
0,04375
|
1,683
|
1,4853103
|
20
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
0
|
204,33
|
0,04209
|
1,707
|
2,6946915
|
След провеждане на планирания експеримент се отчитат стойностите на четирите характеристики на ЗМ, които също са представени в табл. 3.
Емпиричните данни са получени посредством линейни статични анализи, с използване на метода на крайните елементи, съгласно плана на експеримента.
МАТЕМАТИЧЕН МОДЕЛ НА ОСНОВНИТЕ ХАРАКТЕРИСТИКИ НА ЦАНГОВИЯ ПАТРОННИК
Емпиричните данни са обработени с помощта на софтуерния продукт STATISTICA. В резултат са получени уравненията на регресия за основните характеристики на затягащия механизъм (ЗМ). Като автоматизирана процедура за избор на коефициенти от началния модел е използвана стъпкова регресия. Уравненията са представени в неявен вид и съдържат само значимите коефициенти на регресия:
y1 = 199,2482 + 7,1449x1 – 3,5451x2 – 10,0785x3 + 2,3807x22 , (2)
y2 = 0,043992 – 0,002205x1 – 0,005069x3 – 0,001012 x1 x2 , (3)
y3 = 1,702561 – 0,006888x1 – 0,005968x22 , (4)
y4 = (1,943488 – 0,306953x1x2 – 0,30349x2x3 + 0,3721872x32)106 , (5)
където с y1, y2, y3 и y4 са означени целевите функции (параметрите на оптимизация) съответно - екв. max, x, Ky и jр, а с x1, x2 и x3 – входните фактори в кодиран вид.
Стойностите за критерия на Фишер F, коефициента на детерминираност (определеност) R2 и средната процентна грешка при равнище на значимост α = 0,1 за уравненията (2), (3), (4) и (5) са поместени в табл. 4:
Таблица 4. Условия за адекватност
целева функция
|
критерий на Фишер
|
коеф. на определеност
|
процентна грешка
|
изчислена стойност
|
таблична стойност
|
y1
|
F(6,13)= 24,7983
|
Fтаб = 2,283
|
R2 = 0,9196
|
3,68 %
|
y2
|
F(5,14)= 39,404
|
Fтаб = 2,307
|
R2 = 0,9337
|
0,00145 %
|
y3
|
F(4,15)= 5,541
|
Fтаб = 2,362
|
R2 = 0,5964
|
0,00905 %
|
y4
|
F(5,14)= 7,2155
|
Fтаб = 2,307
|
R2 = 0,6257
|
0,406 %
|
Определянето на степента на влияние на отделните фактори върху всяка от характеристиките на ЗМ се извършва чрез последователно изключване на факторите. По този начин се следи промяната на коефициента на детерминираност, като на факторът оказващ най-голямо влияние съответства най малката стойност на R2. Резултатите от сравнението са представени в табл. 5, като факторите с най-силно влияние са подчертани:
Таблица 5. Определяне степента на влияние на входните фактори
целева функция
|
коефициент на определеност R2
|
при изключен фактор x1
|
при изключен фактор x2
|
при изключен фактор x3
|
y1
|
0,72
|
0,8727
|
0,5203
|
y2
|
0,8467
|
0,9392
|
0,4011
|
y3
|
0,4601
|
0,6312
|
0,6342
|
y4
|
0,6271
|
0,6511
|
0,6222
|
Използвайки получените резултати за факторите с най-силно влияние са построени линиите на еднакъв отклик за всяка характеристика (фиг. 4):
y1 = f (x1, x3)
|
y2 = f (x1, x3)
|
y3 = f (x1, x2)
|
y4 = f (x1, x3)
|
Фигура 4. Линии на еднакъв отклик
ЗАКЛЮЧЕНИЕ
В резултат на проведеното многофакторно изследване работата на ЗЦ в неявен вид са изведени уравненията на регресия за разгледаните основни характеристики, отнасящи се за конкретния типоразмер ВЦП.
Фактор оказващ най-голямо влияние върху напрежението в перата на МЦ и върху осовото и радиалното премествания на заготовката е ширината t, а височината на перото b влияе най-силно върху коефициента на усилване на патронника.
Получените математични модели съдържат само значимите фактори и взаимодействия, като за получаването им е използвана стъпкова регресия.
За валидиране на резултатите от многофакторното изследване е необходимо провеждане на реален експеримент.
ЛИТЕРАТУРА
[1] Кузнецов Ю. Н., В. Н. Волошин, П. М. Неделчева, Ф. В. Эль-Дахаби, Зажимные механизмы для высокопроизводительной и высокоточной обработки резанием, Университетско издателство „Васил Априлов”, Габрово 2010, 723 с.
[2] Неделчева П. М., Р. Ст. Атанасов, Динамични характеристики на цангови патронници с напречни пера, Сборник доклади, XIX ННТК „АДП-2010” Созопол, с. 340 – 346
[3] Александрова И. Ст., Основи на инженерните изследвания, УИ „Васил Априлов”, Габрово 2003, 236 с.
[4] Божанов Е., И. Вучков, Статистически методи за моделиране и оптимизиране на статистически обекти, ДИ „Техника”, София 1973, 529 с.
[5] SolidWorks Simulation Professional 2009, Dassault Systemes SolidWorks Corporation, Concord, Massachusetts, USA 2009, 506 p.
Изследванията са подкрепени по договор № BG051PO001-3.3.04/28, „Подкрепа за развитие на научните кадри в областта на инженерните научни изследвания и иновациите”. Проектът се осъществява с финансовата подкрепа на Оперативна програма „Развитие на човешките ресурси” 2007-2013, съфинансирана от Европейския социален фонд на Европейския съюз“.
Сподели с приятели: |