Oʻpka toʻqimasini virtual jarrohlik muhitida real vaqt rejimida deformatsiyalash ikki xil masshtabni bir vaqtning oʻzida yechishni talab qiladi. Jarrohlik asbobi tegadigan kichik hududlarda nozik, yuqori chastotali mahalliy deformatsiyalar yuzaga keladi, oʻpkaning katta qismining siqilishi, choʻzilishi yoki nafas olishga oʻxshash harakatlari esa past chastotali global deformatsiyani hosil qiladi. Extended Position-Based Dynamics (XPBD) yuqori sonli barqarorlik va past hisoblash xarajati tufayli interaktiv yumshoq toʻqima simulyatsiyalarida qulay; biroq klassik cheklovdan-cheklovga mahalliy yechim tuzilmasi katta masshtabdagi xato haqidagi axborotning model boʻylab tarqalishini sekinlashtiradi.

Ushbu tadqiqot XPBD-ni smoothed aggregation algebraic multigrid (SA-AMG) bilan birlashtiradigan koʻp masshtabli yechuvchini taklif qiladi. SA-AMG ochiq geometrik dagʻal toʻrlarni yaratish oʻrniga XPBD-ning algebraik tizim matritsasidagi bogʻlanishlardan avtomatik ravishda yechim iyerarxiyasini tuzadi. Nozik darajalar mahalliy yuqori chastotali xatolarni bartaraf etsa, dagʻal darajalar oʻpkaning umumiy shakli bilan bogʻliq past chastotali xatolarni qisqaroq axborot yoʻllari orqali tuzatadi.

Ikkinchi yangilik — adaptive spectrum switching (ASS) — yechim jarayonida qoldiq xatoning qanchalik yaxshi kamayishini oʻlchab, mahalliy Jacobi silliqlashi bilan global multigrid tuzatishi oʻrtasida dinamik tanlov qiladi. Shu tariqa tizim har bir iteratsiyada zaruratsiz bir xil yechim strategiyasini qoʻllash oʻrniga ayni paytda ustun boʻlgan xato xususiyatiga qarab harakat qiladi.

Haqiqiy bemor CT tasvirlaridan yaratilgan tetraedral oʻpka modellari ustida oʻtkazilgan sinovlarda taklif etilgan usul klassik XPBD-ga nisbatan global deformatsiyalarni NFEM etaloniga yaqinroq hosil qilgan; kuch–siljish sinovida yuqori aniqlikdagi EKF-ML etaloniga nisbatan oʻrtacha taxminan 0,05 N kuch farqini koʻrsatgan. Sinovdan oʻtkazilgan apparatda taxminan 22.000 tetraedral elementgacha ishning 30 Hz real-vaqt mezoniga yaqin ishlashi saqlangan. Bundan tashqari, 15 jarroh ishtirokidagi vazifaga asoslangan baholashda barcha foydalanuvchi baholash oʻlchovlarining oʻrtachasi 85 balldan yuqori boʻlgan.

Oʻpka simulyatsiyasidagi muammo faqat “yumshoqlik” emas

Oʻpka toʻqimasi bitta masshtabda deformatsiyalanmaydi. Jarrohlik forsepsi ushlab turgan bir necha millimetrlik hudud juda tez va mahalliy ravishda shaklini oʻzgartirar ekan, ayni kuchning mexanik ta’siri uzoqroq toʻqima hududlariga ham tarqaladi. Nafas olish, keng hududli siqilish yoki katta deformatsiyalarda bu global javob ayniqsa muhimdir.

Klassik XPBD elastik energiyani cheklovlar shaklida ifodalab, har bir cheklovni mahalliy ravishda yechadi. Bu yondashuv yaqin qoʻshnilar orasidagi yuqori chastotali xatolarni tez kamaytiradi; biroq past chastotali, keng hududli xatolarning model boʻylab uzatilishi uchun koʻp iteratsiya talab etiladi.

Klassik XPBD Nega Oʻpkaning Global Deformatsiyasida Sekin Yaqinlashadi?

Chunki klassik XPBD global tenglamalar tizimining barcha bogʻlanishlarini toʻgʻridan-toʻgʻri yechish oʻrniga, asosan faqat mahalliy cheklov yangilanishlarini bajaradi. Bu jarayon qoʻshni tugunlar orasidagi kichik masshtabli xatolarni tez kamaytirar ekan, butun oʻpkaga tarqaladigan past chastotali xato komponentlari iteratsiyadan iteratsiyaga sekin uzatiladi. Natijada mahalliy kontakt hududi barqaror koʻrinsa ham, modelning uzoq qismlari toʻgʻri mexanik javobga kech yetishi mumkin.

XPBD-ning asosiy tenglamasi

Yumshoq toʻqima harakati tadqiqotda Nyuton tenglamasi bilan boshlanadi:

\[ \mathbf{M}\ddot{\mathbf{x}} = -\nabla U(\mathbf{x})^{T} + \mathbf{f}_{ext} \]

Bu yerda M massa matritsasi, x zarracha pozitsiyalari, U(x) elastik potensial energiya va f_ext tashqi kuchlardir.

XPBD-da potensial energiya cheklov funksiyalari bilan:

\[ U(\mathbf{x}) = \frac{1}{2} \mathbf{C}(\mathbf{x})^{T} \boldsymbol{\alpha}^{-1} \mathbf{C}(\mathbf{x}) \]

koʻrinishida ifodalanadi. Lagrange koʻpaytiruvchilari aniqlangach va Newton-Raphson chiziqli yaqinlashuvi qoʻllangach, yechilishi kerak boʻlgan asosiy Schur-komplement tizimi:

\[ \mathbf{A}\Delta\boldsymbol{\lambda} = \mathbf{b} \]

koʻrinishiga keltiriladi.

Klassik XPBD toʻliq A matritsasini yechish oʻrniga, asosan matritsaning diagonal tuzilmasidan foydalanadigan cheklovga asoslangan taxminiy yangilanishlarni bajaradi. Buning afzalligi tezlik va barqarorlik, kamchiligi esa uzoq erkinlik darajalari oʻrtasida mexanik axborotning sekin tarqalishidir.

Toʻliq global yechim nega ishlatilmaydi?

Manbaga koʻra, toʻliq tizimni bevosita yechish ham mahalliy, ham global xatolarni tez kamaytirishi mumkin; biroq katta tizimlarda taxminan O(m³) xarajat tufayli real vaqt jarrohlik ilovalari uchun amaliy emas. Tadqiqotning koʻp masshtabli yondashuvi shu ikki chegara oʻrtasida samaraliroq yechim yoʻlini yaratishga intiladi.

SA-AMG va ASS Birgalikda Nima Qiladi?

SA-AMG XPBD tizim matritsasidagi kuchli algebraik bogʻlanishlardan dagʻal yechim darajalarini hosil qilib, global past chastotali xatolarni shu darajalarda tuzatadi; ASS esa qoldiq xatoning spektral xulqini kuzatib, ayni iteratsiyada mahalliy silliqlash yoki multigrid tuzatishning qaysi biri samaraliroq boʻlishini tanlaydi. Shu tariqa mahalliy va global deformatsiya bir xil yechim siklida adaptiv tarzda koʻrib chiqiladi.

Geometrik emas, algebraik iyerarxiya

An’anaviy geometric multigrid usullarida dagʻal darajalar uchun haqiqiy geometriyaning soddalashtirilgan yangi toʻrlari yaratiladi. Murakkab oʻpka geometriyalarida va katta deformatsiyalarda XPBD cheklov gradientlarini bu toʻrlar orasida izchil uzatish qiyinlashishi mumkin.

SA-AMG esa tizim matritsasi elementlaridan foydalanib, ikki erkinlik darajasining qanchalik kuchli bogʻlanganini hisoblaydi:

\[ S_{ij} = \begin{cases} 1, & i \neq j \text{ ve } |A_{ij}| \geq \theta_{SoC}\sqrt{A_{ii}A_{jj}} \\ 0, & \text{aksi durumda} \end{cases} \]

Kuchli bogʻlangan cheklovlar parallel greedy algoritmi bilan guruhlanadi va har bir guruh dagʻalroq darajada bitta makro-erkinlik darajasiga aylanadi.

Smoothed aggregation nega kerak?

Oddiy aggregation usulida dagʻal darajadan nozik darajaga uzatish boʻlakma-boʻlak doimiy boʻlishi mumkin. Tadqiqot bu holat murakkab anatomik chegaralar va geterogen material tuzilmalarida keskin oʻtishlarga sabab boʻlishi mumkinligini bildiradi.

Shu sababli tentative prolongation matritsasi weighted Jacobi bilan silliqlanadi:

\[ \mathbf{P} = \left( \mathbf{I} - \omega \mathbf{D}^{-1}\mathbf{A} \right) \mathbf{P}_{tent} \]

va simmetrik tizim uchun:

\[ \mathbf{R}=\mathbf{P}^{T} \]

deb tanlanadi.

Dagʻal tizim Galerkin proyeksiyasi bilan yaratiladi:

\[ \mathbf{A}_{l+1} = \mathbf{P}_{l}^{T} \mathbf{A}_{l} \mathbf{P}_{l} \]

Shunday qilib, nozik darajadagi past chastotali deformatsiya xatolarini oʻlchami kichikroq dagʻal tizimda osonroq yechish mumkin boʻladi.

V-cycle yechim oqimi

Har bir darajadagi asosiy jarayon quyidagicha:

1. nozik darajada pre-smoothing,
2. qoldiq xatoni hisoblash,
3. qoldiq xatoni restriction orqali dagʻal darajaga uzatish,
4. dagʻal xato tenglamasini yechish,
5. xato tuzatishini prolongation orqali nozik darajaga olib oʻtish,
6. post-smoothing.

ASS-ning qaror oʻzgaruvchisi: residual smoothness factor

Tadqiqot bitta mahalliy smoothing qadami qoldiq xatoni qanchalik kamaytirishini:

\[ \eta = \frac{\|\mathbf{r}_{after}\|_2} {\|\mathbf{r}_{before}\|_2} \]

nisbati bilan oʻlchaydi.

Agar:

\[ \eta < \theta_{smooth} \]

boʻlsa, qoldiq xatoda yuqori chastotali komponentlar ustun va mahalliy weighted-Jacobi yangilanishi samarali deb hisoblanadi.

Aksincha:

\[ \eta \geq \theta_{smooth} \]

boʻlganda mahalliy smoothing qoldiq xatoni yetarlicha kamaytira olmaydi; tizim SA-AMG dagʻal-grid tuzatishiga oʻtadi.

Prolongation parametri ham oʻzgarmas qoldirilmaydi

Tadqiqot har bir iyerarxiya tuzilishida oldindan shartlangan matritsaning spektral radiusini bir necha power iteration bilan taxminan baholab:

\[ \omega_{opt} \approx \frac{2} {\rho(\mathbf{D}^{-1}\mathbf{A})} \]

bogʻlanishi orqali weighted prolongation parametrini yangilaydi.

Taklif Etilgan Usul Real Vaqt Virtual Jarrohlik Uchun Qanchalik Tez va Aniq?

Manbaning oʻz tajribalarida usul taxminan 22.000 tetraedral elementgacha 30 Hz real-vaqt mezoni atrofida qolgan, klassik va takomillashtirilgan ayrim fizikaga asoslangan muqobillarga nisbatan yuqoriroq yangilanish tezligini koʻrsatgan hamda EKF-ML kuch–siljish etaloniga nisbatan oʻrtacha taxminan 0,05 N kuch farqini hosil qilgan. Biroq mutlaq aniqlik yuqori aniqlikdagi FEM darajasiga toʻliq yetmaydi va bu natijalar haqiqiy inson toʻqimasidan bevosita oʻlchangan klinik ground truth ga qarshi validatsiya emas.

Real vaqt sinovi

Tadqiqot besh yechim yondashuvini taqqoslagan:

• Transformer-KAN,
• 4D-XPBD,
• taklif etilgan koʻp masshtabli XPBD,
• KF-ADMM,
• WSC-integrated.

Oʻpka modellari 5.000 dan 25.000 gacha tetraedral elementni oʻz ichiga oladi va model hajmi 1.000 tadan element bilan oshiriladi.

Mualliflar Transformer-KAN ning tezlik ustunligiga qarshi oʻz usullarining oʻquv ma’lumotlarini talab qilmaydigan, toʻliq fizikaga asoslangan tuzilishini asosiy afzallik deb baholaydilar.

Kuch–siljish aniqligi

Aniqlik tajribasida ayni oʻpka modeliga 0–2 N oraligʻida ortib boruvchi tashqi kuch qoʻllanib, tanlangan tugunlarning siljishi oʻlchangan.

Ushbu taqqoslashda EKF-ML hisoblash etaloni sifatida tanlangan. Shu sababli natija haqiqiy tirik yoki ex vivo oʻpka toʻqimasining mutlaq oʻlchash xatosi emas.

100 kadrli barqarorlik tajribasi

Manbaning 12-rasmida klassik XPBD va taklif etilgan yechuvchining har bir kadr uchun ishlov berish vaqtlari taqqoslanadi. Taklif etilgan usulda eng yuqori va eng past ishlov berish vaqti orasidagi umumiy diapazon 10 ms dan kam boʻlib qoladi, klassik XPBD da esa ketma-ket kadrlar orasidagi sakrash ayrim nuqtalarda 14 ms gacha chiqadi.

Tadqiqot bu vaqt barqarorligini ASS ning mahalliy va global xato tuzatishlari oʻrtasidagi adaptiv oʻtishi bilan bogʻlaydi.

Jarroh baholashining dizayni

15 jarroh PHANTOM Omni haptik boshqaruv qurilmasi bilan virtual forsepsdan foydalangan. Vazifa oʻpka toʻqimasini A boshlangʻich nuqtadan B maqsad hududiga olib borishda C xavf hududiga tegmaslikdan iborat.

Oltita baholash oʻlchovi qoʻllangan:

• deformatsiya realizmi,
• oʻzaro ta’sir izchilligi,
• idrokka oid qayta aloqa integratsiyasi,
• fizik izchillik,
• tizimning javob berish qobiliyati,
• barqarorlik va chidamlilik.

Ballar 100 dan baholangan va toʻrt ishlash darajasi aniqlangan: Poor (<60), Fair (60–75), Good (75–90) va Excellent (90–100).

Barcha oʻlchovlarning oʻrtacha bali 85 dan yuqori, fizik izchillikning oʻrtachasi esa 92,60 boʻlgan. Manba katta tajribaga ega jarrohlar biroz pastroq ball berishga moyil boʻlganini, ammo umumiy baho barcha tajriba guruhlarida ijobiy qolganini bildiradi.

Kasallik holatlariga moslashtirish

Manba 14-rasmda toʻrtta patologik holatni koʻrsatadi:

• oʻpka fibrozi,
• oʻpka saratoni,
• sil,
• emfizema.

Tadqiqotda asosiy fizik parametrlar oʻzgartirilib, masalan, saraton tugunida mahalliy qattiqlashish yoki fibrotik toʻqimada oʻzgargan elastiklik kabi xatti-harakatlar simulyatsiya qilingani bildiriladi. Biroq bular har bir kasallik uchun mustaqil mexanik validatsiya tadqiqotlari emas.

Tadqiqot Usuli va Natijalari

Dasturiy ta’minot va apparat

Neo-Hookean modeli

Oʻpkaning giperelastik javobi:

\[ \Psi_{Neo}(\mathbf{F}) = \frac{\lambda_{lame}}{2} (\det(\mathbf{F})-1)^2 + \frac{\mu_{lame}}{2} \left( \operatorname{tr}(\mathbf{F}^{T}\mathbf{F})-3 \right) \]

energiya zichligi bilan modellashtirilgan.

Qoʻllangan material parametrlari:

• E = 1,96 kPa,
• ν = 0,49.

Energiya keyin tetraedral elementlar uchun hajmiy va deviatorik XPBD cheklovlariga aylantirilgan:

\[ C^{H}(\mathbf{x}) = \det(\mathbf{F})-1 \]

\[ C^{D}(\mathbf{x}) = \sqrt{ \operatorname{tr}(\mathbf{F}^{T}\mathbf{F}) } \]

Mahalliy deformatsiya validatsiyasi

Ayni tugunlarga ayni mahalliy siqish kuchi qoʻllanib NFEM, taklif etilgan usul va klassik XPBD taqqoslangan. Manbaning 8-rasmida taklif etilgan yechuvchining deformatsiya hududi NFEM etaloniga vizual jihatdan juda yaqin ekani, klassik XPBD esa torroq deformatsiya hududi va qattiqroq javob hosil qilishi koʻrsatiladi.

Global deformatsiya validatsiyasi

Manbaning 9-rasmida uch vaqt qadamida global siqish koʻrsatiladi. Taklif etilgan usul modelning umumiy morfologiyasini NFEM ga yaqin saqlagan holda, klassik XPBD da ayniqsa pastki chegaralar va model chetlarida yaqqol ajralishlar yuzaga keladi.

Tadqiqot talqiniga koʻra, buning sababi klassik XPBD ning past chastotali global xatolarni yetarlicha tez kamaytira olmasligidir.

SA-AMG ning hisoblashdagi roli

Nozik daraja:

• mahalliy deformatsiya,
• yuqori chastotali xato,
• kichik masshtabli kontakt javoblari.

Dagʻal darajalar:

• global siqilish va kengayish,
• past chastotali xato,
• uzoq tugunlar orasida kuch uzatilishi.

ASS ning hisoblashdagi roli

ASS barcha vaqt qadamlarida oʻzgarmas sondagi V-cycle ni bajarish oʻrniga residual smoothness factor qiymatiga qaraydi. Mahalliy smoothing residualni tez kamaytirsa, mahalliy yechim davom etadi; kamayish sekinlashsa, global coarse-grid correction ishga tushadi.

Usulning tadqiqot bevosita qoʻllab-quvvatlaydigan afzalliklari

• Klassik XPBD ga nisbatan global past chastotali deformatsiya xatolarini tezroq kamaytiradi.
• NFEM etaloniga mosroq global va mahalliy deformatsiya shakllarini hosil qiladi.
• N-XPBD va N-MSM ga qaraganda kichikroq kuch farqini koʻrsatadi.
• Sinovdan oʻtkazilgan apparat va mesh oʻlchamlarida real-vaqt interaktivligiga mos yangilanish tezliklarini ta’minlaydi.
• Har bir kadr uchun hisoblash vaqtlaridagi oʻzgaruvchanlikni klassik XPBD ga nisbatan kamaytiradi.
• Geometrik dagʻal toʻr yaratish zaruratisiz algebraik multilevel tuzilma quradi.
• Ma’lumot bilan oʻqitishni talab qilmaydigan fizikaga asoslangan yechim doirasini saqlaydi.

Tadqiqot bevosita qoʻllab-quvvatlamaydigan natijalar

• Usul haqiqiy operatsiyada bemor natijalarini yaxshilashi koʻrsatilmagan.
• 0,05 N farq haqiqiy inson oʻpka toʻqimasiga nisbatan oʻlchangan mutlaq klinik xato emas.
• 15 jarroh ishtirokidagi foydalanish qulayligi tadqiqoti klinik samaradorlik yoki bemor xavfsizligi tadqiqoti emas.
• Haqiqiy CT tasviridan foydalanishning oʻzi bemorga xos biomexanik validatsiya degani emas.
• Fibroz, saraton, sil va emfizema tasvirlari bu kasalliklar uchun keng qamrovli klinik mexanik tekshiruvni tashkil etmaydi.
• Model hozirgi holida kesish, yirtilish va tikish kabi dinamik topologiya oʻzgarishlarini yechmaydi.
• Yuz minglab tetraedral elementli modellarda xuddi shunday real-vaqt ishlashi koʻrsatilmagan.
• XPBD ning mutlaq mexanik aniqligi yuqori aniqlikdagi FEM usullarini butunlay almashtira oladigan darajaga yetgani koʻrsatilmagan.

Masshtablanuvchanlik chegarasi

Manba usul oʻrta masshtabli modellarda yaxshi ishlashini, biroq yuz minglab tetraedral elementni oʻz ichiga olgan juda yuqori aniqlikdagi tuzilmalarda sezilarli ishlash yoʻqotilishi kutilishini tan oladi.

Kelajakdagi ish sifatida dual-space AMG dan foydalanish taklif qilinadi.

Dinamik topologiya chegarasi

Virtual jarrohlikning muhim amallari boʻlgan kesish, yirtilish va tikish geometriyaning bogʻlanish tuzilmasini doimiy ravishda oʻzgartiradi. Mavjud SA-AMG iyerarxiyasini bu oʻzgarishlarga real vaqtda qayta qurish tadqiqotda hal qilinmagan.

Mualliflar kelajakda real-time hierarchy update va dynamic error correction usullarini oʻrganishni rejalashtirganliklarini bildiradilar.

Manba va Usul Haqida Izoh

Asl sarlavha:A multiscale XPBD solver with residual-adaptive algebraic multigrid for lung deformation simulation

Mualliflar: Xiaorui Zhang; Shuaitong Li; Wei Sun; Jiayi Chang; Qifeng Zhou; Sunil Jha.

Mas’ul muallif: Shuaitong Li.

Muassasalar: College of Computer and Information Engineering, Nanjing Tech University; School of Artificial Intelligence, Beijing Institute of Technology; School of Computer Science and Software va School of Automation, Nanjing University of Information Science and Technology; Adani University.

Manba turi: Hisoblash biomexanikasi / virtual jarrohlik tadqiqot preprinti.

Platforma: SSRN.

SSRN qayd raqami: 7367058.

Taqriz holati: Hujjat ochiq ravishda taqrizdan oʻtmagan va asosiy matn “Anonymized for review: Preprint submitted to Elsevier” tarzida belgilangan.

DOI: Yuklangan hujjatda DOI berilmagan.

Ma’lumotlarga kirish: “Data will be made available on request.”

Moliyalashtirish: National Natural Science Foundation of China, Grants 62272236, 62376128 va 62306139; Natural Science Foundation of Jiangsu Province, Grants BK20201136 va BK20191401.

Manfaatlar toʻqnashuvi: Mualliflar ma’lum moliyaviy manfaat yoki shaxsiy munosabat mavjud emasligini bildirgan.

Muallif hissalari: Xiaorui Zhang moliyalashtirish, vizualizatsiya, validatsiya, nazorat, resurslar va ma’lumotlarni boshqarish; Shuaitong Li usul, tadqiqot, formal tahlil, konseptuallashtirish va yozish; Wei Sun nazorat, resurslar va konseptuallashtirish; Jiayi Chang va Qifeng Zhou dasturiy ta’minot va tadqiqot; Sunil Jha yozish va nazoratga hissa qoʻshgan.

Klinik ma’lumot izohi: Tadqiqotda haqiqiy bemor CT tasvirlaridan foydalanilgani va 15 professional jarroh foydalanuvchi baholashida qatnashgani bildiriladi. Biroq yuklangan matnda CT holatlari soni, bemor tanlovi, anonimlashtirish jarayoni, etika qoʻmitasi/rozilik raqami yoki jarrohlar ishtirokiga oid etik jarayon tafsilotlari aniq koʻrsatilmagan. Shu sababli maqola bu unsurlar boʻyicha manbada boʻlmagan ma’lumotni qoʻshmaydi.

Validatsiya chegarasi: NFEM va EKF-ML tadqiqotda ground-truth/reference sifatida ishlatiladi; ular bevosita ex vivo yoki in vivo biomexanik oʻlchov emas. Shu sababli natijalar hisoblash asosidagi taqqoslama validatsiyadir.

Material modeli chegarasi: Asosiy baholashda Neo-Hookean model, E = 1,96 kPa va ν = 0,49 qoʻllanadi. Haqiqiy oʻpka toʻqimasining bemor, kasallik, hudud, yoʻnalish va yuklash tezligiga bogʻliq murakkabroq xatti-harakatlari ushbu yagona parametrlar toʻplami bilan toʻliq ifodalangan deb qabul qilinmasligi kerak.

Ishlash chegarasi: Tajribalar Intel i7-10700 va RTX 2080 Ti apparatida bajarilgan. Bildirilgan Hz qiymatlarini boshqa apparatga bevosita koʻchirish mumkin emas.

Ekspert baholashi chegarasi: 15 jarrohning ballari tadqiqotning oʻz vazifa va baholash doirasiga asoslanadi; bu mustaqil koʻp markazli foydalanuvchi tadqiqoti emas.

Kelajakdagi ish: Dual-space AMG, yuz minglab tetraedral elementli modellar, real-time hierarchy reconstruction va cutting/tearing/suturing kabi dinamik topologiya oʻzgarishlari.