12列作为子矩阵,Р解方程得(用Matlab求解): Р(7)求约束反力和单元应力Р将节点位移矢量u代入总体平衡方程即得各节点的受力。Р Р可见,1节点的水平支反力为1N,方向向右;竖直支反力为0.8846N,方向向上;2节点的水平支反力为3N,方向向左;竖直支反力为0.1154N,方向向上。Р将各个单元节点位移代入方程,可得各单元上的应力分布。Р二、与ANSYS对比Р(1)结果Р在ANSYS 14.5中计算此模型,采用PLANE182单元(没有中节点),节点编号和边界条件如图所示。计算结果如下表所示,与以上结果完全吻合。Р图3 ANSYS中的网格划分Р图4 节点位移解Р图5 1、2两节点的支反力解Р(2)结果对比Р节点位移:(保留4位小数)Р单位Р(mm)Р节点3Р节点4Р节点5Р节点6Рu3Рv3Рu4Рv4Рu5Рv5Рu6Рv6Р手算结果Р-0.0150Р-0.2172Р0.2131Р-0.2710Р0.0542Р-0.5968Р0.3724Р-0.6976РANSYS结果Р-0.0150Р-0.2172Р0.2131Р-0.2710Р0.0542Р-0.5968Р0.3724Р-0.6976РRatioР1.0000Р1.0000Р1.0000Р1.0000Р1.0000Р1.0000Р1.0000Р1.0000Р节点反力:(保留4位小数)Р单位Р(N)Р节点1Р节点2РfxРfyРfxРfyР手算结果Р1.0000Р0.8846Р-3.0000Р0.1154РANSYS结果Р1.0000Р0.8847Р-3.0000Р0.1154РRatioР1.0000Р1.0001Р1.0000Р1.0000Р结论:Р通过结果比较可以发现,在同样的精度下,手算结果和ANSYS计算结果几乎完全吻合,可以肯定手算结果是正确的。误差主要来源于计算机的舍入误差和取结果时的有效位数。
全文预览
