编者注：这篇博客最初发布于 2016 年 7 月 4 日。现已更新，以反映软件的后处理新功能。

2 种 RF 仿真方法简介

在今天的博客中，我们将讨论 2 种仿真方法：渐近波形估计（AWE）法和频域模态（FDM）方法。这两种方法都旨在帮助用户解决采用超精细频率分辨率，或通过常规的 频域 研究运行超宽带仿真时，求解时间过长这一常见问题。当涉及描述具有单一谐振或无谐振的平滑频率响应时，AWE 方法非常有效。而 FDM 方法则适用于快速分析多级滤波器，或者目标通带内有多个谐振的大量元件的滤波器。接下来，我们将讨论这两种方法的典型设置和应用场景。

需要说明的是，AWE 和 FDM 方法几乎都不依赖于所选的频率步长。您可以自由地减小频率步长的值，获得分辨率良好的结果绘图，而不会出现明显的速度减慢或额外的内存消耗。不过，这种做法也存在缺点：降低频率步长值可能会影响最终解中保存的数据量。在本文末尾的数据管理部分，我们将给出能显著减小输出文件大小的建议。

请注意，在使用精细分辨率进行 AWE 或 FDM 计算之前，最好先使用较粗的频率分辨率运行一个初步的 特征频率 和常规的 频域 仿真。这可以帮助您快速地估算谐振位置，总体了解系统的频率趋势，包括实际通带和所需的频率分辨率。

AWE 方法促进降阶模拟

AWE 是一种先进的降阶模拟技术，由于其数值特征和数学算法技术性太强，我们在此不过多赘述，只演示如何在 RF 模块中使用此方法。 自 COMSOL Multiphysics 6.2 版本开始，软件新增了有一个专门的 自适应频率扫描 研究步骤，可以实现 AWE 方法。使用此功能时，需要指定目标输出的频率范围，并选择一个表达式用于对AWE 算法进行误差估计。该方法求解器可执行快速频率自适应扫描，默认情况下，使用 Padé 近似。

自适应频率扫描 研究设置。使用查看默认的 渐近波形估计（ AWE） 表达式。

在模拟谐振电路，尤其是包含许多频率点的带通滤波器类器件时，AWE方法非常有用。例如，COMSOL 案例库中的消失模圆柱腔滤波器模型先运行了一个常规的 频域研究，以 5 MHz 的频率步长对 3.45 GHz 到 3.61 GHz 之间的仿真频率进行初始扫频。

消失模圆柱腔滤波器教程模型 (左) 及其离散频率扫描结果 (右)。谐振频率附近的 S 参数图看起来并不平滑。

假设以更高的频率分辨率再次运行仿真，例如使用100 kHz的频率步长进行扫描，即分辨率提高 50 倍。可以预计，完成仿真所需的时长将提升50倍。但是，在此特定的示例中，使用 自适应频率 扫描研究完成仿真所需的时长几乎与频率扫描仿真所需的相同，并且能获得以100 kHz 频率步长计算的所有因变量解。

求解时间可能会在一定程度上受用户输入的 AWE 表达式的影响。任何模型变量都可以作为 AWE 表达式，只要能生成一个平滑的结果图，如高斯脉冲或平滑的曲线作为频率的函数，但最明显和典型的选择是基于s参数的全局表达式。例如，对于双端口带通滤波器，将S21 的绝对值 （abs(comp1.emw.S21)） 作为 AWE 表达式的输入可以获得非常好的结果。如果 AWE 表达式的频率响应包含无限梯度，例如在单个频率点具有良好阻抗匹配的天线的S11值，完成仿真则需要更长的时间。如果天线损耗可以忽略不计，使用 sqrt(1-abs(comp1.emw.S11)^2) 表达式可能结果更好，且能缩短计算时间。上述表达式是默认的 物理场控制 选项的 渐近波形评估 (AWE) 表达式。作为合理的检查方式，我们始终可以先以较粗的分辨率运行一个频域 扫描，绘制表达式，并选择最平滑的结果。

开始运行 自适应频率 扫描前，需要在研究中设置所需的更精细的频率步长。当仿真完成后，模拟所需时长几乎与离散频率扫描的相同。对比求解出的S参数。由于AWE求解器运行的频率扫描分辨率提高了50倍，因此其频率响应（S参数）结果绘图看起来也更平滑。使用这种方法，不仅可以节省宝贵的时间，还可以获得准确且出色的结果，而且谐振频率的定位也更加准确，如下图所示。如果您对此感兴趣，可以使用相同的分辨率运行一次常规扫描来进行验证，并检查结果是否一致。

分别运行 自适应频率 (AWE) 和离散 频域仿真 计算出的 S 参数图。使用 AWE 方法模拟的结果绘图的分辨率提高了 50 倍。

使用频域模态方法模拟电路谐振

无源电路的带通频率响应由多个谐振组合而成，因此 FDM 方法是提升其仿真效率的最佳选择。此方法通常包含两个步骤。特征频率分析是模拟任意形状器件共振频率的关键。从特征频率分析中获得所有必要信息后，就可以在频域模态研究中重复使用。当需要更精细的频率分辨率来更准确地描述频率响应时，这样做能够提升仿真效率，如AWE方法所示。

为了无缝执行 FDM 分析，有几个方面需要注意。一方面，需要过滤掉 特征频率 解中可能存在的所有不需要的非物理低频结果。另一方面，需要考虑目标频率范围内所有可能影响器件性能的物理模式，以获得正确的结果。要满足这两个要求，需要对 特征频率 研究设置进行一些调整（如下方屏幕截图所示）。首先，建议选择 更大的实部 作为 围绕偏移量的搜索方法 的设置。其次，对于 搜索特征频率 设置，最低通带频率可作为一个参考值。最后，必须调整 所需特征频率数 设置（例如根据初步测试），以包含必要的模式数量。

为模型添加包含两个步骤的 频域模态 研究，图中突出显示了 特征频率 设置。

让我们来看一个 耦合线带通滤波器教程模型，尝试使用 FDM 方法进行仿真分析。先以 50 MHz 的频率步长对 3.00 GHz 和 4.20 GHz 之间的仿真频率进行扫描，运行一个常规的 频域研究 。

耦合线带通滤波器教程模型（左）及以 50 MHz 频率分辨率进行离散频率扫描的仿真结果（右）。整个目标通带内的 S 参数图看起来并不平滑。

接下来，运行 频域模态 研究，并按照上文所述设置每个研究步骤。使用分辨率提高了50 倍的频率步长运行研究，并检查结果是否有更优。与 AWE 方法一样，FDM 方法研究得到的 S 参数图看起来更平滑，包含的信息量更大。例如，它显示了最初缺失的所有 S11 参数波纹。如果您对此方法的验证感兴趣，可以运行相同分辨率的常规扫描，并检查结果是否一致。

请注意，特征频率分析包含一个集总端口，作为额外的载荷因子影响仿真，因此计算出的 S 参数的相位与常规频率模型计算的相位不同。计算结果仅和不包含相位的 S 参数值兼容，如以dB为单位的值、绝对值、反射率或透射率。

分别由 频域，模态（FDM）和离散 频域 仿真计算出的 S 参数图。使用 FDM 方法计算的结果精度提高了 50 倍。

这与本文的主题并无直接关系，但在最后一张图中，可以看到到特殊的图形标记，它可以突出显示 S11 参数图的所有局部极小值以及 S21 参数图的通带。COMSOL Multiphysics 最近对结果评估功能进行了另一项改进，即从图表中提取交互式结果，将结果的信息性和交互式价值提升到了一个新的水平。

使用精细频率分辨率时的数据管理

如前所述， AWE 或 FDM 方法对如何细化频率扫描并没有真正的限制。不过，如果分辨率真的很高，解中将包含大量数据，从而导致在保存模型文件时，文件大小将大幅增加。在大多数无源RF和微波器件设计中，我们通常只对 S 参数感兴趣，因此在这种情况下，没有必要存储所有的场解。在研究的 在输出中存储 部分选择适当的选项，就可以控制模型保存计算的解。例如，可以只添加包含 S 参数计算边界的一个或多个选择。指定这些边界为 端口 或 集总端口，与整个仿真域相比，这些选择通常很小，因此总文件大小可以大幅减小。

请注意，在设置端口时，可以在指定选区后单击 边界选择 部分中的 创建选择 图标添加这种显式选择。然后，可以在相关研究步骤的 在输出中存储 部分添加由端口创建的所需的显式选择。

选择了2个 集总 端口的频域，模态 研究步骤的 在输出中存储 部分。您可以在 图形 窗口中查看这些选择的位置。

案例库中的教程模型

这篇博客介绍了COMSOL 软件中的 2 个强大的仿真功能，可以帮助您更快、更高效地模拟无源 RF 和微波器件。下列为 COMSOL 案例库中的一些相关教程模型，您可以浏览下载，进一步了解这些方法的使用：

• AWE 方法
  • 消失模圆柱腔滤波器
  • 电路量子电动力学 CPW 谐振器
  • 微带贴片天线
  • 采用SMA 连接的Wilkinson功分器
• FDM 方法
  • 级联矩形腔滤波器
  • 耦合线带通滤波器
  • 共面波导带通滤波器
  • 波导虹膜带通滤波器

文中演示的方法和研究是 通用的，不仅适用于RF仿真，在涉及如声学、力学、微机电系统(MEMS)和波动光学计算中，这些方法也尤其有用。

下一步

了解 COMSOL 中可用于RF和微波仿真的其他专用功能:

本文是博客有限元模型中的奇异现象：如何处理模型中的红点的后续内容，该博客介绍了结构力学模型中出现奇异应力的时间和原因，并对奇异现象进行了一般性介绍。如果您是第一次了解这个主题，建议先阅读这篇博客。有关如何处理奇异应力场的详细信息，请阅读本文。

进一步了解奇异应力场

首先，我们来详细分析一下奇异应力场及其与应力集中的关系。二者的相似之处是应力集中都出现在几何不连续处。应力集中与奇点的区别在于：前者的最大应力是有边界的，可以通过在有限元模型中使用足够精细的网格获得精确解。

通常，机械设计人员会通过引入一个半径尽可能大的圆角来减少应力集中。应力集中处的峰值应力通常用应力集中系数 与适当选定的名义应力的乘积描述。对于圆角，有时可以通过下列表达式获得 ：

式中, 是圆角半径， 是圆角缺口处的特征长度。

该方程的背景是求解一个大平板中的椭圆孔处应力集中的解析解，其中 是椭圆较大的半轴的长度。

含一个椭圆孔的大平板。

对于大多数缺口，该表达式只能用于粗略计算 ,因为很难推导出特征长度。但事实上，对于小缺口，峰值应力基本上是随圆角半径平方根的倒数变化。相信任何尝试过减小局部应力集中的工程师都可能为这一事实而苦恼过，因为适度地增大圆角半径会使峰值应力相应地减小。

极限应力集中发生在缺口半径无限小的裂纹尖端处。众所周知，在弹性固体中，裂纹尖端附近的应力场和应变场的解，与到裂纹尖端的距离 的平方根成反比。应力场通常用下式表示

式中，， 和 分别是模式 I（开口）、模式 II（剪切）和模式 III（撕裂）的应力强度因子。函数 , , 和 由裂纹尖端周围的极角 的三角函数组成。(更详细的定义请参见此处）。

由此得出结论，只要到裂纹尖端的距离足够近，裂纹尖端周围的应力场看起来都是一样的，与裂纹的实际形状以及存在裂纹的成分无关。

在线性弹性断裂力学的假设条件下，模式 的断裂准则为 ，其中 是材料参数（称为断裂韧性）。这样，就可以在不明确使用无限应力的情况下研究这种具有特殊奇异性的几何形状。下文，我们将对这一思路进行推广应用。

现在，考虑一种几何形状几乎是奇点的情况：一个圆角或一个圆角半径很小的裂纹，本文将重点讨论这种情况。在距离较远处，我们无法真正区分缺口和奇点。接下来，我们将使用一个例子来解释这句话。

以一个处于拉伸状态的含缺口的长条几何的二维模型为例。通过在这个模型中沿左侧垂直边添加对称条件，同一模型也可用于研究内部缝隙。

含缺口的长条几何结构的应力分布。该模型的参数与缺口深度 () 和缺口半径 () 有关。

对于尖锐的裂纹，这种几何形状的应力强度因子可写成

是裂纹长度； 是外加应力（此处为 1 Pa）； 是长条宽度。函数 有多种表示方法。在此，我们使用以下表达式

本文将这个表达式称作裂纹解。

沿韧带（从切口尖端向 x 方向延伸）的应力分布图，适用于短切口和几种不同的切口半径。由于对称性，只有一个应力分量 c 不为零。

不同缺口半径下沿韧带的垂直应力与到缺口尖端的距离的函数关系。虚线表示相同深度的裂纹的理论值。

一个有趣的现象，在特定情况下，应力场与裂纹解析解的应力场非常相似，即应力-距离对数图中的直线。在靠近缺口的地方，应力是有边界的，因为它是一个缺口，而不是裂纹。峰值应力与 成正比。

在距离尖端较远处，裂纹的局部应力场解在任何情况下都是无效的，不管它是裂纹还是缺口。但在非常近和非常远之间的区域，无论是从观察的角度，还是从物理学和数学的角度来看，都无法真正推断出缺口尖端的真实形状。

为什么这一点很重要？如果知道缺口的形状，那么只要观察一定距离外的应力，就可以确定那里的应力。稍后我们将详细探讨这一想法。

下一步，我们将在同一个图表中绘制大量具有不同缺口半径和切口长度的应力图。现在，通过缺口半径 对横轴进行归一化处理。

不同缺口深度和半径下，沿韧带的垂直应力与到缺口尖端的距离的函数关系。

从图中可以看出，在到缺口尖端的距离小于尖端半径 0.7 时，进入恒定斜率区域。从我们的角度来看，这已经相当接近要求解的问题了。那么这个区域会延伸多长呢？这不是由缺口细节控制的，而是由几何尺寸控制的。通过另一个归一化曲线图，韧带长度 ()，可以得知这一信息。

与上图相同，但通过韧带长度对距离进行归一化处理。

因此得出以下结论，这种情况下的恒定斜率区域延伸到韧带的 10% 左右。再远一些，应力场就不再受裂纹解的控制，而是受更多全局属性控制。对于特定的几何，这一区域的大小取决于该几何特有的长度尺度。

接下来，我们来研究是否可以用裂纹解中的应力场预测缺口尖端的峰值应力。先回到大平板上的椭圆孔。椭圆孔（宽度为，缺口半径为 ）的峰值应力与裂纹（长度为 ）的应力强度因子之比为

假定 ，则峰值应力可用应力强度系数表示为

这样，当计算出应力强度因子时，就可以用以下表达式确定圆形裂纹尖端的应力了。

其中，系数 是一个与配置相关的数量级为1的数字。我们可以在上面的例子中尝试这一假设。

下图显示的表达式

为缺口半径与缺口深度的函数关系。使用两种不同的几何形状进行计算：边缘缺口和中央狭缝。后一种情况是通过在模型中添加对称条件实现。

使用有边缘缺口的几何形状得到的系数 。

使用有中心狭缝的几何形状得到的系数 。

可以看出，只要缺口半径较小，两种情况下假定的乘数 的实际值都接近 1.2。缺口半径大、长度小的情况下，与裂纹的相似度就会降低。使用 进行简化是无效的。

为了绘制这些图，我们使用了 的解析值。在实际情况中，如果不知道这个值，可以使用到缺口一定距离的解，通过数值计算确定 。

事实上，任何尖角都有一个应力场衰减为 的区域，其中 是到尖角的距离。到目前为止，我们已经看到理想的裂纹 。不同开口角度下的 值如下图所示。

不同开口角度下应力奇点衰减的幂次。突出显示了 45^°、90^°和 135^° 的值。

这条曲线是通过求解超越方程绘制

, 为开口角度。

为了完整起见，我们可以在含内角的长条几何拉伸有限元模型中检验超越方程的解。该模型使用拐角的开口角度作为参数。

开口角度为 90^°时，含内角的长条几何中的 von Mises 应力 。

沿韧带的垂直应力。到尖角的距离通过韧带长度进行了归一化处理。虚线表示根据上述 p 值得出的理论解。

可以看出，应力-距离图中有一些几乎是直线的区域，这些区域在拐角附近，与理论斜率非常接近。

另一种奇点是由材料不连续性引起的，在实践中通常与几何奇点同时出现。在此，我们仅研究长条在拉力作用下的纯材料不连续性。

一个下部比上部硬的长条几何模型，其中绘制了载荷方向的应力。

从这幅图中的第一个图就已经可以看出一些通用特性：

• 自由表面出现奇点。
• 硬质材料的应力高于相应位置的软质材料应力。

为了更深入地研究这个问题，我们可以绘制显示应力衰减与材料界面距离的函数关系图。

沿自由边界加载方向上的应力与界面距离的函数关系图。实线表示软质材料的应力结果，虚线表示硬质材料的应力结果。参数 是软质材料与硬质材料的杨氏模量比值。

同样，可以在应力-距离对数图中发现一些直线，这表明应力随距离的变化而变化，如 。在两种材料中，幂 ‘’是相同的（用相同颜色表示的实线和虚线是平行）。奇点的强度受两种弹性模量的比值和泊松比值的控制。

观察上述（放大的）表面图中的变形形状，可以从物理角度对此进行解释：在相同载荷下，较软的材料比较硬的材料延伸得更多。也就是说，在软质材料中，载荷方向上的应变变大。这也意味着，在两种材料具有相同泊松比的情况下，横向收缩也会相应增大。不过，这种收缩会在两种材料的界面处受到抑制，从而产生局部应力奇点。

如果选择

这个奇点就会完全消失。

得出的结论是：在大多数情况下，材料变化会产生奇点。此外，在这种情况下，不连续性附近将存在一个应力随幂律衰减的区域。

至此，我们已经研究了有限元仿真中最常见的奇点类型，并发现它们有一个共同的特性：奇点附近的应力与距离呈幂律关系。

焊缝评估

对焊缝进行设计以使其能够安全地应对失效是工程领域的一项重要功课。虽然在一般情况下无法进行精确的应力评估，但已有大量研究提供了预测失效的系统方法。对于这种情况，造成问题的主要原因是焊缝的实际几何形状未知。根据确切的局部几何形状，判断焊缝是否需要引入应力奇点。更为复杂的是，焊缝中往往存在隐性缺陷。除了需要高质量焊缝的情况外，在这种情况下可以对焊缝进行打磨，并采用某种无损检测方法进行检查。但大多数情况下，对焊缝进行详细的局部应力分析意义不大。

有三种不同局部几何形状的圆角焊缝。

COMSOL博客：如何预测焊缝的疲劳寿命中介绍了焊缝的应力评估。

与其深入了解焊缝分析的细节，不如探讨更有趣的焊缝设计理念：

• 计算指定位置的应力，而不是焊趾本身的应力。
• 确定该应力的允许值，这通常必须通过实验来完成。
• 允许的应力值取决于您同意评估应力的方式和位置，因此它不是真正的材料属性。

在使用纸和笔计算的年代，所有的局部效应都被忽略了。允许的应力值不得不考虑到这一点，因此往往偏低。现代基于有限元的方法考虑了部分应力集中（由整体几何形状引起的部分，但不包括局部焊缝几何形状），因此允许的应力值更高，但仍远远低于纯材料测试所显示的应力值。

在使用有限元计算时，壳模型通常会返回求解所需的应力，而固体力学模型则会计算应力细节，而这在进行焊接疲劳分析时是不需要的。

推荐方法

有限元模型可能包含奇点的原因很多，并且本质不同。例如:

• 本文开头提到的博客中讨论的，边界条件会引起奇点。如果这种奇点给分析带来问题，可以通过完善边界条件来解决。
• 引入尖角的原因是局部几何形状的尺度较小，在全局尺度上建立圆角模型并不合理。在这种情况下，并不存在真正的奇点，而是一个明确定义的应力集中。最准确的方法是建立子模型来确定局部应力状态。在全局模型中，可以利用应力集中附近幂律衰减应力场的振幅了解应力集中的位置。另一种方法是将近场应力场知识与应力场与局部应力集中相关的知识结合起来，从而得出局部应力集中的估计值。

您可以参照焊接评估的方法，但要结合实际情况进行调整。要做到这一点，需要有大量的经验基础。以前的设计哪些失败了，哪些成功了？然后需要对设计进行分析，并尝试找到与经验相关的评估方法。

首先为这些设计建立有限元模型，并尝试确定一个区域，在该区域内的应力或应变场既不受局部缺口几何形状的控制，也不受整体几何形状的控制。至少在制定标准时，可能需要使用子模型。

使用什么标准通常并不明显。由于只是要进行相对比较，而不是将计算出的数字与任何物理强度值联系起来，因此有许多可能的选择。例如：

• 数量应易于计算。
• 数量不应对分析中的不确定性过于灵敏。
• 如果可能，数量应与物理场相关。例如，如果材料是脆性材料，那么查看最大主应力或主应变可能比使用 von Mises 等效应力准则更好。
• 如果疲劳是一个问题，数量必须对逆载荷反灵敏。
• 如果可能，请选择应变准则而不是应力准则。因为应变是直接根据位移计算得出的。应力则是通过应变的组合计算得出的。这表示应变张量中一个不准确的分量将传播到应力张量中的所有元素。
• 在 COMSOL Multiphysics^® 软件中，可以使用 安全 功能来评估大量不同的标准，包括用户定义的标准。

一般情况下，奇点的幂次是未知的。但我们知道，在某一区域，相关量的变化为

和 的值可以通过最小二乘法拟合或简单地使用应力应变对数-图中直线部分上的两点值来获得。由于 必须被看作某类奇点的恒定属性，因此计算出的值可用作方法有效性的检验。

在 已知的情况下，应将 的值与允许值进行比较。这与断裂力学中处理裂纹的方法类似。

百分比法

另一种获得允许应力水平的方法是将参考应力定义为在参考体积的给定部分（例如 5%）中超出的值。如果该参考应力低于允许值，则该设计被接受。使用这种方法，可以避免计算接近奇点的问题。只需计算出超过参考应力的体积即可，而该体积的边界在到奇点一定距离处，此处解可以很好地收敛。

这种方法看似简单，但应用起来却需要一定的标准。其中一个问题就是如何确定参考体积。如果使用结构的总体积，那么只需在低应力区域添加更多材料就可以降低参考应力，这当然是不合理的。参考体积必须与奇点周围特定区域的大小等因素相关。另一个缺点是，优化方法可能会选择重新定位应力，从而使参考应力减小，而峰值应力增大。

同样，我们也只能对类似结构进行比较来确定。

现在，我们来讨论如何使用百分比法计算应力值。在 COMSOL Multiphysics^® 中，无法直接计算 5% 的应力值。下面介绍 3 种替代方法。

方法 1

如果只需要一次求值，最快的方法通常是手动迭代几次。您只需创建一个积分算子（例如，intop1），然后对 intop1(solid.mises>sRef)/intop1(1) 这样的表达式进行求值。通过多次改变参考应力 sRef，很快就能找到与给定百分比相对应的值。

方法 2

使用模型方法自动执行方法 1。

方法 3

可以设置一个额外的方程，求解应力值，下文将对此进行说明。

待解方程如下：

是计算应力。它可能是如 von Mises 的等效应力、第一主应力或其他应力。当然，也可以使用相同的方法来计算应变或能量标准。用 表示参考体积， 为百分比。积分内的布尔表达式假定为 1 时为真，定义为 0 时为假。

为了在计算时更容易处理缩放，最好将方程改写为

可以像第一种方法那样，使用积分算子计算积分。在 COMSOL Multiphysics^® 中使用 全局方程 节点实现该方程的直接方法如下所示：

遗憾的是，这种方法行不通。不等式是不可微的，因此无法形成雅可比矩阵。刚度矩阵将只包含该方程的零点。不过，可以通过手动引入有限差分导数来规避这个问题。该表达式较长，需要您对 COMSOL^® 中基于方程的建模有一定的了解，下面的附加信息部分将给出详细解释。

下图所示是一个修改后的全局方程，它可以解决如何找到能给出预期百分比的应力的问题。

这里，用户定义的参数 dS 是应力增量，在附加信息部分用 表示。

我们以上文的缺口板示例来说明这种方法。由于参考体积应与板的尺寸无关，可以在缺口周围选择一个圆。在这种情况下，圆的半径可以根据结构的以下特征长度来选择：

• 宽度: 1 m
• 最小裂缝长度: 0.1 m
• 最小韧带宽度: 0.3 m
• 最大切口半径: 0.01 m

基于缺口尖端周围半径为 0.05 m 的圆确定的参考体积将远离结构的边界，但也将远离缺口本身的细节。

在不同的狭缝深度和切口半径下，5% 的参考体积所超出的应力水平。

对于所有狭缝深度值，5% 应力基本上与切口半径无关。它只对切口深度敏感。这与基本原理是一致的：避免对局部（可能是奇异的）应力场细节的灵敏性。无论使用尖角还是圆角，都能得到相同的结果。从本质上讲，这种方法提供的信息与应力强度因子相同：它测量的是奇点的强度。如果对具有相同圆角类型的结构进行比较，这种方法可以成为应力奇点处理的标准。

附加信息

该表达式包含两个项：第一个项产生残差，第二个项产生雅可比矩阵。这是一种通常可用于高级建模的解决方案。例如，如果创建精确的雅可比函数成本较高，则可以使用类似的表达式将正确的残差与近似的雅可比函数结合起来。

多个地方使用了 nojac(expr) 算子，用于确保不产生给定表达式的雅可比贡献。

雅可比项乘以系数 (sRef-nojac(sRef))。由于这个表达式的求值总是为零，因此这部分表达式不会产生残差。 sRef 相对于自身的导数就是 1，而表达式的剩余部分就是导数的对称有限差分表达式。

式中， 是应力的有限变化，应选择尽可能小的变化量，同时还要保证 计算出的体积与 有明显差异。一个好的水平是，当一个单元上的应力变化在等值面附近为时，将得出参考应力 。

结束语

理论上，索然无法计算奇点处的梯度和通量（应变和应力），但还是有一些系统的方法可以解决这个问题。不过，这些方法需要有足够的实验数据来解释所选择的临界量。

点击下面的按钮，下载文中使用的模型。

编者注:原博客最初由 Pawan Soami 撰写，发布于 2018 年 12 月 6 日。现已更新以反映最新版本软件的特征与功能。

内容简介

1. 什么是复合材料？
2. 细观力学分析
3. 宏观力学分析
4. 经典层压板理论和物理场接口
5. 材料模型
6. 复合材料仿真的结果计算工具
7. 复合层压板的多物理场分析
8. 复合层压板的优化
9. 多尺度分析

什么是复合材料？

由于复合材料具有特定的力、热、电和磁性能，因此在不同领域有着许多潜在的应用。例如，一些行业正在开发具有传感、驱动、计算、通信和其他功能的“智能”复合材料。在结构工程中，复合材料比传统的整体式材料更坚固、更轻，因此得到了广泛的应用。在使用这些材料设计复合结构之前，工程师必须充分了解它们的性能。

使用复合材料的优势和面临的挑战

与传统材料相比，复合材料具有多项优势，例如：

• 高强度重量比
• 耐冲击性强
• 高抗疲劳性和抗腐蚀性
• 摩擦性和磨损性增强
• 低导热系数和低热膨胀系数
• 耐高温

由于复合材料由多种材料混合而成，因此在使用这些材料时也会遇到一些挑战，包括：

• 各向异性特征
• 复杂的损伤和失效模式
• 原材料和加工成本高昂
• 难以重复利用和处置
• 不同组件的连接性差

复合材料的应用领域

由于复合材料具有以上优点，因此被广泛应用于以下领域：

• 航空航天工程（如卫星的机翼、机身和结构板）
• 国防安全（例如，坦克和潜艇）
• 风力发电机(例如，叶片)
• 建筑和施工 (例如，门、面板、框架和桥梁)
• 化学工程（例如，压力容器、储存罐、管道和反应堆）
• 汽车和运输工具（例如，自行车和汽车零部件）
• 海洋和铁路运输（例如，船体和铁路部件）
• 消费品和体育用品（例如，网球拍和高尔夫球杆）
• 电子产品（例如，配电支柱和接线盒）
• 矫形辅助工具
• 安全设备

复合材料的类型及其分类

复合材料的分类方法有多种，其中的一种方法是根据构成类型（即基体和增强材料）进行分类。根据基体材料的类型，可以将复合材料分为以下几类：

• 聚合物基复合材料 (PMC)
• 金属基复合材料 (MMC)
• 陶瓷基复合材料 (CMC)
• 水泥基复合材料 (CeMC)

根据增强类型，可以将复合材料分为以下几类:

• 纤维复合材料
• 晶须复合材料
• 颗粒复合材料

纤维、晶须和颗粒复合材料示例

纤维增强复合材料

相较于其他层压复合材料，纤维增强聚合物是当今非常流行的一种复合材料。这些材料通常由作为主要承载元件的长纤维和周围用于支撑纤维并传递载荷的基体组成。纤维以指定的方向排列在材料的每一层（或薄层）。许多这样的薄层铺设在一起就形成了可用于构建结构部件的层压复合材料。工业用纤维通常由碳、玻璃、芳纶或硼制成。根据纤维材料的类型，目前业界最常用的两种纤维增强聚合物是碳纤维增强聚合物（CFRP）和玻璃纤维增强塑料（GFRP），也称为玻璃纤维。

虽然我们可以使用 复合材料模块分析任何各向异性层压复合材料，但在这篇博客中，我们将重点讨论单向纤维增强聚合物。

层压板类型

复合层压材料是指由两个或多个单向层/层/薄片按照指定的方式，以一致或变化的纤维取向铺设而成。薄片可以由相同或不同的材料制成，并且可以具有各自的厚度。铺设序列由相对于层坐标系第一个轴的每层纤维的取向定义。

反对称平衡层压板的铺设顺序（0/45/90/-45/0）。

根据铺设顺序，复合材料层压板可以分为以下几种类型:

• 斜角层压板 (例如， 45/30/-45/-30)
• 交叉层压板 (例如， 0/90/0/90)
• 对称层压板 (例如， 45/30/30/45)
• 反对称层压板 (例如， 45/30/-30/-45)

由于纤维、板层和层压板的几何比例完全不同，因此分析复合材料层压板面临很多困难。这也是我们要在细观力学、宏观力学，以及两种（或多种）不同尺度上执行分析的原因。

细观力学

细观力学分析侧重于复合材料的组成层水平。它考虑了组成材料、材料界面以及材料的内部排列。细观力学分析不仅可以计算均质化的材料特性，还有助于了解细观层面的应力、应变、非线性、失效和损伤等。基于细观力学的均质化分析方法可分为两大类：

1. 分析法（例如，混合规则）
2. 数值方法（例如，使用代表性体积单元 (RVE) 或重复单元 (RUC) 进行有限元分析）

在模型开发器树中的材料 节点下，多相材料 和有效材料 节点有多个用于分析计算有效性的混合规则。有效材料 节点内置于复合材料模块，具有以下混合规则：

1. 体积平均
2. 质量平均
3. 谐波体积平均
4. 谐波质量平均
5. 幂律
6. Heaviside 函数
7. Voigt–Reuss 模型
8. 修正的 Voigt–Reuss 模型
9. Chamis 模型
10. Halpin–Tsai 模型
11. Halpin–Tsai–Nielsen 模型
12. Hashin–Rosen 模型

显示了 混合规则选项的 有效材料特征设置窗口

要使用有限元方法数值计算均质材料特性，需要使用代表性体积单元或重复单元。对于周期性材料，代表性体积单元与重复单元相同，但对于非周期性材料，重复单元的概念无效，因此必须使用代表性体积单元材料子体积。

60％ 纤维体积分数的纤维复合材料层的晶胞。

在 COMSOL Multiphysics^® 中，使用 固体力学 接口中的 单元周期性 节点进行基于细观力学的均质化。该接口有两种不同的边界条件：周期性 和均质。周期性 边界条件适用于周期性材料，需要使用重复单元材料子体积。对于非周期性材料，可以通过代表性体积单元材料子卷应用均质 边界条件。在这篇博客中，我们将重点讨论周期性单向纤维复合材料的均质材料特性。

我们从一个包含纤维和基体的晶胞几何结构开始分析。首先需要给出纤维和基体的材料属性。然后，可以使用单元周期性 节点中的操作按钮设置所需的模型节点和研究。自动创建的研究将计算均质材料的材料数据。

6 种不同载荷下，晶胞中的 von Mises 应力分布和变形。

你可以查看纤维复合材料的细观力学模型和复合材料气瓶的细观力学和应力分析案例模型，了解更多内容。

宏观力学分析

宏观力学分析基于均质材料确定复合结构的响应。层压板的均质材料特性可通过细观力学分析或实验方法获得。宏观力学分析的目的是计算层状结构在各种载荷和边界条件下的整体响应。宏观力学分析包括以下几个不同步骤。

复合材料仿真的预处理方法

模拟复合层压板, 需要指定以下几个特性:

• 层数
• 每一层的均质材料特性
• 层压板主要材料方向的定向
• 每一层厚度
• 铺设顺序

复合材料层压板的横截面显示了每一层的纤维厚度和取向。

要定义层压材料的属性，需要使用 多层材料 节点。在该节点中，可以添加所需的层数，输入内容可以直接输入表格，也可以从文本文件中加载。指定输入后，就可以预览层压材料的横截面和铺设顺序。您可以将包含层压板定义的多层材料保存在材料库中，方便后续加载使用。

多层材料节点示例。

使用 多层材料 节点定义层压材料后，就可通过 多层材料链接 或 多层材料堆叠 节点将其连接到几何边界。在此过程中，层压材料坐标系以及几何表面相对于层压材料的位置也会被定义。层压坐标系还能进一步用于解释铺设顺序，并创建多层局部坐标系。多层材料链接 和 多层材料堆叠节点还有更多的选项，可以将多层材料转换为对称、非对称或重复层材料。还包括模拟厚度在空间上变化的模型选项。多层材料堆叠 节点可用于区域建模，在不同的几何选择中，复合材料的铺设顺序会有所不同。

多层材料链接 和 多层材料堆叠 特征的应用示例。

请注意 单层材料 特征是为单层材料设计的特殊的 多层材料 特征。

经典层压板理论和物理场接口

现在，我们已经定义了层压板并将其添加到几何边界上。接下来，我们来介绍经典层压板理论。通常，我们会使用下列三种理论之一分析层压复合壳：

1. 等效单层（ESL）理论
   • 经典层压板理论（CLPT）
   • 一阶剪切变形层压板理论（FSDT）
   • 高阶剪切变形层压板理论
2. 三维弹性理论
   • 三维弹性理论
   • 分层理论
3. 多模型方法

一阶剪切变形等效单层理论: 壳接口

在一阶剪切变形等效单层理论中，计算整个层压板的均质材料特性，并仅在中面上求解方程。该理论采用类壳公式，自由度（DOF）为网格边界上的三个位移和三个旋转。该理论适用于薄至中等厚度的层压板，可用于计算总挠度、特征频率、临界屈曲载荷和面内应力等全局响应。与分层理论相比，一阶剪切变形等效单层理论计算成本较低；但对于较厚的层压板，它需要一个剪切修正系数。

一阶剪切变形等效单层理论中的自由度节点。

在 COMSOL Multiphysics^® 中，壳 接口的 线弹性材料，多层；超弹性材料，多层 和 压电材料，多层 等多层材料特性都是基于 一阶剪切变形等效单层理论。此外，膜 接口中的线弹性材料，多层功能也是基于等效单层理论，可用于对弯曲刚度忽略不计的极薄复合薄膜进行建模。

在风力发电机复合材料叶片的应力和模态分析案例模型中，风力涡轮机复合叶片是使用 壳 接口模拟的。目标是找出在重力和离心力作用下叶片的表层和隔板的应力分布情况。

风力涡轮机复合叶片示例。叶片的表层和隔板的应力分布情况。

您可以查看以下示例，了解更多内容:

• 复合材料轮辋的应力和模态分析
• 复合材料气瓶的屈曲分析

分层理论: 多层壳接口

在这个理论中，方程也在厚度方向上求解。因此，它可用于非常厚的层压板，包括分层区域。该理论采用类似固体的表述方式，其中自由度以三个位移的形式分布在厚度方向上。该理论适用于中等厚度到较厚的层压板，可用于预测正确的层间应力和分层，并进行详细的损伤分析。与一阶剪切变形等效单层理论理论相比，它支持非线性材料模型，并且不需要剪切校正因子。

分层理论中的自由度节点。

从公式的角度来看，分层理论与三维弹性理论非常相似。但是，与后一种理论相比，它具有以下优点：

• 层压板坐标系和层局部坐标系容易定义
• 面内和面外形函数可以具有不同的阶次
• 无需构建具有许多薄层的三维几何结构
• 面内有限元网格划分独立于面外网格划分
• 分层和界面数据容易处理

在 COMSOL Multiphysics^® 中，多层壳 接口基于分层理论。简支复合材料层压板的弯曲案例模型中使用 多层壳 接口和 壳 接口对简支复合材料板进行了弯曲分析，目标是将两种接口得到的厚度应力与给定基准的三维弹性解进行比较。

简支复合板示例。左图：使用 多层壳 接口模拟的板中的 von Mises 应力分布。右图：厚度横向剪应力对比图。

你还可以查看复合材料层压板的强迫振动分析案例模型，了解另一个示例。

多模型法: 壳接口与多层壳接口耦合

多模型法是将等效单层理论与分层理论相结合，应用于复合材料几何结构的不同部位或不同层，以获得可接受的结果，并优化利用计算资源。除了 多层壳 和 壳 接口外，还需要使用 多层壳-壳连接 多物理场耦合节点将这两个不同的物理场接口在厚度方向上进行耦合。

使用多模型方法分析复合材料叶片案例模型中通过耦合 多层壳 和 壳 接口模拟了一个复合材料叶片，目标是对比不同方法的求解时间。

使用不同方法计算的复合材料叶片中的 von Mises 应力分布。

选择合适的层压板理论

基于上述方法，你可以选择合适的层压板理论。一个简单的经验法则是选择基于层压长宽比，即层压板长度与层压板厚度的比值的层压板理论。

基于层压长宽比的两种层压理论的有效性范围。

材料模型

下表列出了在不同物理场接口中用于分析复合材料的材料模型和非弹性效应。

你也可以查看正交材料压力容器 – 壳版本和含压电材料的多层壳案例模型，了解更多内容。

损伤、脱层和首层失效理论

许多复合材料都是准脆性材料，在达到临界应力或应变水平后，初始弹性阶段随后进入非线性断裂阶段。当达到该临界值时，裂纹会逐渐扩展，直至材料断裂。裂纹增长导致的材料刚度下降可以通过多层壳和壳接口中的损伤功能进行模拟。目前有两种损伤模型可供选择：标量损伤模型和 Mazars 混凝土损伤模型。此外，还有几种应变软化损伤演变定律可供选择。为避免网格敏感性，可以选择裂纹带 或隐式梯度选项来使用空间正则化方法。

脱层或层间分离是层压复合材料的一种常见失效模式。包括载荷、材料缺陷和环境条件在内的各种因素都可能引发层间分离的发生和传播。要模拟脱层现象，可以使用 多层壳 接口中的 脱层 功能。脱层理论以内聚力模型（CZM）为基础，并包含多个牵引分离定律。要了解更多信息，请查看 COMSOL 案例库中的 复合材料层压板的混合模式脱层和层压壳中的渐进脱层案例模型。

多层壳 接口和 壳 接口的 安全 功能中提供了多种首层失效理论。具体来说，像 Tsai-Wu、Tsai-Hill、Hoffman、Hashin、Hashin-Rotem、Puck 和 LaRC03 等理论在复合材料仿真中非常有用。要了解更多信息，请参阅层压复合壳的失效预测案例模型。

屈曲

使用这两种层压理论中的任何一种都可能产生线性屈曲；不过，与分层理论相比，一阶剪切变形等效单层理论在寻找临界屈曲载荷系数方面更有效。它可以优化层叠结构，以使临界屈曲载荷最大化。更多信息，请参阅复合材料气瓶的屈曲分析案例模型。

多层材料连续性

如果在多层壳 接口所选的几何结构上激活了一个以上的单层材料、多层材料链接 或 多层材料堆叠，那么默认情况下，这些不同多层材料之间的 DOF 是断开的。使用 连续性 功能可以连接相邻的两个层状材料。利用该功能，你可以对层叠脱落情况进行建模。相关建模示例，请参阅 COMSOL 案例库中的复合板的削层案例模型。

当使用壳 或膜 接口时，DOF 只存在于中面上，因此它们在分层材料之间总是连接的。

在并排放置的两块层压板之间设置连续性的不同方法。

A, B, D 矩阵计算

标准刚度和柔度矩阵可通过 壳 接口中的 线性弹性材料，多层 节点进行计算。可用的四个刚度矩阵包括拉伸刚度矩阵 (A)、弯曲-拉伸刚度矩阵 (B)、弯曲刚度矩阵(D)和剪切刚度矩阵(As)。更多详情，请查看层压复合壳的材料特性案例模型。

有时，复合材料层压板的材料特性由 A、B 和 D 矩阵提供。在这种情况下，可以使用 壳 接口中的 截面刚度 材料特征。

复合材料仿真的结果计算工具

在进行宏观力学分析时，COMSOL Multiphysics® 中有多种功能可用于结果计算。下面我们将讨论其中的一些功能。

多层材料数据

由于几何结构只包含表面，多层材料 数据集用于显示有限厚度几何结构的模拟结果。使用该数据集，可以在法线方向上缩放层压板厚度，这对薄层压板非常有用。 多层材料 数据集还提供在以下位置进行计算的选项：

• 网格节点
• 界面
• 层中面

多层材料 数据集包括选择和取消选择多层材料链接或多层材料堆叠中的不同层的选项。其他一些数据集，如镜像、数组、三维截线，三维截点 和 旋转，也可以与 多层材料 数据集一起使用。

体绘图和表面绘图

多层材料 数据集可直接使用不同的体图、表面图、切面图等。

使用 多层材料 数据集绘制的各种结果图。

多层材料切面图

对于复合材料层压板，多层材料切面 绘图在制作切面时提供了更大的自由度。一些有用实例包括创建切面绘图：

• 通过一个或两个层
• 穿过多个（或所有）层（请注意，不需要在厚度方向上放置切面）
• 在层的特定位置，但不在中面上

使用 多层材料切面绘图创建的层压板每一层中面上的 Von Mises 应力。

全厚度图

该绘图用于确定不同量通过层压板厚度上的变化。你可以在边界上选择一个或多个几何点，也可以选择创建切截点数据集或直接输入点坐标。

层压板上某一点横向剪切应力在厚度方向的变化。

线图或点图

要创建特定变量的线图，需要使用基于 多层材料 数据集的 三维截线 数据集。同样，要创建特定变量的点图，也需要使用基于 多层材料 数据集的 三维截点 数据集。另一种解决方案，可以将包含特殊算子的变量与 多层材料 或 解 数据集一起使用。

复合层压材料的多物理场分析

结构连接

在许多情况下，系统的结构分析需要使用不同的单元类型或物理场接口。下表列出了可用于连接不同结构物理场接口的多物理场耦合。

请参阅多层壳与实体和壳的连接案例模型，查看连接壳和结构单元的示例。

热膨胀

可以使用下列物理场接口模拟复合结构中的热膨胀：

• 壳传热
• 壳 或者多层壳

不同物理场之间的耦合通过下列多物理场耦合节点来定义:

• 热膨胀，多层

有关建模示例，请参阅案例库中的层压复合壳的热膨胀案例模型。

焦耳热和热膨胀

复合结构中的焦耳热和热膨胀可以使用以下物理场接口模拟：

• 多层壳中的电流
• 壳传热
• 多层壳

不同物理场之间的耦合通过下列多物理场耦合节点来定义:

• 电磁热，多层壳
• 热膨胀，多层

声 – 复合材料的相互作用

声–复合材料的相互作用可以通过以下物理场接口模拟:

• 压力声学
• 壳 或者 多层壳

声-结构边界 多物理场耦合节点用于定义这两个物理场接口之间的相互作用。

流体–复合材料的相互作用

流体–复合材料的相互作用可以通过以下物理场接口模拟：

• 层流
• 壳 或者多层壳

流-固耦合 多物理场节点用于定义这两个物理场之间的相互作用。

压电 – 复合材料的相互作用

压电 – 复合材料的相互作用可使用下列物理场接口模拟：

• 多层壳中的电流
• 壳 > 压电材料，多层 或者 多层壳压电材料

压电，多层 多物理场耦合节点用于定义这两个物理场接口的耦合。了解更多内容，请参见含压电材料的多层壳教程模型。

压阻 – 复合材料的相互作用

压阻 – 复合材料的相互作用可以使用以下物理场接口模拟:

• 多层壳中的电流压阻壳
• 多层壳

压阻，多层 多物理场耦合节点用于定义这两个物理场之间的相互作用。

集总机械系统 – 复合材料的相互作用

集总机械系统 – 复合材料的相互作用可使用下列物理场接口模拟：

• 集总机械系统
• 多层壳

集总结构连接 多物理场耦合节点用于定义这两个物理场接口之间的相互作用。

复合层压板的优化

复合层压板是一种合成结构，总是有可能在每层材料、每层厚度和铺层顺序等方面对设计进行优化。利用 优化模块 的功能，可以对复合层压板的不同要素进行优化。要了解此类优化，请查看铺层顺序的优化案例模型，其中根据 Hashin 失效准则对复合层压板的铺层顺序进行了优化。

优化后的复合层压板示例。原始布局（线框）和优化布局（实体）的位移。

多尺度分析

复合材料既可以在宏观尺度上进行分析，也可以在细观尺度上进行分析，无论哪种分析都有其优点和局限性。通过宏观和细观分析，可以深入了解复合材料结构及其成分对宏观加载荷的响应。完整的多尺度分析包括宏观分析和每个材料点的细观分析，计算成本高昂。如果我们将分析限制为只包括几个关键材料点，就可以通过使用 固体力学 接口中的 单元周期性 功能和 多层壳 接口进行多尺度分析。

要查看多尺度分析的实际效果，请参阅失效的细观力学：复合材料结构的多尺度分析案例模型。在这个示例中，首先进行细观力学分析以获得均质材料属性，然后使用分层理论进行宏观力学分析以获得全局响应。最后一步是进行细观力学分析，计算局部应力场和应变场以及基于全局平均应变的失效风险。

多尺度分析示例。左图：基于宏观力学分析的复合材料圆柱体应力。右图：使用细观力学分析法测量不同材料点的应力。

下一步

使用复合材料模块，您可以设计、分析和优化由线性或非线性材料组成的多层复合材料结构。要了解有关复合材料模块的更多信息，请点击以下按钮联系 COMSOL。

理解平面周期性结构的衍射

本文，我们将以入射到具有无限周期的平面结构上的平面光波为例来说明。该平面上方和下方的介质可以具有不同的折射率，并假定为无损耗和无限域。在这些介质的交界面，可以存在材料性质和形状等复杂的周期性结构。入射到周期性结构上的光线至少会发生镜面反射；也会发生折射（也称为镜面透射），通常还会有一些损耗，因为电磁能会转化为热能。我们可以通过斯涅尔定律计算出反射角和折射角，但入射光在周期性结构中的反射、透射或损耗的部分需要通过数值分析计算。

以一定角度入射到平面周期性结构上的平面波。突出显示了周期性结构的一个基本单元。

如前所述，也存在高阶衍射的可能性。当周期性结构散射的光被相长干涉到不同的方向时，就会出现这种情况。下面展示了这种结果的一个示例。

入射到周期性基本单元的线性偏振平面波（黄色）示意图。在反射（红色）和透射（蓝色）中，入射光被散射成几个不同强度和偏振的衍射级。

要确定进入这些其他相似方向的光的比例，同样需要建立一个数值模型，但要了解光会散射到哪些方向，可以通过一种纯几何方法来实现，这种方法被称为埃瓦尔德球结构。在开始数值分析之前，熟悉这种方法是很有帮助的，这也是我们将在这篇文章中介绍的内容。埃瓦尔德球几何结构既可用于单向周期性平面结构，也可用于平面内的双向周期性结构。

单向周期性结构

例如光栅等某些平面周期性结构仅在一个方向上具有周期性变化，即该结构沿三维方向没有变化。当入射光在三维空间的法线平面上传播时，可以被简化模拟为沿一个方向具有周期性的二维平面。

以一定角度入射到单向周期性结构上的平面波，在结构或场中沿面外方向没有变化。突出显示了一个基本单元。

对于这些结构，我们只需考虑基本单元间距 ，并首先在 倒易空间中绘制一组晶格点，因此下图中的尺寸单位为逆长度。这些晶格点的连线对应于周期性结构的界面平面。晶格点之间的间距为 ，晶格点的索引从第四个晶格点 开始，可将其视看作位于基本单元的中间。然后，在晶格点连线的上方和下方绘制两个半圆。反射侧的半径为 ，透射侧的半径为 ，两侧的折射率分别为 和 ， 为自由空间波长。对于与法线夹角为 的入射光，这些圆的公共中心与晶格的第零个点偏移了 。位于这些半圆内的晶格点对应于可能的衍射级。

用于确定单向周期性平面结构的衍射级的几何结构，该结构受到以一定角度入射的平面波的照射。请注意半圆（白点）的中心是如何偏离第零晶格点的。

这种结构还可用于确定衍射的方向，并为每个方向分配一个索引。从半圆中心投影到晶格点的矢量对应于每个衍射级的 矢量。这些晶格点的索引在两侧的符号相反。指向第零个晶格点的箭头始终存在，代表镜面反射和透射。其他衍射级的存在取决于波长、折射率、间距和入射角度。COMSOL案例库中包含了两个建立此类模型的案例：使用 RF 模块的表面等离激元线光栅（RF）和使用波动光学模块的表面等离激元线光栅分析仪（波动光学）。

单向周期性平面结构各种衍射级的波矢量。请注意反射衍射级与透射衍射级之间索引符号的转换。

双向周期性结构

现在，我们来看看在两个方向上具有周期性的平面结构的衍射情况。下图显示了构建平面的矩形、菱形和六边形基本单元。这些单元由两个单元矢量定义： 和 ，它们从一个点开始，沿着相邻的边到达下一个顶点。虽然我们可以自由使用任何坐标和方向，但在本文中，我们将始终选择 向量与全局笛卡尔 x 轴对齐，并始终从光照方向俯视基本单元。此外，还有两个基向量 和 ，描述了基本单元在平面上的移动方式，用于构建平面。也就是说，要构建整个平面，需要在 的 基础上复制基本单元，而 和 的值可以是任意整数。这两个矢量的叉积大小可用于计算基本单元的面积：。

矩形、菱形和六边形基本单元构成了二维平面。单元矢量与单元的两条边相对应，而基矢量则描述了如何移动单元来构建平面。

这些基矢量用于定义两个倒易空间衍射矢量： 和 ，其中 是周期性平面的法向量，即 +z 轴。这些衍射矢量与基矢量垂直，并通过取整数和在周期性平面上创建衍射晶格：，晶格中的每个点对应于 和 方向上 和 的索引对。在基本单元的传输侧，点的位置相同，但索引对调，且符号相反。

在倒易空间绘制的衍射矢量和晶格点。

现在，我们可以在三维空间的周期性平面上将这些衍射点可视化，并在平面上方和下方添加一个半径等于材料中波矢量的半球。通过半球，我们可以得知在反射和透射中存在哪些衍射级。刚开始，我们以点 为半球中心，代表法向入射光线。

平面波光（黄色箭头）通常入射到一个周期性六边形单元上。衍射点绘制在周期性平面上，位于反射半球和透射半球内的突出点表示将出现的衍射级。

接下来，我们来看看入射仰角和入射方位角变化时的情况。考虑到我们习惯上选择保持 向量与球坐标的 +x 轴对齐，增大入射仰角意味着入射波矢量首先绕 –y 轴旋转；然后，入射方位角增大，入射波矢量随之绕 +z 轴旋转。因此，入射仰角从 开始，入射方位角从 开始，如下图所示。入射波矢量和周期性平面的法线定义了入射平面。当光从法线入射：，入射平面被定义为 xz 平面。

入射仰角和入射方位角表示入射波矢量（黄色）的一系列连续旋转，先是绕 –y 轴旋转，然后绕 +z 轴旋转。图中也显式了入射平面。

入射角的变化改变了半球中心的位置。从半球中心到 点的倒易空间距离为 ，该位置在平面内的移动量为 和 ，如下图所示。因此，仰角和方位角的变化往往会导致出现不同的衍射级。

以非零仰角和方位角入射的平面波光会移动半球的中心，从而产生不同的衍射级。

通过这些半球，我们还引导每个衍射级的波矢量。将衍射级点投影到半球上，会得到另一组点，而每个衍射级的波矢量等于从半球中心到这些投影点的矢量。

最后，通过这些矢量，我们还可以知道偏振状态。对于每个衍射级，偏振状态都会根据琼斯矢量的面内和面外分量表示。每个衍射级的平面都是波矢量和周期性平面的法矢量所描述的平面。对于所有衍射级，琼斯矢量的面外分量对应于电场平行于周期平面的波。

结论

综上所述，我们可以得出以下结论：使用埃瓦尔德球的几何构造可以理解平面性周期结构衍射，并且能够获知在反射和透射中会出现哪些较高的衍射级。我们还可以获知波矢量以及用于定义琼斯矢量方向的平面集。在求解数值模型时，会自动得到这些信息，因此这种几何构造并不是必须的，但它有助于我们建立理解和直觉。

进阶学习

如果您想开始学习高阶衍射建模，下面的示例模型是很好的起点，这些模型可以用 RF 模块或波动光模块建立。

• 使用RF模块建立的模型:
  • 表面等离激元线光栅
  • 菲涅尔方程
  • 电波暗室的锥形吸波器建模
  • 频率选择表面周期性互补开口谐振环
  • 六边形光栅
• 使用波动光学模块建立的模型:
  • 表面等离激元线光栅分析仪
  • 菲涅尔方程
  • 六边形光栅
  • 六边形等离激元滤色镜
  • 频率选择表面周期性互补开口谐振环
  • 石墨烯超材料完美吸收器

编者注: 这篇博客最初发表于 2016 年 2 月 2 日，为呈现软件的新特征和新功能，对其进行了更新。

批处理扫描：大多数仿真分析的有效解决方案

在上一篇博客中，我们探讨了批处理扫描中任务并行的附加值。这篇文章，我们将讨论批处理扫描的另一个重要功能：在求解过程中检索参数化扫描的解。使用批处理扫描可以随时检索部分解，即使在求解器无法求解某些参数时，也可以这样做。当所研究的问题中每个参数的解都独立于其他所有参数的解时，就适合采用批处理扫描。

在求解器扫描过程中，希望查看部分解的情况很多，例如：

• 当基于表格中的输入数据扫描时，你发现求解表格中列出的一部分数值所需的时间格外长，但无法判断具体是哪些值。这时，你希望能查看尽可能多参数的解，来判断是否需要中止求解过程，或者在所有扫描完成之前开始分析结果。
• 你用数学表达式描述某一材料属性或边界条件，结果却得到某些参数的非物理解。
• 你用 C 语言或 FORTRAN 语言编写了一个有关定义复杂材料的外部函数，但未检查所有的输入数据，结果某些参数值的输出数据出现了异常。
• 你正运行参数化扫描，其中一个参数为几何尺寸，但参数尺寸的范围定义得过大。当你意识到这一点时，求解已经运行了很长时间，所以不希望停止。

如果你在以上任何一种情况中运行批处理扫描，每个参数都能在一个单独的进程中求解，这个进程可以独立开始和停止。求解完成的参数值都被存储为一个.mph 文件，你可以在求解过程中打开任何数量的文件进行查看。在下面的示例中，你将了解如何将探针结果保存到文本文件中。

静电传感示例

为了演示批处理扫描，我们使用 COMSOL Multiphysics^® 案例库中的 电传感器示例模型来展示批处理扫描的操作。此模型展示了如何通过在盒的边界施加电位差，获得盒内的图像信息。模型计算的表面电荷结果取决于盒内物体的介电常数，该结果也可用于判断例如盒内是否存在物体，甚至可以帮助确定物体的形状。电阻抗断层成像案例中也采用了类似的技术进行医疗诊断，但使用的是交流电。

盒外的表面电荷。

盒内待检测的物体。

原示例模型中将物体厚度拉伸为 0.8 米。本文示例中，物体厚度被拉伸为 0.5 米。

参数化扫描的输入数据

假设有一组包含十个样本物体，这些样本形状相同但由具有不同介电常数的材料制成。我们要探究的是检测物体及其形状的能力，不管物体由什么材料制成。十个样本的介电常数值已存储在内插表 int1 中。实际上我们不会使用这个表格进行任何插入操作，而是把它作为一种存储和调用表格值的简便方法。（在这里，我们可以将内插方法从默认的 线性内插 改为 最相邻内插，但在本示例中结果不变。）

介电常数的表格。

现在用全局参数 sn（样本编号的英文缩写）对这个样本集添加索引。

作为全局参数的样本编号参数 sn 。

此模型在材料节点下定义了三种材料。材料3表示盒内的星形物体，材料 2表示其他物体。材料 1对应盒内物体周围的空气。我们通过适当调用材料3材料属性中定义的内插表 int1，来改变星形物体的介电常数，使用的语法为 int1(sn)。请注意，这里无需考虑单位，因为介电常数是无量纲。

通过参数 sn 在表格中索引得到的星形物体的介电常数。

探测表面电荷

为了监测星形物体上方的表面电荷，我们在盒子上表面的 x = 1.5，y = 1.5，z = 1.0 位置处定义了一个边界点探针。为了快速定义此位置，我们可以单击物体上表面的一个大致位置，然后在 设置 窗口中手动调整坐标值。

研究物体上方的边界点探针

至于点探针表达式，我们选择用变量名 es.nD 和单位 nC/m^2 来表示表面电荷的大小（模值）。

单位为 nC/m^2 的表面电荷的 点探针表达式

稍后，我们会将此探针的输出写入文本文件。为此，单击 更新结果。此操作会生成 派生值 和 表格项。

探针的 更新结果 按钮

在 COMSOL Multiphysics^® 中设置批处理扫描

默认情况下，软件不会显示批处理扫描选项。可以在模型开发器的工具栏选择 显示更多选项 来启用该选项。在 显示更多选项 对话框选择 研究>批处理和集群。

在模型开发器中的 显示更多选项 菜单激活 的高级研究选项

如果右键单击 研究 节点，可以在列表中选择 批处理扫描 (单次运行，选择 批处理，而不是 批处理扫描)。

从 研究 菜单中选择 批处理扫描

在如何定义扫描范围方面，批处理扫描与参数化扫描相似。在这里，我们使用整数 1 到 10 来定义参数 sn 扫描时的步长。

批处理扫描的参数范围。

批处理扫描 模型树 设置 窗口的 批处理设置 栏中包含了用于控制结果的设置。

控制批处理扫描的设置。

默认文件名为 batchmodel.mph，用于自动创建一个 MPH 文件序列，其中的每个文件对应于扫描中的每个参数。可以更改文件名，不过此处演示的示例仍保留此文件名，生成的 MPH 文件名将分别为 batchmodel_sn_1.mph，batchmodel_sn_2.mph，……，batchmodel_sn_10.mph。稍后，我们会探讨这其中的每个文件为何是一个可以单独打开并执行后处理的完整的模型。

扫掠前 部分有两个选项：清除网格和清除解。我们清除了这两个复选框，因为只有在执行远程集群计算或云计算时，这两个选项才真正有用，并希望通过网络传递的文件尽可能小（只有使用网络浮动许可证时才能这样做）。

在扫描过程中 部分有两个选项：同步解 和 同步累积探针表。同步解 将存储的所有文件中的结果从开始仿真时收集到“主模型” MPH 文件中，最终结果与运行参数化扫描的结果极其相似。在这个示例中，选择清除这个复选框是因为我们假设只对查看单个结果感兴趣。此外，如果扫描的范围过大，收集全部结果可能会消耗大量时间和内存。保持选中 同步累积探针表 状态，表示对探针的输出结果进行累积，与完整的求解信息相比，收集这种信息非常轻松。

在 扫掠后部分中，我们勾选输出模型到文件复选框，确保自动保存的 MPH 文件包含（每个参数的）解，同时会增加磁盘空间。

在下一节，将 群集计算 设置更改为 用户控制。然后，将 目录 更改为 C:\COMSOL。这就是存储所有批处理文件的位置。

如果你有网络浮动许可证，并且希望使用远程安装的 COMSOL 软件执行这一计算，可以使用 指定服务器路径 选项。本文我们只讨论包括单用户许可证在内的 COMSOL 所有许可证类型都能使用的功能。

将探针数据保存到文件

为了将探针数据保存到文件，我们不能使用常规的探针表，而是使用一种特殊的 累积探针表。第一步是手动添加一个探针表，即右键单击 结果 下的表格 节点，将表格名称改为 保存的探针表。稍后，我们会通过批处理扫描在其中写入累积探针表的数据。当需要将异步生成的数据在表格中按照一定的顺序排列时，就可以使用累积探针表。请记住，一般来说，我们不清楚不同的参数对应的解应该按照什么顺序排列，而累积探针表帮我们解决这一难题。

手动添加保存的探针表。

在保存的探针表设置中，将 存储表 设置改为 模型和文件中。在本示例中，我们将探针结果保存到 C:\COMSOL\results.txt。如果扫描的数量过于庞大，则可能需要提高 最大行数 的值（这通常适用于需要大量扫描的情况或瞬态解）。

保存的探针表设置。

在 批处理扫描 设置的 求解时输出 一栏中，将 输出表 改为 保存的探针表。这就是我们的累积探针表。

将保存的探针表 用作 累积探针表。

运行批处理扫描

现在，可以单击 计算 启动批处理扫描。求解过程中，图形窗口下的信息窗口显示一个 外部进程 窗口。

外部进程窗口。

通过外部进程窗口，我们可以了解正在运行的批处理进程及其状态的概况。在本文示例中，共有十个这样的进程，分别对应扫描中的每一个参数。

在这个阶段，大部分的用户界面不能交互，但可以单击进度条左侧的脱离作业，重新获得对用户界面的控制。

脱离作业 选项。

使用脱离作业后，就看不到进程状态的实时更新了，不过这时可以对各个进程执行操作。

脱离作业后的 外部进程 窗口。

此时，我们可以停止所有进程，或者单击某一行来停止该进程，甚至可以单击表中状态为 完成 的进程，在另一个 COMSOL Desktop 窗口中打开它。在新打开的窗口中，我们可以执行任何可视化操作或后处理任务，这些任务通常需要在模型中运行。

选择一个完成的外部进程，单击 打开文件，即可在一个新的 COMSOL Desktop 窗口中显示结果。

单击 外部进程 窗口左上角的 附加作业，就可以重新显示进程状态的实时更新。

对 外部进程 窗口附加作业。

探针结果文件

我们可以使用 Notepad 等文本编辑器打开包含探针输出的 results.txt 文件。

包含边界点探针结果的 results.txt 文本文件。

我们也可以在求解过程中查看这个文件，从而在全部参数都完成收敛之前快速浏览结果。不仅是批处理扫描可以这样做，对于常规的参数化扫描，在表格保存为文件后也可以在求解过程中查看。

使用批处理扫描功能时需要考虑的其他因素

性能

如果扫描中每次计算速度都很快，并且不需要很大的内存，如上文中的示例所示，批处理扫描所需的计算时间会比常规的参数化扫描的时间长。因为采用批处理扫描时，每个参数都会启动一个 COMSOL Multiphysics 进程。对于每一个计算要求都较高的参数扫描，这个额外启动进程的时间相对较短，可忽略不计。

参数错误处理

假设我们调用的是用户定义的外部 C 语言函数，而不是内插表；同时假设我们在 C 代码中犯了一个编程错误，导致一个或几个参数在扫描时出现分段错误。这种情况下，我们仍然可以浏览 C:\COMSOL 目录（或者指定保存 .mph 文件的目录），打开其中的文件执行可视化及后处理，甚至进一步的计算工作。

在上文示例中，演示此类错误的一个简单方法是将其中的一个介电常数值设为零，然后启动批处理作业。介电常数为零意味着我们向求解器提供了一个定义模糊（甚至是异常）的方程（类似于 0 = 1），此时会显示一条错误消息，提示求解器找不到出错参数的对应解。同样地，仍可以打开保存的 .mph 文件并执行后处理。如果更正了表格，将介电常数改成非零，就可以重新运行批处理作业并且不会收到错误消息。请注意，如果错误一直存在，则对 外部进程 窗口附加作业时始终会提示错误。此时，仍会更新累积探针表中与其他参数对应的项。更正错误后继续运行扫描，对外部进程窗口附加作业后就可以再次访问该窗口了。

管理多核处理

如果你有一台多核计算机，这是现代电脑的常见配置，就可以更改 批处理 的相关设置，控制批处理扫描中可运行的并发进程数，以及每个进程可用的内核数。假设你使用的是六核机器，那么可以将 并发作业数改为三，将 内核数 改为二。这样设置后，可以并行求解三个参数，每个求解器进程可使用两个内核，所以对于上面的示例，仿真时间减少了一半。

如果所运行的仿真中每个参数的计算量都很少，可以增加并发作业的数量，使其与计算机的内核数相同。对于计算量较大的问题，应该保持并发作业数为一，以充分利用求解器的多核处理能力。

批处理设置中的并发作业数和内核数。

请注意，也可以在批处理扫描 节点研究扩展 一栏下控制并发作业数。这种情况下，软件会用实际内核数除以并发作业数，自动计算出“内核数”。（要实现自动计算，请不要勾选 批处理 设置下的 内核数 复选框，就不会使用以灰色显示的数字 1。）

批处理扫描和 App 开发器

借助 App开发器，我们可以在用于构建 App 的模型中使用 批处理扫描 节点。此时，App 充当批处理扫描的“驱动”，且无法在 外部进程窗口中看到相关信息。也无法在该窗口中看到保存的 .mph 文件，就像此 App 根本没有用户界面。运行结束后，我们可以打开这些文件进行传统的模型后处理工作。如果需要更灵活地创建可使用 批处理扫描 的App，可以使用一个内置的方法来包含 批处理扫描 ，详细请参考 Programming Reference Manual。（注意，在这里使用“录制代码”工具的作用很有限，因为在生成录制的代码时会与实时运行的批处理命令冲突。）

批处理扫描 vs 集群扫描

COMSOL Multiphysics 的所有许可证类型都可以使用 批处理扫描 功能。如果你有网络浮动许可证，则还可以使用 集群扫描 这一附加功能。这两个扫描功能很相似，不过 集群扫描 选项多了远程计算和集群配置等设置。集群扫描功能可以将大量扫描分散到一个（很可能是大型的）集群上。这种做法可以显著提高效率，因为独立扫描（也称作易并行计算）所涵盖的范围通常很广。如果你已熟练掌握了批处理扫描，那么运行集群扫描也不是什么难事。

下一步

想亲自动手尝试模拟这篇博客中讨论的模型吗?欢迎进入 COMSOL 案例库，下载示例模型:

• 电传感器
• 电传感器批处理

拓展阅读

阅读 COMSOL 博客，了解更多关于批处理扫描和集群扫描的文章：

• 混合计算系列博客
• 如何使用集群扫描节点 COMSOL Multiphysics^®
• 如何在 COMSOL Desktop^® 环境中运行集群
• 批处理扫描中任务并行的附加值

目录

1. 简介
2. 终端功率激励
3. 规定电流激励
4. 为什么损耗会随频率变化?
5. 电压激励
6. 关于绝缘的重要说明
7. 连接到电路接口
8. 将固定功率输入模型
9. 总结

简介

如下图所示，有一个同轴电缆模型插入金属腔体内，该腔体内包含相对介电常数为 50 ，电导率为 30mS/m 的有损电介质。

同轴电缆插入一个包含有损介电材料的金属腔范例中。

在上一篇博客中，我们已经对这个模型进行了深入开发，并介绍了如何使用几种不同的建模方法来求解这个模型，所有这些方法在低频范围内几乎给出了相同的结果，但在更高的频率下可能会有差异。基于这些观察结果，我们将在 1-100 MHz 的频率范围内运行模型，并比较电流 物理场接口与电磁波 物理场接口的求解结果。就像我们之前观察到的结果一样，在这个频率范围内，我们应该能得到非常相近的解。在这两个接口中，我们将求解有损电介质和绝缘体内部的场，并假设金属同轴导体内部的集肤深度足够小，以便可以通过边界条件对其进行建模。

模型示意图。

上图突出显示了边界条件。施加了激励的内导体用红色虚线表示，该导体的所有边界都是电连接的。外导体和导电腔壁以及它们之间的绝缘体边界也是电连接的。这里的目标是进一步证明可以使用电流 物理场接口进行这类模拟。 不过，要正确使用这个接口，我们确实需要花点时间来了解各种激励选项。

终端功率激励

首先，我们用终端 类型的终端 边界条件激励模型，这使我们能够指定功率和特征阻抗。这个边界条件模拟了通过已知特征阻抗的无限无损传输线(如 50 欧姆同轴电缆)连接到电源的情况。无限远处的电源提供了一个具有恒定的周期平均功率的正弦信号。该信号沿着无损传输线传播，进入电流 模型，部分反射回传输线，部分在模型内消散。反射分量沿着无限无损传输线传播回来，不与电源相互作用。指定的功率必须始终大于或等于零，并且该功率将在模型内耗散，反射回被终止的终端，或传输到另一个边界条件。

终端 条件必须与另一个规定模型中某处电势的边界条件结合使用，通常是接地 边界条件。接地的概念，我们已经讨论过了，这里需要强调的是，任何一组电连接的边界都可以定义为接地边界条件。在这个模型中，我们将选择将接地定义为外导体和腔壁以及两者之间的边界。也就是说，外部导体和腔壁是电连接的，并且都保持在零电位。

显示了 终止类型的 电流终端条件的屏幕截图，输入了传输线和电源的阻抗。

考虑来自沿传输线的电源信号：部分信号将反射回电源，电源包含一个绝缘体或环行器以保护设备。理想情况下，这种反射最小，这样我们就可以将尽可能多的能量储存到被加热的材料中。反射的幅度由通过该边界条件自动计算的 S 参数量化。

电源类型的 同轴型电磁波，集总端口的屏幕截图，电缆阻抗可以直接规定，也可以从几何结构中计算得到。

我们将其与电磁波 公式中的等效边界条件进行比较。这里，我们考虑的是电连接域的边界集合。在频域模型中，有两组这样的边界：与腔壁电连接的外导体和如图所示的内导体。它们在同轴电缆内部的环形间隙中被分开，在这个环形边界上，我们应用了同轴 类型的集总端口，并再次指定了电源和电缆阻抗。物理解释和之前一样，这个边界条件也会计算 s 参数。

比较使用带功率激励的电流接口和 电磁波接口计算得到的损耗。

我们可以根据上图中损耗材料内沉积的热量来比较这两种方法。图片显示了在大部分频率范围内，两种方法的计算结果具有极好的一致性，但在更高的频率开始出现分歧，这是因为电流 接口仅考虑了电阻和电容效应，而电磁波 接口还额外考虑了感应效应。

理解电流 接口中功率 类型的终端 条件和电磁波 接口中集总端口 条件的物理解释很重要。使用这些条件是为了模拟与无限无损传输线的连接。在无限远处有一个电源沿着这条线发送波传播信号，该信号一部分在模型内以热量的形式消散，一部分反射回无限传输线，继续传播且不会干扰电源。当然，在现实中，反射信号会产生干扰，而高频系统在阻抗不匹配的情况下会有明显的反射，因此源将配置绝缘体或环行器，用于防止信号反射损坏电源。在较低的频率下，特征波长比整个系统的尺寸长得多，我们不能再从波的角度考虑激励，因此选择其他激励条件。

规定电流激励

选择用电流 类型的终端 激励电流 模型，对应用了已知电流的源进行模拟。尽管在物理上仍然有一条传输线连接信号源和模型，但我们不再考虑任何信号被反射回信号源。请注意，在频域模型中输入电流的负数或复数值是合理的。一个复数值，或者一个负数，仅仅意味着具有相移的激励，这篇文章中对此进行了详细描述。这个电流将流过模型，然后流过接地 边界条件或任何其他规定电势的边界条件。

电流类型的 电流终端条件屏幕截图。

电磁波 接口内的等效条件是电流类型的集总端口，如下面的截图所示。

带 电流激励的 同轴类型的 电磁波，集总端口的屏幕截图。

通过电流激励的结果比较。

两种方法的比较如上图所示，可以观察到结果非常一致。与之前使用功率激励的方法明显不同的是，损耗随频率的变化更加明显，这一点值得特别注意。

为什么损耗会随频率变化?

对于一个电流驱动的问题，峰值电流，我们可以从等式计算整个模型内的周期平均损耗

式中，从终端的角度来看， 是系统的阻抗。阻抗是终端，电流 类型和电磁波，集总端口 边界条件自动计算的输出之一。下表显示了频率范围内的值:

表格中列出了由两种公式计算的作为频率的函数的阻抗值。

可以从这些数据中观察到，阻抗的实部随着频率的增长而下降，而虚部变得更小。也就是说，该零件在更高的频率下损耗更少，电容更大。

我们应该寻找一个简单的说法来解释这种趋势，所以让我们忽略所有的几何复杂性，把这个模型单纯地看作是位于两个电极之间的单一材料来考虑。最简单的，我们可以将这种情况看作一个并联 RC 电路，电阻由材料电导率决定，电容由介电常数决定。通过这两条平行路径的电流，即材料电导率引起的传导电流和介电常数引起的位移电流，将随频率变化。在更高的频率下，相对于传导电流，位移电流会更多，因此阻抗的复部与阻抗的实部的幅值比增大。

通过对阻抗及其随频率变化的了解，我们还可以解释为什么功率激励在频率范围内表现不同。相反，功率激励代表一部分反射出模型，一部分吸收到模型中的输入信号。反射由 S 参数给出。对于单端口系统，电缆阻抗可以通过 和 计算:

并且，对于单端口系统，系统吸收的功率为

对于该系统，由于电缆与被加热材料样品之间存在高阻抗失配，大量入射功率将被反射回来。 如果改为通过规定电流源激励，就不用考虑由于阻抗失配引起的反射。所有的电流都将以传导电流或位移电流的形式流过材料，只有传导电流会产生热损耗。

考虑到这些方程和材料属性，我们现在应该理解为什么不同的激励条件会表现出不同的与频率成函数关系的热损耗。

电压激励

接下来,我们看看电压 类型的电流终端。该边界条件在直流和非常低的频率下具有明确的含义，例如，家用电源插座中的 60Hz 线路电压。然而，在更高的频率下，我们需要更仔细地考虑如何从电流 接口和电磁波 接口获得相同的结果。

屏幕截图显示了 电压类型的 电流终端条件。

电磁波物理场接口没有直接等效在整个表面上规定电压的边界条件；您只能将电场设置在两个表面之间的边界上，这在较高频率下没有非常好的物理解释。我们将先研究这种方法，然后再介绍第二种方法。

与电磁波接口中的电流、电压边界条件最接近的模拟是电场 边界条件，应用在与之前方法中的集总端口 相同的环形边界。在这个电场 边界条件下，我们应用这个表达式来表示产生相同电位差的电场分布。对于同轴电缆(与 z 轴对齐)，电介质内的解析电场为

这个表达式在有限元网格上进行评估，并且将随着网格的细化变得更加准确。边界条件的设置和示例网格如下面的屏幕截图所示。

电磁波接口中的 电场边界条件，以及可视化网格的屏幕截图。

虽然这种在电磁波 物理场接口规定电场的方法，是电流 物理场接口的电压 条件的最直接的类比，但很难在物理上解释，因此考虑替代方案可能会有所帮助。

电压条件的近似值也可以通过电压源的诺顿等效和使用集总端口，电流 类型，以及一个沿同轴电缆放置在较短的距离处的一个电阻相对较小的集总单元。然后调整施加的电流，使该集总单元上的电势近似等于我们想要的施加电压。

电磁波接口的 集总单元边界条件的屏幕截图。

使用电压激励时的结果比较。

如上图所示，我们可以评估这些方法的损耗并进行比较。同样，在很宽的频率范围内，两个结果也具有很好的一致性。即使在 100MHz，结果的差异也小于 3。 我们还可以根据阻抗和施加的峰值电压计算损耗:

关于绝缘的重要说明

在上一篇文章中，我们注意到，在功率激励的情况下，同轴导体内部及其周围的绝缘体材料可以从电流 模型中完全省略。对于所有不同的激励类型，值得检查这个结果。我们可以重新求解前面的所有方法，只在有损电介质中求解电流，然后比较结果，如下图所示。可以看出，不管是否建立绝缘模型，功率激励和电压激励显示非常相似的结果。然而，当通过电流激励时，由于从模型中省略了位移电流路径，在更高频率下，结果明显有差异。

求解绝缘体和电介质中的 电流 与只求解电介质中的 电流 时损耗的比较。结果的差异在电流激励下很明显，但在电压或功率激励下可以忽略不计。

连接到电流物理场接口

到目前为止，我们看到的电压和电流源条件是理想化的，因为我们忽略了源和模拟的系统之间的所有中间结构。但通常情况下，我们确实需要考虑一些重要的电气结构。我们可以在模型中包含这些结构的详细细节，但也可以简单的通过电流 物理场接口来实现，并且可以与电流 或电磁波 物理接口相耦合。

例如，假设我们要模拟一个简单的过滤器，如下图所示。正在模拟的是一个与系统串联的电容器、电感器和电阻连接的电压源。电流接口是建立在节点上，并在每个节点上计算电势。这些节点之间的电流是根据节点之间的阻抗计算的，这可以是电路 接口的任何特征，或者可以通过另一个物理场接口中的外部耦合 特征计算。在这个示例中，我们希望将模型放置在节点 3 和 4 之间。

连接到系统模型的电压源、电容器、电感器和电阻的电路图。

电流 接口到电路 接口的耦合是通过外部 I 终端 特征定义的。这将把电路 接口内的电势节点与电流 接口内的一组表面电势耦合起来，正如通过电流 接口的终端 边界条件所定义的那样。使用这种方法，模型中应该只有一个接地 条件，并且应该在电路 接口内。

电路类型的 终端条件。终端选定表面的电势将与电路的一个节点 电流耦合。

外部 I 终端特征指定了 电路节点与 终端电势之间的耦合。

在用电磁波 接口连接电路 时，耦合是通过外部 I vs U 特征定义的。这个特征将两个电路 节点之间的电势差与电磁波 接口的电路 类型的集总端口 之间的电位差相耦合。

电路类型的 集总端口。这个电势差将与电路接口的两个节点之间的电势差相耦合。

外部 I vs. U 特征将电势差与 电路的两个节点之间的 集总端口的电势差 相耦合。

添加电路 接口需要额外的计算成本，因为它通常需要使用直接 求解器来计算全耦合 的方程组(详细请参阅我们的知识库中的文章“理解全耦合与分离求解方法，及直接与迭代线性求解器之间的区别”)。这增加了内存需求和计算时间，结果如下图所示，两种方法之间的差异可以忽略不计。注意损耗与频率图中的峰值是由集总电感引起的。

当电容器、电感器和电阻被串联在电压源和模型之间时，损耗与频率的关系。

将固定功率输入模型

最后，假设您要对一个包含反馈环路的源建模，该反馈环路可以调整输入来输入固定功率。有几种方法可以实现这一点。我们先简要介绍两个内置选项并讨论它们的局限性，然后介绍第三种更通用的方法，该方法需要更多设置但具有更大的灵活性。

在电流，终端 边界条件，功率 选项下，通过在模型中引入一个额外的耦合方程来实现反馈，该方程将调整施加的电流，从而使指定的峰值至峰值功率输入模型中。也就是说，对于频域模型，终端内指定的任何功率都将是域内耗散功率的两倍。这个额外的、完全耦合的方程会产生一个非线性方程组，它是非对称和非线性的，这意味着收敛取决于初始值、缩放比例和求解器设置，并且需要完全耦合的直接求解器。 因此，它的计算量很大，应该仅用于模型中存在非常强的非线性的情况，例如等离子体建模。

在电磁波 物理场接口中，有一个终端边界条件和启用主动端口反馈 选项。这个端口边界条件可用于任何类型的波导或传输线。 这个边界条件对计算的要求较低。但是这两种方法都有一个缺点：在所有可能引入损耗的域和边界条件上，它们将指定的功率存储在整个模型中。

假设我们只想将输入的功率指定到模型的一个域（或边界）中。在这种情况下，我们需要使用不同的方法，引入一个全局方程，迭代地调整输入来获得理想的输出，就像这篇博客中所讨论的。

简单来说，无论要求解的是什么物理问题，使用这种方法都需要三个步骤:

1. 使用一个通过全局方程 定义的变量指定激励。
2. 为全局方程 写一个更新的表达式，它将在求解过程中迭代更新激励的大小。
3. 确保使用分离式求解器，以使全局方程 作为一个单独的步骤进行求解。

我们将详细研究在电流 接口中施加电流激励的例子，因为其他情况在设置上非常相似。

首先，施加的电流不是常数，而是通过变量指定，如下面的屏幕截图所示:

屏幕截图显示了一个指定为变量的电流。

接下来，通过全局方程 接口指定施加的电流的更新方程，如下图所示。

通过一个额外的 全局方程执行反馈方程。请注意，单位设置为安培。

迭代更新方程为:

I_applied-nojac(if(intopSample(ec.Qh)==0,10[mA],I_applied*sqrt(50[mW]/intopSample(ec.Qh)))).

这个等式使用了一个名为 intopSample() 的非局部耦合，计算了样本域内损耗的积分。该等式首先检查积分损耗是否为零（初始状态），如果处于初始状态，则该等式将电流设置为 10mA，否则将计算的耗散与 50 mW 的目标耗散进行比较，并对施加的电流进行缩放。因此，这个表达式既包含电流表达式的初始化，也包含更新的表达式。 因此，有必要确保使用分离式求解器，并且来自全局 ODE 和 DAE 接口的变量先被求解，如下面的截图所示。第二个分离步骤，求解电流 接口，然后可以使用计算效率高的迭代求解器。

分离式求解器设置。首先求解初始化和反馈方程，然后使用默认的迭代求解器求解物理场。

使用这种方法，电介质中的损耗在所有频率下均为 50mW，我们可以绘制施加的电压或电流，来获得耗散。

在样本域内将产生指定耗散的终端电流和终端电压。

结束语

我们对用于频域电流建模的几种激励方式进行了介绍，并将它们与电磁波 物理场接口的类似激励进行了比较。文章重点讨论了在系统的 1–100MHz 频率范围计算的结果，但我们的结论可以扩展到低于或高于该频率范围，具体取决于所讨论的情况。重点是理解激励选项之间的差异以及如何使用它们。点击下方的按钮可以获得相关模型文件。

带铜线和电容的电路板模型设置

假设你有一块 1.5mm 厚的印刷电路板（PCB），它类似于一个三明治结构，其底部是导电的基底层，顶部为 200µm 厚的覆铜层，中间是电介质。在这块电路板顶层的导线之间，焊接了几个表面贴装电容，如下图所示。其中一条铜走线被激励，由于三个电容器的电容不同，信号将在输出端之间被分割。在 100kHz 左右的工作频率下，集肤深度将与铜的厚度相当，因此我们应该对铜导体的体积以及电介质和空气的体积进行建模，以正确捕捉传输线的行为。(关于传输线建模的更多细节，请参阅我们学习中心的文章《TEM 和准 TEM 传输线建模》）。我们不想建立包含表面贴片电容的真实几何形状的模型，因此只在传输线间引入附加的耦合电容，以观测施加信号的变化。

图为一块具有 200um 厚的接地底板、铜迹线和连接铜迹线的三块表面贴片电容的 PCB 板模型。PCB 板上的其他结构不在建模范围内，铜迹线是阻抗已知的信号传输线。

由于导体内部存在明显的集肤效应和邻近效应，因此适用于磁场 接口。为了在这个接口中对电容器进行建模，我们使用了集总元件 功能，并将类型设置为用户定义。这个功能应该施加在一个矩形面上，为我们想引入集总元件的区域内架起一个桥梁。这个面需要在几何序列中被勾画出来。

我们需要在集总元件 功能中指定三个基于几何形状的输入端口，集总元件的输入方向必须是一与导体平行的矢量（如下图所示）。也就是说定义器件中流过的电流方向与这个矢量方向相同。集总端口的高度为导电域之间间隙的长度，可用进行电场积分以及计算电压降。集总端口宽度方向为矩形面的正交方向。

图中显示了选中用户定义选项并指定集总端口高度、宽度和方向的集总元件功能的界面。

除了几何集总元件的几何特征外，我们还需要指定元件随频率变化的阻抗。我们可以在用户定义、电容、电感器、平行 LC、串联 LC、平行 RLC 和串联 RLC 等效阻抗等内置选项之间切换。用户定义 选项允许输入一个与频率相关的复数表达式，这样做可以将一个任意的等效电路添加到一个频域模型中。

按照这个模型的激励，我们可以认为是在模拟 4 条阻抗已知的、连接接地平板和铜迹线的信号传输线。为此，我们可以使用集总端口 功能，并将其类型设置为用户定义。在使用和功能上，集总端口 功能与集总元件 功能几乎相同，唯一的不同是它允许施加激励并监测S参数。

集总端口的类型也可以被设置为电路，这样就可以通过电路 接口将任意复杂的集总电路元件组合引入模型中。对于频域模型，使用电路 类型在功能上与使用包含用户定义 的频率相关阻抗的集总元件 功能相同。另一方面，对于时域模型，需要使用电路 接口添加集总电容或电感。

在设置方面，需要注意模型的边界条件。我们只对电路板和传输线的一小部分进行建模，并将假设周围环境对建模域不产生影响；也就是说，我们将忽略周围结构的任何串扰。我们选择将模型放置在一个更大的空气域内，沿其外侧使用理想磁导体边界条件，模拟一个绝缘的外壳。关于这一点的深入讨论，请参阅我们的博客 “如何为线圈建模选择边界条件“。

评估结果

在 100kHz 的工作频率下求解这个模型后，我们可以评估 S 参数并绘制导体中的电流，如下图所示。观察集肤效应和由集总元件引入的电容耦合的电流的拆分情况。由模型结果可知，我们已经在模型表面利用集总元件实现了电流通路。

铜导线中的电流分布图。

在更高的频率下求解

如果将工作频率提高到 10MHz 会发生什么。在这个频率下，集肤深度大约是迹线厚度的 1/10，所以不再需要对铜迹线的内部体积进行建模。我们可以在铜迹线所有边界上使用过渡边界条件。这样做是合理的，具体原因请参阅博客 “如何模拟时变磁场中的导体“。通过只求解空气和电介质内的磁场，我们以较少的总自由度解决了问题。现在有可能将集总端口和集总元件从用户定义 类型切换到均匀 类型。由于施加这些边界的面现在两边都是导电边界，均匀 类型的设置将自动确定端口的宽度、高度和方向。

在更高的频率下，不需要对导体的内部进行建模，可以使用 过渡边界条件，表征电流在导体的表面流动。

求解更薄的导体层

在讨论的最后，让我们来看看如何对较薄的铜线进行建模。随着铜线厚度的减小，用于划分网格的结果单元变得更小，这就增加了计算成本。可以通过使用过渡边界条件 避免建立铜迹厚度的几何模型，可以将过渡边界条件 施加在内部边界上，也就是说，边界上的场是在两侧求解的。这使得我们可以将迹线建模为一个几何厚度为零的边界，这进一步减少了模型的自由度，尽管我们可能希望在迹线上有稍微精细的网格。在这样的模型中，不包括由于迹线的有限高度而产生的电容效应，但我们可以合理地假设它们很小。

过渡边界条件的用户界面截图，它根据材料属性和厚度计算内部边界的损耗。

使用 过渡性边界条件的结果。电流在几何厚度为零的表面上流动。

结束语

这篇博客，我们介绍了如何使用集总单元 功能结合固体导体建模，建立一个等效电路元件。然后，讨论了另外两种情况：当集肤深度非常小时，通过过渡边界条件 对固体导体建模，以及通过过渡边界条件 对非常薄的固体建模，所有这些都与集总单元功能相结合使用。这些建模技术对任何模拟电路板或想在其电磁模型中包含集总电路元件的应用都很有用。值得一提的是，以上技术不仅可以在 AC/DC 模块的磁场 公式中使用，也适用于 RF 模块的电磁波、频域 公式。

编者注：在之前的博客“屈曲，当结构突然倒塌时”中，我们介绍了屈曲分析的各个方面，并讨论了一些更专业的屈曲分析技术。这篇博客提供的信息适用于 COMSOL 软件 6.0 版本以前的用户。

包括几何缺陷

众所周知，由于存在几何缺陷，一些结构对屈曲的承载能力会显著降低。在现实生活中，结构中总会存在一些制造缺陷。另外，施工过程中的安装不完善，或者结构可能因服役载荷效应而发生变形。因此，考虑缺陷很重要。

对此，我们有不同的方法。对于一个简单的结构，例如单个支柱，可以将它假设为其他相似的几何形状。如果要将缺陷包含在模型中，可以直接创建包含该缺陷的几何结构，也可以先构建一个完美的几何结构，然后添加包含指定变形 节点的变形几何 接口。

对于更复杂的结构，一般很难创建包含缺陷的几何结构。一种解决方案是进行线性屈曲分析，然后使用一种或多种屈曲模式作为缺陷。这样做的前提是，结构对屈曲模式本身很敏感。

请注意，在研究缺陷的影响时，通常不能使用线性屈曲分析。需要逐步增加载荷，直到结构失效。变形是渐进的，因此失效准则通常是最大允许位移或应力。

从 COMSOL Multiphysics 6.0 版本开始，有一种可以自动设置基于屈曲模式缺陷研究的方法。让我们来看看如何设置这种研究。

基于缺陷建立研究的 7 个步骤

第 1 步：从线性屈曲分析开始

我们先从构建一个理想结构的普通线性屈曲分析开始。如果要在缺陷形状中包含比第一个屈曲模式更多的屈曲模式，就需要更改待求的屈曲模态数(默认值为 1)。

更改待求的屈曲模态的数量。

步骤 2：添加屈曲缺陷节点

现在，屈曲模式和相应的临界载荷因子已经设置好了，下一步是在定义 下添加一个屈曲缺陷 节点。

左图显示了如何添加 屈曲缺陷节点，右图显示了新添加节点的 设置窗口。

步骤 3：输入模式编号及它们的比例因子

现在需要调整添加的屈曲缺陷 节点的设置，包括输入模式编号和比例因子。要确定合适的比例因子，首先需要知道屈曲模式的最大挠度。屈曲模式的大小是任意的，因此求解器在求解时会成比例地缩小屈曲模式的位移大小，以方便计算。默认情况下，模式会被缩小，使最大位移为几何图形对角边界框长度的 10^-6 倍。

实际的缺陷大小应反映真实结构的几何质量。或者，缺陷的大小可以通过一些设计规范给出。假设真实的几何形状尺寸可能与理想形状相差 2mm，该几何尺寸为 1m。如果使用单屈曲模式，比例因子将为 2000。如果使用多种模式，那么为所有的模式都分配 2000 的比例因子可能是保守的，因为这些值会相加。这是不可忽略的，因为某些模式可能会在相反的方向上起作用。我们可能需要检查模式，甚至为其中一些模式分配负比例因子来获得预期的形状。

从上到下，前三个图显示了 Euler 2 悬臂梁的前三种屈曲模式，底部的图显示了这些模式的纯叠加。标记显示了最大垂直位移。（您可以通过从案例下载页面下载模型的 MPH 文件 euler2_buckling_imperfection.mph，尝试自己建模。）

步骤 4：配置变形几何体

接下来，我们来设置变形几何。首先，单击屈曲缺陷设置中的变形几何 上部的配置 按钮。

创建变形几何结构。

这样，特殊配置了指定变形 的变形几何 节点将被添加到模型开发器树中。

新添加的 指定变节点。

如果要更改包含的模式、比例因子或重新运行线性屈曲研究，则无需重复这个步骤。

步骤 5：创建加载参数

添加一个将充当载荷的乘数的参数，然后将该参数插入到用于线性屈曲分析的所有载荷特征中。

添加了一个新的参数 lf，并用作所有载荷的乘数。

接下来，在载荷参数 拉菜单中选择 If（载荷因子） 参数。

载荷因子已经选择好了。

步骤 6：配置新的屈曲研究

要为包含缺陷的增量分析创建研究，请单击屈曲缺陷分析中的非线性屈曲研究 上部的配置 按钮。

创建非线性屈曲研究。

现在，我们已经创建好了一个新研究，并进行了一些特殊设置：

• 包括几何非线性。
• 使用选定的载荷参数进行辅助扫描。
• 基于计算屈曲模式的最低临界载荷因子建立参数值列表。使用这些设置后，最大载荷比线性屈曲研究的预测值高出约 10%。

新研究的设置窗口。

步骤 7：运行研究

该研究在较高的载荷水平下可能无法收敛。这通常意味着到目前为止已经远远超过极限负载了，但中间步骤仍会被存储并可用于计算。

对于非线性屈曲分析，没有唯一的临界载荷。通常，需要根据载荷参数绘制相关的位移和应力曲线，并基于结构的物理特性使用失效准则。

下列示例与用于表示上述三种屈曲模式（长度为 2m 的悬臂梁）的 Euler 2 示例相同。假设最大允许的尖端位移为 10mm，最大允许的 von Mises 等效应力为 400 MPa。

显示载荷参数与位移和应力关系的曲线图。屈曲分析产生的临界载荷因子对水平轴进行了归一化处理。图中的标记点分别显示了失效准则的位移为 10mm，应力为 400 MPa。

下表显示了选择不同屈曲模式缩放时的临界载荷值的比较。

可以看出，中等水平的缺陷对实际选择生成模式的敏感性是有限的。

这个例子很简单，而且 Euler 支柱对缺陷不是很敏感。壳结构在这方面通常问题更大，所以接下来，我们来探讨一个壳示例。

一个更高级的示例

在这个示例中，考虑一个直径为 1.5m、高度为 2m 的钢圆柱体，其壁厚度为 10mm，并有两个厚度为 20mm 的加强环。使用轴对称壳单元进行分析。（您可以通过 COMSOL 案例下载页面下载相关的 MPH 文件 cylindrical_shell_buckling_cleared.mph，来探索这个模型。）

几何形状和载荷。圆柱体的下端是固定的。

许多屈曲模式具有相似的临界荷载因子。

前 5 种屈曲模式对应的临界荷载因子。

在这种情况下，检查不同的缺陷模式很重要。可以通过在模式编号上添加外部参数扫描来自动进行这项操作，如下图所示。

屈曲缺陷节点中的 模式选择表格。

关于这个表格，需要考虑以下信息：

•  将类型（mode==1）的布尔表达式作为过滤器，来检查哪种模式将使用缺陷。 mode 是扫描参数。
• 设计比例因子，使每种模式的峰值都变为 1 mm。
• 使用最大算子 maxop1（在非局部耦合 → 定义 下添加）获取每个屈曲模式在径向方向上的最大位移。
• withsol 算子的目的是从特定解中选取结果。在这个示例中，它被用来检索单个屈曲模式。（阅读 COMSOL 学习中心的文章 “withsol 算子的示例”，了解关于这个算子的更多信息。

现在，非线性分析包括将前 4 种模式的外部扫描作为缺陷，以及每种模式的载荷都会增加的相同的内部辅助扫描。由于任何非线性解都可能无法在较高载荷下收敛，因此必须设置参数化扫描，使整个分析不会因此而失败。

非线性研究和 参数化扫描节点。

在 参数化节点中， 必须将误差选项设置为 存储空解。

现在，我们可以绘制所有情况下的最大应力和位移的变化。

显示了位移和应力归一化处理后施加载荷的曲线图。

在这个示例中，允许的应力在所有情况下都是在相当低的外部载荷下被超过的。在线性区域已经能够看到这种情况。最终的结果是，这种结构将在屈曲发生之前由于塑性坍塌而失效。有可能通过在模型中加入塑性来完善分析。

从这个示例中，我们可以看到，使用非线性分析的变形很大程度上与实际的变形一致。

不同的缺陷在接近失效载荷时的变形（缩小了 5 倍）。颜色表示径向位移。

必须注意的是，不能完全假设轴对称壳的所有屈曲模式也是轴对称的。真正的第一屈曲模式不是轴对称的，看起来像这样：

使用完整的三维公式时的第一种屈曲模式。

活荷载和静荷载的组合

线性屈曲分析中的载荷因子可以被认为是相对于所施加载荷的一种安全系数。有时，只有某一组载荷可以变化，其他载荷具有明确定义的值，比如自重就是一个典型的例子。如果假设结构不会仅因重力而失效，那么需要回答的问题是：如果考虑部分承载能力已经被自重利用，服役载荷的安全因子是多少？

本文中将不变的荷载称为静荷载，会发生变化的载荷称为活荷载。 COMSOL Multiphysics 6.1 版本新增了分析活载荷与静荷载共同作用的功能。

需要注意的是，不可能首先计算静载荷利用了多少承载能力，然后相应地减少允许的活载荷。假设有两个独立的荷载， 和 ，并且对于每个荷载都有一个单独的屈曲临界值 和 。如果施加的载荷 是这两个载荷的线性组合，那么当 时，临界状态通常不会发生。

在仅具有活荷载的普通线性屈曲分析中，我们需要有任意活荷载水平的静荷载工况，通过特征值分析计算临界荷载因子和相应的振型。这非常简单，当添加线性屈曲研究时，研究序列就会生成。如果要进行也包括静载荷的类似分析，需要两个稳态的研究步骤 — 从这些步骤中得出的结果需要由特征值求解器以不同的方式加以考虑。这样的研究可以通过不同的方式建立。接下来，您将看到我们建议的工作流程的步骤。

设置一个包含静载荷研究的 7 个步骤

第 1 步：添加研究

正常添加线性屈曲 研究。

步骤 2：定义活荷载

像往常一样，用一个任意的荷载水平创建活荷载。

步骤 3：定义静载荷

使用载荷的实际值创建静载荷。最好同时选中载荷特征的线性屈曲 部分中的视为静载 复选框。严格来说，只有当一个载荷是随动型的时候才需要这样做（取决于变形）。如果要在载荷特征中显示线性屈曲 部分，请确保选中显示更多选项 下的高级物理场选项。

添加一个静载荷。

第 4 步：添加额外的研究步骤

在线性屈曲 步骤之前再插入一个稳态研究 步骤。

添加第二个稳态研究步骤。

第 5 步：在研究步骤中停用活荷载

在两个稳态研究步骤中的其中一个步骤中，我们需要同时分析两个载荷的影响。这与模型树的设置一样。

在另一个步骤中，应该只分析静载荷。这意味着我们必须禁用描述活载荷的所有载荷特征，进行这项任务最方便的是使用修改研究步骤的模型配置 选项。

图中 载荷1 已被禁用。

步骤 6：选择两个稳态解进行屈曲分析

在线性屈曲 研究步骤设置中，需要指定两个解。为此，首先需要为研究运行显示默认求解器。线性点是仅具有静载荷的研究，活荷载的解包含活荷载和静荷载的总和。

选择两个解。

第 7 步：运行全部研究

您可能想知道为什么两组载荷的静态解不是单独计算的。这是因为如果静态解是非线性的，那么给定的方法在静荷载解给出的应力状态下为屈曲分析提供了一个更加精确的线性点。

有关这类分析的示例，请查看静载作用下的桁架塔线性屈曲分析模型以及相关的应用程序文件。您可以在下文中找到这个模型和相关设置：

第一种屈曲模式。

在这个示例中，有两个效应会影响静载荷。除了自重外，支撑线中还有一个预紧力，会在塔的下部引起压应力。

结束语

最新版本的 COMSOL Multiphysics 能够建立高级屈曲研究，并且操作简单。单击下面的链接，进入 COMSOL 案例下载页面，尝试自己模拟文中讨论的模型：

• 静载作用下的桁架塔线性屈曲分析
• 具有初始缺陷的屈曲

背景知识

卫星热设计是一个具有挑战性的问题。一颗卫星是由许多对温度敏感的部件组成的。传感器、照相机、收音机、电子器件、电池、姿态控制系统和太阳能电池板都需要保持在一定的温度范围内，甚至卫星结构本身也必须保持在一定的温度范围内，以防止过度的热变形。许多组件都会散发热量，卫星也会经受来自环境的多种不同的红外(IR)热载荷。设计一颗卫星需要了解如何最好地将所有这些热量辐射出去，并使卫星保持在理想的工作条件下。

各种电子元件产生的热量通常很容易定义，但环境热载荷可能出奇地复杂。首先，在面向太阳的任意表面上，有直接入射的准直太阳光通量。其次，对于近地轨道上的卫星，入射到地球日光侧的太阳光通量会被漫反射到卫星朝向地球一侧的表面。这些反射的大小取决于地球的局部表面特性以及不断变化的大气条件。总的来说，漫反射太阳光通量大约是直接太阳光通量的三分之一，被称为反照率通量。当卫星进入日食时，这些直接的太阳光通量和反照率载荷下降为零，但有一个第三环境热源始终存在。地球是温暖的，相当于一个漫射器，其红外辐射的大小是纬度和经度的函数。

知道这些随时间变化的环境通量，以及它们在卫星表面的分布，是计算卫星温度所需要的输入，这涉及到求解固体部分的热传导和所有暴露表面的辐射。通常将这些环境通量分为两个波段：太阳波段和环境波段。这样做的原因是，太阳温度在 5780K 左右，主要发出短波长的辐射，而卫星和地球温度都在 300K 左右，主要发出长波长的红外辐射。这种划分是很重要的，因为为了热管理，卫星外部涂层的表面吸收性能往往是专门定制的，是波长的函数。例如，为了使卫星的工作温度尽可能低，一种方法是使用在太阳波段具有低吸收率(发射率)但在环境波段具有高发射率的表面涂层。

卫星在轨道上经历的辐射传热取决于它相对于地球和太阳的位置和方向。地球图片来源：可视地球和美国国家航空航天局。

一颗质量恒定的卫星围绕着一颗质量大得多的球形行星运行时，会沿着一条椭圆的路径，这条路径可以用开普勒轨道根数来描述，它描述了一个周期性的轨道。轨道根数被用来计算（通过中心方程）卫星在赤道坐标系（ECS）中的坐标。

在了解卫星随时间变化的位置的同时，也有必要了解它的方向性。这首先要通过一组卫星轴来定义卫星坐标系。根据任务参数，这些卫星轴是朝向特定方向的，如朝向地球、太阳、前进方向，或天球上的一个固定点。有时候也会想要改变这些轴的定义和它们的方向，例如为了使仪器朝向一个地理位置。卫星也可能围绕一个或多个轴缓慢旋转或相对快速地翻滚。方向的变化会影响热载荷以及卫星的阴影。例如，卫星上的一个突出物，如太阳能电池板或仪器，将对它后面的表面产生阴影。如果有旋转的太阳能电池板或其他铰接元件，这也将改变阴影和热载荷。另一方面，如果卫星在快速翻滚，这意味着环境热载荷是均匀的。

一旦获取了所有的任务参数，就有可能计算出所有的环境热载荷，然后就可以很直接地计算出卫星在一段时间内的温度曲线。现在，让我们看一下轨道热载荷 的用户界面，探索在设置卫星热分析时它是如何帮助输入数据的。

轨道热载荷用户界面概述

轨道热载荷 接口的工作方式与 COMSOL Multiphysics 产品库中的其他接口类似，并使用一致的工作流程。您可以从以下两个方面入手：

• 使用核心 CAD 建模功能或设计模块在软件中创建自己的结构 CAD 描述
• 导入一个 CAD 文件，例如 Parasolid^®、ACIS^® 或 STEP 文件

• 使用 LiveLink for CAD 产品之一，将 COMSOL^® 与您的 CAD 平台双向连接。

从这里开始，您将使用与其他工作流程类似的 CAD 简化和清理以及网格划分、求解和结果评价功能。因此，如果你已经是 COMSOL 用户，你将很快掌握这个新接口。

作为典型工作流程的一部分，轨道热载荷 接口分三个步骤使用，对应于三种不同的研究类型。首先，定义轨道根数、卫星轴和方向，并使用轨道计算 研究步骤解决一个或几个轨道周期的问题。这让我们在通过轨道热载荷 研究步骤计算所有环境来源的辐照度之前，快速验证任务参数。一旦这些辐照度被求解并存储，它们就被用来计算结构的温度和所有暴露表面之间随着时间推移的表面对表面辐射，使用轨道温度 研究步骤。对于每个轨道的环境热载荷相同的情况，可以只计算一个轨道的热载荷，并在热分析中及时重复计算。

轨道热载荷 接口可以单独用于计算环境热载荷，但更多的是与固体传热 和事件 接口一起使用。固体传热接口计算卫星固体结构内的温度分布，而事件 接口则跟踪日食和调整方向，以及任何其他状态的瞬时变化，如仪器的开启和关闭。

截图显示了轨道热载荷接口，以及三种相关的研究类型。

该接口的设置如上图所示，其与表面对表面辐射 接口的设置类似，但默认使用太阳光和环境 双波段模型。辐射总是存在于所有暴露的表面之间，包括卫星的外部和内部。

现在，让我们看一下接口内的默认功能。

首先，太阳属性 功能定义了入射的太阳光在 ECS 中的方向，以及太阳光通量。太阳方向 有四个预定义的选项。这些选项是根据地球的二分时刻和至点来定义的，这将控制太阳光矢量和太阳光通量。你也可以输入你自己的太阳光矢量，使它成为一个随时间变化的表达式。如果使用多波段光谱模型，太阳可以被看作是一个黑体发射器，在每个波段内有一个确定的通量，或者通过通量的分布作为波长的函数。通常把太阳当作黑体发射体就足够了，这是默认行为。

太阳属性特征指定了太阳光矢量的方向和太阳光通量。

接下来，行星属性 功能指定了计算轨道和日食所需的几个参数。开始时间的行星经度 使卫星下的行星定向，当以地面上的位置指定轨道机动时，或者当行星属性是纬度和经度的函数时，这很重要。辐射属性 部分可以用来启用反照率和行星红外载荷，反照率和行星红外通量 都可以成为经纬度的函数。这些数据可以从电子表格或图像数据中读入。行星是通过离散化处理的，即行星的可见部分被细分为一组具有相等角系数的斑块。当反照率或行星的红外特性变化很大时，以及对于低空轨道，可能需要更精细的离散化。

行星属性功能在启动时将行星定位在卫星下方，并描述行星的辐射属性。

下一个默认功能，轨道参数，提出了在六个开普勒轨道要素方面定义轨道的设置：半长轴、离心率、倾角、升交点经度、近心点角和起始时间的真近点角。圆形、赤道椭圆形 和赤道圆形 轨道也可以用一组减少的参数来定义。

轨道参数功能是用来输入轨道单元的。

航天器轴 特征定义了卫星坐标系的轴方向。这些轴可以在 CAD 坐标中指定，也可以选择卫星的一个面，在这种情况下，使用面的法线方向。当把仪器指向一个特定的方向时，这很有用。选择的次要指向方向不需要与主要方向垂直；采用次要矢量对主要方向的法线平面的投影。第三个轴完成了一个右手坐标系。可以定义不同的坐标系，与航天器方向 特征一起使用。

航天器轴 特征定义了主要和次要的指向性方向。

航天器方向 特征定义了航天器的主轴和次轴方向，以及围绕主轴、次轴和第三个方向的旋转（如果有的话）。方向可以是天顶/天底、太阳/与太阳相反、速度/与速度方向相反、轨道法向/与轨道法向相反、天体点或地面站中的任何一个。卫星将精确地朝向主轴方向并旋转，以使航天器的第二指向方向指向第二个方位方向。

航天器方向特征与 航天器轴特征一起，定义了卫星在一段时间内的定位方式。

如果只有一个航天器轴 和航天器方向 特征，那么这些特征将在整个分析中使用。可以有这些特征的多个实例，并可以在它们之间进行切换，以引入各种轨道机动。为了在各种定义的组合之间进行切换，使用了生成事件接口 特征，它允许轨道机动的序列，例如在卫星处于射程内时指向一个特定的地理位置。

事件接口 还有另一个目的，在所有用例中都是持续的：追踪日食。卫星进入和离开日食的时间(如果它们发生的话)被用来向求解器表明热载荷正在发生变化。

生成事件接口特征将填充事件接口中的 轨道热载荷事件节点。始终考虑日食。

地面指向特征可用于设置额外的事件，根据不同的条件将卫星指向某个地理位置。

接口中的其余特征，即漫反射表面、初始值 和不透明度 特征，都与各种建模表面的发射和反射的建模有关。

从热建模的角度来看，一旦环境热载荷被计算出来，工作流程与任何其他涉及热传导和辐射的模型是相同的。轨道热载荷 接口解决了辐射传热问题，并与固体传热 接口相耦合，后者考虑了卫星固体结构内的热传导，也允许你定义体积或表面内的热载荷，这些热载荷可以是恒定的或时间变化的。除此之外，你还拥有传热模块的完整功能集，包括：

• 薄壁部件的传导传热
• 交界面处的接触热阻
• 相变材料
• 流体中的对流传热
• 集总热连接和组件

一旦完全定义了问题并解决了轨道热载荷以及温度问题，将得到一组默认的图表，也可以创建任何数量的其他数据的可视化图表。让我们看一下其中的一些…

显示围绕地球的轨道、太阳光矢量和卫星方向的图。地球图片来源：Visible Earth 和 NASA。

所有来自环境源的总辐照量，以及太阳矢量随时间变化和天底方向的图。

显示一个小卫星内几个部件的温度的图。一些外部表面被隐藏起来。

在几个轨道上收敛到周期性状态的几个组件的温度图。

结束语

有了新的轨道热载荷 接口，现在就可以快速建立一个在轨卫星的热模型，用于预测运行性能。该接口对于从事卫星规划和设计的工程师来说是一个很好的工具。对于那些想马上开始使用传热模块的用户，请看下面的例子：

• 轨道计算
• 轨道热载荷
• 航天器分析

LiveLink 是 COMSOL AB 的商标。Parasolid 是 Siemens Product Lifecycle Management Software Inc. 或其在美国和其他国家的子公司的商标或注册商标。ACIS 是 Spatial Corporation 的注册商标。

在这篇博客中，我们将介绍表达式运算符 的工作原理，以及在使用 COMSOL Multiphysics^® 软件时如何在你的建模项目中使用这项功能。

表达式运算符

表达式运算符 节点是自 COMSOL Multiphysics 5.5 版本起引入的功能，右击模型开发器 中的 全局定义 或定义 节点，可以从变量实用程序 子菜单(通过选择显示更多选项 对话框中常规 下的变量实用程序 复选框激活)中获得。(当添加到全局定义 下时，表达式运算符 在整个模型中有效)。表达式运算符 节点的设置包括一个定义 部分，可以在其中定义运算符。下面是具体的操作步骤。

步骤 1：输入默认表达式

在表达式 字段中，输入默认表达式，该表达式在其选择上定义了表达式运算符（或者在全局定义 下添加时全局定义）。表达式可以包含任何常量、变量、函数和运算符，它们在要计算表达式运算符的上下文中是有效的。此外，它应该至少使用一个输入参数。请注意，默认表达式可以被子选择上的运算符贡献 节点所覆盖。

步骤 2：定义输入参数

在表达式 字段下面的表中，定义输入参数。在参数 列中键入参数的名称，并在维度 列中键入用于定义其维度 的单位;在参数类型 列中，从下拉列表中选择表达式 选项(默认设置)或值 选项。通过值传递输入参数意味着在将它作为数字插入运算符表达式之前，在调用表达式运算符的上下文中对它进行计算。由表达式传递的参数直接被替换到运算符表达式中。如果将参数作为参数传递给在不同上下文中计算其参数的耦合运算符，这可能会有所不同。例如，在一个参数中传递空间坐标变量时，必须小心，该参数被转发给作为表达式运算符定义一部分的积分 耦合。

步骤 3：输入名称

在名称 字段中，为表达式运算符输入一个合适的名称。它将是你在 COMSOL 模型中引用表达式运算符时使用的名称。

表达式运算符的优点

表达式运算符 对各种建模项目和方案都有好处。下面列出了一些具体的内容：

• 一旦定义表达式运算符并为它确定了描述性名称后，就可以在 COMSOL 模型中使用它来隐藏复杂的数学表达式。这意味着可以用速记的方式书写的表达式代替复杂的表达式。

• 表达式运算符 节点可以使用运算符贡献 子节点在模型几何图形的不同部分定义不同的表达式。可以将 运算符贡献 子节点添加到组件 节点下定义的表达式运算符 节点。对于每个这样的子节点，可以定义几何实体级别和选择使用它的几何实体。然后，可以在运算符贡献 子节点的表达式 字段中输入新的运算符 表达式。例如，表达式运算符 节点这一功能是解决多层壳接口后处理困难的好工具。对于这类后处理，能够使用单个运算符是有好处的，该运算符可以根据不同的选择扩展为不同的运算符。然后，它作为一个包装器，将合并不同域或其他几何实体上的不同运算符和表达式。

• 表达式运算符 还可以用来创建混合函数，这些函数混合了本地计算的参数和变量。

表达式运算符与分析函数

你可以把表达式运算符看作是一个参数化的变量，而不是一个分析函数。与分析函数相反，表达式运算符可以定义为仅在某些域中处于活动状态。根据定义，表达式运算符可以使用任何被定义的能计算出结果的表达式。

相反，分析函数应该是其参数的纯函数。由于它的全局性质，分析函数可用于定义模型参数;而表达式运算符不能。

使用表达式运算符的示例

这个示例演示了如何将表达式运算符 节点用于以下两个目的：

1. 根据表达式运算符的计算位置，将表达式作为参数替换为不同的运算符。在本示例中，它被用于绘制双镜对称。
2. 通过混合函数混合参数和变量，例如本示例中的运算符 phase(kx,ky)。

在这个示例中，固体力学 接口仅在模型几何图形的第一个象限(右上角)，也就是单位圆中定义。两个广义拉伸 运算符， yAxisMirror 以及 xAxisMirror，分别为 –x 和 –y 坐标，x 和 –y 坐标提供从源到目标的镜像操作。然后使用三个运算符贡献 子节点定义一个名为 doubleMirror 的表达式运算符，每个子节点分别在第二、第三和第四圆象限中提供双镜像转换 yAxisMirror(expr) 、xAxisMirror(yAxisMirror(expr)) 和 xAxisMirror(expr)。下图显示了第三象限的运算符贡献 子节点的设置：

第三象限中的 表达式运算符的定义。

下图显示了有效 von Mises 应力超过某个极限的区域：

超过应力极限的区域绘图。表达式运算符提供的双镜包括整个圆内的解决方案。

由于使用了表达式运算符，我们只需一个简洁的表达式，第一象限计算的解就可以在整个圆中被可视化为双镜图。

对于表达式运算符 的第二种使用方式，在整个几何体中定义了一个叫做 phase 的表达式运算符，并将参数和变量混合在一起，进行局部计算。表达式运算符 被定义为 exp(-i(kx*x+ky*y))，其中，kx 和  ky 是两个输入参数，前者代表波向量的 x 分量，后者代表 y 分量。它们也被定义为表达式，如下图中的设置窗口所示。

包含两个表达式参数的 phase 表达式运算符的 设置窗口。

相位的表达式 等效于使用 COMSOL Multiphysics 语法的数学表达式 。

你可以在曲面图中使用 phase 表达式运算符 来显示具有不同参数值的相位。下图显示相位 kx 设置为 5 ，相位 ky 设置为 20：

曲面绘图显示了使用 表达式运算符 的相位。

后续步骤

这篇博客解释了什么是表达式运算符 节点，以及如何在 COMSOL Multiphysics 中建模时利用它来发挥自己的优势。我们鼓励你练习使用这个功能来简化和改进你的 COMSOL 模型。如果您对这个主题有任何疑问，请通过以下按钮联系 COMSOL。

延伸阅读

想要了解有关 表达式运算符 节点的更多信息？查看我们的博客文章获取统计信息：计算标准偏差和其他统计数量，其中表达式运算符 用于简化在多个位置中使用的标准偏差的表达式。