第2讲_CAE技术基本求解过程.ppt

* 三、非线性分类及有限元法基本流程 状态变化； 几何非线性； 材料非线性。 非线性结构 几何非线性 * 固体力学问题，从本质上讲是非线性的，线性假设仅是实际问题中的一种简化 假设节点位移无限小； 材料的应力与应变关系满足虎克定律： 加载时边界条件的性质保持不变。 线性弹性体系 非线性问题 否则 三、非线性分类及有限元法基本流程 * 非线性问题 材料非线性 几何非线性 由于材料的应力与应变关系的非线性引起的 结构的位移使体系的受力状态发生了显著变化，以至不能采用线性体系的分析方法 ①大位移小应变问题; ②大位移大应变问题; ③结构的变位引起外载荷大小、方向或边界支承条件的变化等问题。 * 有限单元法分析非线性问题 三个基本步骤组 1．单元分析 2．整体组集 3．非线性方程组的求解 材料非线性问题，使用材料的非线性本构关系 几何非线性问题，计算位移高阶导数的应变－位移矩阵 整体刚度方程写成增量形式 增量法 迭代法 混合法 * 增量逐步解法的基本思想 假定 t 时刻的解已知，?t为选择的时间增量， 在t+?t时刻，有： t时刻的解为已知 外载荷节点力向量 单元应力引起的节点力向量 * 节点力增量向量 材料和几何条件的切向刚度矩阵 节点位移增量向量 * 修正Newton法的迭代公式 迭代步数 初始值 t时刻的解 迭代所用的初始值正是t时刻的解 满足精度指标时 迭代终止 * Area=b*h=12.5mm2 Izz=b*h3 /12=26.0417 E=70GPa, u=0.0 5mm 2.5mm b h P=5N 150mm 悬臂梁静力分析 教材P50： 例题问题描述：一侧固定的方板如图2.22所示，长宽均为1m，厚度为5cm，方板的右侧受到均布拉力q=200MPa的作用。材料的弹性模量为E=2.1×105MPa，泊松比为0.3。求孔板内的应力分布。 * * * * * * * * /PREP7 !* ET,1,PLANE2 !* KEYOPT,1,3,3 !* R,1,0.1, !* UIMP,1,EX, , ,2.1e11, UIMP,1,NUXY, , ,0, n,1 ,,,, n,2 ,6,,, n,3 ,6,6,, n,4 ,0,6,, !* !* E,1,3,4 e,1,2,3 FINISH /SOLU D,1,all D,2,ux ,0 d,3,ux,0 F,3,FY,-300, F,4,FY,-300, solve FINISH /POST1 * 1.ANSYS 2.ADINA 3.MARC 4.ABAQUS 5.ADMAS 6.PRO/E - CAD 7.Phoenics CFD 8. Fluent6.3 CFD 我校购置的正版工程软件 * 第2章 CAE技术基本求解过程 一、有限单元法基本知识 二、线性分析有限元法的基本计算步骤 三、非线性分类及有限元法基本流程 第2讲 * 一、有限单元法的基本概念 有限元分析（FEA） 利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 * F 一、有限单元法的基本概念 轮齿有限元模型 网格 节点 边界 单元 有限单元 节点数是有限的 单元数目也是有限的 * 按有限元规则，下面单元编号正确与否? k i j 3节点单元 6节点单元 k i j l m n 单元编号那些是正确的？ 1. k, j, i 2. i, j, k 3. i, k, j 4. j, k, i 单元编号那些是正确的？ 1. k, j, i, n, m , l 2. i, j, k, l, m, n 3. i, n, j, m, l, k 4. j, k, l, m, n, i 5. i, k, m, j, l, n * 下面单元编号正确与否? 5 3 12 10 2 1 6 7 5 4 3 8 单元编号那些是正确的？ 1. 5, 3, 12, 10 2. 3, 12, 10, 5 3. 10, 12, 3, 5 4. 10, 5, 3, 12 单元编号那些是正确的？ 1. 1, 2, 3, 4, 6, 7, 8, 5 2. 1, 2, 3, 4, 8, 5, 6, 7 3. 3, 4, 1, 2, 7, 8, 5, 6 4. 2, 3, 4, 1, 6, 7, 8