结构几何非线性有限元分析与应用.ppt

结构几何非线性有限元分析与应用 丁泉顺 副研 同济大学桥梁系 一、非线性有限元方法 1、概述： 结构几何非线性问题是连续介质力学中固体力学的一个分支 几何线性与非线性问题的区别：线性问题假设结构在变形前后的受力特征是一样，而非线性问题则考虑结构在变形之后的受力特征发生变化 非线性问题的分类：大变位（大位移大转动）小应变理论，大应变理论 从理论方法上来说，结构几何非线性问题的研究已经非常完善 杆系结构几何非线性问题 2、非线性分析方法 按结构参考形态： 完全的拉格郎日方法（T.L. Total Lagrange） 更新的拉格郎日方法（U.L. Updated Lagrange） 欧拉方法（结构分析中不常用） 引入随旋坐标系： CR－TL（Corotational－TL） CR－UL（Corotational－TL） 2、分析方法(cont’l) 2、分析方法(cont’l) 3、非线性问题的有限元求解 Newton-Raphson迭代 方法 修正Newton-Raphson迭代 方法 3、非线性问题的有限元求解(Cont’l) Newton-Raphson迭代 方法 修正Newton-Raphson迭代 方法 3、非线性问题的有限元求解(Cont’l) 4、小应变时单元内力增量计算 平面梁单元的变形 4、小应变时单元内力增量计算(Cont’l) 空间梁单元的变形（详细情况见论文） 4、小应变时单元内力增量计算(Cont’l) 有限转动公式(见数学手册): 4、小应变时单元内力增量计算(Cont’l) 二、数值算例 1、坦拱的大挠度分析 1、坦拱的大挠度分析（Cont’l） 单元划分：将拱的一半分18个单元； 每步加载量为 ?P=2.0； 2、45?弯梁空间弯扭大位移分析 2、45?弯梁空间弯扭大位移分析（Cont’l） 单元划分：将梁划分为8个单元 ； 每步加载量为 ?P=10.0； 2、45?弯梁空间弯扭大位移分析（Cont’l） * * (N为轴向拉力) (3) ，， 图1 平面梁单元变形图 图2 各坐标系统示意图 1 x3 2 x31 x1 x2 x31 j j i i ( (j (i (0 图4.1 Newton-Raphson法计算示意(第一加载步) 图4.2 修正的Newton-Raphson法计算示意(第一加载步) 图4.1 Newton-Raphson法计算示意(第一加载步) 图1 平面梁单元变形图 图2 各坐标系统示意图 1 x3 2 x31 x1 x2 x31 j j i i ( (j (i (0 图4.5 随旋坐标系下从ti 时刻到ti＋(t 时刻的变形情况 图1 平面梁单元变形图 图2 各坐标系统示意图 1 x3 2 x31 x1 x2 x31 j j i i ( (j (i (0 图4(a) t+(t时刻端截面(节点1)位置图 图4(b) 单元随转坐标系的确定 2 1 端截面 1 图3(a) 矢量转动总图 图3(b) 垂直于转动轴的平面 O N P Q V ( Q P N V ( 图5 坦拱在集中荷载下的大挠度 图6 450弯梁大位移分析 R R 450 u v w P r s t 固定端 R=100.0 b=1.0 h=1.0 (=0.0 E=107 ( P H h R w0 L R=133.114 h=0.1875 b=1.0 L=34.0 H=1.09 (=7.33970 A=0.1875 I=0.0005493 E=107 (=0.2 表1 坦拱无量纲拱顶位移w/H的比较 Tab.1 Comparison of nondimensional displacements at the top of arch 加载步 5 10 15 17 18 20 Ansys 0.03522 0.08655 0.1839 0.2954 1.6073 1.6303 ADINA 0.03525 0.08677 0.1853 0.3056 1.6081 1.6316 本文解 0.03525 0.08675 0.1852 0.3044 1.6077 1.6308 图5 坦拱在集中荷载下的大挠度