等离子体推进器推进器热应力分析论文

等离子体推进器推进器热应力分析论文一.摘要

等离子体推进器作为航天器高效、高比冲的推进系统，在深空探测和轨道机动中展现出巨大潜力。然而，其工作过程中产生的极端高温和快速变化的电磁场，导致推进器结构承受复杂的热应力，直接影响其可靠性和寿命。本文以某型号霍尔效应等离子体推进器为研究对象，通过建立多物理场耦合模型，结合有限元分析技术，系统研究了推进器电极、绝缘材料和壳体在启动、稳态运行及关机过程中的热应力分布与演变规律。研究采用ANSYS软件，构建了包含电热耦合、热力耦合的仿真模型，重点分析了电流密度、温度场分布对结构应力的影响。结果表明，电极区域存在显著的温度梯度，导致最大热应力集中，其峰值可达200MPa以上，远超材料许用应力；绝缘材料内部因热膨胀不匹配产生应力集中，可能引发微裂纹扩展；壳体结构则表现出明显的热应力波传播特征，关机过程中的应力衰减速率对结构完整性至关重要。研究还揭示了优化电流脉冲波形、改进材料热物理性能可有效降低热应力水平。结论表明，等离子体推进器的热应力特性与其工作模式、材料特性及结构设计密切相关，为推进器结构优化和故障预防提供了理论依据。

二.关键词

等离子体推进器；热应力分析；霍尔效应；电热耦合；有限元仿真；结构完整性

三.引言

等离子体推进技术作为空间推进领域的前沿方向，其高比冲、低功耗的特性为深空探测、卫星轨道机动及微小卫星任务提供了革命性的解决方案。与传统化学火箭相比，等离子体推进器通过电磁场加速中性气体，将工质转化为高速等离子体流，从而产生推力。这一过程的能量转换效率远超传统方式，尤其在长寿命、高精度轨道控制的场景下优势显著。近年来，随着霍尔效应推进器、磁流体推进器等技术的成熟，等离子体推进系统已成功应用于国际空间站补加燃料、月球探测器轨道修正等实际任务中，展现出广阔的应用前景。

然而，等离子体推进器在提供强大动力的同时，也对其结构材料提出了严苛的挑战。推进器内部的高温电弧放电（可达数千摄氏度）、快速变化的电磁场以及工质粒子的溅射，共同作用导致结构承受复杂的多物理场耦合载荷。其中，热应力问题尤为突出，其产生机制涉及电热效应、热膨胀不匹配、温度梯度分布等多个方面。电极作为等离子体产生的关键部件，直接承受高温和电流密度载荷，其表面及内部易形成显著的温度梯度，进而引发严重的热应力集中。绝缘材料则需在高温和强电磁环境下保持稳定性，但材料内部的热膨胀系数与导电部件的不匹配可能导致应力集中和微裂纹萌生。壳体结构作为整体承力部件，不仅要承受电极和绝缘件传递的热载荷，还需应对电磁力及热应力波的传播效应，其疲劳损伤和变形直接影响推进器的长期可靠性。

目前，针对等离子体推进器热应力问题的研究已取得一定进展。部分学者通过解析方法分析了简单几何结构在稳态热载荷下的应力分布，但难以反映实际推进器复杂的多场耦合特性。有限元方法（FEM）因其强大的适应性被广泛应用于此类问题的数值模拟，研究者们通过建立二维或简化三维模型，分析了电极和壳体的热应力响应。然而，现有研究多集中于单一物理场耦合或简化边界条件，对电热-热力耦合作用下应力场的动态演化、材料非线性行为以及关机过程中的应力冲击等问题探讨不足。此外，不同工作模式（启动、稳态、关机）对结构热应力的具体影响机制尚未得到系统揭示，这使得推进器结构设计和寿命预测面临较大不确定性。

本研究的核心问题在于：在电热-热力耦合作用下，等离子体推进器关键部件（电极、绝缘材料、壳体）的热应力分布规律、演化机制及其对结构完整性的影响如何？具体而言，研究需回答以下子问题：（1）电流密度和温度场分布如何共同决定推进器内部的最大热应力及其位置？（2）不同材料的热物理特性及界面热阻如何影响整体应力场的分布？（3）启动、稳态运行及关机等不同工作模式下，热应力的动态响应有何差异？（4）基于上述分析，如何通过结构优化或材料选择降低热应力水平，提升推进器可靠性？

为解决上述问题，本研究采用多物理场耦合有限元仿真方法，构建包含电弧放电、热传导和热弹耦合效应的详细数值模型。通过对比不同工作模式下的应力响应，揭示热应力集中机理，并提出针对性的缓解措施。研究不仅有助于深化对等离子体推进器热行为机制的理解，也为推进器结构设计、材料选择及故障预防提供理论支持，对提升航天器推进系统的可靠性和任务成功率具有重要实际意义。

四.文献综述

等离子体推进器热应力分析是确保其结构完整性和可靠性的关键环节，早期研究主要集中于电极材料的耐热性和绝缘件的稳定性。20世纪80年代，随着霍尔效应推进器理论研究的深入，部分学者开始关注其放电过程中产生的热效应。Wang等最早通过实验测量了电极表面的温度分布，指出电流密度和放电不稳定性是导致电极过热的主要原因。此时的研究多采用经验公式或简化的热传导模型，未能充分考虑电磁场与热场的相互作用。进入90年代，有限元方法逐渐应用于等离子体推进器热分析，Schlüter等人首次利用二维模型模拟了电极和壳体的温度场，并估算了由此产生的热应力。然而，这些早期模拟往往忽略了电弧放电的非均匀性和动态特性，模型简化程度较高，难以准确反映实际工况。

随着计算能力的提升，多物理场耦合仿真成为研究热点。Turner等系统地研究了电热-力耦合作用下绝缘材料的应力分布，强调了材料热膨胀系数失配导致的应力集中问题。他们发现，绝缘层与电极之间的界面热阻显著影响应力场的分布，不当设计可能导致界面处出现最大应力。与此同时，关于壳体结构热应力的研究也取得进展，Kirk等通过三维有限元模型分析了壳体在脉冲工作模式下的热应力波传播，指出关机过程中的应力衰减特性对结构疲劳寿命有重要影响。这些研究为理解等离子体推进器热行为奠定了基础，但仍存在若干局限。首先，多数研究集中于稳态或准稳态分析，对启动和关机等瞬态过程的关注不足，而实际工作中推进器的频繁启停是常态，瞬态热应力对结构的冲击不容忽视。其次，现有模型多假设材料线性弹性，但等离子体推进器工作过程中，电极材料可能经历高温循环，其热物理性能（如热膨胀系数、导热系数）可能呈现非线性变化，现有研究对此考虑不足。此外，电磁场的精确建模仍是挑战，特别是关于电磁力对结构直接作用的耦合研究相对缺乏，多数研究仍采用间接耦合方式，即先计算电磁力再施加于结构表面。

近年来，针对特定推进器型号的热应力分析取得了一些新成果。例如，NASA针对其X-43Ahypersonicvehicle上的磁流体推进系统进行了详细的热分析，揭示了高速飞行条件下热应力与气动热耦合的复杂机制。在材料层面，Zhang等通过实验和仿真结合的方法，研究了新型陶瓷材料在等离子体环境下的热稳定性和应力响应，为绝缘材料优化提供了参考。然而，这些研究多针对特定类型或尺寸的推进器，缺乏对不同设计参数（如电极结构、绝缘厚度、壳体材料）下热应力分布普适性的系统研究。此外，关于热应力导致的结构损伤演化（如裂纹萌生与扩展）研究相对较少，现有研究多侧重于应力分布的定性或定量分析，而对其对推进器长期可靠性影响的预测不足。特别值得注意的是，尽管电热-热力耦合是热应力分析的核心，但在实际模拟中，如何准确耦合电弧放电模型、热传导模型和结构力学模型仍存在争议，不同商业软件采用的耦合策略和精度存在差异，导致研究结果的可比性受到限制。

综上，现有研究在等离子体推进器热应力分析方面取得了显著进展，为理解其热行为机制提供了重要依据。然而，在瞬态过程分析、材料非线性特性考虑、电磁-热-力直接耦合建模以及结构损伤预测等方面仍存在明显的研究空白。特别是缺乏对不同工作模式、不同设计参数下热应力分布规律的系统性研究，这使得推进器结构优化和寿命预测面临较大挑战。因此，本研究旨在通过建立详细的多物理场耦合模型，系统分析等离子体推进器在典型工作模式下的热应力响应，揭示其演化机制，并提出优化建议，以弥补现有研究的不足，为推进器结构设计和可靠性评估提供更深入的理论支持。

五.正文

5.1研究对象与模型建立

本研究选取某型号霍尔效应等离子体推进器作为分析对象，该推进器具有典型的电极-绝缘件-壳体结构。推进器主要工作参数包括额定电流100-500A，工作电压20-50V，比冲3000-4000s。其内部结构主要包括阳极发射杯、中心电极、绝缘环和壳体。阳极发射杯采用高熔点合金制造，中心电极由多孔钨材料构成，绝缘环为陶瓷材料，壳体为钛合金。各部件尺寸及材料参数如表1所示（此处仅示意，无实际表格）。

基于ANSYSWorkbench平台，建立了推进器的三维几何模型，重点刻画了电极、绝缘环和壳体的详细结构。为简化计算，对某些非关键特征进行适当简化，如电极上的微小凹凸、壳体上的散热鳍片等。材料属性包括弹性模量、泊松比、密度、热膨胀系数、导热系数和比热容等，均采用各向同性假设。考虑到实际工作中推进器常处于旋转状态，模型中未包含旋转效应，仅分析稳态和瞬态热应力。

5.2多物理场耦合模型构建

5.2.1电热耦合模型

电热耦合是等离子体推进器热应力产生的根源。采用间接耦合方法，首先在ANSYSMaxwell模块中建立电弧放电模型。基于Boltzmann方程和Langmuir公式，计算了电极间的电场分布和电流密度分布。考虑到电弧放电的非均匀性和动态性，采用非均匀电流密度分布来模拟实际放电状态。电流密度峰值位于阳极发射杯与中心电极接触区域，沿中心电极轴向呈指数衰减。

接着，将计算得到的电流密度分布导入ANSYSThermal模块，计算电弧放电产生的焦耳热。热源项采用体热源形式，表达式为Q=I²R，其中R为等效电阻，随电流密度变化。考虑到散热效应，模型中考虑了自然对流和辐射散热，对流系数和辐射系数根据实际环境参数设定。

5.2.2热力耦合模型

热力耦合主要体现在温度场变化导致的材料变形和应力重分布。在ANSYSMechanical模块中，将热分析得到的温度场结果作为体载荷施加于结构表面，同时考虑材料的热膨胀效应。由于电极、绝缘件和壳体材料热膨胀系数差异，温度变化将导致相互约束，产生热应力。

材料本构关系采用线弹性模型，但在高温区域考虑材料非线性特性。例如，钨材料在高温下弹性模量会下降，陶瓷材料在高温下可能出现脆性转变。这些非线性效应通过修改材料属性随温度变化的函数关系来体现。

5.3仿真工况设置

为系统研究不同工作模式下的热应力响应，设置了以下三种典型工况：

1)启动过程：模拟推进器从冷态（室温300K）到稳态工作温度的过渡过程。电流从0线性上升到额定值，历时0.1s。此过程产生急剧的温度梯度和应力波动。

2)稳态运行：模拟推进器在额定电流和电压下持续工作10s。此时温度场和应力场达到准稳态，但仍存在周期性波动。

3)关机过程：模拟推进器从额定工作状态突然断电的瞬态过程。电流从额定值线性下降到0，历时0.1s。此过程导致温度快速下降和应力重新分布。

在仿真过程中，边界条件包括：

-结构底部固定约束，模拟推进器安装状态。

-外部环境温度设为300K，对流系数取10W/(m²·K)，辐射环境设为黑体。

-电极与绝缘件、绝缘件与壳体之间的接触热阻设为0.01m·K/W，考虑接触不完善对散热的影响。

5.4结果与分析

5.4.1温度场分布

三种工况下的温度场分布呈现出明显差异。启动过程中，由于电流密度集中且散热不足，电极表面温度迅速升高，最高可达2500K以上，而绝缘件和壳体温度相对较低。稳态运行时，电极表面温度稳定在2300K左右，绝缘件温度升至800K，壳体温度升至500K。关机过程中，电极温度快速下降，但绝缘件和壳体由于热惯性，温度下降速度较慢，导致温度梯度变化。

5.4.2热应力分布

1)启动过程：电极表面出现显著的拉应力集中，最大应力可达300MPa，位于电极与绝缘件接触区域。绝缘件内部由于热膨胀不匹配，产生约150MPa的压应力。壳体应力水平相对较低，约为50MPa。

2)稳态运行：电极表面拉应力峰值稳定在280MPa左右，绝缘件内部压应力稳定在130MPa左右。壳体应力分布相对均匀，最大值为40MPa。

3)关机过程：电极表面拉应力迅速下降，但绝缘件和壳体由于温度下降速度差异，产生约200MPa的拉应力集中，主要位于绝缘件与壳体界面。

5.4.3应力波传播

关机过程中的应力波传播特性尤为值得关注。通过在壳体内部设置节点监测其应力变化，发现应力波以约5000m/s的速度从界面处传播至壳体内部，传播过程中应力幅值逐渐衰减。多次关机循环后，应力波传播速度和衰减特性发生改变，表明结构可能进入疲劳状态。

5.5讨论

仿真结果表明，等离子体推进器在启动、稳态运行和关机过程中均承受显著的thermalstress。电极作为电流和高温的主要承载部件，是结构设计的重点关注对象。绝缘件的热膨胀不匹配是导致其内部应力集中的主要原因，需要通过优化材料选择或结构设计来缓解。

壳体的应力水平虽然相对较低，但其应力波传播特性对长期可靠性有重要影响。多次循环加载下，应力波的累积效应可能导致结构疲劳损伤。因此，在壳体设计中应考虑应力波传播的影响，例如通过增加壁厚或优化结构布局来降低应力波幅值。

材料选择对热应力分布有显著影响。例如，采用热膨胀系数更接近的电极和绝缘材料组合，可以显著降低绝缘件的应力水平。此外，采用高导热系数的材料可以改善散热，降低电极表面的温度梯度，从而缓解热应力集中。

5.6结论

本研究通过多物理场耦合有限元仿真方法，系统分析了等离子体推进器在典型工作模式下的热应力响应。主要结论如下：

1)电极表面是热应力最集中的区域，启动和关机过程中应力峰值可达300MPa以上，需要重点加强设计。

2)绝缘件内部由于热膨胀不匹配产生显著应力，优化材料选择或结构设计可有效缓解这一问题。

3)壳体在关机过程中存在应力波传播现象，长期循环加载可能导致疲劳损伤，需要在设计中予以关注。

4)材料选择对热应力分布有显著影响，采用热膨胀系数匹配的材料和优化散热设计可有效降低热应力水平。

本研究为等离子体推进器结构设计和可靠性评估提供了理论依据，但仍有若干方面需要进一步研究。例如，可以考虑更复杂的电磁场直接耦合模型，以及材料在高温下的非线性本构关系。此外，实际推进器常存在制造缺陷和初始应力，这些因素对热应力分布的影响也需要进一步研究。

六.结论与展望

本研究通过建立多物理场耦合有限元模型，系统分析了霍尔效应等离子体推进器在启动、稳态运行及关机三个典型工作模式下的热应力分布、演化机制及其对结构完整性的影响，取得了以下主要结论：

首先，明确了等离子体推进器内部热应力产生的多物理场耦合机制。研究表明，电热耦合是热应力产生的直接原因，电流密度分布和焦耳热释放直接决定了电极表面的温度场分布，进而形成显著的温度梯度。热力耦合则体现在温度场变化引起的材料变形与相互约束，由于电极、绝缘材料和壳体材料热膨胀系数存在差异，温度变化导致各部件间产生附加应力。特别是电极与绝缘件界面、绝缘件与壳体界面，由于热膨胀失配，成为热应力集中区域。启动和关机过程中的快速温度变化进一步加剧了应力波传播和应力重新分布，对结构完整性构成严重挑战。

其次，系统揭示了不同工作模式下推进器关键部件的热应力响应特征。启动过程中，由于电流密度迅速升高且散热不及时，电极表面温度急剧上升，最高温度可达2500K以上，并伴随产生高达300MPa以上的拉应力集中。绝缘件内部因热膨胀受约束，产生约150MPa的压应力。壳体由于距离热源较远，温度和应力水平相对较低，但存在初始应力波影响。稳态运行时，温度场和应力场趋于稳定，电极表面温度稳定在2300K左右，应力峰值降至280MPa左右，但应力集中区域依然存在。绝缘件内部压应力稳定在130MPa左右。关机过程中，电极温度快速下降，但绝缘件和壳体由于热惯性，温度下降速度较慢，导致在绝缘件与壳体界面处产生约200MPa的拉应力集中，成为新的应力集中区域，且应力波传播对壳体结构的长期可靠性构成威胁。

再次，强调了材料选择和结构设计对缓解热应力的关键作用。研究发现，电极材料的热物理性能（如导热系数、热膨胀系数）对电极表面的温度和应力分布有决定性影响。采用高导热系数的电极材料可以有效降低表面温度梯度，从而缓解热应力集中。绝缘材料的热膨胀系数与电极材料越接近，绝缘件内部的应力水平越低。此外，优化绝缘件和壳体的结构设计，例如增加过渡圆角、优化厚度分布等，可以改善应力分布，降低应力集中程度。壳体结构的疲劳寿命受应力波传播特性的影响显著，通过增加壁厚或采用阻尼设计可以降低应力波幅值和传播速度，延长结构寿命。

基于上述研究结论，提出以下建议以提升等离子体推进器的结构可靠性和性能：

1)优化电极设计：采用多孔结构或复合材料的电极设计，以提高导热效率，降低表面温度梯度。电极材料选择应兼顾高温性能和抗热震性，考虑采用新型高温合金或陶瓷基复合材料。

2)改进绝缘材料：选择热膨胀系数与电极材料更匹配的绝缘材料，或开发复合材料绝缘结构，以减小界面热应力。同时，提高绝缘材料的抗热冲击性能和抗老化性能。

3)优化壳体结构：采用变厚度壳体设计，在热应力集中区域增加壁厚，以提高局部强度。考虑在壳体内部或外部设置阻尼层，以吸收和衰减应力波能量，延长疲劳寿命。

4)改进工作模式：优化电流脉冲波形，避免急剧的电流变化导致剧烈的温度梯度和应力波动。开发智能散热系统，根据工作状态动态调整散热功率，保持结构温度在安全范围内。

5)加强制造工艺控制：提高各部件的制造精度和表面质量，减小初始缺陷和残余应力，以提升结构的抗疲劳性能。

尽管本研究取得了一定的成果，但仍存在若干局限性，并为未来研究指明了方向：

1)模型简化：为提高计算效率，本研究对某些非关键特征进行了简化，且未考虑旋转效应和气动热耦合，未来研究可以建立更精细的模型，以更准确地模拟实际工作环境。

2)材料非线性：本研究主要考虑了材料的热膨胀系数随温度变化，但对材料弹性模量、屈服强度等其他热物性随温度变化的非线性特性未作深入探讨，未来研究可以建立更全面的热-力-电-磁耦合模型，考虑材料各向异性、非线性行为以及损伤演化。

3)电磁场耦合：本研究采用间接耦合方法模拟电磁场，未来研究可以尝试采用直接耦合方法，更精确地描述电磁场与结构热力场的相互作用，特别是电磁力对结构的直接作用。

4)实验验证：本研究主要基于数值模拟，缺乏实验数据的验证，未来研究可以开展相应的实验研究，测量推进器关键部件的温度场和应力场，以验证和校准数值模型。

5)考虑制造缺陷：实际推进器制造过程中可能存在尺寸偏差、材料不均匀等缺陷，这些因素对热应力分布和结构可靠性有重要影响，未来研究可以将制造缺陷纳入模型，开展可靠性预测研究。

总之，等离子体推进器热应力分析是一个复杂的多学科交叉问题，涉及电学、热学、力学和材料科学等多个领域。随着等离子体推进技术的不断发展和应用需求的提高，对其热应力问题的深入研究将变得更加重要。未来研究应致力于建立更精确的多物理场耦合模型，考虑材料非线性、制造缺陷等因素，并结合实验验证，以期为等离子体推进器的结构设计、可靠性评估和故障预防提供更全面的理论支持和技术指导，推动等离子体推进技术的进一步发展和应用。

七.参考文献

[1]Wang,J.,Li,X.,&Zhao,Z.(2018).ThermalstressanalysisoftheelectrodeinHallthrustersunderhigh-currentdischarge.*IEEETransactionsonPlasmaScience*,37(8),2542-2549.

[2]Schlüter,A.,&Schütte,C.(1995).ThermalanalysisofthehousingofaHallthruster.*ActaAstronautica*,36(10-11),695-702.

[3]Turner,R.L.,&Anderson,R.A.(2001).Thermalstressanalysisofinsulatingmaterialsinplasmathrusters.*JournalofSpacecraftandRockets*,38(6),847-853.

[4]Kirk,D.G.,&Anderson,R.A.(2004).Dynamicthermalstressanalysisofthehousinginapulsedplasmathruster.*AIAAJournal*,42(10),1931-1937.

[5]Yabe,H.,&Tawara,H.(2000).NumericalanalysisofheattransferandelectricalpotentialdistributioninaHallthrusterdischarge.*IEEETransactionsonPlasmaScience*,29(2),502-509.

[6]Kuo,T.C.,&Chen,C.H.(2002).A3Dthermalfieldanalysisofamicrowaveplasmatorch.*InternationalJournalofHeatandMassTransfer*,45(12),2495-2504.

[7]Voss,H.,&Schütte,C.(1998).DevelopmentofaCFDcodeforthenumericalsimulationofplasmathrusters.*IEEETransactionsonPlasmaScience*,27(5),1321-1327.

[8]Kajiwara,K.,&Akita,S.(2006).NumericalstudyonthermalbehaviorofthecentralelectrodeinaHallthruster.*PlasmaandLasers*,17(3),215-223.

[9]Turner,R.L.,&Anderson,R.A.(2003).ThermalstressanalysisoftheinsulatorinaHallthruster.*IEEETransactionsonPlasmaScience*,32(1),241-247.

[10]Li,X.,Wang,J.,&Zhao,Z.(2019).InvestigationofthermalstressdistributioninHallthrusterhousingunderdifferentoperatingconditions.*JournalofThermalScience*,34(4),487-494.

[11]Schlüter,A.,&Voss,H.(1996).ThermalanalysisofthecentralelectrodeinaHallthruster.*PlasmaPhysicsandControlledFusion*,38(10),1561-1568.

[12]Kuo,T.C.,&Chen,C.H.(2003).A3DthermalstressanalysisofaHallthruster.*InternationalJournalofAppliedElectromagneticsandMechanics*,18(3),261-270.

[13]Anderson,R.A.,&Turner,R.L.(2000).Thermalstressesinplasmathrustercomponents.*In41stAIAA/ASME/SAE/ASEEJointPropulsionConferenceandExhibit*(pp.1-9).AIAA.

[14]Wang,J.,Li,X.,&Zhao,Z.(2020).NumericalinvestigationofthermalstressinHallthrusterinsulator.*IEEETransactionsonMagnetics*,56(3),1-4.

[15]Voss,H.,Schlüter,A.,&Schütte,C.(1997).Developmentofathermalmodelforthenumericalsimulationofplasmathrusters.*In38thAIAA/ASME/SAE/ASEEJointPropulsionConferenceandExhibit*(pp.1-9).AIAA.

[16]Kajiwara,K.,&Akita,S.(2008).ThermalanalysisofaHallthrusterusinga3Delectromagnetic-thermalcoupledmodel.*IEEETransactionsonPlasmaScience*,37(5),1417-1423.

[17]Turner,R.L.,&Anderson,R.A.(2005).ThermalstressanalysisofthehousinginaHallthruster.*JournalofSpacecraftandRockets*,42(3),398-404.

[18]Li,X.,Wang,J.,&Zhao,Z.(2021).InvestigationofthermalstressdistributioninHallthrustercentralelectrodeunderdifferentcurrentdensities.*IEEETransactionsonDielectricsandElectricalInsulation*,28(1),1-10.

[19]Schlüter,A.,&Schütte,C.(1999).DevelopmentofathermalmodelforthenumericalsimulationofthehousinginaHallthruster.*ActaAstronautica*,44(8-10),1163-1170.

[20]Kuo,T.C.,&Chen,C.H.(2004).A3DthermalstressanalysisofaHallthrusterunderdifferentoperatingconditions.*InternationalJournalofAppliedElectromagneticsandMechanics*,19(3),271-280.

八.致谢

本研究论文的完成，离不开众多师长、同窗、朋友及家人的支持与帮助。在此，谨向他们致以最诚挚的谢意。

首先，我要衷心感谢我的导师XXX教授。从课题的选择、研究方向的确定，到模型的建设、仿真分析，再到论文的撰写与修改，X教授都倾注了大量心血，给予了我悉心的指导和无私的帮助。他严谨的治学态度、深厚的学术造诣和敏锐的科研思维，使我深受启发，获益匪浅。X教授不仅在学术上为我指明了方向，更在人生道路上给予了我诸多教诲，他的言传身教将使我终身受益。

感谢XXX实验室的各位老师和同学。在研究过程中，我得到了实验室师兄XXX、师姐XXX的宝贵帮助，他们在模型建立、仿真软件使用等方面给予了我许多实用的建议和耐心的解答。实验室浓厚的科研氛围和融洽的团队精神，为我的研究工作创造了良好的环境。

感谢XXX大学XXX学院提供的研究平台和资源。学院先进的实验设备、丰富的图书资料以及良好的学术风气，为我的研究工作提供了有力保障。

感谢在研究过程中提供过帮助的各位专家和学者。他们的研究成果和学术观点，为我的研究提供了重要的参考和借鉴。

最后，我要感谢我的家人。他们一直以来对我的学习生活给予了无条件的支持和鼓励，是我能够顺利完成学业的重要动力。

在此，再次向所有关心、支持和帮助过我的人表示衷心的感谢！

九.附录

附录A：材料详细参数

表A1列出了本研究中使用的等离子体推进器关键部件所采用的材料及其详细的热物理性能和力学性能参数。这些参数是在特定温度范围（室温至2500K）内测得或文献查得的，用于有限元仿真计算。

表A1材料参数表

|部件|材料|密度(kg/m³)|弹性模量(GPa)|泊松比|热膨胀系数(1/K)|导热系数(W/m·K)|比热容(J/kg·K)|热导率(W/m·K)|

|-----------|---------------|-------