简介：这份资源提供了一套用MATLAB读取Patran Bdf文件数据的m文件工具，面向从事有限元前后处理、需要对Nastran Bdf模型进行批量解析与二次开发的工程师和研究人员。其核心功能包括：将Bdf文件逐行读入cell数组保存，自动剔除空白行与注释行，并按每列8字符的固定格式切分为10列，进而识别各类卡片内容并归类存入名为output的结构体，便于后续提取节点、单元等信息。压缩包为rar格式，共2个文件，包含1个m脚本与1个bdf示例文件，整体约18KB，体积轻量，可直接在MATLAB中运行调试。目前已有928人学习关注，说明该工具在Bdf数据处理场景中具有一定实用参考价值。借助这份脚本，读者可以省去从零编写解析逻辑的时间，快速获得可复用的读取框架与分类思路，并对照示例文件验证结果，适合作为Bdf数据自动化处理的基础工具。
1. 当 Patran 的 BDF 文件在 MATLAB 里变成一堆乱码时
如果你做过有限元前后处理，大概率遇到过这种场景：Patran 导出的 BDF 文件用文本编辑器打开看着挺规整，一拖进 MATLAB 用readmatrix或者importdata读，要么报错，要么读出来一堆 NaN 和错位的字符串。这不是你操作有问题，而是 BDF 这种格式本身就不是给通用表格读取函数设计的——它是 Nastran 的卡片式输入格式，固定字段宽度、大小写敏感、注释和续行规则都很讲究。这个 MATLAB 读取 Patran BDF 文件的 m 文件项目，解决的就是这个痛点：把 BDF 里的节点坐标、单元连接关系、材料属性、边界条件这些数据，按卡片类型解析成 MATLAB 能直接用的矩阵和结构体。适合正在做有限元二次开发、需要批量处理多组 BDF 模型、或者想把 Patran 前处理结果接入自研 MATLAB 求解器的开发者。下面按实际拆项目的顺序，从文件结构讲到跑通再到改造。
2. 先看清 BDF 的文件骨架：卡片、字段与续行规则
2.1 BDF 不是 CSV，别用读表格的思路去读它
Patran 导出的 BDF 文件本质上是 Nastran 的 Bulk Data File，每一行是一条“卡片”，卡片有固定格式。最常见的几种卡片包括GRID（节点）、CTETRA/CHEXA/CQUAD4（单元）、PSOLID/PSHELL（属性）、MAT1（材料）、SPC1/FORCE（边界与载荷）。每条卡片的前 8 个字符是字段名，后面按 8 字符一个字段切分，但 Nastran 允许“短字段”和“长字段”两种写法，短字段就是每 8 字符一个，长字段是逗号分隔。Patran 默认导出短字段格式，所以你会看到类似这样的内容：
GRID 1 0.000000 0.000000 0.000000 GRID 2 10.000000 0.000000 0.000000 CTETRA 1 1 2 3 4注意GRID后面第一个字段是节点 ID，接着是坐标系编号（通常 0），然后是三个坐标值。CTETRA是四节点四面体单元，后面跟的是四个节点 ID。这些字段在文件里是右对齐的，但 MATLAB 的textscan如果按空格分割，遇到连续空格或者负号紧贴数字的情况就会错位。所以这个 m 文件项目的第一步，就是自己写一个按固定宽度切分的解析器，而不是依赖textscan的默认分隔符。
2.2 目录结构：一个主入口加若干卡片解析函数
拿到这个资源后，我建议先看目录结构。典型的组织方式是这样：
bdf_reader/ ├── read_bdf.m % 主入口，读文件、分发卡片 ├── parse_grid.m % 解析 GRID 卡片 ├── parse_ctetra.m % 解析 CTETRA 卡片 ├── parse_chexa.m % 解析 CHEXA 卡片 ├── parse_cquad4.m % 解析 CQUAD4 卡片 ├── parse_psolid.m % 解析 PSOLID 属性 ├── parse_mat1.m % 解析 MAT1 材料 ├── parse_spc1.m % 解析 SPC1 约束 ├── utils/ │ ├── split_fixed_width.m % 固定宽度切分 │ └── strip_comment.m % 去掉 $ 注释行 └── demo/ └── demo_read.m % 示例脚本主入口read_bdf.m负责打开文件、逐行读取、判断卡片类型，然后调用对应的parse_xxx函数。每个解析函数只关心自己那类卡片的字段含义，返回结构体或矩阵。这种拆分方式的好处是：你要新增一种卡片支持，只需要加一个parse_xxx.m并在主入口注册，不用动其他代码。常见做法是在主入口里用一个containers.Map或者switch语句做分发，我一般会用switch，因为卡片类型不多，可读性更好。
2.3 固定宽度切分的实现细节
固定宽度切分是核心工具函数，写不好后面全错。下面这段代码是我从项目里摘出来的简化版：
function fields = split_fixed_width(line, width) % SPLIT_FIXED_WIDTH 按固定宽度切分一行文本 % line : 输入字符串，通常是一行 BDF 卡片 % width : 每个字段的字符宽度，Nastran 短字段为 8 % fields: 切分后的字符串元胞数组 line = strtrim(line); n = numel(line); nf = ceil(n / width); fields = cell(1, nf); for k = 1:nf idx1 = (k-1)*width + 1; idx2 = min(k*width, n); fields{k} = strtrim(line(idx1:idx2)); end end逻辑说明：先去掉行首尾空格，然后按width计算需要切几个字段。每个字段用strtrim去掉前后空格，这样GRID后面的节点 ID 和坐标值就能干净地取出来。参数width默认传 8，但如果遇到长字段格式（逗号分隔），这个函数就不适用了，需要在主入口里先判断格式。注意strtrim会把空字段变成空字符串，后面用str2double转数值时，空字符串会变成 NaN，所以解析函数里要判断字段是否为空，空就跳过或者给默认值。
3. 跑通第一个模型：从 demo 脚本到数据落地
3.1 准备一个最小 BDF 文件
在跑 demo 之前，你需要一个 BDF 文件。如果手头没有 Patran，可以自己写一个最小模型，比如两个四面体单元组成的立方体。下面是一个合法的 BDF 片段：
BEGIN BULK GRID 1 0.000000 0.000000 0.000000 GRID 2 1.000000 0.000000 0.000000 GRID 3 0.000000 1.000000 0.000000 GRID 4 0.000000 0.000000 1.000000 GRID 5 1.000000 1.000000 1.000000 CTETRA 1 1 2 3 4 CTETRA 2 2 3 4 5 MAT1 1 2.100000E+11 3.000000E-01 7.800000E+03 PSOLID 1 1 ENDDATA注意BEGIN BULK和ENDDATA之间的内容才是有效数据，$开头的行是注释。Patran 导出的文件通常还有BEGIN BULK之前的头部信息，解析时要跳过。
3.2 demo 脚本逐段拆解
项目里的demo_read.m一般长这样：
% demo_read.m % 读取 BDF 文件并打印节点和单元信息 filename = fullfile('data', 'simple.bdf'); model = read_bdf(filename); fprintf('节点数: %d\n', size(model.grid, 1)); fprintf('单元数: %d\n', size(model.ctetra, 1)); fprintf('材料数: %d\n', numel(model.mat1)); % 打印前 5 个节点坐标 disp('前 5 个节点坐标:'); disp(model.grid(1:min(5, end), 2:4)); % 打印第一个单元的节点连接 if ~isempty(model.ctetra) fprintf('第一个 CTETRA 单元节点: %s\n', ... mat2str(model.ctetra(1, 2:5))); end逻辑说明：read_bdf返回一个结构体model，里面按卡片类型分字段。model.grid是一个 N×4 矩阵，第一列是节点 ID，后三列是 x、y、z 坐标。model.ctetra是 M×5 矩阵，第一列是单元 ID，后四列是节点 ID。model.mat1是结构体数组，每个元素包含材料 ID、弹性模量、泊松比、密度。参数方面，filename可以用绝对路径或相对路径，建议用fullfile拼路径，避免 Windows 和 Linux 分隔符差异。
3.3 运行结果与常见报错
跑通后你应该看到类似输出：
节点数: 5 单元数: 2 材料数: 1 前 5 个节点坐标: 0 0 0 1 0 0 0 1 0 0 0 1 1 1 1 第一个 CTETRA 单元节点: [1 2 3 4]如果报错Undefined function 'read_bdf'，说明当前工作目录不在项目根目录，用addpath(genpath(pwd))把子目录加进来。如果报错Index exceeds matrix dimensions，大概率是 BDF 文件里有空行或者注释行没跳过，检查read_bdf.m里有没有if isempty(strtrim(line)) || startsWith(line, '$')这样的过滤。如果读出来的坐标全是 NaN，检查split_fixed_width的width参数是不是 8，有些 Patran 版本导出的是长字段格式，需要改成按逗号分割。
4. 核心解析函数拆解：GRID、CTETRA 与 MAT1 怎么写
4.1 GRID 解析：坐标字段的容错处理
parse_grid.m是使用频率最高的函数，因为节点数据量大。下面是一个健壮版本：
function grid = parse_grid(fields) % PARSE_GRID 解析 GRID 卡片 % fields: split_fixed_width 切分后的字段元胞 % grid : 1×4 向量 [节点ID, x, y, z] grid = zeros(1, 4); grid(1) = str2double(fields{2}); % 节点 ID % fields{3} 是坐标系编号，通常忽略 for k = 1:3 val = str2double(fields{3 + k}); if isnan(val) val = 0.0; % 空字段按 0 处理 end grid(1 + k) = val; end end逻辑说明：fields{1}是'GRID'，fields{2}是节点 ID，fields{3}是坐标系编号，fields{4:6}是坐标。这里对坐标做了 NaN 容错，因为有些 BDF 文件里坐标值为空表示 0。参数方面，如果遇到大模型，逐行调用parse_grid会有性能问题，常见优化是预分配矩阵，每读一行往里面填，而不是每次[grid; new_row]拼接。
4.2 CTETRA 解析：单元节点数不固定怎么办
CTETRA是四节点四面体，但 Nastran 还有CTETRA的变体，比如带中间节点的十节点四面体。这个项目里通常只处理四节点版本，但解析函数要能识别字段数量：
function ctetra = parse_ctetra(fields) % PARSE_CTETRA 解析 CTETRA 卡片 % 返回 [单元ID, n1, n2, n3, n4] ctetra = zeros(1, 5); ctetra(1) = str2double(fields{2}); for k = 1:4 ctetra(1 + k) = str2double(fields{2 + k}); end end逻辑说明：fields{2}是单元 ID，fields{3:6}是四个节点 ID。如果遇到十节点版本，字段会多出 6 个，这时候要么扩展返回向量，要么在read_bdf里根据字段数量分发到不同函数。我一般会在主入口里判断numel(fields)，大于 7 就调用parse_ctetra10。参数方面，节点 ID 必须是正整数，如果str2double返回 NaN，说明字段切分错了，回头检查split_fixed_width。
4.3 MAT1 解析：材料参数的字段映射
MAT1卡片字段比较多，常见的是弹性模量、泊松比、密度，但 Nastran 允许只写部分字段，后面的用空字段占位。解析时要按位置取：
function mat = parse_mat1(fields) % PARSE_MAT1 解析 MAT1 卡片 % 返回结构体，包含 id, E, nu, rho mat.id = str2double(fields{2}); mat.E = str2double(fields{3}); mat.nu = str2double(fields{4}); mat.rho = str2double(fields{5}); if isnan(mat.E), mat.E = 0; end if isnan(mat.nu), mat.nu = 0; end if isnan(mat.rho), mat.rho = 0; end end逻辑说明：fields{2}是材料 ID，fields{3}是弹性模量，fields{4}是泊松比，fields{5}是密度。注意有些 BDF 文件里MAT1后面还有剪切模量、热膨胀系数等，这个项目只取最常用的三个。参数方面，如果材料数据在后续计算中要用到，建议把单位统一成国际单位制，Patran 导出的单位取决于建模时的设置，解析函数不做单位转换，需要使用者自己保证一致。
5. 避坑与排查：BDF 读取时最容易翻车的几个地方
5.1 现象：读出来的节点坐标整体偏移了一位
原因：split_fixed_width的width参数设成了 8，但某些 Patran 版本导出的 BDF 用的是长字段格式，字段之间用逗号分隔，不是固定宽度。这时候按 8 切分，第一个字段会包含GRID,1,这样的内容，str2double直接返回 NaN。
解决：在read_bdf.m里加一个格式判断，如果某行包含逗号，就改用strsplit(line, ',')分割，然后再strtrim每个字段。判断逻辑可以放在读第一行有效卡片时做一次，后面统一用同一种切分方式。
5.2 现象：单元连接关系里出现 0 或者负数
原因：BDF 文件里有些卡片有续行，续行以+或*开头，如果没处理续行，第二行的字段会被当成新卡片解析，导致节点 ID 错位。另外，Nastran 允许节点 ID 从 0 开始，但 MATLAB 矩阵索引从 1 开始，如果直接拿节点 ID 当索引就会报错。
解决：在read_bdf.m里维护一个“当前卡片”缓冲区，遇到续行就追加到缓冲区，直到遇到新卡片再解析。节点 ID 到矩阵索引的映射用一个containers.Map或者sparse数组做转换，不要直接拿 ID 当索引。
5.3 现象：大模型读取速度极慢，内存爆掉
原因：逐行调用parse_grid并动态扩展矩阵，每次[grid; new_row]都会重新分配内存，节点数上万后时间复杂度是 O(n²)。另外，把所有卡片都存成结构体数组，字段开销也很大。
解决：预分配矩阵。先扫一遍文件统计各类卡片数量，然后zeros(n, 4)预分配，读的时候按行号填入。如果模型特别大，可以考虑用memmapfile或者分块读取。我一般会在read_bdf里加一个prealloc选项，默认开启。
5.4 现象：注释行和空行导致解析中断
原因：BDF 文件里$开头的注释行和空行很常见，如果read_bdf没有过滤，split_fixed_width会把注释内容也切分，str2double返回 NaN，后续逻辑就乱了。
解决：在逐行读取的循环里，第一件事就是line = strtrim(line); if isempty(line) || line(1) == '$', continue; end。注意$可能出现在行首，也可能出现在行中间表示行内注释，但 Patran 导出的文件通常$只在行首。
5.5 现象：MATLAB 版本差异导致函数不兼容
原因：项目里可能用了startsWith、endsWith、string类型这些较新的函数，在 R2016b 之前的版本会报错。另外，containers.Map在旧版本里行为也有差异。
解决：如果需要在旧版本运行，把startsWith(line, '$')换成strncmp(line, '$', 1)，把string类型换成char。我一般会在项目根目录放一个check_version.m，启动时检查 MATLAB 版本并给出提示。
6. 二次开发：把 BDF 读取接入你自己的求解流程
6.1 扩展新卡片：以 CQUAD4 和 PSHELL 为例
假设你要处理壳单元模型，需要支持CQUAD4和PSHELL。步骤很直接：在parse_cquad4.m里按字段取单元 ID 和四个节点 ID，在parse_pshell.m里取属性 ID、材料 ID、厚度。然后在read_bdf.m的switch里加两个case。注意CQUAD4的字段顺序和CTETRA不同，fields{2}是单元 ID，fields{3:6}是节点，但PSHELL的厚度在fields{4}，不是fields{3}。这种字段映射最好在函数头部用注释写清楚，方便后面改。
6.2 数据接口：从结构体到求解器输入
读进来的model结构体是中间格式，你的求解器可能需要不同的数据结构。常见做法是写一个model2solver.m转换函数，把model.grid的坐标矩阵和model.ctetra的连接矩阵转成求解器需要的nodes和elements。如果求解器用稀疏矩阵组装刚度阵，节点 ID 要先映射成 1 到 N 的连续索引。下面是一个映射示例：
function [nodes, elements, id_map] = model2solver(model) % MODEL2SOLVER 把 BDF 模型转成求解器输入 % nodes : N×3 坐标矩阵 % elements: M×4 单元连接，索引从 1 开始 % id_map : 原节点 ID 到新索引的映射 original_ids = model.grid(:, 1); id_map = containers.Map(original_ids, 1:numel(original_ids)); nodes = model.grid(:, 2:4); elements = zeros(size(model.ctetra, 1), 4); for k = 1:size(model.ctetra, 1) for j = 1:4 elements(k, j) = id_map(model.ctetra(k, 1 + j)); end end end逻辑说明：id_map把原始节点 ID 映射到 1 到 N 的连续整数，elements里存的是新索引。这样求解器组装刚度阵时可以直接用elements(k, :)去索引nodes。参数方面，如果模型里有孤立节点（没有被任何单元引用），id_map里仍然会有映射，但elements不会用到，不影响计算。
6.3 验证方法：用已知模型对拍
改完解析函数后，怎么确认没读错？我一般会做两件事：一是用 Patran 或者 Nastran 自带的工具把同一个 BDF 文件转成文本报告，对比节点数和单元数；二是构造一个简单模型，手算刚度阵或者质量阵，和 MATLAB 读进来的数据算出来的结果对拍。比如两个四面体单元的立方体，总节点数 5，单元数 2，材料弹性模量 2.1e11，泊松比 0.3，密度 7800，这些数值在 BDF 里写死，读进来后打印出来核对。如果对不上，就逐字段检查split_fixed_width的切分结果。
6.4 一个我踩过的坑：单位不一致导致结果差 1000 倍
有一次读一个 Patran 导出的 BDF，节点坐标是毫米，材料弹性模量是 MPa，我默认按米和帕斯卡算，结果位移小了 1000 倍。后来在read_bdf里加了一个unit参数，默认'mm'，在model2solver里统一转成米和帕斯卡。从那以后我每次读新模型都强制走一遍单位检查：先打印几个节点坐标和材料参数，确认量级对不对，再往下算。希望这个习惯能帮到你，少走一段弯路。
本文还有配套的精品资源，点击获取
