单元分布力 =>>单元节点力（曾攀有限元分析）

最新推荐文章于 2023-03-24 09:48:34 发布

王.伟

最新推荐文章于 2023-03-24 09:48:34 发布

阅读量2.8k

点赞数

分类专栏：有限元理论文章标签：线性代数矩阵几何学

本文链接：https://blog.csdn.net/weixin_43717845/article/details/105569662

版权

有限元理论专栏收录该内容

7 篇文章 11 订阅

订阅专栏

在这里插入图片描述

$\qquad分布力\bar{p}(x)所做的外力功，使用形函数N(x)与节点位移q^e代替位移函数v(x),\\ 分布力\bar{p}(x)，所以P^e_{(4\times1)}等于形函数与均布力的积分。$
$P^e_{(4\times1)}=\int_{l}N_{(x)}^T\cdot\bar{p}(x)\,d(x)$
$下面以三角形均布荷载左端 A 点弯矩为例，计算其等效节点弯矩$
$\begin{aligned} 左端弯矩对应的形函数为N_{(x)}中的第二项(\xi-2\xi^2+\xi^3)l\\ \xi=\frac{x}{l}\\ \\ N_{(x)}=(\frac{x}{l}-2(\frac{x}{l})^2+(\frac{x}{l})^3)l\quad\qquad\quad\qquad\\ \end{aligned}$
$\bar{p}(x)=\begin{cases} 2p_0\cdot\frac{x}{l}, & 0\leq x<\frac{L}{2}\\ 2p_0\cdot(1-\frac{x}{l}), & \frac{L}{2}\leq x\leq L \end{cases}$

$代入积分，左端 A 点弯矩：$
$\begin{aligned} M_A=\int_0^{{L}}N_{(x)}\cdot\bar{p}(x)\,d(x)\quad\qquad\quad\quad\qquad\qquad\qquad\\ =\int_0^{\frac{L}{2}}(\frac{x}{L}-2(\frac{x}{L})^2+(\frac{x}{L})^3)L\cdot2p_0\cdot\frac{x}{L}\,\ d(x)+\qquad\\ \int_{\frac{L}{2}}^{L}(\frac{x}{L}-2(\frac{x}{L})^2+(\frac{x}{L})^3)L\cdot2p_0\cdot(1-\frac{x}{L})\,\ d(x)\\ \\ \end{aligned}$
$计算得$
$M_A=\frac{5}{96}p_0L^2$