Het huidige werk is volledig gebaseerd op analytische modellering, symbolische berekening en numerieke simulaties om gekoppelde roterende magneto-foto-thermo-elastische interacties in anisotrope halfgeleidermedia te onderzoeken. Er waren geen menselijke deelnemers, dierproeven, klinische gegevens of biologische monsters betrokken bij deze studie. Daarom waren ethische goedkeuring en geïnformeerde toestemming niet vereist.
Wiskundige formulering van het magneto-foto-thermo-elastische probleem in een roterende, vezelversterkte anisotrope halfgeleiderhalfruimte
In de huidige studie wordt een tweedimensionale roterende, vezelversterkte anisotrope halfgeleiderhalfruimte onderzocht onder optische excitatie en een toegepast magnetisch veld. Het medium bezet het semi-oneindige gebied x ≥ 0, waarbij de rand bij x = 0 het blootgestelde oppervlak vertegenwoordigt dat aan externe optische belasting is blootgesteld. Het coördinatensysteem wordt zo gekozen dat de x-as in het medium uitloopt, terwijl de y-as langs het oppervlak ligt, wat het in-vlak gedrag van de structuur weergeeft. Het toegepaste magnetisch veld en de hoeksnelheidsvector worden beide genomen langs de z-as. Het halfgeleidermedium wordt verondersteld homogeen en lineair elastisch te zijn, terwijl anisotropie wordt geïntroduceerd door uitgelijnde versterkende vezels die in de x-richting zijn ingebed. Optische absorptie aan de grens produceert lokale verwarming en overtollige ladingsdragers, wat resulteert in gekoppelde thermische, mechanische en draaggolfinteracties binnen het medium. Daarnaast introduceert het magnetisch veld elektromagnetische koppelingseffecten, terwijl de rotatiebeweging inertiële effecten veroorzaakt die de voortplanting van thermo-elastische golven en de algehele fysieke respons aanzienlijk beïnvloeden. Dienovereenkomstig wordt de fysieke toestand van het medium weergegeven door het temperatuurveld T(x, y, t)(K), de draagkrachtdichtheid N(x, y, t)(m-3) en de verplaatsingscomponenten u(x, y,t)(m) en v(x,y,t)(m), onder de aanname van kleine vervormingen. Figuur 1 illustreert de geometrie van het probleem, inclusief het coördinatensysteem, optische excitatie, magnetisch veld, rotatie-effect en vezeloriëntatie. De huidige formulering is toepasbaar op homogene anisotrope vezelversterkte halfgeleidermedia die opereren binnen het kleindeformatieregime en het kader van lineaire thermoelasticiteit. Het model gaat uit van een vaste vezeloriëntatie en constante materiaaleigenschappen door het hele medium. Daarom worden niet-lineair materiaalgedrag, grote vervormingen, materiaalschade en ruimtelijke variaties in materiaaleigenschappen niet meegenomen in de huidige studie. Daarom is het voorgestelde model bedoeld voor matige belastingsomstandigheden waarbij de respons binnen het lineaire bereik blijft. In de huidige studie wordt de optische excitatie gemodelleerd met behulp van voorgeschreven randvoorwaarden voor oppervlaktetemperatuur en fotogegenereerde draagdrijvendheid. Het gedetailleerde laser-materie interactieproces, inclusief optische absorptie, penetratiediepte en intensiteitsverdeling, wordt niet expliciet behandeld. In plaats daarvan wordt het netto-effect weergegeven door de randamplitudes θ0 en N0, die de thermische en draaggolfexcitaties karakteriseren die door het invallende optische veld worden geïnduceerd.

Figuur 1: Schematische weergave van de roterende anisotrope, vezelversterkte halfgeleiderhalfruimte die wordt blootgesteld aan optische excitatie en een extern magnetisch veld. De figuur illustreert de fysieke configuratie van het probleem, inclusief het coördinatensysteem, optische excitatie, toegepast magnetisch veld, vezeloriëntatie en rotatie-effecten die in de huidige formulering worden beschouwd. Klik hier om een grotere versie van deze figuur te bekijken.
De constitutieve relatie voor de spanningstensor in een vezelversterkt anisotrop thermo-elastisch halfgeleidermedium kan in de gegeneraliseerde vorm als volgt worden uitgedrukt: 12,15,16.
(1)
Hier duidt σij de componenten van de spanningstensor aan, Cijkl zijn de elastische stijfheidscoëfficiëntententen, ekl de rektensor, T is de temperatuurtoename ten opzichte van de referentietemperatuur T0, en N duidt de overtollige draagkrachtdichtheid aan. De tensoren βij en ηij komen overeen met respectievelijk de thermo-elastische en dragerkoppelingscoëfficiënten. Dienovereenkomstig kan de constitutieve relatie inclusief de expliciete invloed van vezelversterking worden geschreven als 12,15:
(2)
In deze formulering zijn λ en μT de elastische Lamé-constanten, terwijl μL de longitudinale schuifmodulus langs de vezelrichting aanduidt. De parameters α en β beschrijven de versterkingseffecten die samenhangen met de ingebedde vezels. De grootheid δij is het Kronecker-deltasymbool, en ai zijn de componenten van de eenheidsvector die de vezeloriëntatie definieert. Voor het huidige model zijn de versterkingsvezels uitgelijnd langs de x -richting zodanig dat a = (1,0). De termen die βijθ en ηijN bevatten, vertegenwoordigen respectievelijk thermische en draaggolfkoppelingseffecten. Voor de huidige tweedimensionale configuratie verminderen de regelende spanningscomponenten tot de volgende vormen12,15:
. (3)
. (4)
. (5)
Hier
duiden en
de verplaatsingscomponenten respectievelijk aan langs de x- en y-richtingen , terwijl Aij de effectieve elastische coëfficiënten zijn van het anisotrope vezelversterkte medium. De thermo-elastische en dragerkoppelingscoëfficiënten worden als volgt gedefinieerd
,
,
,
.
In bovenstaande relaties vertegenwoordigen αij de thermische expansiecoëfficiënten, terwijl ξij de draaggolfexpansiecoëfficiënten aanduiden die geassocieerd zijn met het halfgeleidermedium. De effectieve elastische coëfficiënten van het vezelversterkte anisotrope medium worden gegeven door
,
,
, . 
Deze coëfficiënten karakteriseren de anisotrope elastische respons van het versterkte halfgeleidermateriaal en beschrijven hoe vezeloriëntatie het gekoppelde thermo-elastische gedrag beïnvloedt. Om rekening te houden met de invloed van elektromagnetische interacties in de huidige magneto-foto-thermo-elastische formulering, wordt aangenomen dat er een uniform magnetisch veld wordt aangelegd langs de z-richting, die normaal is op het x - y vervormingsvlak. Dienovereenkomstig wordt de magnetische veldvector beschouwd in de vorm23,25 ,
, waarbij H0 de constante magnetische veldintensiteit aanduidt. Aangezien de huidige formulering beperkt is tot tweedimensionale vervormingen, wordt het verplaatsingsveld van het medium genomen als
, waarbij
en
respectievelijk de verplaatsingscomponenten langs de x- en y- richting vertegenwoordigen.
Onder de aannames van kleine vervormingen en een langzaam bewegend elektrisch geleidend halfgeleidermedium, genereert de interactie tussen de deeltjessnelheid en het aangelegde magnetisch veld een geïnduceerd elektrisch veld.
Op basis van Maxwells elektromagnetische relaties voor bewegende geleidende media kan de geïnduceerde elektrische veldvector worden uitgedrukt als19,23
. (6)
waarbij μ0 de magnetische permeabiliteit aanduidt en
de deeltjessnelheidsvector is. Door de uitdrukkingen van en
in de bovenstaande relatie te plaatsen, krijg je
. (7)
wat resulteert in
. (8)
Door de tijdsafgeleide van het geïnduceerde elektrische veld te nemen, verkrijgen we
. (9)
De magnetische perturbatievector die wordt gegenereerd door de vervorming van het halfgeleidermedium wordt gedefinieerd als23˒25
. (10)
De bovenstaande uitdrukking voldoet automatisch aan Maxwells divergentievoorwaarde voor het magnetische perturbatieveld, namelijk
. (11)
Om de elektrische stroomdichtheid te bepalen, wordt eerst de rotatie van de magnetische perturbatievector geëvalueerd. In determinantvorm kan de rotatie-operator worden geschreven als23
. (12)
Het uitbreiden van de determinant leidt tot
. (13)
De elektrische stroomdichtheidsvector wordt dan verkregen uit Maxwells elektromagnetische vergelijking23
. (14)
waarbij ε0 de elektrische permittiviteit van het medium aanduidt.
De elektromagnetische lichaamskracht die op het halfgeleidermedium werkt, wordt bepaald met behulp van de Lorentzkrachtrelatie23
. (15)
Het vectorproduct
kan in determinantvorm worden geëvalueerd als
. (16)
Door de voorgaande uitdrukkingen in vergelijking (15) in te vullen, worden de componenten van de elektromagnetische lichaamkrachtvector
. (17)
. (18)
Deze relaties geven duidelijk aan dat het toegepaste magnetisch veld extra koppelingsmechanismen toevoegt aan de beheersende vergelijkingen via zowel stijfheidsachtige elektromagnetische termen als aangepaste traagheidstermen evenredig aan
. Daardoor beïnvloedt het magnetisch veld aanzienlijk de voortplantingseigenschappen van thermo-elastische golven en het algehele dynamische gedrag van het vezelversterkte anisotrope halfgeleidermedium. Naast de elektromagnetische effecten wordt de invloed van rotatie in de huidige formulering opgenomen om de dynamische respons van het medium te beschrijven wanneer deze vanuit een roterend referentiekader wordt waargenomen. Het vezelversterkte halfgeleidermedium wordt verondersteld een uniforme starre rotatie te ondergaan met een constante hoeksnelheidsvector gegeven door33
, waarbij
de constante hoeksnelheid rond de z-as aanduidt. Aangezien de rotatieas normaal is op het x-y-vlak , ontstaan er bij het formuleren van bewegingsvergelijkingen in een roterend coördinatensysteem extra traagheidsversnellingen door de niet-traagheidsaard van het roterende frame. Deze versnellingen bestaan voornamelijk uit Coriolis- en centrifugale versnellingen. De Coriolisversnelling is geassocieerd met het deeltjessnelheidsveld en wordt uitgedrukt als31˒33
. (19)
Door de uitdrukkingen van
en
te substitueren, verkrijgen we
. (20)
Daarom wordt de Coriolisversnelling
. (21)
De centrifugale versnelling hangt direct af van het verplaatsingsveld zelf en wordt weergegeven door 31,35
. (22)
Eerst wordt het vectorproduct
geëvalueerd als
. (23)
Vervolgens levert het invullen van het verkregen resultaat in de centrifugale versnellingsrelatie op
. (24)
Dienovereenkomstig kan de totale rotatiebijdrage die in de besturende vergelijkingen voorkomt, worden uitgedrukt als
. (25)
Daarom worden de componenten van rotatieversnelling in de x- en y- richting
, (26)
. (27)
De bovenstaande uitdrukkingen tonen aan dat rotatiebeweging extra koppeling introduceert tussen de verplaatsingscomponenten via de Coriolisversnelling, naast verplaatsingsafhankelijke traagheidseffecten door de centrifugale versnelling. Daarom veroorzaakt de gecombineerde werking van het magnetisch veld en de rotatie aanzienlijke aanpassingen aan de dynamische respons en golfvoortplantingskenmerken van het roterende, vezelversterkte anisotrope halfgeleidermedium. Om de gecombineerde invloed van elektromagnetische interacties en rotatiebeweging te integreren, worden de bewegingsvergelijkingen voor het vezelversterkte anisotrope halfgeleidermedium gegeneraliseerd om zowel de elektromagnetische lichaamskracht als de extra traagheidsversnellingen die ontstaan in een roterend referentiekader te omvatten. Dienovereenkomstig kan de algemene bewegingsvergelijking voor een vervormbaar roterend continuüm als volgt worden uitgedrukt:31˒32
. (28)
Hier duidt ρ de massadichtheid aan, en Fi vertegenwoordigt de componenten van elektromagnetische lichaamskracht. Bovendien duidt de
hoeksnelheidsvector van het roterende frame aan. De gelijkstellingen. (22) en (24) komen overeen met respectievelijk de Coriolis- en centrifugale versnellingen. De bewegingsvergelijkingen in de x- en y- richting kunnen worden geschreven als
. (29)
. (30)
Door de eerder verkregen elektromagnetische lichaamskrachtcomponenten in bovenstaande vergelijkingen in te vullen, krijg je
. (31)
. (32)
Vervolgens wordt de constitutieve relaties die overeenkomen met het vezelversterkte anisotrope halfgeleidermedium in bovenstaande vergelijkingen worden vervangen, wat de gekoppelde bewegingsvergelijkingen oplevert in termen van de verplaatsingscomponenten, temperatuurveld en draaggolfdichtheid32
. (33)
. (34)
Ten slotte kunnen, met behulp van de uitdrukking van het magnetische perturbatieveld zoals eerder gegeven in vergelijking (15), en deze in bovenstaande vergelijkingen invullen, de bepalende bewegingsvergelijkingen in hun definitieve gekoppelde vorm worden geschreven als
. (35)
. (36)
Deze vergelijkingen tonen de gekoppelde invloed van rotatie- en magnetische veldeffecten op de thermo-elastische respons van het medium. De rotatietermen verklaren zowel Coriolis- als centrifugale bijdragen, terwijl het magnetisch veld extra elektromagnetische koppeling introduceert en het dynamische gedrag van het systeem verandert. Daardoor creëren de beheersende vergelijkingen een uniform kader voor het analyseren van golfvoortplanting en multifysische interacties in roterende magneto-foto-thermo-elastische vezelversterkte halfgeleiders. Bij optische excitatie wordt het thermisch gedrag van het halfgeleidermedium aanzienlijk beïnvloed door de wisselwerking tussen warmtegeleiding, draaggolftransport en mechanische vervorming, wat leidt tot een sterk gekoppeld thermo-fotoelastisch proces. In tegenstelling tot het klassieke warmtegeleidingsmodel wordt de temperatuurverdeling in halfgeleidermaterialen niet alleen beïnvloed door thermische diffusie, maar ook door dragerrecombinatie en thermo-elastische koppeling. Daarom kan de gegeneraliseerde warmtegeleidingsvergelijking voor het anisotrope vezelversterkte halfgeleidermedium als volgt worden uitgedrukt7˒23.
. (37)
De bovenstaande vergelijking toont duidelijk aan dat het thermisch veld binnen het roterende magneto-foto-thermo-elastische halfgeleidermedium wordt bepaald door de gecombineerde invloed van anisotrope warmtegeleiding, dragerrecombinatieprocessen en thermo-elastische interacties. De term
beschrijft de thermische energie die wordt opgewekt door dragerrecombinatie onder optische excitatie, terwijl de koppelingstermen met
en
de invloed van tijdsafhankelijke mechanische vervorming op de thermische respons van het medium aangeven. Daardoor wordt het temperatuurveld sterk gekoppeld aan zowel de draaggolfdichtheid als het elastische veld, wat een belangrijke rol speelt in de voortplantingskenmerken van thermo-elastische golven in vezelversterkte halfgeleidermaterialen. In de huidige formulering wordt de evolutie van de dragerconcentratie N(x, y, t) binnen het halfgeleidermedium bepaald door de gecombineerde effecten van draaggolfdiffuutifsie, recombinatieprocessen en thermische activatie die worden opgewekt door optische excitatie. Dienovereenkomstig kan de dragertransportvergelijking die de niet-evenwichtsdragerdynamica beschrijft als volgt worden geschreven 4,23
. (38)
Hier duidt DE de draaggolfdiffusiecoëfficiënt aan, terwijl
de tweedimensionale Laplaceiaanse operator in het x -y vlak vertegenwoordigt. De term
correspondeert met het dragerrecombinatie-effect dat geassocieerd is met de draaggolflevensduur τ. Bovendien is κ de thermodragerkoppelingsparameter gedefinieerd door
, waarbij N0 de evenwichtsdragerconcentratie aanduidt. De koppelingsterm κT beschrijft de invloed van het temperatuurveld op de generatie van overtollige dragers binnen het halfgeleidermedium. De bovenstaande vergelijking toont aan dat de dragerconcentratie sterk gekoppeld is aan het thermisch veld via thermisch geactiveerde draagdragergeneratiemechanismen. Daardoor worden de draagkrachtdynamica sterk afhankelijk van zowel thermische diffusie- als recombinatie-effecten, die de gekoppelde foto-thermo-elastische respons van het roterende, vezelversterkte anisotrope halfgeleidermedium aanzienlijk beïnvloeden. De beheersende vergelijkingen en wiskundige formulering van het gekoppelde magneto-foto-thermo-elastische halfgeleidersysteem zijn nu volledig vastgesteld. De fysische en materiaalparameters voor het siliciummedium worden samengevat in Tabel 2, samen met hun numerieke waarden, eenheden en bijbehorende referenties. Deze parameters worden vervolgens gebruikt in de numerieke berekeningen en de nondimensionalisatieprocedure.
Niet-dimensionale formulering van het roterende magneto-foto-thermo-elastische vezelversterkte halfgeleidermodel
Om de beheersvergelijkingen te vereenvoudigen en een compacte wiskundige representatie van het gekoppelde roterende magneto-foto-thermo-elastische systeem te verkrijgen, worden passende karakteristieke schalen geïntroduceerd om de fysische variabelen niet-dimensionaal te maken. Deze niet-dimensionalisatieprocedure vermindert het aantal bepalende materiaalparameters en vergemakkelijkt de analytische en numerieke behandeling van de gekoppelde vergelijkingen. De geselecteerde karakteristieke grootheden worden consistent gekozen met de thermo-elastische, elektromagnetische, rotatie- en draagdragertransporteigenschappen van het halfgeleidermedium16,21 Dienovereenkomstig worden de volgende dimensieloze variabelen geïntroduceerd:
, , , , ,
,
, 



, , . 


Hier duidt CT de karakteristieke elastische golfsnelheid aan, terwijl t* de karakteristieke thermische relaxatietijd vertegenwoordigt die geassocieerd is met het gekoppelde thermo-elastische proces. Bovendien definieert de parameter de niet-dimensionale rotatieparameter
die de invloed van rotatiebeweging op het dynamische gedrag van het medium karakteriseert. Door bovenstaande dimensieloze grootheden in de eerder afgeleide beheersvergelijkingen te substitueren, wordt het gekoppelde systeem genormaliseerd in een genormaliseerde vorm. Deze transformatie vereenvoudigt de wiskundige structuur van de vergelijkingen aanzienlijk en biedt een geschikt kader om de gecombineerde effecten van magnetisch veld, optische excitatie, rotatie, anisotropie en drager-transportinteracties te onderzoeken. Voor de eenvoud wordt de priemnotatie die bij de dimensieloze grootheden hoort, in de volgende analyse weggelaten. Dienovereenkomstig kunnen de beheersvergelijkingen van het gekoppelde roterende magneto-fotothermo-elastische systeem worden geschreven in de volgende niet-dimensionale vorm16,20:
, (39)
, (40)
, (41)
. (42)
De overeenkomstige niet-dimensionale spanningscomponenten van het roterende anisotrope vezelversterkte halfgeleidermedium worden als volgt verkregen:
, (43)
, (44)
. (45)
De niet-dimensionale coëfficiënten ai (i = 1,2,...,18) vertegenwoordigen combinaties van de fysische, thermische, elektromagnetische, draaggolf- en rotatieparameters van het gekoppelde halfgeleidermedium. Deze coëfficiënten karakteriseren de invloed van anisotropie, vezelversterking, magnetisch veld, thermo-elastische koppeling, draagdragertransport en rotatiebeweging op het algehele gedrag van het systeem. Daardoor bieden de verkregen niet-dimensionale regelvergelijkingen een compact en efficiënt wiskundig model voor het analyseren van de gekoppelde golfpropagatiefenomenen en multifysische interacties in roterende magneto-foto-thermo-elastische, vezelversterkte halfgeleidermedia. De niet-dimensionale parameters ai, γi en δi werden geïntroduceerd om compacte combinaties van de fysische en materiaaleigenschappen te representeren die het gekoppelde roterende magneto-foto-thermo-elastische gedrag van het vezelversterkte anisotrope halfgeleidermedium bepalen. Elke coëfficiënt weerspiegelt een specifiek interactiemechanisme binnen het gekoppelde systeem en geeft inzicht in de relatieve invloed van de onderliggende fysieke processen. Voor duidelijkheid worden de dimensieloze parameters en hun bijbehorende definities samengevat in Tabel 1.
Tabel 1: Definities en fysieke interpretaties van de dimensieloze parameters die in de huidige formulering worden gebruikt. De tabel vat de dimensieloze parameters samen die voorkomen in de bepalende vergelijkingen, samen met hun fysieke betekenissen en hun rollen bij het beschrijven van de gekoppelde thermo-elastische, elektromagnetische, dragerdichtheids- en rotatieinteracties. Klik hier om deze tabel te downloaden.
Analytische oplossing met behulp van de normaalmodustechniek
Om de analytische oplossing van het gekoppelde roterende magneto-fotothermo-elastische systeem af te leiden, wordt de normaalmodustechniek toegepast. Deze methode wordt veel gebruikt in gegeneraliseerde thermoelasticiteit en halfgeleidertheorieën vanwege de effectiviteit ervan in het transformeren van gekoppelde partiële differentiaalvergelijkingen tot een gereduceerd systeem van gewone differentiaalvergelijkingen. Een dergelijke benadering is bijzonder nuttig bij de analyse van golfvoortplanting, demping en multifysische interacties in anisotrope halfgeleidermedia. Na de normalmodusanalyse wordt aangenomen dat alle fysische veldgrootheden harmonisch variëren ten opzichte van de tijd en de transversale ruimtelijke coördinaat . Dienovereenkomstig worden het temperatuurveld, de draagkrachtdichtheid, de verplaatsingscomponenten en de spanningsgrootheden weergegeven in de exponentiële vorm24,27
. (46)
Hier duidt ω de complexe frequentieparameter aan die de temporele variatie van de fysische velden bepaalt. Het reële deel van ω wordt geassocieerd met de temporele demping (of groei) van de golfamplitude, terwijl het imaginaire deel het oscillerend gedrag van de voortplantende modus vertegenwoordigt. Deze interpretaties zijn consistent met de conventionele normaalmodusanalyse die in deze studie is toegepast, terwijl a het golfgetal vertegenwoordigt dat geassocieerd is met de ruimtelijke variatie langs de y-richting. De grootheden
, en
komen overeen met de veldamplitudes die alleen afhankelijk zijn van de ruimtelijke coördinaat x. Door bovenstaande normaalmodusrepresentaties te vervangen in de eerder verkregen niet-dimensionale bestuursvergelijkingen en de resulterende uitdrukkingen te vereenvoudigen, wordt het oorspronkelijke gekoppelde partiële differentiaalsysteem getransformeerd in een set gewone differentiaalvergelijkingen ten opzichte van de ruimtelijke coördinaat x. Bijgevolg nemen de beheersende vergelijkingen in het getransformeerde domein de volgende vorm aan:
, (47)
, (48)
, (49)
. (50)
Bovendien worden de overeenkomstige getransformeerde spanningscomponenten verkregen als
, (51)
, (52)
. (53)
Hier duidt
, de differentiaaloperator aan ten opzichte van de ruimtelijke coördinaat . Het verkregen getransformeerde systeem vormt de wiskundige basis voor het construeren van de karakteristieke vergelijking en het afleiden van de volledige analytische oplossing van het gekoppelde roterende magneto-foto-thermo-elastische probleem. De coëfficiënten die voorkomen in de getransformeerde vergelijkingen worden als volgt gedefinieerd:
, , , , ,
,
,
, 



, 


Deze coëfficiënten bevatten de gecombineerde bijdragen van anisotrope elasticiteit, interactie tussen magnetisch veld, thermische koppeling, dragertransport en rotatie-effecten. Daarom biedt het getransformeerde systeem een compacte weergave die geschikt is om de karakteristieke wortels te verkrijgen en het voortplantingsgedrag van gekoppelde golven binnen het roterende vezelversterkte halfgeleidermedium te onderzoeken.
Aanvullend bestand 1: Analytische oplossing met behulp van matrixformulier. Dit bestand bevat de gedetailleerde matrixformulering, de eigenwaardeoplossingsprocedure, de afleiding van karakteristieke vergelijkingen en tussenliggende analytische stappen die zijn gebruikt om de algemene oplossing van het gekoppelde magneto-foto-thermo-elastische model te verkrijgen. Klik hier om dit bestand te downloaden.
De gedetailleerde matrixformulering, eigenwaardeoplossingsprocedure en karakteristiekvergelijkingsafleiding zijn opgenomen in Aanvullend Bestand 1.
Randvoorwaarden en bepaling van de onbekende constanten
Om de analytische formulering te voltooien, werden de verkregen algemene oplossingen vervangen door de voorgeschreven randvoorwaarden die zijn opgelegd aan het oppervlak x = 0. Deze substitutie genereerde een gekoppeld algebraïsch systeem met de onbekende amplitudeconstanten
. Elke randvereiste die verband houdt met temperatuur, draagkrachtdichtheid, mechanische verplaatsing en spanningsbeperkingen werd uitgedrukt in termen van de toelaatbare eigenmodi, wat resulteerde in een lineaire set vergelijkingen die de coëfficiënten
met elkaar relateerden. Voor het gemak werd het resulterende algebraïsche systeem herschreven in compacte matrixvorm als BC = D, waarbij B de coëfficiëntmatrix aanduidt die is geconstrueerd uit de eigenvectorcomponenten geëvalueerd op het randoppervlak,
de vector van onbekende constanten vertegenwoordigt, en D overeenkomt met de vector die wordt gegenereerd uit de opgelegde randvoorwaarden, inclusief de thermische belastingsparameter θ0, de draaggolfexcitatieterm N0, en de voorgeschreven verplaatsingsvoorwaarden. Na het evalueren van deze constanten werden ze teruggezet in de algemene uitdrukkingen van de veldvariabelen om de volledige analytische oplossingen te verkrijgen. Deze uitdrukkingen werden vervolgens gebruikt in de numerieke berekeningen en grafische visualisatie van de thermo-elastische, draagkrachtdichtheids- en verplaatsingsvelden binnen het roterende, vezelversterkte anisotrope halfgeleidermedium. De aangenomen randvoorwaarden vertegenwoordigen een optisch verlicht halfgeleideroppervlak dat gelijktijdig wordt onderhevig aan thermische en draaggolfexcitaties. De voorgeschreven temperatuurconditie modelleert de thermische belasting die wordt opgewekt door het invallende optische veld, terwijl de dragerdichtheidsconditie rekening houdt met de fotogegenereerde overtollige dragers die door optische verlichting worden geproduceerd. Daarnaast vertegenwoordigt de transversale verplaatsingsbeperking de mechanische opsluiting van het oppervlak in de -richting, terwijl de nul-schuifspanningsvoorwaarde overeenkomt met een tangentieel tractievrije grens. Daardoor bieden de geselecteerde gemengde thermische, elektronische en mechanische randvoorwaarden een fysisch consistente weergave van gekoppelde foto-thermo-elastische interacties aan het halfgeleideroppervlak en bieden ze de noodzakelijke beperkingen om de onbekende constanten van de oplossing te bepalen. De opgelegde randvoorwaarden worden als volgt gegeven:
Temperatuurbeperking:
. (78)
Deze voorwaarde vertegenwoordigt een harmonisch variërende oppervlaktetemperatuur die wordt geïnduceerd door periodieke optische verwarming. Het fungeert als de primaire thermische excitatie die de gekoppelde thermo-elastische en dragertransportprocessen binnen het medium aandrijft. De amplitude θ0 karakteriseert de intensiteit van de toegepaste thermische belasting.
Dragerdichtheidsbeperking:
. (79)
Deze randvoorwaarde beschrijft de foto-gegenereerde draagwaterdichtheid die voortkomt uit optische verlichting. Het weerspiegelt de elektronische excitatie door fotonabsorptie en de harmonische modulatie ervan die consistent is met het invallende optische veld.
Verplaatsingsbeperking:
. (80)
Deze voorwaarde geeft aan dat de grens mechanisch beperkt is in de transversale richting. Daarom treedt er geen verplaatsing op langs de -richting aan het oppervlak.
Afschuifspanningsbeperking:
. (81)
Deze voorwaarde komt overeen met een trekvrije grens met betrekking tot afschuifspanning. Het zorgt ervoor dat er geen tangentiële krachten op het oppervlak werken, wat consistent is met een mechanisch vrije grens in de tangentiële richting. Naast de randvoorwaarden bij x = 0 werd de fysieke eis bij oneindig opgelegd als:
het waarborgen van begrensde fysieke oplossingen binnen het semi-oneindige domein. Om een duidelijk overzicht te geven van de analytische en computationele procedure die in deze studie is toegepast, worden de belangrijkste stappen van de oplossingsmethodologie samengevat in Figuur 2.

Figuur 2: Stroomdiagram van de analytische oplossingsprocedure die in deze studie is toegepast, inclusief de formulering van de leidende vergelijkingen, normalmodusanalyse, eigenwaardeoplossing, toepassing van randvoorwaarden en evaluatie van de fysische veldvariabelen. Het diagram vat de belangrijkste rekenkundige stappen samen om de analytische oplossing te verkrijgen en vervolgens de gekoppelde magneto-foto-thermo-elastische respons numeriek te evalueren. Klik hier om een grotere versie van deze figuur te bekijken.