为何PCB光板过回流也会发生翘曲？
在SMT产线上，许多工程师都曾遇到过这样一个“反直觉”的现象：一块没有任何元器件贴装的PCB光板，空板过回流炉，竟然也会出现明显的翘曲，翘曲量最高可达0.5%~1.5%，严重时甚至导致后续贴装对位偏移、连接器装配困难等问题。
为什么没有搭载BGA芯片的单纯PCB光板也会“笑脸”或“哭脸”呢？（搭载了BGA芯片的PCBA回流翘曲参见上一篇：PCBA回流翘曲仿真)
答案的核心就隐藏在PCB内部——层间铜分布不均匀与不对称。多层板中各铜层的铺铜面积、走线密度、铜箔厚度往往差异明显，一面大面积整铜平面，另一面走线稀疏，这种结构上的“先天不对称”，在回流高温环境下形成了显著的热膨胀系数失配和内建弯矩，就像板子内部藏着一组不断弯曲的“隐形双金属片”。
传统解决方式只能靠反复烧板、换板材、加治具等方式试错，不仅周期长、成本高，而且难以找到根本原因。而有限元仿真（FEA）为我们提供了一条高效路径：它可以在电脑里完整复现回流焊的整个过程，定量计算出铜分布不对称究竟产生了多大的翘曲量、哪几层铜是主要贡献者，以及不同优化方案下的改善效果。
本文将聚焦于“层间铜分布不均匀”这一核心诱因，从物理机理出发，带你一步步建立PCB光板回流翘曲的有限元模型，输入典型材料参数和回流温度，揭示翘曲的量化规律。
首先需要获取PCB的设计ECAD文件，此文件包含PCB的叠层信息，然后根据这个文件建立PCB的详细叠层几何。在此以一块简单的双面板为例，PCB的T面和B面铜层分布如下图所示。
PCB几何特性高度匹配
PCB本质是一块规则的矩形多层平板结构，长宽方向尺寸远大于厚度，层间界面相互平行。笛卡尔网格与板子的正交几何特征完全贴合，能最大程度减少网格扭曲和畸变，保证单元质量极高。 便于精确描述铜分布
划分笛卡尔网格后，可以根据单个单元内的材料分布，赋予不同的材料参数，通过将网格细化，尽可能准确的描述铜层的分布。 计算效率与收敛性显著优势
相比四面体非结构网格，六面体笛卡尔网格单元数量更少、雅可比矩阵更规整、求解速度更快。在同样的计算资源下，能划分出更细的网格，从而更准确地捕捉Z向层间剪切应力。 工程验证经验支持
行业内大量案例表明，对于PCB这类薄板大面积结构的翘曲预测，高质量的笛卡尔网格在精度和效率上通常优于自由网格，尤其在评估整体“笑脸/哭脸”形态和最大翘曲量时表现突出。
当然，笛卡尔网格也存在一定局限（如对复杂走线细节无法完全建模），但是可以通过将网格划分得更精细来降低其不利影响。总之，将PCB划分为结构化笛卡尔网格既可以极大简化PCB的六面体网格划分工作量，又保证了计算精度。设置平均网格尺寸为0.1mm后，划分的PCB结构笛卡尔网格如下图所示。
模型中共有790万一阶六面体单元，在划分网格的同时，自动计算每个单元的材料，并分配到不同的body里面。铜层的模量为110GPa，CTE为16.4ppm。而FR4的材料则需要考虑材料的各向异性，如下图所示。
完成网格划分和材料赋予后，下一步就是施加温度载荷，由于本文是分析铜层分布不均导致的翘曲，所以不考虑回流过程的温度分布不均匀，施加的是全局均匀载荷，初始温度为20℃，最终温度为220℃。并且为了避免出现刚体 位移，设置惯性释放。最终计算的翘曲结果如下图所示。
在A/C处，T面铜层存在由于线路导致的间隙，而在B面是完整的铜箔，考虑到本文设置的铜箔CTE大于FR4的CTE，在升温过程中铜箔膨胀程度更大，导致变形呈现T面凹B面凸的变形模式。而在B处，则正好相反，T面是完整铜层，B面存在线路，导致在这个局部呈现T面凸B面凹的变形模式。在整体层面，由于T面的残铜率较B面低10%，所以整体变形模式为T面凹B面凸的变形模式，在Z向上最大位移值为0.36mm，最小位移值为-0.24mm，整体Z向翘曲量为0.6mm。
至此，我们完整走完了从“为何PCB光板过回流也会发生翘曲？”到“如何用有限元仿真量化翘曲”的全过程。
核心结论非常清晰：PCB光板回流翘曲的根本原因，并非来自外部元件，而是板子内部层间铜分布的不均匀与不对称。而有限元仿真为我们提供了一把利器——它能在投产前就定量预测翘曲量、判断主因层位，并帮助我们通过优化铜平衡、调整叠层、选择合适板材等方式，从源头上把问题压制住。
仿真无法完全替代产线验证，但它能让我们把大量试错工作从实物板转移到电脑里，大幅缩短开发周期、降低报废成本。这正是仿真真正的价值所在。在设计之初就把铜分布对称性做好，远比后期烧板调治具来得高效和经济。
各位工程师在实际工作中有遇到由于PCB本身翘曲导致的可靠性问题吗？欢迎留言分享。我们下一期再见！（下一期我们分享：划分精细的PCB笛卡尔网格数量太多，计算量太大，有没有降低计算量又准确计算翘曲的方法？）
