Utafiti wa kitaaluma, lugha inayoeleweka

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

27 Septemba 2026, Jumapili
VERİANLAUchapishaji huru wa sayansi
Fungua au funga menyu
...
Home / Sayansi Tumizi / Uhandisi / Utoaji wa Trajektoria ya Rola ya Uundaji kwa Njia ya Projection na Fidia ya Hitilafu kwa Regression ya Mchakato wa Gaussian
Uhandisi

Utoaji wa Trajektoria ya Rola ya Uundaji kwa Njia ya Projection na Fidia ya Hitilafu kwa Regression ya Mchakato wa Gaussian

Utafiti huu unatengeneza mbinu ya kupima na kupunguza mkengeuko wa kimfumo unaotokea wakati trajektoria ya mwendo wa rola inayopatikana kutoka simulizi ya vipengele finiti inapohamishwa kwenda kwenye mashine halisi ya uundaji wa metali kwa rola.

11/08/2026  Veri Anla Imetazamwa mara 21
Utoaji wa Trajektoria ya Rola ya Uundaji kwa Njia ya Projection na Fidia ya Hitilafu kwa Regression ya Mchakato wa Gaussian

Utafiti huu unatengeneza mbinu ya kupima na kupunguza mkengeuko wa kimfumo unaotokea wakati trajektoria ya mwendo wa rola inayopatikana kutoka simulizi ya vipengele finiti inapohamishwa kwenda kwenye mashine halisi ya uundaji wa metali kwa rola. Katika simulizi ya ABAQUS, ganda dogo la 1 × 1 × 1 mm liliunganishwa kwenye rola iliyofafanuliwa kama mwili mgumu wa kianalitika; koordinati za nodi za ganda hili zilitolewa fremu kwa fremu kwa kutumia Python na trajektoria ya kituo cha rola ikahesabiwa. Trajektoria ya simulizi ilibadilishwa kwenda kwenye mfumo wa koordinati wa mashine halisi; mwendo halisi wa mashine ya uundaji yenye rola moja ulipimwa kwa kihisi cha uhamishaji cha waya wa kuvuta. Katika majaribio matatu yasiyo na fidia, hitilafu ya wastani ya projection ilikuwa 0,503 mm, ilhali baada ya fidia ya mwelekeo wa normal iliyoundwa kwa njia ya projection na Gaussian Process Regression, hitilafu ya wastani katika uendeshaji mpya wa mashine iliripotiwa kuwa 0,229 mm. Chanzo kinatoa mabadiliko haya kama punguzo la hitilafu la %54,47. Hata hivyo, utafiti hauonyeshi kwamba mbinu ya fidia iliboresha unene wa ukuta, ubora wa uso au usahihi wa vipimo wa sehemu iliyoundwa kwa kiwango hicho hicho; uthibitishaji wa majaribio ulilenga hasa trajektoria ya mwendo wa rola.

Lengo la utafiti si kutafuta “trajektoria bora zaidi ya uundaji” kwa maana ya kawaida. Utafiti unachunguza kwa nini trajektoria iliyofafanuliwa kwenye modeli ya vipengele finiti haifuati kabisa njia ileile ya kijiometri inapohamishwa kwenda kwenye mwendo halisi wa CNC/mashine, na jinsi tofauti hiyo kati ya simulizi na mashine inaweza kufidiwa. Badala ya kuunda modeli tofauti za tofauti za X au Z, waandishi hutafuta projection iliyo karibu zaidi ya kila nukta iliyopimwa kwenye mviringo bora na kufafanua hitilafu katika mwelekeo wa normal wa ndani wa mviringo huo.

Gaussian Process Regression ilitumiwa kuunda modeli ya jinsi hitilafu ya mwelekeo wa normal inavyobadilika kwenye trajektoria. Sampuli 370 za projection kutoka katika uendeshaji mitatu ya mashine bila fidia ziliunganishwa, kisha zikapunguzwa hadi sampuli 102 za mafunzo kwa kuchukua wastani ndani ya visanduku vya upana wa 1,0 mm kwenye urefu wa arc. Katika modeli ya GPR, kernel ya Matérn (\(\nu=1{,}5\)), constant kernel na white-noise term zilitumiwa pamoja; hitilafu ya kimfumo iliyotabiriwa na modeli ilitumika katika mwelekeo wa normal ulio kinyume na trajektoria ya simulizi ili kuunda trajektoria mpya yenye fidia.

Uchambuzi wa kitakwimu wa matokeo ya fidia haukutegemea tu thamani ya wastani. Chanzo kililinganisha mgawanyo wa hitilafu kamili za projection kabla na baada ya fidia kwa Welch t-test na Mann–Whitney U test. Thamani za p zilizoripotiwa ni \(2{,}442\times10^{-4}\) na \(2{,}040\times10^{-5}\) mtawalia, ambazo ziko chini ya kizingiti cha umuhimu \(\alpha=0{,}05\) kilichotumika katika utafiti. Matokeo haya yanaunga mkono kwamba katika mashine na hali ya trajektoria zilizochunguzwa, punguzo la hitilafu lililopimwa ni vigumu kuelezwa kwa mabadiliko ya nasibu ndani ya data ya majaribio; punguzo lilelile halijaonyeshwa bado kwa mashine au hali nyingine za uundaji.

Kwa upande wa Türkiye, utafiti huu ni muhimu kwa kuonyesha kwamba simulizi ya vipengele finiti inaweza kutumika si kama zana ya kubuni pekee, bali pia kuunda mnyororo wa uthibitishaji uliofungwa unaotegemea vipimo wakati matokeo yanahamishwa kwenda kwenye vifaa halisi vya CNC/uundaji wa metali. Katika matumizi ya metal spinning na precision forming nchini Türkiye, mbinu sawa inaweza kubadilishwa ili kubaini mkengeuko wa kimfumo kati ya trajektoria ya simulizi na mwendo halisi wa mhimili. Hata hivyo, utafiti huu ulifanywa kwa mashine moja tu ya spinning iliyotengenezwa na watafiti wenyewe, familia moja ya trajektoria ya koni mbili na idadi ndogo ya marudio ya majaribio; pia inaelezwa wazi kuwa jaribio la kulinganisha trajektoria lilifanywa bila workpiece. Kwa hiyo, matokeo ya 0,229 mm hayapaswi kuhamishwa moja kwa moja kwenda kwenye mashine nyingine au toleransi za sehemu za mwisho.

Tatizo kuu la utafiti ni nini?

Katika finite element analysis, mwendo wa rola unaweza kufafanuliwa kihisabati juu ya mviringo sahihi. Mashine halisi, hata hivyo, hutekeleza trajektoria hiyo kupitia CNC interpolation, mabadiliko ya koordinati, mifumo ya uhamishaji na mwendo halisi wa mihimili. Kwa hiyo, tofauti za kimfumo zinaweza kujitokeza kati ya njia ya kijiometri iliyoundwa kwenye simulizi na mwendo halisi wa rola.

Pengo la fasihi linalolengwa na watafiti ni tofauti na kuboresha umbo la trajektoria katika metal spinning. Swali kuu ni jinsi ya kuhamisha kwa kuaminika trajektoria ya kituo cha rola kutoka simulizi ya vipengele finiti kwenda kwenye mashine halisi na jinsi ya kufidia hitilafu ya kijiometri ya eneo kati ya trajektoria inayotekelezwa na mashine na trajektoria ya simulizi.

Chanzo kinasisitiza hasa kwamba mbinu za kawaida za hitilafu kwa vijenzi vya koordinati zinaweza kutowakilisha kikamilifu mkengeuko wa kijiometri wa eneo kwenye trajektoria zilizopinda kwa kuendelea. Kwa hiyo, hitilafu ilifafanuliwa kama hitilafu moja ya kijiometri yenye ishara katika mwelekeo wa normal wa eneo wa mviringo, badala ya tofauti mbili tofauti za koordinati kwenye mihimili X na Z.

Modeli ya vipengele finiti inawakilisha sehemu gani?

Katika mazingira ya ABAQUS 2022, modeli ya spinning yenye rola moja ya sehemu ya koni mbili iliundwa. Modeli ina vipengele vitatu vikuu:

  • mandrel/core die,
  • blank ya mviringo ya aloi ya aluminium 6061,
  • rola ya spinning ya aina α.

Kipenyo cha mwisho mdogo wa mandrel ni 30 mm, kipenyo cha mwisho mkubwa 160 mm, nusu-pembe ya koni ya kwanza 30°, nusu-pembe ya koni ya pili 45°, urefu wa jumla 78,39 mm na radius ya arc ya mpito kati ya sehemu mbili za koni 50 mm.

Blank ya mviringo ina kipenyo cha ndani 15 mm, kipenyo cha nje 150 mm na unene 2,2 mm.

Kwa jiometri ya rola, sehemu husika ya modeli inatoa working angle ya 30°, corner radius ya 5 mm na kipenyo cha rola cha 200 mm.

Mesh ya vipengele finiti iliundwaje?

Blank ya aluminium 6061 iliwekewa mesh kwa vipengele vya hexahedral vya nodi nane C3D8R vyenye reduced integration. Modeli ya mwisho ina:

  • vipengele 5850,
  • nodi 11.960

Njia ya adaptive mesh ilitumika kuboresha msongamano wa mesh kwenye uso wa nje wa blank. Kwa kuwa rola na mandrel zilichukuliwa kuwa analytical rigid bodies, mesh ya finite element haikuundwa kwa vipengele hivi viwili.

Masharti ya contact yalifafanuliwaje?

Jozi mbili kuu za surface-to-surface contact ziliundwa: rola–blank na mandrel–blank. “Hard contact” ilitumika katika mwelekeo wa normal, na penalty friction model katika mwelekeo wa tangential.

Jozi ya contactMgawo wa msuguano
Rola – blank0,05
Mandrel – blank0,10

Mandrel inapozunguka kwa angular velocity thabiti, rola hufanya radial infeed na axial feed kando ya trajektoria iliyofafanuliwa mapema.

Ni thamani gani za nyenzo zilitumika kwa aluminium 6061?

Sifa ya nyenzoThamani iliyotolewa kwenye chanzo
Moduli ya elasticity68,9 × 109 Pa
Msongamano2730 kg/m³
Poisson ratio0,33
Brinell hardness95 HB

Kwa nini kuna kutokulingana ndani ya chanzo kuhusu vigezo vya simulizi?

Maandishi yanayoeleza vigezo vya simulizi yanataja roller installation angle ya 30°, corner radius ya 5 mm, spindle speed ya 120 r/min na feed ratio ya 0,4 mm/r. Hata hivyo, Table 3 katika sehemu hiyo hiyo inatoa vigezo viwili vya kwanza kama 25° na 6 mm mtawalia.

KigezoMaandishi ya maelezoTable 3
Roller installation angle30°25°
Roller corner radius5 mm6 mm
Spindle speed120 r/min120 r/min
Feed ratio0,4 mm/r0,4 mm/r

Chanzo hakielezi ni seti ipi kati ya hizi mbili za angle/radius iliyotumika katika simulizi ya mwisho. Kwa hiyo, katika maelezo ya Verianla hakuna seti iliyochaguliwa kimya kimya kama “thamani sahihi”.

Kituo cha analytical rigid roller kilitolewaje?

Kwa kuwa trajektoria endelevu ya kituo cha kijiometri cha rola kilichofafanuliwa kama analytical rigid body katika ABAQUS haikuweza kuchukuliwa moja kwa moja kutoka output ya kawaida, watafiti walitumia ganda dogo la 1 × 1 × 1 mm kama alama msaidizi.

Ganda hili dogo liliunganishwa na rola kwa tie constraint na kusogezwa rigidly pamoja na rola. Koordinati za nodi moja ya kona yenye sifa maalum kwenye ganda zilirekodiwa katika kila fremu ya uchambuzi.

COORD field output ilisomwa fremu kwa fremu kutoka ABAQUS ODB database kwa script ya Python 3.13; node number na three-dimensional coordinates ziliandikwa kwenye text file.

Nodi yenye sifa:

\[ A_i=(x_{A_i},y_{A_i},z_{A_i}) \]

na kituo cha rola:

\[ P_i=(x_{P_i},y_{P_i},z_{P_i}) \]

zinapofafanuliwa, chanzo kinatoa uhusiano mkuu kama:

\[ P_i=A_i+(\Delta x,\Delta y,\Delta z) \]

kwa namna ilivyoainishwa.

Geometric offsets zilizochapishwa ni:

\[ \Delta x=\frac{a}{2}+\frac{d}{2}\cos\theta \]

\[ \Delta y=\frac{a}{2} \]

\[ \Delta z=\frac{d}{2}\sin\theta+h \]

Hapa \(a\) ni edge length ya discrete shell, \(d\) ni unene wa rola, \(\theta\) ni roller installation angle na \(h\) ni additional geometric offset.

Baada ya kutoa \(a=1\) mm, \(d=30\) mm na \(\theta=30^\circ\), chanzo kinaandika transformation ya kituo cha rola kama \(P_1=(x_A+13{,}433,\ y_A-22{,}267,\ z_A+0{,}5)\). Kwa kuwa offset hii ya nambari haiwezi kuzalishwa moja kwa moja kwenye mihimili yote kutoka Equations (2)–(7) zilizochapishwa, inapaswa kutazamwa kama utata wa ndani wa chanzo kuhusu coordinate transformation.

Trajektoria ya simulizi ilihamishwaje kwenda kwenye mashine halisi?

Koordinati za kituo cha rola zilizotolewa kwanza zilibadilishwa kutoka mfumo wa koordinati wa ABAQUS kwenda kwenye mfumo wa koordinati wa mashine. Watafiti wanaeleza kwamba hili lilifanywa kwa kulinganisha origins za koordinati na kufanya spatial positioning inayolingana na jiometri ya mashine.

Table 5 inaonyesha kuwa trajektoria iliyobadilishwa ilikuwa na pointi 401 za data. Kisha, kutokana na hitaji la incremental interpolation la mashine ya CNC, miondoko kati ya pointi jirani ilichakatwa kwa first-order differencing.

Hapa ni muhimu kuzingatia uwasilishaji wa unit. Mihimili ya Figure 10 imeandikwa “m”, wakati katika majedwali na grafu za majaribio zinazofuata koordinati za X–Z zinaonyeshwa kwa mm. Chanzo hakielezi kwa kina scale transformation kati ya uwasilishaji hizi mbili.

Trajektoria halisi ya rola ilipimwaje?

Watafiti walitumia suluhisho la kipimo la nje lisiloingilia mfumo ili kutotegemea sensors za ndani za mashine. Absolute draw-wire displacement sensor iliyochaguliwa ina output ya RS485.

Sifa ya sensorThamani iliyoripotiwa katika chanzo
Measurement range1000 mm
Resolution iliyotajwa katika chanzo1024
Displacement resolution0,098 mm
Linear accuracy±%0,1
Repeatability±%0,01
Protection classIP54
CommunicationRS485 / Modbus-RTU

Kwa kuwa chanzo hakitoi ufafanuzi kamili wa kimwili au kidijitali wa thamani “1024”, thamani hiyo haijapatanishwa na 0,098 mm resolution kupitia formula kutoka nje ya chanzo.

Mfumo wa virtual–real mapping unafanya nini?

Mfumo wa virtual–real mapping unaotumia Unity3D 2022.3019f1c1 uliundwa kwa ufuatiliaji wa trajektoria wa wakati halisi. Data za sensor zilipokelewa kupitia Modbus-RTU protocol na kuonyeshwa kwa XChart.

Kwenye interface, target trajectory kutoka finite element model na measured trajectory ya rola ya mashine halisi huonyeshwa kwenye grafu moja. Picha katika kurasa 13–16 zinaonyesha kuwa computer interface na actual machine-axis positions zilifuatiliwa pamoja.

Jaribio la trajektoria lilifanywa kwenye sehemu halisi?

Hapana; chanzo kinaeleza wazi kwamba jaribio la trajectory fitting lilifanywa bila workpiece. Lengo katika hatua hii lilikuwa kupima kama rola inaweza kufuata simulation trajectory kwenye mashine. Red frame inayoonekana kwenye picha husika pia inawakilisha scale reference, si sehemu.

Kwa hiyo, maana ya kisayansi ya matokeo ya jaribio imezuiliwa kwenye swali “mwendo wa axis/rola wa mashine halisi unakaribia kwa kiasi gani njia ya kijiometri ya simulizi?” Unene wa ukuta, surface roughness au geometric tolerance ya final formed part si outcome variables katika jaribio hili.

Majaribio matatu kabla ya fidia yalionyesha nini?

Trajektoria ya mashine iliendeshwa mara tatu na ilionekana kwamba actual trajectory na simulation trajectory zilikuwa na tofauti ndogo lakini za kimfumo, hasa kwenye eneo la koni. Kwenye grafu, mwendo halisi unaonekana stepped, wakati simulation trajectory ni smooth curve.

Jaribio lisilo na fidiaHitilafu ya wastani ya projection (mm)
Jaribio la 10,458
Jaribio la 20,504
Jaribio la 30,546
Wastani wa jumla0,503

Ukaribu lakini kutokuwa sawa kabisa kwa wastani wa majaribio matatu unaonyesha kwamba trajectory deviation ina systematic repeatable component pamoja na components zinazobadilika kati ya majaribio.

Kwa nini projection method ilitumika?

Kwa sababu actual machine trajectory ni stepped, kuchukua simple vertical difference katika X coordinate ileile au kutathmini X na Z errors kivyake kunaweza kutowakilisha kikamilifu local geometry ya curve. Kwa hiyo, kwa kila measured point, point iliyo karibu zaidi kwenye simulation curve ilihesabiwa.

Measured point:

\[ P_i=(x_i,y_i) \]

na simulated parametric curve:

\[ C(u)=(x(u),y(u)) \]

zinapofafanuliwa, squared distance ni:

\[ D_i(u)=\|P_i-C(u)\|^2 \]

kwa namna hiyo.

Nearest-point condition ilipatikana kama:

\[ (x_i-x(u))x'(u)+(y_i-y(u))y'(u)=0 \]

na nonlinear equation ikatatuliwa iteratively kwa Newton method:

\[ u_{k+1}=u_k-\frac{f(u_k)}{f'(u_k)} \]

Mchakato huu huamua projection point ya kila measured point kwenye ideal curve.

Hitilafu ya mwelekeo wa normal ilifafanuliwaje?

Error vector relative to projection point ilifafanuliwa kama:

\[ \mathbf e=P-C(u^*) \]

ilifafanuliwa kwa namna hiyo.

Unit tangent vector ya curve ni:

\[ \mathbf t(u)= \frac{C'(u)}{\|C'(u)\|} \]

na two-dimensional unit normal vector ni:

\[ \mathbf n(u)=(-t_y,t_x) \]

hivyo signed normal error huhesabiwa kama:

\[ e(u)=\left(P-C(u)\right)\cdot\mathbf n(u) \]

huhesabiwa kwa namna hiyo.

Ufafanuzi huu unalenga kuunda modeli ya geometric deviation kwenda nje au ndani kutoka ideal trajectory, badala ya tangential component ya hitilafu.

Data set ya GPR iliundwaje?

Pointi zote za projection error kutoka majaribio matatu yasiyo na fidia ziliunganishwa. Chanzo kinatoa initial data set kama:

sampuli 370 za projection

.

Ili kupunguza athari ya stepped structure na repeated positions katika vipimo halisi kwenye GPR, arc length iligawanywa katika intervals za upana wa 1,0 mm na error values ndani ya kila interval zikapewa wastani. Baada ya hatua hii:

sampuli 102 za mafunzo ya GPR

zilibaki.

Chanzo kinaeleza wazi kwamba train/test split haikufanywa kwa sababu idadi ya data ilikuwa ndogo. Majaribio yote matatu yalitumika katika mafunzo ya combined GPR error model moja.

GPR inatumia kernel gani?

Kwa dhana kwamba trajectory error ni continuous kwa ujumla lakini inaweza kuwa na local irregularities, Matérn kernel ilitumika. Smoothness parameter ilichaguliwa kuwa:

\[ \nu=1{,}5 \]

.

Matérn kernel iliyotumika kwenye chanzo ni:

\[ k(u,u')= \sigma^2 \left( 1+\frac{\sqrt3 r}{l} \right) \exp \left( -\frac{\sqrt3 r}{l} \right) \]

na:

\[ r=|u-u'| \]

.

Hapa \(\sigma^2\) ni signal variance na \(l\) ni characteristic length scale.

White-noise term pia iliongezwa ili kuzingatia measurement noise:

\[ k_{total}(u,u')= k_{matern}(u,u') + \sigma_n^2\delta(u,u') \]

Implementation inajumuisha constant kernel + Matérn kernel + white-noise kernel. Signal variance, length scale na noise level ziliboreshwa kwa maximizing log marginal likelihood. Hata hivyo, chanzo hakiripoti numerical values za final optimized hyperparameters kwenye main text.

Trajektoria yenye fidia iliundwaje?

Ikiwa predicted mean normal error ya GPR katika arc length fulani ni \(\hat e(u)\), compensated trajectory iliundwa kama:

\[ C_{comp}(u) = C(u) - \hat e(u)\mathbf n(u) \]

.

Kwa maneno mengine, simulation trajectory ilisogezwa katika mwelekeo ulio kinyume na systematic deviation tendency ya mashine na sambamba na curve normal.

Figure 24 inaonyesha compensation amount kama curve inayobadilika kwenye trajektoria, na average compensation magnitude ya 0,3873 mm imewekwa kwenye grafu. Hii haimaanishi kwamba constant offset ya 0,3873 mm ilitumika kwenye kila point ya trajektoria; compensation inabadilika kwa nafasi.

Verianla Live: Mnyororo wa mchakato kutoka simulizi hadi trajektoria ya mashine yenye fidia

Mchakato huu unaorodhesha hatua halisi zilizotumika katika utafiti kuhamisha mwendo wa rola kutoka finite element simulation kwenda mashine halisi, kuupima na kuurekebisha kwa msingi wa data.

HatuaMaelezoChanzo
1. Finite element modelModeli ya spinning ya rola moja ya sehemu ya koni mbili ya aluminium 6061 huundwa katika ABAQUS 2022.Sehemu 2.1
2. Discrete shellGanda dogo la 1 × 1 × 1 mm huunganishwa kwenye analytical rigid roller.Figure 7 na Sehemu 2.2.1
3. Utoaji wa koordinatiData za COORD hutolewa fremu kwa fremu kutoka ODB file kwa Python 3.13.Table 4
4. Trajektoria ya kituo cha rolaKoordinati za kituo huhesabiwa kwa kutumia geometric offset kati ya characteristic node na roller center.Equations 1–7
5. Machine coordinate transformationSimulation trajectory hubadilishwa kwenda machine coordinate system na kuchakatwa kwa CNC.Tables 5–6
6. Kipimo cha trajektoria halisiMwendo halisi wa rola hurekodiwa kwa draw-wire displacement sensor na Modbus-RTU.Sehemu 3.1–3.2
7. Nearest projectionKila measured point hu-projected kwenye point iliyo karibu zaidi ya simulation curve kwa Newton method.Equations 8–16
8. Normal errorTofauti ya measurement–simulation huhesabiwa kama signed error katika local normal direction ya curve.Equations 17–19
9. GPR error modelError samples kutoka majaribio matatu huunganishwa; GPR hufunzwa kwa Matérn 1,5 na white-noise kernel.Sehemu 4.1.2
10. Normal-direction compensationGPR prediction hutumika katika mwelekeo wa normal ulio kinyume kwenye simulation trajectory ili kuunda trajektoria mpya.Equation 26
11. Uendeshaji mpya wa mashineCompensated trajectory huendeshwa tena kwenye mashine halisi na trajektoria mpya ya mwendo kupimwa.Sehemu 4.2
12. Statistical validationError distributions kabla na baada ya fidia hulinganishwa kwa Welch t-test na Mann–Whitney U test.Figure 27 na Sehemu 4.2
 

Verianla Live: Muda ni wa uchezaji wa interface pekee; hauwakilishi muda wa simulizi au jaribio la kimwili. Scientific source-of-truth ni jedwali la mbinu linaloonekana hapo juu.

Hitilafu ilipungua kwa kiasi gani baada ya fidia?

Compensated trajectory iliwekwa tena kwenye mashine na trajektoria mpya ya mwendo halisi ikapimwa. Chanzo kinaripoti average deviation baada ya fidia kuwa 0,229 mm. Combined average ya majaribio matatu kabla ya fidia ni 0,503 mm.

Verianla Live: Hitilafu ya wastani ya trajektoria kabla na baada ya GPR compensation

Ulinganisho unatumia average projection errors zilizoripotiwa moja kwa moja katika utafiti. Thamani ndogo ina maana geometric deviation ndogo kati ya simulation na measured machine trajectory.

HaliHitilafu ya wastani ya trajektoria (mm)KigezoChanzo
Kabla ya GPR compensation0,503Wastani wa jumla wa majaribio matatu yasiyo na fidiaSehemu 3.3
Baada ya GPR compensation0,229Uendeshaji mpya wa compensated trajectorySehemu 4.2
 

Verianla Live: Chanzo kinaripoti mabadiliko haya kama %54,47 average error reduction. Grafu hutumia tu source values zilizo kwenye jedwali linaloonekana.

Statistical test ilionyesha nini?

Data set ya kabla ya fidia iliundwa kutoka absolute projection errors zote za majaribio matatu yasiyo na fidia, na data set ya baada ya fidia kutoka projection errors za compensated run.

Chanzo kilitumia tests mbili tofauti:

Statistical testAlternative hypothesisp valueSignificance threshold
Welch t-testErrors baada ya fidia ni ndogo2,442 × 10−40,05
Mann–Whitney UErrors baada ya fidia ni ndogo2,040 × 10−50,05

p values zote mbili ziko chini ya 0,05. Watafiti wanatafsiri hili kama ushahidi kwamba kupungua kwa hitilafu baada ya fidia ni statistically significant ndani ya data set yao ya majaribio.

Uchambuzi huu haujaribu generalizability ya mbinu kwa mashine au sehemu nyingine. Statistical significance inaunga mkono kwamba katika experimental setup hii, error distributions kabla na baada ya fidia ni tofauti.

Ni reference gani ilitumika katika statistical comparison?

Chanzo kinatumia hapa maelezo yanayohitaji tahadhari. Sehemu ya kwanza ya Section 4.2 inafafanua 0,229 mm kama average deviation kati ya compensated real trajectory na compensated simulation trajectory. Lakini sehemu ya takwimu inayofuata inasema kwamba both uncompensated and compensated measurements zili-projected kwenye “original finite element simulated trajectory”.

Chanzo hakielezi kwa kina jinsi definitions hizi mbili za evaluation zinavyolinganishwa. Kwa hiyo, 0,229 mm na p-values zimehifadhiwa kwa definitions zake za chanzo bila kufanya reinterpretation ya kimya ya reference trajectory.

Nguvu ya GPR model ni nini?

Mbinu inayopendekezwa hailazimishi polynomial form maalum mapema kwa error. GPR inaweza kutoa kwa pamoja mean value ya systematic normal error na prediction uncertainty kulingana na trajectory position.

Katika Figure 22, red points ni training samples, blue curve ni GPR prediction mean, na shaded area ni prediction band iliyoandikwa kama %95 confidence interval. Hii inaonyesha wazi kwamba error model si constant offset kwenye trajektoria, bali ni local-varying function.

Kikomo cha GPR model ni nini?

Data kutoka uendeshaji mitatu tu wa mashine bila fidia zilitumika kufundisha modeli. Chanzo hakina independent train/test split au cross-validation kwa tofauti za operating conditions.

Kwamba raw projection samples 370 zote zimetokana na uendeshaji huo huo mitatu ni muhimu: sampuli 370 hazimaanishi majaribio huru 370. Sampuli hizo ni pointi zinazohusiana zilizokusanywa kwenye trajektoria ya mfumo huo huo wa kimwili.

Uendeshaji mpya baada ya fidia unaonyesha practical effect ya mbinu katika mashine ileile na general trajectory condition ileile; generalization kwa part geometry tofauti, rotational speed tofauti, feed ratio tofauti au mashine tofauti haijathibitishwa na data hizi.

Matokeo yanayoungwa mkono na utafiti

  • Trajektoria ya mwendo wa analytical rigid roller iliweza kutolewa kutoka ABAQUS ODB data kupitia discrete shell ndogo iliyounganishwa nayo.
  • Simulation trajectory ilibadilishwa kwenda machine coordinate system na kuendeshwa kwenye mashine halisi ya spinning yenye rola moja.
  • Mwendo halisi wa mashine ulipimwa kwa external displacement sensor na kulinganishwa na simulation trajectory.
  • Average projection errors za uendeshaji mitatu bila fidia zilikuwa 0,458, 0,504 na 0,546 mm; overall average ilikuwa 0,503 mm.
  • Projection method ilitumika kufafanua trajectory error katika local normal direction.
  • GPR model iliundwa kwa normal-direction error samples kutoka majaribio matatu.
  • Kwa kutumia GPR prediction, position-dependent normal-direction compensation ilitumika kwenye simulation trajectory.
  • Katika machine run mpya baada ya fidia, average error ilipimwa kuwa 0,229 mm.
  • Chanzo kinaripoti error reduction ya %54,47.
  • Welch t-test na Mann–Whitney U test ziliona error reduction baada ya fidia kuwa statistically significant ndani ya data set iliyopo.

Matokeo ambayo utafiti hauungi mkono au haujajaribu

  • Utafiti hauonyeshi kwamba dimensional error ya final spun part ilipungua kwa %54,47.
  • Katika sehemu kuu ya trajectory fitting experiment, physical workpiece haikutumika.
  • Wall-thickness distribution, surface roughness au formed-part geometry hazikulinganishwa kabla/baada ya fidia.
  • GPR haijathibitishwa kwa kujitegemea kwenye machine models tofauti.
  • Performance ileile haijaonyeshwa katika rotational speeds tofauti, feed ratios tofauti au trajectory geometries tofauti.
  • Dynamic trajectory errors zinazotokana na machine dynamics, vibration na control response hazikujumuishwa katika main scope ya modeli.
  • Independent test set haikutumika kwa 102 GPR training samples.
  • Utafiti hauonyeshi online real-time adaptive GPR compensation; compensation model huundwa kutoka measured data na kisha huendeshwa tena kama trajektoria mpya.
  • Haijaonyeshwa kwamba 0,229 mm ni accuracy limit inayoweza kupatikana kwa mashine zote za metal spinning.

Mbinu na Matokeo ya Utafiti

Vigezo vya msingi vya kijiometri vya modeli ya koni mbili

KipengeleKigezoThamani
MandrelKipenyo cha mwisho mdogo30 mm
MandrelKipenyo cha mwisho mkubwa160 mm
MandrelNusu-pembe ya koni ya kwanza30°
MandrelNusu-pembe ya koni ya pili45°
MandrelUrefu78,39 mm
MandrelRadius ya arc ya mpito50 mm
BlankKipenyo cha ndani15 mm
BlankKipenyo cha nje150 mm
BlankUnene2,2 mm
RolaKipenyo200 mm
RolaWorking angle iliyotolewa katika sehemu ya modeli30°
RolaCorner radius iliyotolewa katika sehemu ya modeli5 mm

Muhtasari wa kiufundi wa finite element solution

KipengeleTaarifa iliyotolewa katika chanzo
SoftwareABAQUS 2022
Blank element typeC3D8R
Idadi ya elements5850
Idadi ya nodes11.960
RolaAnalytical rigid body
MandrelAnalytical rigid body
Roller–blank friction coefficient0,05
Mandrel–blank friction coefficient0,10
Normal contactHard contact
Blank material6061 aluminium alloy
Spindle speed120 r/min
Feed ratio0,4 mm/r

Kwa kuwa kuna tofauti ya 30°/5 mm dhidi ya 25°/6 mm kati ya source text na Table 3 kuhusu roller angle na corner radius, vigezo hivi viwili havijaunganishwa katika jedwali la kiufundi hapo juu kama final process values moja.

Mnyororo wa data wa trajectory extraction

HatuaMaelezo ya kiufundi
Auxiliary geometric marker1 × 1 × 1 mm discrete shell
Mechanical connectiontie constraint kwa synchronous movement na body
ABAQUS field outputCOORD
Data extractionPython 3.13 / ODB interface
Characteristic pointCorner node ya discrete shell katika mwelekeo wa roller installation
Transformed trajectory401 coordinate points katika Table 5
CNC preparationFirst-order differencing kati ya pointi jirani

Kipimo cha trajektoria halisi ya mashine

KipengeleTaarifa iliyoripotiwa katika chanzo
MashineMashine ya spinning yenye rola moja iliyotengenezwa na Changchun University of Science and Technology
External sensorAbsolute draw-wire displacement sensor
Measurement range1000 mm
Displacement resolution0,098 mm
Linear accuracy±%0,1
Repeatability±%0,01
CommunicationRS485 / Modbus-RTU
Visualization environmentUnity3D 2022.3019f1c1
Graph pluginXChart

Idadi za majaribio kwa projection analysis

Kundi la dataIdadi ya projection pointsHitilafu ya wastani (mm)
Jaribio la 1 lisilo na fidia1180,458
Jaribio la 2 lisilo na fidia1270,504
Jaribio la 3 lisilo na fidia1250,546
Combined uncompensated data3700,503 overall average

118 + 127 + 125 = 370 projection samples ziliunda initial error data set ya GPR model inayofuata.

GPR data preparation na model parameters

Model elementTaarifa iliyoripotiwa katika chanzo
Initial projection samples370
Arc-length bin width1,0 mm
Training samples baada ya binning/averaging102
Train/test splitHaikufanywa
Main kernelMatérn
Matérn smoothness parameterν = 1,5
Additional kernelConstant kernel
Noise modelWhite noise kernel
Hyperparameter optimizationMaximizing log marginal likelihood
Final numerical hyperparameter valuesHazijatolewa katika main text

Matokeo kabla na baada ya fidia

Performance indicatorKabla ya fidiaBaada ya fidia
Average trajectory error0,503 mm0,229 mm
Average error reduction iliyoripotiwa katika chanzo—%54,47
Average compensation magnitude—0,3873 mm katika Figure 24

Ni muhimu kutofautisha compensation amount ya 0,3873 mm na residual error ya 0,229 mm baada ya fidia. Ya kwanza ni average correction magnitude iliyotumika kwenye trajektoria, ya pili ni average geometric deviation iliyobaki baada ya fidia.

Statistical validation

Testp valueαInterpretation ya chanzo
Welch t-test2,442 × 10−40,05Error baada ya fidia ni ndogo
Mann–Whitney U2,040 × 10−50,05Punguzo linaungwa mkono pia katika non-parametric test

Kutumia tests mbili tofauti kunazuia chanzo kutegemea result moja ya takwimu inayotokana tu na dhana ya normal distribution. Hata hivyo, uncompensated group ina runs tatu na compensated group ina run moja; kuwa na pointi nyingi kwenye trajektoria hakuongezi idadi ya independent physical experimental repetitions.

Figures zinaungaje mkono matokeo?

Figure 17 inaonyesha kwamba katika uendeshaji mitatu bila fidia target na actual trajectories zinafanana kwa umbo la jumla, lakini tofauti inaonekana kwenye enlarged conical region. Figure 18 inaonyesha moja ya sababu muhimu za tofauti hii: actual machine trajectory ni stepped wakati simulation curve ni continuous.

Katika Figures 19–21, black simulation curve, green actual machine trajectory na red projection points zinaonyeshwa pamoja. Visualizations hizi zinaeleza kwa nini projection method ilionekana kufaa zaidi kuliko coordinate-component difference.

Figure 22 inaonyesha GPR ikimodeli normal error inayobadilika na prediction uncertainty kwenye trajektoria; Figure 23 inaonyesha compensated trajectory ikisogezwa kuelekea general direction ya measurement sets tatu; Figure 24 inaonyesha compensation magnitude si constant dhidi ya arc length.

Figure 25 inaonyesha ukaribu wa compensated trajectory na newly measured machine trajectory; Figure 26 inaonyesha distribution ya signed na absolute errors baada ya fidia kwa arc length; Figure 27 inaonyesha box plots za uncompensated data groups tatu na compensated data group pamoja na p values za statistical tests mbili.

Muhtasari wa technical uncertainties ndani ya chanzo

MadaTofauti inayoonekana kwenye chanzoMbinu ya Verianla
Simulation roller parametersText 30° / 5 mm; Table 3 25° / 6 mmSeti mbili za thamani zimehifadhiwa tofauti.
Roller-center offsetPrinted equations na numerical P1 offset hazipatani moja kwa moja.Source value haikuhesabiwa upya.
Trajectory coordinate unitFigure 10 “m”; later experimental graphs/tables zinatumia mm.Source-internal unit ambiguity imeachwa wazi.
Statistical referenceSehemu moja inarejea compensated reference, maandishi yanayofuata original FE.Definitions mbili hazijaunganishwa kimya kimya.
Sensor resolution“1024” na 0,098 mm zinatolewa pamoja; unit ya 1024 haijaelezwa.Hakuna conversion kutoka nje ya chanzo.

Mipaka ya kazi za baadaye iliyotajwa na watafiti

Sehemu ya conclusion inaeleza wazi kwamba modeli ya sasa inalenga hasa static trajectory deviations. Dynamic errors zinazotokana na machine dynamics, vibration na control response hazikujumuishwa katika scope ya utafiti huu.

Pia error model imejengwa kwa limited experimental data. Watafiti wanapendekeza kuongeza dynamic error modeling methods na kujaribu generalizability ya mbinu kwa data set pana zaidi chini ya operating conditions tofauti katika kazi za baadaye.

Maelezo ya Chanzo na Mbinu

Jina kamili la utafiti asilia: Roller Trajectory Extraction and Error Compensation Based on Projection Method and Gaussian Process Regression

Waandishi: Wei Liang, Shi-Yuan Kong, Xin Zhao, Pi-Yao Liu.

Mpangilio wa waandishi: Mpangilio wa chanzo umehifadhiwa kama ulivyo.

Mwandishi wa mawasiliano: Shi-Yuan Kong.

Equal contribution/co-first author: Hakuna taarifa kama hiyo kwenye chanzo.

Taasisi: School of Mechanical and Electrical Engineering, Changchun University of Science and Technology, Changchun, China; Zhongshan Institute, Changchun University of Science and Technology, Zhongshan, China.

Jarida: Machines.

Mchapishaji: MDPI.

Bibliographic record: Machines 2026, 14(6), 676.

DOI: 10.3390/machines14060676.

Kiungo rasmi cha chapisho:https://doi.org/10.3390/machines14060676

Tarehe ya uchapishaji: 10 Juni 2026.

Aina ya chanzo na hali ya peer review: Ni original peer-reviewed research article. Utafiti unachanganya finite element simulation, physical machine trajectory measurement na data-driven error compensation.

Leseni: Creative Commons Attribution (CC BY).

Ufadhili: Utafiti uliungwa mkono chini ya introduced innovative research team project ya Zhongshan Research Institute of Changchun University of Science and Technology, kwa project number CXTD2023006.

Upatikanaji wa data: Waandishi wamesema raw data zinazounga mkono matokeo ya utafiti zitatolewa kwa ombi. Pia transformed trajectory coordinates, processed coordinates, projection data sets tatu, compensated trajectory na measured coordinates baada ya fidia zimeorodheshwa katika supplementary material kama Tables S1–S7.

Conflict of interest: Waandishi wametangaza kuwa hakuna conflict of interest.

Author contributions: Wei Liang alishiriki katika conceptualization, methodology, resources, review/editing, supervision na funding acquisition; Shi-Yuan Kong katika software na validation; Xin Zhao katika original draft writing; Pi-Yao Liu katika data curation. Chanzo kinasema waandishi wote walisoma na kuidhinisha published version.

Finite element model, geometries, material na contact values, trajectory extraction method, sensor system, projection equations, GPR structure, training data counts, compensation values, experimental results na statistical analyses katika makala hii ya Verianla zinategemea utafiti uliokaguliwa. Hakuna experimental finding mpya au performance result mpya iliyoongezwa kutoka vyanzo vya nje zaidi ya bibliographic identity verification.

Basic experimental limit: Chanzo kinaeleza wazi kwamba workpiece haikutumika katika roller trajectory fitting experiment. Jaribio la baada ya fidia pia linathibitisha machine trajectory motion; chanzo hakitoi before/after compensation comparison ya dimensional accuracy, wall thickness au surface quality ya final double-cone part. Kwa hiyo, %54,47 error reduction ni ya roller trajectory matching na haiwezi kutafsiriwa kama part-quality improvement kwa kiwango kilekile.

GPR validation limit: Sampuli zote 370 za projection kutoka majaribio matatu yasiyo na fidia zilitumika katika model development. Baada ya 1,0 mm arc-length binning, GPR ilifundishwa kwa sampuli 102 na independent random train/test split haikufanywa. Kuendesha tena compensated trajectory kwenye mashine kunathibitisha athari ya mbinu katika mfumo huo huo, lakini si ushahidi wa generalizability kwa mashine tofauti na process conditions tofauti.

Source-internal process parameter note: Roller installation angle na corner radius zimetolewa kama 30° na 5 mm katika explanatory text, lakini 25° na 6 mm katika Table 3. Chanzo hakielezi kwa maana moja ni pair ipi ya values iliyotumika kwenye final simulation.

Source-internal coordinate note: Katika discrete-shell–roller-center transformation, si components zote za numerical offset \(P_1=(x_A+13{,}433,\ y_A-22{,}267,\ z_A+0{,}5)\) zilizoandikwa baadaye kwenye chanzo zinapatana moja kwa moja na Equations (2)–(7). Pia Figure 10 inaandika machine coordinate axes kwa meters, wakati experimental trajectory data zinazofuata zimewasilishwa kwa scale ya mm. Taarifa hizi zimehifadhiwa kama zilivyo katika chanzo.

Statistical reference note: Maelezo ya kwanza ya Section 4.2 yanafafanua 0,229 mm kama deviation kati ya compensated real na compensated simulation trajectory, wakati maelezo ya takwimu yanayofuata yanasema data sets zote, pamoja na compensated measurement, zili-projected kwenye original finite element trajectory. Chanzo hakipatanishi definitions hizi mbili za reference.

Model scope: Watafiti wanaeleza kwamba utafiti huu unalenga static trajectory deviations, haujachunguza dynamic errors kutoka machine dynamics, vibration na control response, na unahitaji additional validation kwa experimental data sets pana zaidi.


Shiriki:

Maoni huchapishwa baada ya kukaguliwa.Maoni yako yatapitia mchakato wa idhini na yataonekana yakikubaliwa.

Acha maoni

Anwani yako ya barua pepe haitachapishwa. Sehemu za lazima zimewekewa alama ya *

Your experience on this site will be improved by allowing cookies Cookie Policy