组装各单元刚度矩阵,求解整体刚度方程:
$$\mathbf{K}\mathbf{u}= \mathbf{f}$$
单元刚度矩阵(方向余弦 $c,s$):
$$k_e = \frac{AE}{L}\begin{bmatrix}c^2&cs&-c^2&-cs\\cs&s^2&-cs&-s^2\\-c^2&-cs&c^2&cs\\-cs&-s^2&cs&s^2\end{bmatrix}$$
杆件轴力:$N = \frac{AE}{L}(c\Delta u_x + s\Delta u_y)$
基于直接刚度法对二维铰接桁架进行有限元分析。从普拉特、沃伦、悬臂等5种预设构型中选择,实时计算并可视化杆件拉压内力与变形图。
组装各单元刚度矩阵,求解整体刚度方程:
$$\mathbf{K}\mathbf{u}= \mathbf{f}$$
单元刚度矩阵(方向余弦 $c,s$):
$$k_e = \frac{AE}{L}\begin{bmatrix}c^2&cs&-c^2&-cs\\cs&s^2&-cs&-s^2\\-c^2&-cs&c^2&cs\\-cs&-s^2&cs&s^2\end{bmatrix}$$
杆件轴力:$N = \frac{AE}{L}(c\Delta u_x + s\Delta u_y)$
整个分析的核心是求解整体刚度方程,它表达了所有节点上的力与位移之间的关系:
$$\mathbf{K}\mathbf{u}= \mathbf{f}$$其中,$\mathbf{K}$ 是整体刚度矩阵(由所有单元组装而成),$\mathbf{u}$ 是所有节点的位移向量,$\mathbf{f}$ 是所有节点上的外力向量。解这个方程,就能从已知的力 $\mathbf{f}$ 求出未知的位移 $\mathbf{u}$。
整体刚度矩阵 $\mathbf{K}$ 由一个个杆单元的刚度矩阵 $k_e$ 像拼图一样组装起来。每个杆单元的刚度矩阵由其几何和材料属性决定:
$$k_e = \frac{AE}{L}\begin{bmatrix}c^2 & cs & -c^2 & -cs \\ cs & s^2 & -cs & -s^2 \\ -c^2 & -cs & c^2 & cs \\ -cs & -s^2 & cs & s^2 \end{bmatrix}$$这里,$A$是截面积,$E$是弹性模量,$L$是杆长。$c = \cos\theta$ 和 $s = \sin\theta$ 是杆件方向的方向余弦($\theta$是杆件与水平轴的夹角)。这个矩阵建立了杆件两端两个节点的力与位移的联系。
桥梁工程:普拉特桁架和沃伦桁架是铁路和公路桥的经典形式。工程师使用有限元分析精确计算各杆件内力,以确定最优的杆件截面,确保桥梁在车辆荷载下安全可靠,同时避免材料浪费。
屋顶与体育馆结构:大型工业厂房、体育场馆的屋顶常采用轻质桁架结构。分析时需考虑风荷载、雪荷载,通过模拟不同荷载方向和大小,找到最危险的杆件(通常是压力最大的杆件,需防屈曲)。
起重机与施工设备:塔吊的臂架和塔身本质上是大型悬臂桁架。设计时必须分析其在吊起最大重量时的位移和杆件内力,确保作业时臂架顶端不会晃动过大,且关键受压杆件不会失稳。
通讯塔与电力塔架:这些高耸结构主要承受风荷载。通过有限元分析可以优化其斜杆(腹杆)的布置,用最少的材料抵抗来自各个方向的风力,并控制塔顶的位移在允许范围内。
开始使用这个工具时,有几个需要注意的地方。首先切勿忘记“铰接”假设。这个模拟器构建的是一个所有杆件端部都以完全自由旋转的“铰链”连接的世界。因此杆件只会产生轴力(拉伸或压缩)。但实际的焊接或螺栓连接具有一定程度的“刚性连接”,也会传递弯矩。如果不理解这个差异就直接应用于实际设计会很危险。例如,细长杆件受压屈曲的风险可以用这个工具再现,但连接处的局部应力集中是无法显示的。
其次是参数设置的现实感。随意拖动滑块虽然有趣,但例如将“杨氏模量E”从10 GPa(类似橡胶的值)改为200 GPa(钢材)时,变形量会减小到百分之一以下。在实际结构中,改变材料会显著影响成本和重量,所以这种直觉很重要。另外,将“截面积A”设置得过小,受压杆件的变形会异常增大,这其实显示了线性分析的局限。实际上,在发生那种变形之前,结构会因屈曲这种非线性现象而破坏。不要盲目相信工具的结果,要带着怀疑的眼光思考:“这根受压杆件是不是太细了?”
最后是支座条件的解读。固定支座(△)与滚动支座(▽)的区别在于是否约束水平方向的移动。例如,桥梁一端固定、另一端设为滚动支座,是为了给温度变化引起的桥梁伸缩留出余地。如果两端都设为固定,工具上会显示刚度提高、变形减小,但实际上会产生巨大的热应力。仅仅一个支座的设置,就会从根本上改变结构的行为和力的传递路径。
熟悉这个工具后,建议下一步学习“矩阵结构分析”的基础。工具在背后进行的“直接刚度法”的核心,在于将单元刚度$\mathbf{k_e}$组装成整体刚度$\mathbf{K}$的“装配”过程。具体来说,就是将$\mathbf{k_e}$的分量累加到与各单元节点编号对应的自由度位置上。亲手用笔算体验一次这个“累加”过程(用节点数少的简单桁架即可),会立刻驱散对有限元法的黑箱感。
在数学上,这会连接到数值分析领域,即如何高效求解线性方程组$\mathbf{K}\mathbf{u}= \mathbf{f}$。实际工程中的有限元软件,正是用乔列斯基分解或预处理共轭梯度法等算法来求解这个拥有数百万未知数的方程。此外,从桁架扩展到框架结构(承受弯曲的杆项)时,单元刚度矩阵会变成6x6(每个节点有位移和转动共3个自由度),从而能够表现更多的实际结构。
最终,可以尝试用市售的通用CAE软件(例如Ansys、Abaqus、Nastran等)进行类似的桁架分析。届时,选择“桁架(杆)单元”作为单元类型,体验网格划分、材料定义、边界条件设置以及后处理中应力可视化这一系列流程,就能切身理解这个浏览器工具对应实际工作中的哪个环节。首先“用工具直观理解现象”,其次“用理论深挖原理”,最后“用专业软件重现”,这是掌握CAE技术的可靠路径。
深化理论:在本工具的简化模型基础上,进一步研究非线性效应、三维行为和时间依赖现象。阅读专业教材和学术论文,掌握严格的数学推导,是提升工程解题能力的关键。
数值方法:系统学习有限元法(FEM)、有限差分法(FDM)和有限体积法(FVM),理解商业CAE求解器的内部运行机制,这将显著提升您设置有效仿真的能力。
实验验证:理论和仿真结果必须通过实验数据加以验证。养成将计算结果与测量值进行对比的习惯,这正是V&V(验证与确认)的精髓所在。
CAE工具:准备好后,可进一步探索Ansys、Abaqus、OpenFOAM、COMSOL等业界主流工具。通过本模拟器培养的物理直觉,将帮助您更有效地配置和使用这些工具。
以普拉特桁架为例:跨度L=6m、高度h=2m、杆件采用Q345钢(E=210GPa)、上弦杆截面积A=1500mm²、下弦杆A=1800mm²、腹杆A=800mm²。在中点施加竖向荷载P=50kN,求解结果显示:中点竖向位移δ=8.3mm、上弦杆最大拉力N=145kN、竖腹杆最大压力N=-72kN、系统自由度数DOF=24