3. 一维杆件与桁架材料建模实战

一维杆单元是有限元方法中最基础、也最适合用来建立“直觉”的单元类型。它只承受轴向力，节点只有沿杆轴方向的平动自由度，因此刚度矩阵可以完全手算推导，并与软件结果逐项对照。本节以一根两节点杆单元为对象，完成从刚度矩阵推导、材料参数输入、轴向载荷与温度载荷施加，到手算与软件结果比对的完整闭环。

一、杆单元刚度矩阵的推导

取一根长度为 L、截面积为 A、弹性模量为 E 的等截面直杆，两端节点编号为 1、2，节点轴向位移分别为 u₁、u₂。规定沿杆轴向右为正方向，节点力 F₁、F₂ 也以向右为正。杆的轴向应变为

ε = (u₂ − u₁) / L

由胡克定律 σ = Eε，轴力 N = σA = EA(u₂ − u₁)/L。节点 1 受到的力为 −N，节点 2 受到的力为 +N，写成矩阵形式即得到杆单元刚度方程：

{ F1 }   EA {  1  -1 } { u1 }
{    } = -- {         } {    }
{ F2 }    L { -1   1 } { u2 }

记单元刚度矩阵为 k^e = (EA/L)[[1, −1], [−1, 1]]。它有两个重要性质：一是奇异性，行列式为零，说明单元本身存在刚体位移，必须施加边界条件后才能求解；二是对称性，k₁₂ = k₂₁，这是线弹性系统能量守恒的体现。对于桁架结构，各杆单元按节点自由度“对号入座”组装成整体刚度矩阵 K，再引入边界条件求解 Ku = F。

二、材料参数输入

以一维钢杆为例，采用示例参数：弹性模量 E = 2.0×10⁵ MPa，截面积 A = 200 mm²，长度 L = 500 mm，线膨胀系数 α = 1.2×10⁻⁵ /℃。则单元轴向刚度系数为

EA/L = 2.0×10⁵ × 200 / 500 = 8.0×10⁴ N/mm

在通用有限元软件中，杆单元通常对应 LINK180（ANSYS）、T2D2（Abaqus）或 Truss 单元（COMSOL）。材料参数输入时需注意三点：第一，E 必须与长度单位一致，若几何用 mm、力用 N，则 E 用 MPa；第二，A 是实常数而非材料属性，需在截面/实常数中单独定义；第三，温度载荷需要额外输入线膨胀系数 α 和参考温度 T_ref，否则热应变无法计算。

三、轴向载荷与温度载荷

轴向载荷分为节点集中力和沿杆分布力。节点集中力直接加到对应自由度上；分布力 q（N/mm）需按静力等效原则转化为节点力，对均布载荷有 F₁ = F₂ = qL/2。温度载荷则通过初应变方式进入方程：温度变化 ΔT = T − T_ref 产生的自由热应变为 ε_th = αΔT，等效节点力为

F_th = EAαΔT · [−1, 1]^T

若杆两端完全约束，温度升高时杆内产生压应力 σ = −EαΔT；若一端自由，则杆自由伸长 ΔL = αΔT·L，应力为零。这一区别是判断温度载荷结果是否合理的关键依据。

四、手算与软件对照验证

以“一端固定、一端受 20 kN 拉力”的杆为例。固定节点 1，即 u₁ = 0，节点 2 受 F₂ = 20 kN。由刚度方程第二行得

u₂ = F₂·L/(EA) = 20000 × 500 / (2.0×10⁵ × 200) = 0.25 mm

轴向应力 σ = F/A = 20000/200 = 100 MPa，轴向应变 ε = σ/E = 5×10⁻⁴，与 u₂/L = 0.25/500 = 5×10⁻⁴ 一致。软件中读取节点 2 的轴向位移应为 0.25 mm，单元轴向应力应为 100 MPa，反力 R₁ = −20 kN。

再叠加温度载荷：设 ΔT = 40 ℃，两端固定。此时热等效节点力为 EAαΔT = 2.0×10⁵×200×1.2×10⁻⁵×40 = 19200 N = 19.2 kN。由于两端固定，节点位移为零，热应力 σ_th = −EαΔT = −2.0×10⁵×1.2×10⁻⁵×40 = −96 MPa，即压应力。若改为仅节点 1 固定，则节点 2 自由伸长 ΔL = αΔT·L = 1.2×10⁻⁵×40×500 = 0.24 mm，热应力为零。软件中可通过“约束端反力是否为零”和“自由端位移是否等于 αΔT·L”来快速判断温度载荷施加是否正确。

通过本节的一维杆件实战，可以把“刚度矩阵—材料参数—载荷—结果判读”这条链路完整走通。后续桁架建模只是把多个杆单元按节点自由度组装，其单元层面的推导、参数输入和验证逻辑与本节完全一致。