弹塑性力学第八章.ppt

2020 2 7 1 8 1位移法求解 第八章柱体的自由扭转问题 8 2按应力函数求解 8 3薄膜比拟 8 4等截面杆扭转按应力函数举例 8 5薄壁杆的自由扭转 2020 2 7 2 在第五章的最后我们以圆柱形杆的扭转问题为例来说明空间三维问题的求解过程 无体力 对于圆杆扭转 扭矩Mz MT 应力 x y z xy 0 位移分量 u Kyz v Kxz w 0 K为单位长扭转角 2020 2 7 3 对于一般等截面杆扭转w 0称为自由扭转 为了求解一般等截面杆自由扭转 参考圆杆扭转解进行假设 半逆解 8 1位移法求解 对于一般等截面杆自由扭转 可设位移分量 u Kyz v Kxz u v与园杆扭转一致 w K x y w不能为零 为x y函数 而 x y 称为扭曲函数 2020 2 7 4 8 1位移法求解 无体力等截面杆扭转位移表达式已设定 未知量为 K和 x y 工程 应变分量 u Kyz v Kxz w K x y 2020 2 7 5 8 1位移法求解 应力分量 x y z xy 0 所有物理量均由K和 x y 表示 2020 2 7 6 8 1位移法求解 按位移法求解 基本方程为平衡微分方程 三个 或 2 0 两个平衡微分方程自然满足 而第三个方程为 2020 2 7 7 8 1位移法求解 基本方程仅为一个 求解 x y 的方程 由基本方程可见 x y 为一个调合函数 同时在基本方程中不出现K K的确定当然也应通过边界条件来确定 扭曲函数 x y 除了满足 2 0 还需要满足边界条件 2020 2 7 8 8 1位移法求解 首先考察扭杆侧边的边界条件 主要边界 在侧边上方向余弦 l m n l m 0 面力 满足 2020 2 7 9 8 1位移法求解 上式也可以用 边界条件用 x y 的偏微分表示 由于 则 代入侧面边界条件 2020 2 7 10 8 1位移法求解 在扭杆端面 如z 0 法线的方向余弦 l m n 0 0 1 杆端截面法线方向面力 满足 合力为零 合力矩为 而在杆端截面面内的面力分布不清楚 应用圣维南原理 在 x y方向面力分量不清楚 但要求 2020 2 7 11 8 1位移法求解 上式也可以表示为 可以证明当扭曲函数 x y 在主要边界上力边界条件满足时 则和自然满足 见以下 2020 2 7 12 8 1位移法求解 利用格林公式 2 0 2020 2 7 13 8 1位移法求解 而第三个方程为 扭矩MT与K和 x y 的关系 小结 用位移法求解扭转问题归结为求解扭曲函数 x y 和单位扭转角K 2020 2 7 14 2 0在V上在杆侧边上 由 求 x y 8 1位移法求解 当 x y 确定后 利用杆端面条件 求K 令 扭转刚度 当 x y 和K均找到后 则扭杆的位移 应力均可求出 2020 2 7 15 作业 证明扭曲函数能用来求椭圆截面杆的扭转问题 其中a和b为椭圆截面的半轴长度 并且扭矩为 2020 2 7 16 8 2按应力函数求解 按位移法求解扭转问题要求在V内求解调和方程 2 0 其边界条件 x y 的微分形式 但能满足边界条件调合函数 x y 是不易找到的 下面讨论按应力法求解等截面杆扭转问题基本方程以及应力函数法求解等截面杆扭转问题的作法 2020 2 7 17 8 2按应力函数求解 2 1按应力法求解方程 同圆杆扭转类似 设 x y z xy 0 仅存在 zx x y xz和 zy x y yz 两个应力分量 将应力分量代入应力法的基本方程九个 三个平衡和六个相容方程 2020 2 7 18 8 2按应力函数求解 三个平衡方程 前两式自然满足 剩下一个控制方程 无体力相容方程为 由于设 x y z 0 0 2020 2 7 19 8 2按应力函数求解 则相容方程中有四个自然满足 仅剩下两个控制方程 2 zx 0和 2 zy 0 按应力法求解基本方程为三个 2 zx 0 2 zy 0 2020 2 7 20 8 2按应力函数求解 边界条件 在侧边 方向余弦 l m n l m 0 面力 前两个方程满足 第三个力边界条件 l zx m zy 0 在端面 方向余弦 l m n 0 0 1 面力 满足 2020 2 7 21 8 2按应力函数求解 在x y方向面力应用圣维南原理 2020 2 7 22 8 2按应力函数求解 2 2按应力函数 x y 求解 设应力分量与应力函数的关系为 则应力法第一个基本方程 平衡微分方程 自然满足 2020 2 7 23 8 2按应力函数求解 常数C是什么 C和位移法公式中的系数有什么关系 将上式代入应力法的其它两个基本方程 得 2 C 泊松方程 由应力函数法和位移法可知 2020 2 7 24 8 2按应力函数求解 将应力函数 代入杆侧边的边界条件l zx m zy 0 2020 2 7 25 8 2按应力函数求解 而 代入边界条件 得 则应力函数在扭杆侧边应该为常数 s C1 l zx m zy 0 2020 2 7 26 8 2按应力函数求解 对于单连域 可取 s 0 对于复连域 可取一条边界线上 s为零 而其它边界 s为非零常数 s0 0 si Ci 0 i 1 2 3 再将 x y 代入端面上的边界条件 方向余弦 l m n 0 0 1 面力 满足 2020 2 7 27 8 2按应力函数求解 在x y方向面力应用圣维南原理 第一 二方程恒满足 第一个方程 第二个方程 2020 2 7 28 8 2按应力函数求解 在x y方向面力应用圣维南原理 第三个方程 2020 2 7 29 8 2按应力函数求解 当为单连域时 在s上 s 0 当为多连域时 s0 0 si Ci 0 i 1 2 3 2020 2 7 30 8 2按应力函数求解 Ai为si围成的面积 2020 2 7 31 8 2按应力函数求解 总结 按应力函数 x y 求解 x y 须满足 2 2KG C 且 x y 与MT之间满足 单连域 多连域 2020 2 7 32 8 2按应力函数求解 在柱体侧边 s 0 单连域 si Ci 多连域 当k和 x y 由上述方程确定后 可求出 zx zy以及应变和位移 2020 2 7 33 8 3薄膜比拟 对于截面形状比较复杂的柱体 不管采用位移法还是应力法求解扭转问题解答 解析解 是很困难的 而普朗特 Prondtl 在1903年提出了薄膜比拟 它利用薄膜在均匀压力下的垂度与等截面直杆扭转问题中的应力函数在数学上的相似性 用薄膜来比拟扭杆 它可以帮助我们寻找扭转问题的解答 尤其是对截面较复杂的扭转可以避开数学上的困难 而采用实际薄膜比拟实验测定 形象的获得一些有价值的解 2020 2 7 34 8 3薄膜比拟 一均匀薄膜形状同扭杆截面 周边固定 并使薄膜受均匀微小压力q作用 薄膜将微微凸起 而形成曲面z z x y 薄膜仅承受张力 拉力 T 下面来寻求薄膜垂度z z x y 所应满足的方程和边界条件 2020 2 7 35 8 3薄膜比拟 寻求z z x y 应满足的方程 即求解方程是由薄膜微元dxdy的z方向的平衡条件来确定 Fz 0 2020 2 7 36 8 3薄膜比拟 整理后 得 或 z x y 所应满足的方程 2020 2 7 37 8 3薄膜比拟 与扭转问题应力函数 x y 所应满足方程和边界条件相比 2 2KG s 0 与z之间存在比拟关系 薄膜垂度z z x y 所应满足的边界条件 zs 0 单连域 2020 2 7 38 8 3薄膜比拟 薄膜垂度z x y 可由实验测定 再根据上再根据上式可确定 的分布规律 在应力函数解扭转问题时 考虑边界条件还有 由此式确定比例系数 单连域 扭矩MT与薄膜垂度所围成体积的两倍之间也同样存在一致的比拟关系 2020 2 7 39 8 3薄膜比拟 对于多连域 在孔边上应为常数 所以在薄膜比拟试验中 开孔区应用平行于x y平面的无重刚性平板来代替 扭杆剪应力 剪应力分量的大小与该薄膜垂度上对应点沿垂直方向的斜率成正比 2020 2 7 40 8 3薄膜比拟 扭转截面上任意点总剪应力 应力矢量t 数值和方向确定 任意点总剪应力数值 可利用薄膜等高线 平行于x y面的平面与薄膜相截可获得一系列闭合曲线 薄膜等高线 2020 2 7 41 8 3薄膜比拟 在等高线上任意点应力可沿x y方向分解 也可沿n s方向分解 根据剪应力分量与薄膜垂度沿垂直方向斜率成比例 在等高线上 所以在等高线上 2020 2 7 42 8 3薄膜比拟 任意点总剪力 等高线切方向 与垂度等高线的垂直方向斜率成正比 薄膜等高线为扭杆横截面上的剪应力线 发生在薄膜具有最陡斜率的点处 一般在杆边界上 截面上的最大剪应力 2020 2 7 43 8 3薄膜比拟 总结薄膜比拟与杆扭转各物理量之关系 2020 2 7 44 8 4等截面杆扭转按应力函数举例 采用应力函数解法求扭转问题 应力函数 x y 在域内满足方程 2 2KG 1 例题1 椭圆截面杆的扭转 在边界上满足方程 s 0 2 以及 3 2020 2 7 45 8 4等截面杆扭转按应力函数举例 由于椭圆杆截面方程为 因此 可设应力函数 x y 为 则 x y 自然满足方程 s 0 2020 2 7 46 8 4等截面杆扭转按应力函数举例 代回 x y 再代回 3 式 注意 将 x y 代入基本方程 2 2KG 得 2020 2 7 47 8 4等截面杆扭转按应力函数举例 再代回 3 式 注意 得 2020 2 7 48 8 4等截面杆扭转按应力函数举例 应力分量 各点总剪应力 最大剪应力在柱截面边界上 设a b 当y b时 为最大 2020 2 7 49 8 4等截面杆扭转按应力函数举例 应变 x y z xy 0 位移 当不考虑刚体位移时 椭圆杆扭转时 杆纵向发生位移 2020 2 7 50 8 4等截面杆扭转按应力函数举例 例题2 等边三角形截面 高为a 受扭矩Mz作用 求截面剪应力 x a 0 解 对于单连域 应力函数 s 0 考虑此原因 设 时同椭圆杆扭转一样 取三角形截面杆的边界方程为 的因子 2020 2 7 51 8 4等截面杆扭转按应力函数举例 设 则 x y 自然满足方程 s 0 2 2KG得4ma 2KG 将 x y 代入基本方程 2020 2 7 52 8 4等截面杆扭转按应力函数举例 利用 或 得 则 2020 2 7 53 8 4等截面杆扭转按应力函数举例 应力分量为 截面上最大剪应力 2020 2 7 54 8 4等截面杆扭转按应力函数举例 例题3 矩形截面杆的扭转 矩形截面 a b 受扭矩Mz作用 应力函数中要求 s 0 如果假设应力函数为 满足 s 0 但 2 2kG不能满足 2020 2 7 55 8 4等截面杆扭转按应力函数举例 所以直接采用上述 s 0的假设式不能作为扭转的应力函数 利用薄膜比拟 来判断狭矩形截面的应力函数特点 对于矩形截面杆扭转首先考虑a b时的情况 情况1 2020 2 7 56 8 4等截面杆扭转按应力函数举例 同样形状薄膜周边固定受均匀压力作用时 薄膜垂度变化如图 可见垂度曲面沿x方向很长一段为柱面 在此段 只是在狭矩形截面两端部 但区域很小 近似法忽略两端影响 2020 2 7 57 8 4等截面杆扭转按应力函数举例 y 这样狭矩形界面扭转应力函数 也认为 应满足基本方程为 1 2020 2 7 58 8 4等截面杆扭转按应力函数举例 s 0 2 3 2020 2 7 59 8 4等截面杆扭转按应力函数举例 由式 1 积两次分 得 将上式代入式 2 得C1 0 C2 GKb2 4 GK y2 b2 4 则 2020 2 7 60 8 4等截面杆扭转按应力函数举例 最后将 代入式 3 得 GK y2 b2 4 2020 2 7 61 8 4等截面杆扭转按应力函数举例 解得 则 应力分量 截面上最大剪应力 y b 2 2020 2 7 62 8 4等截面杆扭转按应力函数举例 原因是忽略了 zy 近似的 如不忽略 zy 很小 但力臂大 产生另一半Mz 2 按近似计算 偏于保守 实际上 将应力分量对截面形心取矩 得 2020 2 7 63 8 4等截面杆扭转按应力函数举例 情况2 一般矩形截面扭矩 a b且a b 按应力函数求解 则 x y 应满足 2 2KG b 2 0 a 2 0 x y 和K为待定 2020 2 7 64 8 4等截面杆扭转按应力函数举例 1 将求解 x y 的问题转化为求解一个调和函数F x y 的问题 考虑在狭矩形截面的应力函数为 1 GK y2 b2 4 能满足 2 1 2KG和 1 y b 2 0条件 选一般矩形截面的 x y 1 F x y GK y2 b2 4 F x y 2020 2 7 65 8 4等截面杆扭转按应力函数举例 由于 x y 满足 2 2KG s 0 因此F x y 需要满足 2F 0 F y b 2 0 F x a 2 GK y2 b2 4 GK y2 b2 4 F x y 2020 2 7 66 8 4等截面杆扭转按应力函数举例 2 根据F x y 为调合函数以及满足对称边界条件 F x y 亦采用级数形式的分离变量函数 即 Am为待定系数 2020 2 7 67 8 4等截面杆扭转按应力函数举例 3 利用边界条件F a 2 GK y2 b2 4 将GK y2 b2 4 展开为cos m y b 的级数 可将Am用GK表示 2020 2 7 68 8 4等截面杆扭转按应力函数举例 4 最后利用 将GK用Mz表示 并可确定应力分量 zx zy 具体过程参看徐芝纶 上册 P 330 333 2020 2 7 69 8 5薄壁杆的自由扭转 薄壁杆件在工程中经常碰见 它们可分为开口薄壁和闭口薄壁杆件 下面分别讨论它们的计算方法 5 1开口薄壁杆件的自由扭转 开口薄壁杆为单连域 其截面可由曲边等宽狭长矩形截面或由几个直边等宽狭长矩形截面组成 2020 2 7 70 8 5薄壁杆的自由扭转 对于曲边狭长形截面可近似以等宽的直边狭长截面代替进行计算 从薄膜比拟看两者围成的体积和最大斜率不会有多大差别 当两者受相同扭矩时 两个柱体的K和剪应力没有多大差别 2020 2 7 71 8 5薄壁杆的自由扭转 直边狭长截面剪应力计算式 2020 2 7 72 8 5薄壁杆的自由扭转 对于由几个 若干个 同样材料的狭矩形截面组成的薄壁杆 其中第i个狭矩形截面长ai 宽bi 则它应承受扭转为 总的扭转为 2020 2 7 73 8 5薄壁杆的自由扭转 则代回Mi表达式 第i个狭矩形截面上的最大剪应力为 2020 2 7 74 8 5薄壁杆的自由扭转 5 2闭口薄壁杆扭转 闭口薄壁杆为多连域 按应力函数求解时基本方程 2 2KG s0 0 si Ci 0 i 1 2 Ai为si围成的面积 2020 2 7 75 8 5薄壁杆的自由扭转 对于二连域薄壁扭转杆 一个孔洞