1 问题描述

  如下图所示的框架结构,其顶部梁跨中受集中力作用,用有限元方法分析该结构的应力、应变及位移。结构中各个截面的参数都为:E=3.0×1011PaE=3.0\times 10^{11}Pa,I=6.5×10−7m4I=6.5\times 10^{-7}m^4,A=6.8×10−4m2A=6.8\times 10^{-4}m^2。尝试改变各种不同的边界条件及载荷施加位置,进行相应的有限元分析,并对各种工况进行对比分析(至少选择2种不同的工况)。

图1.1 框架结构示意图

图1.1 框架结构示意图(A、B 为固定支座,C、D 为刚接节点,E 为载荷施加位置)

  本文中使用的Abaqus版本是2025。

2 Abaqus有限元模型与计算结果分析

  根据已知数据,假设梁截面为矩形。可以算出梁截面为 0.00635m×0.10710m0.00635m\times 0.10710m——该尺寸与题给的 II、AA 自洽(按 I=bh3/12I=bh^3/12 与 A=bhA=bh 反算,误差均在 0.01% 以内)。

  本文对比两种工况:中间载荷(P 作用于跨中 E 点)与偏置载荷(P 偏离跨中)。每种工况给出模型、应力与位移三张结果图。

2.1 中间载荷

  中间载荷Abaqus源文件点此处下载。

图2.1 框架结构有限元分析模型(中间载荷)

图2.1 框架结构有限元分析模型(中间载荷)

图2.2 框架结构应力云图(中间载荷)

图2.2 框架结构应力云图(中间载荷,Mises 应力峰值 3.362×10⁹ Pa)

图2.3 框架结构位移云图(中间载荷)

图2.3 框架结构位移云图(中间载荷,最大位移 4.422×10⁻¹ m)

2.2 偏置载荷

  偏置载荷Abaqus源文件点此处下载。

图2.4 框架结构有限元分析模型(偏置载荷)

图2.4 框架结构有限元分析模型(偏置载荷)

图2.5 框架结构应力云图(偏置载荷)

图2.5 框架结构应力云图(偏置载荷,Mises 应力峰值 2.498×10⁹ Pa)

图2.6 框架结构位移云图(偏置载荷)

图2.6 框架结构位移云图(偏置载荷,最大位移 3.038×10⁻¹ m)

3 结果对比与讨论

  两种工况的计算结果如下表。

工况 Mises 应力峰值 最大位移
中间载荷(E 点) 3.362×109Pa3.362\times 10^{9}Pa 4.422×10−1m4.422\times 10^{-1}m
偏置载荷 2.498×109Pa2.498\times 10^{9}Pa 3.038×10−1m3.038\times 10^{-1}m

  由结果可见,载荷施加位置显著改变了结构的力学响应:应力峰值相差约 35%,最大位移相差约 46%。载荷偏离跨中后,两侧立柱分担的弯矩不再对称,靠近载荷一侧的柱与梁柱节点承受更大内力,因而应力分布与变形形态都随之改变。

  两种工况下峰值应力都出现在顶部横梁与立柱的刚接节点处(靠载荷一侧),而非载荷施加点正下方。以中间载荷工况为例,红色高应力区集中在横梁右端的节点附近;偏置载荷工况同样如此。

  这与框架的受力特点有关。竖向载荷在横梁内产生弯矩,梁端弯矩需通过刚接节点传递给立柱;节点处既要传递梁端弯矩,又要承受立柱的柱端弯矩,两个方向的弯矩叠加,使该处成为全结构弯矩最大、应力最高的位置。而载荷作用点处虽有集中力,但该处截面承受的主要是局部压力,量值远低于节点处的弯曲应力。

  需要说明的是,本题仅对比了载荷位置的变化,未改变边界条件(A、B 始终保持固定支座)。若进一步固定支座改为铰支,结构将成为机构或改变内力分布形式,是另一个值得单独分析的方向。