Orratzaren alaka geometriak bihurgune anplitudeari eragiten dio ultrasoinu bidez anplifikatutako orratz fineko biopsian

Eskerrik asko Nature.com bisitatzeagatik.Erabiltzen ari zaren arakatzailearen bertsioak CSS laguntza mugatua du.Esperientzia onena lortzeko, eguneratutako arakatzailea erabiltzea gomendatzen dugu (edo Internet Explorer-en bateragarritasun modua desgaitzea).Bitartean, laguntza etengabea bermatzeko, gunea estilorik eta JavaScript gabe errendatuko dugu.
Duela gutxi frogatu da ultrasoinuen erabilerak ehunen errendimendua areagotzen duela ultrasoinu bidez lagundutako orratz finaren aspirazioan (USeFNAB) ohiko orratz finaren aspirazioarekin (FNAB) alderatuta.Orain arte, alakaren geometriaren eta punta-mugimenduaren arteko erlazioa ez da ondo aztertu.Azterketa honetan, orratzaren erresonantziaren eta deflexioaren anplitudearen propietateak ikertu ditugu orratz-alakaren geometria ezberdinetarako.3,9 mm-ko lantza alakadun ohiko bat erabiliz, punta desbideratzeko potentzia-faktorea (DPR) airean eta uretan 220 eta 105 µm/W-koa zen, hurrenez hurren.4 mm-ko punta alakatua baino altuagoa da, 180 eta 80 µm/W DPR airean eta uretan, hurrenez hurren.Azterketa honek alakaren geometriaren tolestura-zurruntasunaren arteko erlazioaren garrantzia azpimarratzen du txertatzeko bide ezberdinen testuinguruan, eta, beraz, zulatze osteko ebaketa-ekintza kontrolatzeko metodoak ezagut ditzake, orratzaren alakaren geometria aldatuz, eta hori garrantzitsua da.USeFNAB aplikazio baterako ezinbestekoa da.
Orratz finko aspirazioko biopsia (FNA) ustezko patologiarako ehun-laginak lortzeko metodo bat da1,2,3 orratz baten bidez.Franseen puntak ohiko lancet4 eta Menghini5 aholkuek baino diagnostiko errendimendu handiagoa ematen duela frogatu da.Aldapa axisimetrikoak (hau da, zirkunferentziazkoak) ere iradokitzen dira ale histopatologikoki egokiak izateko probabilitatea areagotzeko.
Biopsia batean, orratz bat larruazaleko eta ehuneko geruzetatik pasatzen da lesio susmagarrietara sartzeko.Azken ikerketek frogatu dute ultrasoinuek ehun bigunetara sartzeko behar den sartze-indarra murrizten dutela7,8,9,10.Orratz-alakaren geometriak orratzaren elkarreragin-indarren eragina duela frogatu da, adibidez, alak luzeagoek ehunen sartze-indar txikiagoak dituztela11.Orratza ehunaren gainazalean sartu ondoren, hau da, zulatu ondoren, orratzaren ebaketa-indarra ehunarekin duen elkarrekintza-indarraren %75 izan daiteke12.Frogatuta dago zulatu osteko fasean, ultrasoinuak (ultrasoinuak) ehun bigunen biopsiaren eraginkortasuna areagotzen duela.Ehun gogorren laginak hartzeko ultrasoinuen bidez hobetutako beste teknika batzuk garatu dira, baina ez da biopsiaren etekina hobetzen duen emaitzarik jakinarazi.Azterketa ugarik ere baieztatu dute desplazamendu mekanikoa areagotzen dela ultrasoinu-estresak jasaten dituztenean16,17,18.Orratz-ehunen elkarreraginetan indar estatiko axialak (longitudinalak) ikerlan asko dauden arren19,20, FNAB ultrasoinuen azpian orratzaren alakaren denborazko dinamikari eta geometriari buruzko azterketa mugatuak daude (USeFNAB).
Azterketa honen helburua alaka geometria ezberdinek orratz-puntaren mugimenduan duten eragina ikertzea izan da, ultrasoinu-tolesturak bultzatutako orratz batean.Bereziki, injekzio-erdiak orratz-punta desbideratzen duen eragina ikertu dugu orratz-alaketan zulatu ondoren (hau da, USeFNAB orratzak hainbat helburutarako, hala nola aspirazio selektiboa edo ehun bigunak eskuratzeko).
Azterketa honetan alaka geometria desberdinak sartu ziren.(a) Lancet zehaztapenak ISO 7864:201636 arauarekin bat egiten du non \(\alpha\) alaka primarioa den, \(\theta\) alaka sekundarioaren biraketa angelua den eta \(\phi\) alaka sekundarioa den. angelua., biratzean, gradutan (\(^\circ\)).(b) Urrats bakarreko txanfla asimetriko linealak (DIN 13097:201937 "estandarra" izenekoak) eta (c) Urrats bakarreko txanfla lineal asimetrikoak (zirkunferentzialak).
Gure ikuspegia alaketan zehar okertze-uhin-luzeraren aldaketa modelatzen hasten da lantza-geometria konbentzionalen, axisimetrikoen eta etapa bakarreko alakaren geometria asimetrikoetarako.Ondoren, azterketa parametriko bat kalkulatu dugu hodiaren maldak eta luzerak transferentziaren jariakortasun mekanikoan duen eragina aztertzeko.Hori beharrezkoa da orratz prototipo bat egiteko luzera optimoa zehazteko.Simulazioan oinarrituta, orratz-prototipoak egin ziren eta haien portaera erresonatzailea esperimentalki ezaugarritu zen, tentsioaren islapen-koefizienteak neurtuz eta airearen, uraren eta %10eko (w/v) gelatina balistikoko potentzia-transferentziaren eraginkortasuna kalkulatuz, eta hortik abiatuta funtzionamendu-maiztasuna zehaztu zen. .Azkenik, abiadura handiko irudiak orratzaren puntan dagoen tolestura-uhinaren desbideratzea zuzenean neurtzeko erabiltzen da airean eta uretan, baita angelu zeihar bakoitzean emandako potentzia elektrikoa eta deflexio-potentzia erlazioaren geometria kalkulatzeko ere ( DPR) injektatutako medioari..
2a irudian ikusten den bezala, erabili 21 galgako hodi bat (0,80 mm OD, 0,49 mm ID, hodiaren hormaren lodiera 0,155 mm, horma estandarra) orratzaren hodiaren luzera (TL) eta alaka angeluarekin (BL) definitzeko ISO-ren arabera. 9626:201621) 316 altzairu herdoilgaitzean (Young-en modulua 205 \(\text {GN/m}^{2}\), dentsitatea 8070 kg/m\(^{3}\) eta Poisson-en ratioa 0,275).
Makurdura-uhin-luzera eta elementu finituen ereduaren (FEM) sintonizazioa orratz eta muga-baldintzetarako.(a) Alakaren luzera (BL) eta hodiaren luzera (TL) zehaztea.(b) Hiru dimentsioko (3D) elementu finituen eredua (FEM) \(\tilde{F}_y\vec {j}\) puntu-indar harmoniko bat erabiliz, orratza proximalki gidatzeko, puntua desbideratzeko eta abiadura neurtzeko. aholkua (\ ( \tilde {u}_y\vec {j}\), \(\tilde{v}_y\vec {j}\)) jariakortasun mekanikoaren transferentzia kalkulatzeko.\(\lambda _y\) \(\tilde{F}_y\vec {j}\) indar bertikalarekiko tolesturako uhin-luzera bezala definitzen da.(c) Grabitate-zentroaren, A sekzio-eremuaren eta x eta y ardatzen inguruan \(I_{xx}\) eta \(I_{yy}\) inertzi momentuen definizioak, hurrenez hurren.
irudian ikusten den bezala.2b,c, A sekzio-eremua duen eta habearen ebakidura-tamaina baino uhin-luzera handiagoa duen habe infinitu (mugagabea) \( c_{EI }\) fase-abiadura okertua (edo okertua) \( c_{EI }\) 22 arabera zehazten da. :
non E Young-en modulua den (\(\text {N/m}^{2}\)), \(\omega _0 = 2\pi f_0\) kitzikapen-maiztasun angeluarra den (rad/s), non \( f_0 \ ) maiztasun lineala da (1/s edo Hz), I interes ardatzaren inguruko eremuaren inertzi momentua\((\text {m}^{4})\), \(m'=\ rho _0 A\ ) luzera unitateko masa da (kg/m), non \(\rho _0\) dentsitatea\((\text {kg/m}^{3})\) eta A gurutzea den habe-eremuaren sekzioa (xy planoa) (\(\ text {m}^{2}\)).Gure adibidean aplikatzen den indarra y-ardatz bertikalarekiko paraleloa denez, hau da, \(\tilde{F}_y\vec {j}\), x ardatz horizontalaren inguruan dagoen inertzia-momentu erregionala bakarrik interesatzen zaigu. hau da, \(I_{xx}\), beraz:
Elementu finituen eredurako (FEM), desplazamendu harmoniko hutsa (m) suposatzen da, beraz, azelerazioa (\(\text {m/s}^{2}\)) \(\partcial ^2 \vec) gisa adierazten da. {u}/ \ partziala t^2 = -\omega ^2\vec {u}\) \(\vec {u}(x, y, z, t): = u_x\vec {i} + u_y\ vec {j } + u_z\vec {k}\) koordenatu espazialetan emandako hiru dimentsioko desplazamendu-bektorea da.Azken honen ordez, COMSOL Multiphysics software paketean (5.4-5.5 bertsioak, COMSOL Inc., Massachusetts, AEB) inplementatzearen arabera, momentuaren oreka-legearen deformazio finituko forma Lagrangiarra honela ematen da:
non \(\vec {\nabla}:= \frac{\partziala}}{\partziala x}\vec {i} + \frac{\partziala}}{\partziala y}\vec {j} + \frac{ \partial }{\partial z}\vec {k}\) tentsore-dibergentzia-eragilea da, \({\underline{\sigma}}\) bigarren Piola-Kirchhoff tentsore tentsorea (bigarren ordena, \(\ text { N/ m}^{2}\)) eta \(\vec {F_V}:= F_{V_x}\vec {i}+ F_{V_y}\vec {j}+ F_{V_z}\vec {k} \) gorputzaren indar-bektorea (\(\text {N/m}^{3}\)) deformatutako bolumen bakoitzeko, eta \(e^{j\phi }\) fase-angeluaren bektorea\(\phi \) (pozik).Gure kasuan, gorputzaren bolumen-indarra nulua da, gure ereduak linealtasun geometrikoa eta deformazio elastiko huts txiki bat hartzen ditu, hau da, non \({\underline{\varepsilon}}^{el}\) eta \({\azpimarratu {\varepsilon}}\) tentsio elastikoa eta tentsio osoa (bigarren ordena, dimentsiorik gabekoa), hurrenez hurren.Hooke-ren elastikotasun isotropo tentsore eratzailea \(\underline{\underline{C}}\) Young-en E modulua erabiliz kalkulatzen da (\(\text {N/m}^{2}\)) eta Poisson-en ratioa v zehazten da, beraz, hau da \(\azpimarra{\underline{C}}:=\azpimarra{\underline{C}}(E,v)\) (laugarren ordena).Beraz, tentsioaren kalkulua \({\underline{\sigma}} := \underline{\underline{C}}:{\underline{\varepsilon}}\ bihurtzen da.
Kalkuluak 10 nodoko elementu tetraedriko bat erabiltzen du, \(\le\) 8 µm-ko elementu-tamaina duena.Orratza hutsean modelatzen da, eta transferitutako mugikortasun mekanikoaren balioa (ms-1 N-1) honela definitzen da: \(|\tilde{Y}_{v_yF_y}|= |\tilde{v}_y\vec { j}|/ |\ tilde{F}_y\vec {j}|\)24, non \(\tilde{v}_y\vec {j}\) esku piezaren irteerako abiadura konplexua den eta \( \tilde {F}_y\ vec {j }\) hodiaren hurbileko muturrean kokatutako indar eragile konplexu bat da, 2b irudian ikusten den bezala.Itzuli jariakortasun mekanikoa dezibeletan (dB) erreferentzia gisa balio maximoa erabiliz, hau da, \(20\log _{10} (|\tilde{Y}|/ |\tilde{Y}_{max}|) \ ) .FEM azterketa guztiak 29,75 kHz-ko maiztasunean egin ziren.
Orratzaren diseinua (3. irudia) 21 galgako orratz hipodermiko konbentzionalez osatuta dago (Kat. zk. 4665643, Sterican\(^\circledR\), kanpoko diametroa 0,8 mm, luzera 120 mm, AISI 304 kromo-nikel herdoilgaitza. altzairua , B. Braun Melsungen AG, Melsungen, Alemania) polipropilenoz egindako plastikozko Luer Lock mahuka batekin hornitua hurbileko muturrean eta behar bezala aldatuta amaieran.Orratzaren hodia uhin-gidara soldatzen da 3b irudian ikusten den moduan.Uhin-gidak altzairu herdoilgaitzezko 3D inprimagailu batean inprimatu ziren (EOS 316L altzairu herdoilgaitza EOS M 290 3D inprimagailu batean, 3D Formtech Oy, Jyväskylä, Finlandia) eta, ondoren, Langevin sentsoreari lotu ziren M4 torlojuen bidez.Langevin sentsorea 8 eraztun piezoelektrikoz osatuta dago, bi muturretatik bi masekin kargatuta.
Lau punta motak (argazkia), komertzialki eskuragarri dagoen lanceta bat (L) eta fabrikatutako fase bakarreko hiru alaka axisimetrikoek (AX1-3) 4, 1,2 eta 0,5 mm-ko alakaren luzera (BL) izan zuten, hurrenez hurren.(a) Amaitutako orratzaren puntaren lehen planoa.(b) 3D inprimatutako uhin-gidari soldatutako lau pinen goiko ikuspegia eta, ondoren, M4 torlojuekin Langevin sentsorera konektatuta.
Hiru alaka-punta axisimetriko (3. irud.) fabrikatu ziren (TAs Machine Tools Oy) 4,0, 1,2 eta 0,5 mm-ko alaka-luzerak (BL, 2a. irudian definitzen den bezala), \(\approx) 2 \(^ \\-ri dagozkionak). circ\), 7\(^\circ\) eta 18\(^\circ\) hurrenez hurren.Uhin-gidaren eta orratzaren masa 3,4 ± 0,017 g da (batez bestekoa ± sd, n = 4) L eta AX1-3 alaketarako, hurrenez hurren (Quintix\(^\circledR\) 224 Design 2, Sartorius AG, Göttingen, Alemania) .3b irudiko L ​​eta AX1-3 alaketarako, orratzaren puntatik plastikozko mahukaren amaierarainoko luzera osoa 13,7, 13,3, 13,3 eta 13,3 cm-koa izan zen, hurrenez hurren.
Orratz-konfigurazio guztietan, orratzaren muturretik uhin-gidaren puntarainoko luzera (hau da, soldadura-eremuraino) 4,3 cm-koa zen, eta orratzaren hodia ebakia gorantz orientatuta zegoen (hau da, Y ardatzarekiko paraleloa). , irudian ikusten den bezala.c (2. irudia).
Ordenagailu batean exekutatzen den MATLAB-en (R2019a, The MathWorks Inc., Massachusetts, AEB) (Latitude 7490, Dell Inc., Texas, AEB) script pertsonalizatu bat erabili zen 25 eta 35 kHz arteko ekorketa sinusoidal lineal bat sortzeko 7 segundoz. pasatzea Digital-analogiko (DA) bihurgailu batek (Analog Discovery 2, Digilent Inc., Washington, AEB) seinale analogikoa bihurtzen du.Seinale analogikoa \(V_0\) (0,5 Vp-p) irrati-maiztasun (RF) anplifikatzaile dedikatu batekin anplifikatu zen gero (Mariachi Oy, Turku, Finlandia).Irteerako inpedantzia duen RF anplifikadorearen \({V_I}\) erortzen den tentsio anplifikatua 50 ohmioko sarrerako inpedantzia duen orratz-egituran eraikitako transformadore batera elikatzen da.Langevin transduktoreak (aurreko eta atzeko geruza anitzeko transduktore piezoelektriko astunak) uhin mekanikoak sortzeko erabiltzen dira.RF anplifikadore pertsonalizatuak kanal bikoitzeko uhin geldikorreko potentzia-faktorea (SWR) neurgailu batekin hornituta dago, \({V_I}\) eta islatutako tentsio anplifikatu\(V_R\) analogikotik digitalerako (AD) moduan erregistratzen dituena.300 kHz-ko laginketa-abiadurarekin (Discovery 2 analogikoa).Kitzikapen-seinalea hasieran eta amaieran anplitudea modulatzen da, anplifikadorearen sarrera iragankorrekin gainkarga ez dadin.
MATLABen inplementatutako script pertsonalizatu bat erabiliz, maiztasun-erantzunaren funtzioa (FRF), hau da, \(\tilde{H}(f)\), lineaz kanpo kalkulatu zen bi kanaleko ekorketa sinusoidala neurtzeko metodoa erabiliz (4. irudia), eta horrek suposatzen du. linealtasuna denboran.sistema aldaezin.Gainera, 20 eta 40 kHz arteko banda-iragazkia aplikatzen da seinaletik nahi ez diren maiztasunak kentzeko.Transmisio-lineen teoriari erreferentzia eginez, kasu honetan \(\tilde{H}(f)\) tentsioaren islapen koefizientearen baliokidea da, hau da, \(\rho _{V} \equiv {V_R}/{V_I}\ ) \) \({V_R}^ 2 /{V_I}^2\ ) \(|\rho _{V}|^2\) berdintzen du.Potentzia elektrikoaren balio absolutuak behar diren kasuetan, intzidente potentzia \(P_I\) eta islatutako potentzia \(P_R\) potentzia (W) kalkulatzen dira dagokion tentsioaren balio eff (rms) hartuz, adibidez.kitzikapen sinusoidala duen transmisio-lerro baterako \( P = {V}^2/(2Z_0)\)26, non \(Z_0\) 50 \(\Omega\) berdina den.\(P_T\) kargari (hots, txertatutako euskarria) hornitutako potentzia elektrikoa \(|P_I – P_R |\) (W RMS) gisa kalkula daiteke, baita potentzia transferentziaren eraginkortasuna (PTE) eta ehunekoa ( % forma nola ematen den zehaztu daiteke, beraz 27:
Ondoren, maiztasun modal azikularrak \(f_{1-3}\) (kHz) eta haiei dagozkien potentzia transferitzeko faktoreak \(\text {PTE}_{1{-}3} \) kalkulatzen dira FRF erabiliz.FWHM (\(\text {FWHM}_{1{-}3}\), Hz) zuzenean estimatua \(\text {PTE}_{1{-}3}\), 1. taulako A alde bakarretik espektro lineala deskribatutako maiztasun modalean lortzen da \(f_{1-3}\).
Orratz-egituren maiztasun-erantzunaren (AFC) neurketa.Bi kanaleko ekorketa-neurketa sinusoidala25,38 erabiltzen da maiztasun-erantzunaren \(\tilde{H}(f)\) eta bere bulkada-erantzuna H(t) lortzeko.\({\mathcal {F}}\) eta \({\mathcal {F}}^{-1}\) mozketa digitalaren Fourier transformatua eta bere alderantzizkoa adierazten dute, hurrenez hurren.\(\tilde{G}(f)\) maiztasun-domeinuko bi seinaleren biderkadura esan nahi du, adibidez, \(\tilde{G}_{XrX}\) alderantzizko eskaneatu produktua\(\tilde{ X} r (f)\ ) eta jaitsiera-tentsioa \(\tilde{X}(f)\) hurrenez hurren.
5. irudian erakusten den moduan, abiadura handiko kamera (Phantom V1612, Vision Research Inc., NJ, AEB) makro-lente batekin hornituta dago (MP-E 65 mm, \(f\)/2.8, 1-5\).(\times\), Canon Inc., Tokio, Japonia), tolestura-kitzikapenean (maiztasun bakarrekoa, sinusoide jarraitua) 27,5-30 kHz-ko maiztasunetan mutur desbideratzeak erregistratzeko.Itzal-mapa bat sortzeko, orratzaren puntaren atzean intentsitate handiko LED zuri baten elementu hoztu bat jarri zen (pieza zenbakia: 4052899910881, LED zuria, 3000 K, 4150 lm, Osram Opto Semiconductors GmbH, Regensburg, Alemania).
Konfigurazio esperimentalaren aurrealdeko ikuspegia.Sakonera bitartekoaren gainazaletik neurtzen da.Orratz-egitura motordun transferentzia-mahai batean finkatu eta muntatzen da.Erabili abiadura handiko kamera bat handitze handiko lente batekin (5\(\x\)) angelu zeiharkako desbideraketa neurtzeko.Dimentsio guztiak milimetrotan daude.
Orratz alaka mota bakoitzerako, 128 \(\x\) 128 pixeleko abiadura handiko kamera baten 300 fotograma grabatu genituen, bakoitza 1/180 mm-ko bereizmen espaziala (\(\gutxi gorabehera) 5 µm), eta 310.000 fotograma segundoko denborazko bereizmena.6. Irudian ikusten den bezala, fotograma bakoitza (1) mozten da (2), orratzaren punta markoaren azken lerroan (behean) dagoela eta irudiaren (3) histograma kalkulatzen da, beraz, Canny 1 eta 2 atalaseak zehaztu daitezke.Ondoren, aplikatu Canny ertzaren detekzioa 28(4) Sobel operadorearekin 3 \(\times\) 3 eta hipotenusa ez diren pixelen posizioak kalkulatu (\(\mathbf {\times }\)) kabitaziorik gabe 300 denbora-pausotan.Puntako desbideratze-barrutia zehazteko, kalkulatu deribatua (diferentzia zentralaren algoritmoa erabiliz) (6) eta zehaztu desbiderazioaren tokiko muturrak (hau da, gailurra) dituen fotograma (7).Kavitaziorik gabeko ertzaren ikuskapen bisual bat egin ondoren, fotograma pare bat (edo erdiko tartea duten bi fotograma) hautatu ziren (7) eta puntaren desbideratzea neurtu zen (\(\mathbf {\times} gisa adierazten da). \) ).Aurrekoa Python-en (v3.8, Python Software Foundation, python.org) inplementatzen da OpenCV Canny ertzak detektatzeko algoritmoa erabiliz (v4.5.1, kode irekiko ordenagailuaren ikusmenaren liburutegia, opencv.org).Azkenik, desbideratze potentzia-faktorea (DPR, µm/W) gailurtik gailurreko desbideratzeak transmititutako potentzia elektrikoaren \(P_T\) (Wrms) proportzioan kalkulatzen da.
7 urratseko algoritmo bat erabiliz (1-7), mozketa (1-2), Canny ertzaren detekzioa (3-4), kalkulua barne, neurtu punta desbideratzeko ertzaren pixelaren posizioa goi-mailatik hartutako fotograma batzuk erabiliz. abiadura kamera 310 kHz-an ( 5) eta bere denbora-deribatua (6), eta, azkenik, puntako desbideratze-eremua ikusmen egiaztatutako fotograma-pareetan neurtzen da (7).
Airean (22,4-22,9 °C), ur deionizatuan (20,8-21,5 °C) eta % 10 (w/v) gelatina balistiko urtsuan (19,7-23,0 °C), \(\text {Honeywell}^{ \ text { TM}}\) \(\text {Fluka}^{\text {TM}}\) Behi- eta txerri-hezur-gelatina I. motako analisi balistikoetarako, Honeywell International, Ipar Carolina, AEB).Tenperatura K motako termopare anplifikadore batekin (AD595, Analog Devices Inc., MA, AEB) eta K motako termopare batekin (Fluke 80PK-1 Bead Probe No. 3648 mota-K, Fluke Corporation, Washington, AEB) neurtu da.Erabili Z ardatz motorizatu bertikal bat (8MT50-100BS1-XYZ, Standa Ltd., Vilnius, Lituania) komunikabideen gainazaletik (Z ardatzaren jatorri gisa ezarrita) sakonera neurtzeko urrats bakoitzeko 5 µm-ko bereizmenarekin.
Laginaren tamaina txikia zenez (n = 5) eta normaltasuna ezin zenez onartu, bi lagin bi isatseko Wilcoxon rank batura proba (R, v4.0.3, R Foundation for Statistical Computing, r-project.org) erabili zen. hainbat alaketarako orratz-puntaren bariantza-kopurua alderatzeko.Aldapa bakoitzerako hiru konparaketa egin ziren, beraz, Bonferroni zuzenketa aplikatu zen 0,017ko esangura-maila egokituarekin eta %5eko errore-tasa batekin.
Beheko 7. irudiari erreferentzia egiten zaio.29,75 kHz-tan, 21 galgako orratz baten uhin-luzera erdi kurbatua (\(\lambda _y/2\)) \(\gutxi gorabehera) 8 mm-koa da.Makurdura-uhin-luzera maldan zehar txikitzen da puntara hurbildu ahala.Puntan \(\lambda _y/2\) 3, 1 eta 7 mm-ko alaka mailakatuak daude, hurrenez hurren, lanceta arruntetarako (a), asimetrikoetarako (b) eta ardatz-simetrikoetarako (c).Hortaz, honek esan nahi du lantza \(\inguru\) 5 mm-ko desberdina izango dela (lantzaren bi planoek 29,30 puntu bat osatzen dutelako), malda asimetrikoa 7 mm aldatuko dela eta malda simetrikoa. 1 mm-koa.Malda aximetrikoak (grabitate-zentroak berdin jarraitzen du, beraz, hormaren lodiera soilik aldatzen da maldan zehar).
FEM azterketa 29,75 kHz eta ekuazioa aplikatzea.(1) Kalkulatu flexio-uhin-erdi aldaketa (\(\lambda _y/2\)) lanceta (a), asimetrikoa (b) eta ardatz-simetrikoa (c) geometria zeiharretarako (1a,b,c irudietan bezala).).Malda lancet, asimetriko eta axisimetrikoen batez besteko \(\lambda_y/2\) 5,65, 5,17 eta 7,52 mm-koa da, hurrenez hurren.Kontuan izan alaka asimetriko eta axisimetrikoetarako punta-lodiera \(\gutxi gorabehera) 50 µm-ra mugatuta dagoela.
Mugikortasun gailurra \(|\tilde{Y}_{v_yF_y}|\) hodiaren luzera optimoa (TL) eta inklinazio luzera (BL) konbinazioa da (8, 9. irudia).Lanceta konbentzional baterako, bere tamaina finkoa denez, TL optimoa \(\gutxi gorabehera\) 29,1 mm-koa da (8. irud.).Malda asimetriko eta axisimetrikoetarako (9a, b irudiak, hurrenez hurren), FEM azterketak 1 eta 7 mm arteko BL barne hartu zuen, beraz, TL tarte optimoak 26,9 eta 28,7 mm (tartea 1,8 mm) eta 27,9 eta 29,2 mm (tartea) izan dira. 1,3 mm).) ), hurrenez hurren.Aldapa asimetrikoetarako (9a. irudia), TL optimoa linealki handitu zen, BL 4 mm-ko lautada batera iritsiz, eta, ondoren, nabarmen jaitsi zen BL 5-tik 7 mm-ra.Malda axisimetrikoetarako (9b. irudia), TL optimoa linealki handitzen da BL luzapenarekin eta azkenean BLn egonkortzen da 6 eta 7 mm-ra.Malda axisimetrikoen azterketa hedatu batek (9c. irudia) \(\gutxi gorabehera) 35,1-37,1 mm-n kokatutako TL optimoen multzo ezberdin bat erakutsi zuen.BL guztietarako, TL optimoen bi multzoen arteko distantzia \(\gutxi gorabehera\) 8 mm-koa da (\(\lambda _y/2\)-ren baliokidea).
Lancet transmisioaren mugikortasuna 29,75 kHz-tan.Orratz-hodia 29,75 kHz-ko maiztasunean flexionatu da, bibrazioa amaieran neurtu eta transmititutako mugikortasun mekanikoaren (dB balio maximoarekiko) TL 26,5-29,5 mm-rako (0,1 mm-ko urratsa) adierazi zen.
29,75 kHz-ko maiztasuneko FEM-aren azterketa parametrikoek erakusten dutenez, punta axisimetrikoaren transferentzia-mugikortasunak hodiaren luzera aldaketek bere pareko asimetrikoak baino eragin txikiagoa dute.Alakaren luzera (BL) eta hodiaren luzera (TL) azterketak alaka geometria asimetrikoak (a) eta axisimetrikoak (b, c) maiztasun-domeinuko azterketetan FEM erabiliz (muga-baldintzak 2. irudian ageri dira).(a, b) TL 26,5 eta 29,5 mm bitartekoa (0,1 mm urratsa) eta BL 1-7 mm (0,5 mm urratsa).(c) Angelu zeihar simetriko hedatua TL 25-40 mm (0,05 mm-ko urratsa) eta 0,1-7 mm-koa (0,1 mm-ko urratsa) barne hartzen duen azterketa \(\lambda_y/2\) \(\lambda_y/2\) mutur baten muga-baldintza solteak betetzen dira.
Orratz-egiturak hiru maiztasun natural \(f_{1-3}\) eskualde modal baxu, ertain eta altuetan banatuta ditu 1. taulan agertzen den moduan. PTEaren tamaina 10. irudian erakusten da eta, ondoren, 11. irudian aztertzen da. Eremu modalitate bakoitzeko emaitzak:
Bat-bateko potentzia-transferentzia-eraginkortasuna (PTE) erregistratutako anplitude tipikoak 20 mm-ko sakoneran frekuentzia sinusoidala erabiliz lortutako lanceta (L) eta AX1-3 malda axisimetrikoak airean, uretan eta gelatinan.Alde bakarreko espektroa erakusten da.Neurtutako maiztasun-erantzuna (300 kHz-ko lagin-abiadura) behe-pasea iragazi zen eta, ondoren, 200 faktore batez lagintu zen analisi modala egiteko.Seinale-zarata erlazioa \(\le\) 45 dB da.PTE fasea (puntudun marra morea) gradutan (\(^{\circ}\)) erakusten da.
Erantzun modalaren analisia 10. irudian (batez bestekoa ± desbideratze estandarra, n = 5) erakusten da L eta AX1-3 maldetarako airean, uretan eta % 10eko gelatina (20 mm-ko sakonera) hiru eskualde modalen (goian) (baxua). , ertaina, altua).), eta haiei dagozkien maiztasun modalak\(f_{1-3}\) (kHz), (batez besteko) energia-eraginkortasuna\(\text {PTE}_{1{-}3 }\) diseinu-ekuazioak erabiltzen ditu.(4) eta (behean) zabalera osoa dira neurtutako balio maximoaren erdian \(\text {FWHM}_{1{-}3}\) (Hz), hurrenez hurren.Kontuan izan PTE baxua grabatzean, hau da, AX2 malda baten kasuan, banda-zabaleraren neurketa baztertzen dela, \(\text {FWHM}_{1}\).\(f_2\) modua da plano inklinatuen deflexioa alderatzeko egokiena dela, potentzia transferentziaren eraginkortasun mailarik handiena erakusten baitu (\(\text {PTE}_{2}\)), arte. %99.
Lehen eskualde modala: \(f_1\) ez dago txertatutako euskarri-motaren araberakoa, baina alaka geometriaren araberakoa da.\(f_1\) txikiagotzen da alakaren luzera txikitzean (27,1, 26,2 eta 25,9 kHz AX1-3rako, hurrenez hurren, airean).Eskualdeko batez bestekoak \(\text {PTE}_{1}\) eta \(\text {FWHM}_{1}\) % 81 eta 230 Hz dira hurrenez hurren.\(\text {FWHM}_{1}\) izan zen Lancet-eko gelatina gehien (L, 473 Hz).Kontuan izan \(\text {FWHM}_{1}\) gelatinako AX2rako ezin dela estimatu jakinarazitako maiztasun-erantzunen magnitude txikia dela eta.
Bigarren eskualde modala: \(f_2\) itsatsi eta alaka euskarri motaren araberakoa da.Airean, uretan eta gelatinan, batez besteko \(f_2\) balioak 29,1, 27,9 eta 28,5 kHz dira, hurrenez hurren.Eskualde modalitate honetako PTE ere %99ra iritsi zen, neurketa talde guztien artean altuena, eskualdeko batez bestekoa %84koa izanik.Eremuaren batez bestekoa \(\text {FWHM}_{2}\) \(\gutxi gorabehera\) 910 Hz-koa da.
Hirugarren eskualde modala: \(f_3\) Maiztasuna txertatzeko euskarria eta alaka motaren araberakoa da.Batez besteko \(f_3\) balioak 32,0, 31,0 eta 31,3 kHz dira airean, uretan eta gelatinan, hurrenez hurren.\(\text {PTE}_{3}\) \(\gutxi gorabehera\) %74ko eskualdeko batez bestekoa du, edozein eskualdetan baxuena.Eskualdeko batez bestekoa \(\text {FWHM}_{3}\) \(\gutxi gorabehera\) 1085 Hz-koa da, hau da, lehenengo eta bigarren eskualdea baino handiagoa.
Ondorengoak irudiari egiten dio erreferentzia.12 eta 2. taula. Lancetak (L) desbideratu zuen gehien (puntu guztietarako esangura handiarekin, \(p<\) 0,017) bai airean bai uretan (12a. irudia), DPR altuena lortuz (220 µm/ arte). W airean). 12 eta 2. taula. Lancetak (L) desbideratu zuen gehien (puntu guztietarako esangura handiarekin, \(p<\) 0,017) bai airean bai uretan (12a. irudia), DPR altuena lortuz (220 µm/ arte). W airean). Следующее относится к рисунку 12 и таблице 2. Ланцет (L) отклонялся больше всего больше всего всего ех наконечников, \(p<\) 0,017) как в воздухе, так и в воде (рис. 12а), достигая самого высокого высокого . 12. Irudiari eta 2. Taulari honako hau aplikatzen zaio. Lanceta (L) desbideratu zen gehien (puntu guztietarako esangura handia duena, \(p<\) 0,017) bai airean bai uretan (12a. irudia), DPR handiena lortuz.(220 μm/W egin airean).Beheko 12. irudia eta 2. taulara aipatzen dira.柳叶刀(L) 在空气和水中(图12a)中偏转最大(对所有尖端具有高度意义高度意乼,N((对所有尖端具有高度意乎,N\(p0N((对所有最高DPR (空气中高达220 µm/W).柳叶刀(L) du airean eta uretan desbideratze handiena (图12a) (对所述尖端是对尖端是是电影,\(p<\) 0,017), eta DPR altuena lortu du (µm/220 arte). W airean). Ланцет (L) имеет наибольшее отклонение (весьма значимое для всех наконечников, \(p<\) весьма значимое для всех наконечников, \(p<\) весьма значимое для всех наконечников, \(p<\) весьма значимое . 12а), достигая самого высокого DPR (220 мкм/Вт воздухе). Lancetak (L) du desbideratze handiena (oso esanguratsua punta guztietarako, \(p<\) 0,017) airean eta uretan (12a. irud.), DPR altuena lortuz (220 µm/W airean). Airean, AX1 BL handiagoa zuen AX2–3 baino altuago desbideratu zen (esanguratsu, \(p<\) 0,017), AX3 (BL txikiena zuena) AX2 baino gehiago desbideratu zuen 190 µm/W-ko DPR batekin. Airean, AX1 BL handiagoa zuen AX2–3 baino altuago desbideratu zen (esanguratsu, \(p<\) 0,017), AX3 (BL txikiena zuena) AX2 baino gehiago desbideratu zuen 190 µm/W-ko DPR batekin. В воздухе AX1 с более высоким BL отклонялся выше, чем AX2–3 (со значимостью \(p<\) 0,017 тским AX), (со мостью 0,017) м BL) отклонялся больше, чем AX2 с DPR 190 мкм/Вт. Airean, BL handiagoa zuen AX1 AX2–3 baino altuago desbideratu zen (\(p<\) 0,017ko esangurarekin), AX3 (BL baxuena duena) AX2 baino gehiago desbideratu zuen DPR 190 µm/W-rekin.在空气中,具有较高BL 的AX1 偏转高于AX2-3(具有显着性,\(p<\) 0,017),耜弎(具)有显着性,转大于AX2,DPR 为190 µm/W. Airean, BL handiagoa duen AX1-en desbideratzea AX2-3rena baino handiagoa da (nabarmen, \(p<\) 0,017), eta AX3-ren desbideratzea (BL txikiena duena) AX2rena baino handiagoa da, DPR 190 da. µm/W. В воздухе AX1 с более высоким BL имеет большее отклонение, чем AX2-3 (значимо, \(p<\) имеет большее отклонение, чем AX2-3 (значимо, \(p<\) см 0,01AXд), см 0,01AXда зким BL) имеет большее отклонение, чем AX2 с DPR 190 мкм/Вт. Airean, BL handiagoa duen AX1 AX2-3 baino desbideratze handiagoa du (esanguratsua, \(p<\) 0,017), AX3 (BL baxuena duena) AX2 baino desbideratze handiagoa du 190 μm/W-ko DPR-arekin. 20 mm-ko uretan, ez da desberdintasun esanguratsurik aurkitu (\(p>\) 0,017) desbideratzean eta PTEan AX1-3-rako. 20 mm-ko uretan, ez da desberdintasun esanguratsurik aurkitu (\(p>\) 0,017) desbideratzean eta PTEan AX1-3-rako. В воде на глубине 20 мм достоверных различий (\(p>\) 0,017) по прогибу и ФТР для AX1–3 наруж. 20 mm-ko sakoneran dagoen uretan, desbideratzean eta FTR-n desberdintasun esanguratsuak (\(p>\) 0,017) antzeman ziren AX1-3rako.在20 mm 的水中,AX1-3 的挠度和PTE 没有显着差异(\(p>\) 0,017)。 20 mm-ko uretan, ez zegoen AX1-3 eta PTEren arteko alde nabarmenik (\(p>\) 0,017). На глубине 20 мм прогиб eta PTE AX1-3 существенно не отличались (\(p>\) 0,017). 20 mm-ko sakoneran desbideratzea eta PTE AX1-3 ez ziren nabarmen ezberdintzen (\(p>\) 0,017).PTE-mailak uretan (% 90,2-98,4) airean (% 56-77,5) baino handiagoak izan ziren (12c. irudia), eta kabitazio-fenomenoa uretan egindako esperimentuan nabaritu zen (13. irudia, ikus gehigarria ere). informazioa).
Punta tolestura anplitudearen neurketak (batez bestekoa ± desbideratze estandarra, n = 5) L eta AX1-3 txanflan airean eta uretan (20 mm-ko sakonera) txamflanaren geometria aldatzearen eragina agerian utzi zuten.Neurketak maiztasun bakarreko kitzikapen sinusoidala jarraitua erabiliz lortzen dira.(a) Desbideratze gailurra (\(u_y\vec {j}\)) erpinean, (b) dagozkien maiztasun modalen \(f_2\) neurtuta.(c) Potentzia transmisioaren eraginkortasuna (PTE, rms, %) ekuazio gisa.(4) eta (d) Desbideratze-potentzia-faktorea (DPR, µm/W) desbideratze gailurra eta transmisio-potentzia \(P_T\) (Wrms) gisa kalkulatuta.
Abiadura handiko kamera baten itzal grafiko tipikoa, lantza-puntaren (lerro berdeak eta gorriak) desbideratze osoa erakusten duena uretan (L) eta punta axisimetrikoa (AX1-3) uretan (20 mm-ko sakonera), ziklo erdia, gidatzeko maiztasuna. \(f_2\) (maiztasuna 310 kHz-ko laginketa).Hartutako gris-eskalako irudiak 128×128 pixeleko neurriak ditu, \(\gutxi gorabehera) 5 µm-ko pixeleko tamainarekin.Bideoa informazio gehigarrian aurki daiteke.
Horrela, makurtze-uhin-luzeraren aldaketa modelatu dugu (7. irud.) eta transferentziarako mugikortasun mekanikoa kalkulatu dugu hodi-luzeraren eta alakaren konbinazio lantzeolatu, asimetriko eta axial konbentzionalen (8, 9. irudia).Alakatutako geometria simetrikoa.Azken horretan oinarrituta, puntatik soldadurarako distantzia optimoa 43 mm-koa dela estimatu dugu (edo \(\approx\) 2,75\(\lambda_y\) 29,75 kHz-tan) 5. irudian erakusten den moduan, eta hiru alaka axisimetriko fabrikatu ditugu. alakaren luzera desberdinak.Ondoren, haien maiztasun-erantzunak airean, uretan eta %10eko (w/v) gelatina balistikoko lantzeta konbentzionalekin alderatuta (10., 11. irudiak) zehaztu genituen eta inklinazio-desflexio modua alderatzeko kasurik onena zehaztu genuen.Azkenik, puntako desbideratzea neurtu dugu 20 mm-ko sakoneran airean eta uretan olatu okertuz eta injektatutako medioaren potentzia transferitzeko eraginkortasuna (PTE, %) eta desbideratze potentzia-faktorea (DPR, µm/W) kuantifikatu ditugu okertu bakoitzeko.mota (12. irud.).
Emaitzek erakusten dute geometriaren okertze-ardatzak eragiten duela punta-ardatzaren anplitude-desbideratzean.Lancetak kurbadura handiena eta DPR altuena ere izan zuen alaka axisimetrikoarekin alderatuta, eta alaka axisimetrikoak batez besteko desbideratze txikiagoa zuen (12. irudia). Ardatz-simetriko 4 mm-ko alaka (AX1) alaka luzeena duena, estatistikoki esanguratsua den desbideratzerik handiena lortu zuen airean (\(p < 0,017\), 2. taula), beste ardatz-simetriko orratz batzuekin alderatuta (AX2-3), baina ez zen alde nabarmenik antzeman, orratza uretan sartuta. Ardatz-simetriko 4 mm-ko alaka (AX1) alaka luzeena duena, estatistikoki esanguratsua den desbideratzerik handiena lortu zuen airean (\(p < 0,017\), 2. taula), beste ardatz-simetriko orratz batzuekin alderatuta (AX2-3), baina ez zen alde nabarmenik antzeman, orratza uretan sartuta. Осесимметричный скос 4 мм (AX1), имеющий наибольшую длину скоса, достиг статистически зоголь ьна лонения в воздухе (\(p <0,017\), таблица 2) по сравнению с другими осесимметричными иглами (AX2–3ми). Alaka axisimetrikoak 4 mm-ko (AX1), alakaren luzerarik luzeena izanik, estatistikoki esanguratsua den desbideratze handiagoa lortu zuen airean (\(p < 0,017\), 2. taula) beste orratz axisimetriko batzuekin alderatuta (AX2-3).baina ez ziren alde nabarmenik antzeman orratza uretan sartzerakoan.与其他轴对称针(AX2-3) 相比,具有最长斜角长度的轴对称4 mm 斜角(AX1) 在空氡气丟角长度的轴对称着的最高偏转(\(p < 0,017\),表2),但当将针头放入水中时,没有观察到显着差异。 Axialki simetriko beste orratz batzuekin (AX2-3) alderatuta, airean 4 mm-ko angelu zeiharkor luzeena du axialki simetrikoa (AX1), eta estatistikoki esanguratsua den desbideratze maximoa lortu du (\(p < 0,017\), 2. taula) , baina orratza uretan sartuta, ez zen alde nabarmenik antzeman. Осесимметричный скос 4 мм (AX1) с наибольшей длиной скоса обеспечивает статистически значимое максимальное отклонение в воздухе По сравнению с другими ос другими осесимметричными иглами (\2-3) (\ (P <0,017 \), Таблица 2), Но существенной разницы не было. 4 mm-ko luzera handiena duen malda axisimetrikoak (AX1) aireko desbideratze maximo estatistiko esanguratsua eman zuen beste malda axisimetrikoekin alderatuta (AX2-3) (\(p < 0,017\), 2. taula), baina ez zen egon. alde nabarmena.orratza uretan sartzen denean ikusten da.Beraz, alakaren luzera handiagoak ez du abantaila nabarmenik punta puntako desbideratzeari dagokionez.Hori kontuan hartuta, ondorioztatzen da ikerketa honetan ikertzen den malda-geometriak eragin handiagoa duela anplitudearen desbideratzean malda-luzerak baino.Hau tolestura-zurruntasunarekin erlazionatu daiteke, adibidez, tolestuta dagoen materialaren eta eraikuntza-orratzaren lodiera orokorraren arabera.
Azterketa esperimentaletan, islatutako uhin flexuralaren magnitudeak muturraren muga-baldintzek eragiten dute.Orratzaren punta uretan eta gelatinan sartu zenean, \(\text {PTE}_{2}\) %95eko batez bestekoa izan zen \(\gutxi gorabehera\) eta \(\text {PTE}_{2}\) balioen batez bestekoa. ​​% 73 eta % 77 dira (\text {PTE}_{1}\) eta \(\text {PTE}_{3}\), hurrenez hurren (11. irudia).Horrek adierazten du energia akustikoaren gehienezko transferentzia galdaketa-euskarrira (adibidez, ura edo gelatina) \(f_2\\) gertatzen dela.Antzeko portaera ikusi zen aurreko ikerketa batean 41-43 kHz-ko maiztasunetan gailu-egitura sinpleagoak erabiliz, non egileek intercalatutako medioaren modulu mekanikoarekin lotutako tentsio-islaketa-koefizientea frogatu zuten.Sartze-sakonera32 eta ehunaren propietate mekanikoek karga mekanikoa ematen diote orratzari eta, beraz, UZeFNAB-en oihartzun-jokabidean eragina izatea espero da.Hori dela eta, 17, 18, 33 bezalako erresonantzia-jarraipenaren algoritmoak erabil daitezke arkatza bidez emandako soinuaren potentzia optimizatzeko.
Tolestura-uhin-luzeraren modelizazioak (7. irudia) erakusten du ardatz-simetrikoak egitura-zurruntasun handiagoa duela (hau da, tolestura-zurruntasun handiagoa) puntan lancetak eta alaka asimetrikoak baino.(1)tik eratorrita eta ezaguna den abiadura-maiztasuna erlazioa erabiliz, lancetaren tolestura-zurruntasuna, mutur asimetrikoak eta axisimetrikoak \(\gutxi gorabehera) 200, 20 eta 1500 MPa-ko malda gisa estimatzen ditugu, hurrenez hurren.Hau (\lambda _y\) 5,3, 1,7 eta 14,2 mm-i dagokio 29,75 kHz-an, hurrenez hurren (7a-c. irudiak).USeFNAB prozeduraren segurtasun klinikoa kontuan hartuta, geometriak alakaren diseinuaren zurruntasunean duen eragina ebaluatu behar da34.
Alakaren eta hodiaren luzeraren parametroen azterketak (9. irudia) erakutsi zuen asimetrikorako TL tarte optimoa (1,8 mm) alaka axisimetrikorako (1,3 mm) baino handiagoa zela.Gainera, mugikortasun-lautada 4 eta 4,5 mm-tik eta 6-7 mm bitartekoa da inklinazio asimetrikorako eta ardatz-simetrikorako, hurrenez hurren (9a, b irudia).Aurkikuntza honen garrantzia praktikoa fabrikazio-perdoietan adierazten da, adibidez, TL optimoaren tarte txikiagoak luzera zehaztasun handiagoaren beharra suposa dezake.Aldi berean, etekin-plataformak tolerantzia handiagoa ematen du maiztasun jakin batean malda luzera aukeratzeko errendimenduan nabarmen eragin gabe.
Azterketak honako muga hauek ditu.Orratzaren desbideratzea zuzenean neurtzeak ertzak detektatzeko eta abiadura handiko irudiak erabiliz (12. irudia) esan nahi du optikoki gardenak diren euskarrietara mugatzen garela, hala nola airea eta ura.Era berean, adierazi nahi dugu ez dugula esperimenturik erabili transferentzia-mugikortasun simulatua probatzeko eta alderantziz, baizik eta FEM azterketak erabili ditugula fabrikatutako orratzaren luzera optimoa zehazteko.Muga praktikoen ikuspuntutik, lantzek puntatik mahukarainoko luzera beste orratz (AX1-3) baino 0,4 cm luzeagoa da, ikus irud.3b.Horrek egitura azikularraren erantzun modalean eragina izan dezake.Horrez gain, uhin-gidaren berunezko soldaduraren formak eta bolumenak (ikus 3. irudia) pin diseinuaren inpedantzia mekanikoan eragina izan dezake, inpedantzia mekanikoan eta tolestean jokabidean akatsak eraginez.
Azkenik, esperimentalki frogatu dugu alakaren geometriak USeFNAB-en deflexio-kopuruari eragiten diola.Desbideratze-anplitude handiagoak orratzak ehunean duen eraginan eragin positiboa izan dezakeen egoeretan, adibidez, zulatu ondoren ebaketa-eraginkortasuna, USeFNAB-rako ohiko lanceta bat gomenda daiteke, desbideratze-anplitude handiena ematen baitu zurruntasun nahikoa mantenduz. diseinuaren puntan.Horrez gain, azken ikerketa batek frogatu du punta desbideratze handiagoak efektu biologikoak hobetu ditzakeela, hala nola kavitazioa, eta horrek esku-hartze kirurgiko gutxien inbaditzaileak egiteko aplikazioak garatzen lagun dezake.Potentzia akustiko osoa handitzeak USeFNAB13-ren biopsiaren etekina areagotzen duela frogatu denez, laginaren etekinaren eta kalitatearen azterketa kuantitatibo gehiago behar dira aztertutako orratz-geometriaren onura kliniko zehatza ebaluatzeko.
Frable, WJ Orratz finaren aspirazioko biopsia: berrikuspena.Humph.Gaixo.14:9-28.https://doi.org/10.1016/s0046-8177(83)80042-2 (1983).


Argitalpenaren ordua: 2022-10-13
WhatsApp Online Txata!