基于特征值法的各向异性纤维增强硅中耦合热-光弹性场的时空分析

摘要

本研究采用正则模态与特征值方法，分析了各向异性纤维增强硅材料中的热–光弹性耦合场。结果表明，该场具有空间衰减特性和时间依赖性的演化行为，且对各向异性表现出强敏感性。热图展示了场分布与局域化特征。

摘要

本研究探讨了各向异性纤维增强硅半导体介质中热-光弹性耦合系统，旨在揭示热场、载流子场与力学场之间的相互作用。此类各向异性纤维增强材料在现代工程应用中具有关键作用，包括微电子与光电子器件、基于激光的技术、传感器以及先进复合结构，这些领域均要求具备方向性特性及增强的力学性能。在承受热载荷与光载荷的半导体元件设计中，此类材料尤为重要，因为准确预测多场耦合行为对于提高器件可靠性与优化性能至关重要。控制方程基于耦合物理模型建立，并进一步转化为无量纲形式，以简化分析并突出相关参数的相对影响。通过采用正规模态法将问题转化为一阶向量-矩阵微分系统，继而运用特征值方法求解，获得满足半无限域内给定边界条件的解析解。通过数值分析研究了时间变化对所有物理场的影响，结果表明各物理场具有显著的空间衰减特性，其耦合行为受材料各向异性与纤维增强效应主导。采用时空热图可视化各物理场的演化过程与局域化特征，揭示了多物理场相互作用的物理机制，并验证了所采用解析方法的有效性。

引言

热-光弹性学已成为一门重要的多学科领域，用于描述半导体材料在光热激发下热、力学与光学效应之间的相互作用。在涉及激光加热和光激发的现代应用中，这些场之间的耦合效应尤为显著。例如，Saeed¹采用双曲双温度模型研究了半导体中的热-光弹性相互作用，证明了热弛豫效应在准确预测系统行为中的重要性。纤维增强复合材料的热力响应也因其增强的力学强度和各向异性特性而受到广泛关注。Li 和 Lambros²分析了此类复合材料的动态热力行为，突出了其在先进工程应用中的适用性。类似地，Kalkal 等人³研究了旋转功能梯度纤维增强介质在磁场作用下的二维（2D）变形，揭示了各向异性和外部效应对系统响应的显著影响。此外，Pitarresi 等人⁴探讨了宏观非均质性在热弹性行为中的作用，强调了对复合材料进行精确建模的必要性。近年来的研究已将这些工作拓展至具有多物理场耦合效应的半导体介质。Mondal 等人⁵研究了考虑记忆响应和磁场作用的增强型半导体中的波传播，揭示了复杂的耦合相互作用。实验研究，如 Akai 等人⁶的工作，通过热弹性温度变化评估了纤维增强复合材料的疲劳损伤，证实了其实际重要性。此外，Lu 等人⁷展示了采用微力学方法估算有效热弹性性能的应用。外部载荷和环境条件的影响也得到了广泛研究。Barak 和 Dhankhar⁸分析了功能梯度纤维增强介质中的斜向加载，而 Kundu 和 Kalkal⁹研究了重力和移动热载荷作用下的光热相互作用。Chaudhary 等人¹⁰利用双相滞后模型研究了温度依赖性材料性能，Pandit 等人¹¹则应用分数阶应变模型来捕捉非局部变形行为。这些研究凸显了在热弹性分析中考虑真实载荷和材料条件的重要性。

半导体系统中的动态和振动相关现象也得到了研究。Song 等人¹²研究了半导体结构中的光热振动，而 Mondal 和 Sur¹³分析了具有记忆效应的正交各向异性介质中的波传播。Abouelregal 等人¹⁴考虑了粘弹性与微结构效应，Hobiny 等人¹⁵研究了斜坡型加热效应。此外，Zenkour¹⁶发展了诸如三相滞后理论等广义热弹性模型，以提高预测精度。结合电磁效应和微极效应的先进多物理场模型进一步加深了对热–光弹行为的理解。Al-Hazaemh 等人¹⁷研究了旋转半导体介质中的光电磁热弹性激发，而 Nazir 和 Kumar¹⁸分析了微极热弹性相互作用。Song 等人¹⁹还考察了非耗散热弹性相互作用，Nasr 和 Abouelregal²⁰研究了含空腔半导体中的光吸收过程。非局部和分数阶模型在近期发展中发挥了关键作用。Gupta 等人²¹研究了非局部多孔介质中的光热激发，而 Hobiny 和 Abbas²²分析了半导体中分数阶波传播。Hafed 和 Zenkour²³研究了倾斜加载效应，Oliinyk 等人²⁴考察了非稳态热–光弹效应。Awwad 等人²⁵还探讨了功能梯度半导体在激光激发下的行为，而 Gupta 等人²⁶利用依赖记忆的模型研究了热–压电–光电耦合效应。Greene 和 Patterson²⁷以及 Barone 和 Patterson²⁸建立了结合热学与光学技术的经典实验方法，为应力分析提供了可靠手段。此外，Kaur 和 Singh²⁹为半导体谐振器开发了非局部依赖记忆的模型，Abbas 等人³⁰分析了具有变热导率的光热相互作用。Lu 等人³¹对增强复合膜进行了实验与数值研究，而 Purkait 和 Kanori³²研究了旋转纤维增强介质中的记忆响应。Abo-Dahab 等人³³进一步考察了载荷条件下纤维增强热弹性介质中的波反射。最近，先进的分数阶与非局部热弹性模型被提出用于描述复杂材料。Abouelregal 等人³⁴利用分数阶模型研究了生物组织中的热响应，Selvamani 等人^35,36研究了纳米梁中的非局部波传播与振动行为。此外，Abouelregal 等人³⁷将双相滞后和粘弹性模型应用于微结构，同时发展了具有记忆核的分数阶热弹性公式^38,39，以捕捉复杂的耦合现象。

尽管已有大量研究进展，但大多数现有研究主要集中在解析解或数值解上，未能提供物理场的详细时空可视化。在许多实际应用中，特别是在各向异性的纤维增强半导体介质中，系统响应强烈依赖于空间和时间的变化，因此可视化对于准确解释至关重要。此外，尽管数值和实验方法能够提供有价值的见解，但基于正则模态和特征值技术的解析方法在涉及半无限域和耦合多物理场系统的问题中具有显著优势。这些方法能够获得闭式解，从而更深入地理解波的传播、衰减及耦合机制，并为验证数值和实验结果提供可靠的基准。

本研究全面探讨了各向异性纤维增强半导体介质中热-光弹性耦合行为。本研究的创新之处在于，将基于特征值的解析方法与时空热图可视化相结合，以深入揭示物理场的演化过程及其局域化特性。本工作的目的在于分析各向异性和纤维增强对热、力学及载流子场之间相互作用的影响，并展示所提出方法在解释与现代工程应用（如半导体器件、传感器和激光技术）相关的复杂多物理场现象中的适用性。

方案

本研究完全基于理论建模和数值模拟，不涉及人类受试者、动物实验或生物样本。因此，无需伦理审批和知情同意。

纤维增强各向异性介质中光热弹性的数学表述
本研究考虑了受表面光激发作用的二维纤维增强各向异性半导体半空间。该介质占据以下区域 x ≥ 0，其中边界位于 x = 0 表示暴露的表面。坐标系的定义使得 x-轴延伸至培养基中，同时 y-轴沿表面方向排列，用于描述面内行为。由于存在定向排列的增强纤维，材料被假定为均匀但各向异性，这使得其弹性性能和耦合特性具有方向依赖性。表面处的光吸收产生局部加热和过剩载流子，导致热场、力学场与载流子场之间发生完全耦合的相互作用。因此，系统状态由温度描述 θ(x，y，t) (K），载流子密度 N (x, y, t (m^-3）和位移分量 u (x, y, t）和 v (x，y，t)(m），假设形变较小。物理域、坐标系、纤维取向及施加的光学激发示意图如图所示 图1所有符号和数值计算均使用 Wolfram Mathematica（版本 12.0）完成。

图1。半无限光纤增强型半导体介质在边界 x = 0 处受到光激发的示意图。 坐标系（x, y）已标出，其中纤维取向沿 x 方向（a = (1, 0)），展示了介质的几何构型及其方向相关的各向异性特征。请点击此处查看该图的放大版本。

纤维增强各向异性热弹性半导体介质中应力张量的本构关系以一般形式通过公式1^1,5给出。在此表述中，θ表示相对于参考温度T₀的温度增量，而T则表示适用情况下的绝对温度。

.   (1)

其中，C_ijkl 为弹性刚度系数，e_kl 为应变张量，而 β_ij 和 η_ij 分别表示热弹性张量和载流子耦合张量。当存在纤维增强时，材料响应呈现方向依赖性，并由纤维取向矢量 a = (a_i) 所决定，该矢量在弹性项和耦合项中均引入了各向异性贡献。因此，本构关系被扩展为显式地包含纤维增强的影响^2,3：

. (2)

其中，λ 和 μ_τ 为拉梅常数，μ_L 是沿纤维方向的纵向剪切模量。参数 α 表示纤维增强效应，不同于表示热膨胀系数的 α_ij。单位矢量定义了纤维取向，从而在应力-应变响应中引入了方向依赖性。在当前的二维模型中，假设纤维沿 x 轴方向排列，因此取向矢量明确表示为 a = (1, 0)。该设定为纤维方向提供了清晰的参数化表达，并确保各向异性贡献被一致地纳入控制方程中，直接反映了纤维增强引起的方向性行为。对于当前的二维构型，控制应力分量简化为：

，（3）

，（4）

. (5)

这些方程说明了各向异性、纤维增强以及多物理场耦合效应的综合作用。β_ij 和 η_ij 两个系数按材料参数定义如下：

，

，

，

.

此处，系数 A_ij 表示纤维增强各向异性介质的有效弹性常数，其定义如下：

. (6)

其中，λ、μ_L 和 μ_T 为各向异性纤维增强介质的弹性常数，而 α_ij 和 ξ_ij 分别表示热膨胀系数和载流子膨胀系数。热-光弹性半导体介质中弹性波的传播遵循动量守恒原理，该原理构成了动态热弹性分析的基础。在无体力作用的情况下，可变形连续体的运动通式基于^1,15表示如下：

. (7)

其中，ρ 为质量密度，σ_ij 为应力张量。在本研究中，模型被限定为 x - y 平面上的二维构型，位移场由 u(x, y, t) 和 v(x, y, t) 表示。根据热-光弹性介质中的标准公式，二维运动控制方程可写为如下形式：

，（8）

. (9)

将各向异性的纤维增强本构关系代入上述方程后，得到如下耦合的偏微分方程组（DEs）：

，（10）

. (11)

此处，下标表示对空间和时间变量的偏微分。这些方程突出了各向异性、纤维增强、温度梯度和载流子扩散对介质动态响应的耦合影响。在光激发存在的情况下，半导体内部的温度场受到载流子浓度与机械变形相互作用的显著影响，从而导致一个完全耦合的能量输运过程。与经典热传导不同，此类介质中的温度演化受到来自载流子复合和热弹性效应的附加源项所支配，这些源项显著改变了热传播特性。在广义热弹性理论框架下，热传导方程表达如下^16,20：

. (12)

此处，C_E 为恒定应变下的比热容，表示材料的热容量；T₀ 表示介质在平衡状态下的参考绝对温度。对于当前的二维构型，该方程简化为^16,20：

. (13)

该方程表明，温度场不仅受到方向性热导率的影响，还通过项 受到载流子复合的影响，并且还通过热弹性耦合项受到与时间相关的变形影响。该模型捕捉了各向异性纤维增强半导体中热传递所涉及的关键多物理场相互作用，突出了载流子动力学和力学响应在改变系统热行为中的作用。当半导体介质受到光激发时，由于入射辐射的吸收，会产生大量载流子。这些载流子经历包括空间扩散、复合以及热驱动产生在内的输运过程，而这些过程均与材料内部的温度场内在关联。因此，载流子浓度成为决定耦合热–光弹性响应的关键变量之一。

在当前的模型中，载流子浓度的演化通过扩散机制、衰减效应和热激活过程之间的平衡来描述，由此得到以下控制关系^1,5”

. (14)

此处，D_E 表示载流子扩散系数， 是 x - y 平面上的二维拉普拉斯算子。项 考虑了具有弛豫时间 τ 的复合效应，而 k 是热-载流子耦合系数，定义为 ，用于表征平衡载流子浓度 N₀ 对温度变化的敏感性。该关系突显了温度作为载流子生成驱动力的作用，并建立了各向异性纤维增强半导体介质中热场与电子场之间的直接耦合关系。

已建立耦合光热弹载流子系统的控制方程及数学表述。与硅（Si）介质相对应的物理和材料参数及其数值、单位和参考文献列于表1中。这些参数随后被用于数值计算和无量纲化过程。

符号	数值	单位	参考文献
λ	3.64 × 10¹⁰	N/m²	12
μT	5.46 × 10¹⁰	N/m²	12
μL	3.20 × 10¹⁰	N/m²	12
ρ	2330	kg/m³	13
CE	695	J/(kg·K)	30
k11	0.0921 × 10³	W/(m·K)	30
k22	0.0963 × 10³	W/(m·K)	30
DE	2.5 × 10⁻³	m²/s	22
τ	5 × 10⁻⁵	s	15
T₀	300	K	15
Eg	1.11 × 10⁻¹⁹	J	12
α11	3.1 × 10⁻⁶	K⁻¹	30
α22	3.5 × 10⁻⁶	K⁻¹	30
ξ11	−7 × 10⁻³¹	m³	21
ξ22	−9 × 10⁻³¹	m³	21
κ	2.16 × 10²¹	m⁻³·s⁻¹·K⁻¹	21
α	−1.28 × 10¹⁰	N/m²	28
β	220.90 × 10¹⁰	N/m²	28
ω	2.95 + 1i	s⁻¹	12
a	1	— (无量纲)	13
y	0.6	m	13
θ₀	1	— (无量纲)	15
N₀	1	— (无量纲)	15

表1. 各向异性纤维增强半导体介质数值分析中使用的材料属性和参数。所有量均以国际单位制（SI）表示，除非另有说明。无量纲参数已相应标明。所列数值对应于基于硅的材料属性及本计算中采用的模型参数，这些参数来自所引用的参考文献。热-载流子耦合系数 κ 定义为 κ = (∂N₀/∂T)(1/τ)，遵循热-光弹性半导体模型中的标准公式。

耦合各向异性光热弹模型的无量纲化表述
为简化控制方程并获得耦合热-光弹系统的统一无量纲表示形式，引入了适用于空间坐标 x,y、时间 t、位移分量 u, v、温度 T、载流子浓度 N 和应力 σ 的适当特征尺度。这些标度参数的选取基于介质的固有物理性质以及热、力学与载流子场之间的耦合机制，遵循文献中已建立的表述方法^16,21。据此，无量纲变量定义如下：

, , , , , , , , , 。

该变换减少了独立材料参数的数量，并提供了耦合系统的归一化表示。通过将上述无量纲变量代入先前推导的控制方程，系统被重写为无量纲形式。为简便起见，随后省略了与无量纲变量相关的上标符号。该过程得到了一组紧凑的无量纲偏微分方程，可写为如下形式：

，（15）

，（16）

，（17）

. (18)

应用无量纲变换后，该系统的应力分量被表示为以下归一化形式：

，（19）

，（20）

。（21）

引入了无量纲参数 a_i 以表示控制耦合各向异性光-热-弹性行为的物理和材料属性的紧凑组合。每个系数反映了系统内特定的相互作用机制，并揭示了基础物理过程的相对影响。 表示法向耦合刚度与主弹性刚度之间的比值，反映了两个位移分量之间各向异性相互作用的程度。 表征横向变形对法向应力分量的相对贡献。 衡量热弹性耦合的方向性变化，指示热膨胀效应中的各向异性。 描述载流子密度对诱发弹性变形的各向异性影响。 表示归一化的剪切刚度，量化剪切变形相对于法向变形的贡献。 考虑了控制位移方程中法向与剪切变形之间的联合耦合。 表达横向刚度与剪切刚度之间的比值，突出各向异性变形行为。 表示归一化的惯性参数，关联波传播效应与剪切刚度。 表征不同空间方向上位移梯度之间的耦合。 量化热效应在横向位移场中的相对贡献。 衡量载流子诱发变形相对于剪切刚度的影响。：表示不同空间方向上热导率的各向异性。 表征载流子复合对介质内热生成的影响。 表示热效应与时间依赖性弹性变形之间的耦合。 考虑两个空间方向上各向异性热膨胀的联合影响。 表示控制载流子输运速率的归一化扩散参数。 表征载流子复合效应的相对强度。 描述热变化与载流子生成过程之间的耦合。

利用法向模态技术的解析解
为了获得耦合各向异性热-光弹系统的解析解，采用了法向模态技术，因其能有效将控制偏微分方程组简化为更易处理的常微分方程组。该方法广泛应用于波传播现象的分析，包括色散和衰减。因此，假设场变量在时间和横向空间方向上均呈谐波变化^1,12,23。于是，位移分量、温度、载流子密度和应力可表示为如下指数形式：

. (22)

此处，ω 表示控制场量时间行为的复频率，而 a 表示与 y 方向空间变化相关的波数。选择这些参数是为了满足稳定性要求，并确保在半无限区域内获得物理上可接受的有界解。将上述假设形式代入先前推导的无量纲控制方程，并对所得表达式进行简化，可将原始的耦合偏微分方程组约化为关于空间坐标 的常微分方程组，其形式如下：

，（23）

，（24）

，（25）

. (26)

此外，变换域中相应的应力分量可表示如下：

，（27）

, (28)

. (29)

此处，D 表示微分算子 。这些方程代表了控制方程组在法向模态域中的简化形式，为后续推导特征方程并构建通解解析表达式提供了基础。系数定义如下：，，，，，，，，，。

矩阵微分方程的建立与特征值分析
在应用正则模态变换后，由方程23–26给出的控制方程组被简化为一组关于空间坐标的二阶常微分方程。为了便于系统求解，通过引入对应于各场变量一阶导数的辅助变量，将该方程组转化为等价的一阶微分方程组。具体地，定义了以下变量：

, . (30)

根据这些定义，将方程 23–26 重写为以下八个一阶微分方程组：

，（31）

，（32）

，（33）

，（34）

. (35)

上述系统以紧凑的矩阵 A 形式表示如下：

. (36)

状态矢量由以下给出：

. (37)

系统矩阵具有如下显式形式：

. (38)

该公式将原系统转化为一个特征值问题^1,15。特征方程由下式得到

. (39)

得到一个关于特征值的八阶多项式。以简化形式表示，特征多项式可写为

。（40）

其中，Z_i 的系数是系统参数的函数，其具体定义见下文。所得的特征值决定了解的空间行为，包括衰减和传播特性。仅保留满足 Re(m) > 0 的特征值，以确保解在物理上合理，并在 x → ∞ 时呈指数衰减。

。（41）

特征多项式的根定义了本征值 m，这些本征值决定了解的空间行为。这些本征值通过 Mathematica 数值计算得到：首先利用 CharacteristicPolynomial 函数 构造特征多项式，然后使用 NSolve 求解所得的代数方程。由于该问题定义在半无限域（x ≥ 0）上，因此仅考虑在 x → ∞ 时保持有界的物理上可接受的解。据此，仅保留满足 Re(m) > 0 的本征值，以确保当 x → ∞ 时解具有 exp(−mx) 形式的指数衰减行为。其余的根被舍弃，因为它们对应于非衰减或无界的解，不符合模型的物理要求。

对于每个保留的特征值 m，其对应的特征向量由相关的代数方程组求得

，（42）

并以以下形式表达：

. (43)

展开上述矩阵方程，得到以下线性方程组：

，（44）

，（45）

，（46）

，（47）

. (48)

由于本征值问题具有齐次性，本征向量仅能确定至任意乘法常数因子。为了获得唯一且一致的表示，通过固定本征向量中的一个分量施加了归一化条件。在本研究中，选定第一个分量使得 q₁ = 1，其余分量则根据上述方程组依次确定。从计算角度来看，该归一化通过将某一分量赋值为1，并求解由此得到的线性方程组以计算其余分量来实现。此方法为计算每个可接受本征值对应的本征向量提供了一种系统且可重复的途径。

。（49）

其余分量可根据系统关系相应得出。这些特征向量描述了每个模态中温度、载流子密度和位移场的相对贡献。因此，该问题的通解被构造为所有可接受本征模的线性组合，每个本征模均对应一个特征值及其相应的特征向量，从而为半空间介质中耦合各向异性光-热-弹性行为提供了完整的解析描述。该系统的通解因此可写为如下形式：

。（50）

其中，C_i 是由边界条件确定的常数。通过展开上述矢量表达式，得到场变量如下：

，（51）

，（52）

，（53）

. (54)

该图示表明，该解由指数模式的叠加构成，其中每个特征值-特征向量对独立地贡献于整体物理响应。选择允许的特征值，使其实部为正，从而确保解是有界的且具有物理意义。 x → ∞.

边界条件与物理约束
通过将通解代入 x = 0 处给定的边界条件，得到了一组关于常数 C_i 的线性代数方程组。具体而言，每个边界条件（温度、载流子密度和位移约束）均用本征模式展开式表示，从而得到一组关联系数 C_i 的方程。该过程导出一个可写成矩阵形式 BC = D 的线性方程组，其中 B 是由在边界处计算的本征向量分量构成的系数矩阵，C = (C₁, C₂, C₃, C₄)^T 是未知常数向量，其值由施加的边界值（如 θ₀, N₀ 及位移约束）确定。该线性方程组通过 Mathematica 进行数值求解，其中系数矩阵和右端项向量被显式组装，并利用 LinearSolve 例程求得未知常数。随后，将这些常数代回通解中，构建出物理场的完整表达式，并进一步用于结果的数值计算与图形化表示。

施加的边界条件如下所示：

温度限制：

. (55)

该条件表示由周期性光加热引起的谐波变化的表面温度。它作为主要的热激发源，驱动介质内的耦合热弹性与载流子输运过程。振幅 θ₀ 表征所施加热载荷的强度。

载流子密度限制：

。（56）

该边界条件描述了由光学照射引起的光生载流子密度，反映了由于光子吸收导致的电子激发及其与入射光场一致的谐波调制。

位移约束：

. (57)

该条件表明边界在横向受到机械约束，因此在表面沿 v 方向不发生位移。

剪切应力约束：

. （58）

该条件对应于剪切应力意义上的无牵引力边界。它确保表面上无切向力作用，这与切向上的力学自由边界一致。除了在 x = 0 处的边界条件外，还施加了无穷远处的物理条件：，以保证在半无限区域内物理解的有界性。在展示数值结果之前，本研究采用的整体计算流程已在图2中总结。计算中使用的激励参数 θ₀、N₀、复频率 ω 和波数 a 的数值列于表1中。表1所列参数包括用于归一化公式的无量纲参数以及有量纲材料常数。在数值计算中，空间域定义为，横向坐标固定为 y = 0.6，时间域取为。所有数值计算和图形表示均采用上述范围。

图 2．所提出方法的计算工作流程。 该图展示了从数学建模到数值结果的完整步骤序列：控制方程、无量纲化、法向模态技术的应用、转换为一阶系统、矩阵形式表达、特征值与特征向量分析、边界条件的施加、常数的确定，以及数值图的生成。请点击此处查看此图的放大版本。

结果

数值结果
本节通过数值计算分析了各向异性纤维增强半导体介质中耦合热-光弹-载流子系统的行为。所考虑的材料为硅（Si），其物理和材料参数列于表1中。这些材料常数被直接代入控制方程，并在数值计算中用于求解各场变量。所有参数均采用定义的无量纲化尺度进行一致化处理。数值计算在Mathematica中使用默认的精度和收敛设置完成。计算网格为均匀网格，空间和时间方向的步长均为0.01。

首先，将控制方程以无量纲形式建立，并用主要场变量表示。然后在 Mathematica 中采用模态分析法和特征值方法对所得的解析表达式进行符号和数值计算。根据表1设定材料常数，并代入推导出的表达式中。随后，应用边界条件以确定未知常数，并确保在半无限域内获得具有物理意义的解。

通过数值计算对所得表达式进行求解，以获得温度θ*、载流子浓度N*以及位移分量u*和v*的...

讨论

所获得的结果为各向异性纤维增强半导体介质中的热–光弹耦合行为提供了清晰的物理洞察。本研究提出了一种基于特征值的解析框架，用于研究此类介质中热载荷、载流子生成与弹性变形之间的相互作用。所观察到的响应本质上由这些物理过程之间的强耦合关系所决定。边界处对光能的吸收导致局部加热和载流子激发，进而引发介质内部的热弹性应力和位移。该耦合机制解释了温度、载流子浓度和力学场的同步演化过程，其中每个场都主动影响其他场，而非独立行为。

所有物理量的空间衰减是半无限介质内能量耗散和波阻尼的直接结果。这种衰减行为在物理上与所施加的边界条件一致，并反映了系统的稳定性，即扰动始终局限于激励区域附近。相比之下，由于热能与载流子动力学之间的持续相互作用，场的振幅随时间演化而增强，这些相互作用作为内部源随时间推移强化了系统的响应。这些发现通过提供对各向异性半导体系统中耦合多物理场相互作用的清晰解析理解，推动了该领域的发展，而这类系统往往难以仅通过纯数值方法进行解释。

纤维增强引入的各向异性在塑造介质的定向行为方面起着关键作用。与各向同性材料不同，其力学响应在不...

材料

本文使用的材料清单

姓名	公司	目录编号	评论
计算软件（符号与数值分析）	Wolfram Research	使用了 Wolfram Mathematica（版本 12.0）	用于特征值计算、解析解实现及数值评估
数据可视化工具（等高线与热图生成）	Wolfram Research	使用 Wolfram Mathematica（版本 12.0）生成二维等高线图和时空热图	用于生成二维等高线图和时空热图
材料参数数据集（硅半导体特性）	多种文献来源	不适用	计算中使用的物理常数（弹性、热学、载流子相关参数）（见表1）
个人计算机/工作站	HP	不适用	在运行 Windows 操作系统的标准个人计算机上进行计算，具备足够的内存以支持数值模拟
公式编辑器	Microsoft  Word 和 MathType	不适用	用于稿件中数学表达式的排版与呈现
参考文献管理软件	Elsevier	不适用	用于管理参考文献并格式化引文（温哥华格式）