有限单元法第三章研究.pptVIP

* * * * * * * 由于N 足够大，可以近似认为 同时得到 求出位移 之后，即可以求出结构的应力场 。 3.5 约束处理及求解 第三章 用有限单元法计算空间刚架结构，在原理上及推导过程与计算平面刚架结构相同。由于空间的每一结点一般具有六个自由度，故计算较之复杂。 设两杆的杆长和截面尺寸相同， 杆件长 m。 第三章 3.6 计算示例 结构离散化后 将结构划分为4个结点、3个单元 截面积 ，惯性矩 (2) 求结点载荷 首先须求局部坐标系中固定端内力 单元1作为两端固定梁反力示意图 单元2作为两端固定梁反力示意图 第三章 3.6 计算示例 单元1 单元2 在局部坐标系下单元载荷列向量 单元1 单元2 单元3 第三章 3.6 计算示例 为了求出在整体坐标下的载荷列向量，先求单元得坐标转换矩阵 单元1、2 单元3 第三章 3.6 计算示例 求各单元在整体坐标下的等效结点载荷 第三章 3.6 计算示例 求刚架的等效结点载荷 第三章 3.6 计算示例 无结点载荷作用时，总结点载荷即为等效结点载荷。 (3) 求单元刚度矩阵 由于单元1、2、3的尺寸相同，材料弹性模量相同，故 梁单元的局部坐标下的刚度矩阵表达式 第三章 3.6 计算示例 （4）求整体坐标系中的 单元1 单元2 单元3 第三章 3.6 计算示例 （5）求结构整体刚度矩阵 利用刚度集成法 第三章 3.6 计算示例 （6）建立原始平衡方程式 第三章 3.6 计算示例 （7）引入约束条件解方程组 由于1、3、4为固定端， 修改整体刚度矩阵中的1~3，6~12行与列， 以及载荷列向量中的相应的行，约束处理。 建立基本平衡方程 得到 （8）求各杆的杆端力 单元3结点位移列向量 第三章 3.6 计算示例 单元1杆端内力计算 单元2杆端内力计算 单元3杆端力计算 第三章 3.6 计算示例 （9）作内力图 （a） 刚架轴力图 （b） 刚架剪力图 （c） 刚架轴弯矩图 第三章 3.6 计算示例 3.7 ANSYS桁架结构计算示例 三种常用的单元 Link1-二维轴向拉伸-压缩单元，具有两个结点，不计杆件的弯曲和扭转变形。模拟桁架、连杆和弹簧。 Link8是三维轴向的拉伸-压缩杆，具有两个结点，每个结点有三个评议自由度，模拟两端铰接的空间杆件，机械系统的连杆、三维线弹簧等，不考虑弯曲和扭转。 三维杆单元Link10具有双线性刚度矩阵，只能承受单向手拉活受压。当它是单向受拉选项时，只要受压（与松弛的绳索和链条相似），它的刚度变为0，适用于分析静态的电缆、钢索等问题，也可用于动态问题的分析中。计入惯性和阻尼的影响。 第三章 3.7 ANSYS桁架结构计算示例 =1m; =1m; 材料为Q235; （1）选择单元类型 运行PreprocessorElement TypeAdd/Edit/Delete 在结点8上施加竖直向下的集中载荷F＝60000N， 约束为结点1处约束X,Y方向自由度，结点5处约束Y方向自由度。 桁架结构示意图 桁架各单元横截面图 第三章 单元类型库对话框 单元类型库对话框 （2）设置材料属性 运行PreprocessorMaterial PropsMaterial Models 选择材料属性对话框 设置材料1属性对话 （3）设置单元截面形式 选择菜单PreprocessorElementLink8 截面设置对话框 3.7 ANSYS桁架结构计算示例 第三章 （4）定义实常数 运行Real ConstantsAdd/Edit/Delete ，定义面积为4.91E-4m2 设置LINK1单元的实常数 （5）建立模型 首先生成结点，运行主菜单PreprocessorModeling Create Nodes In Active CS; 再生成单元，运行主菜单 PreprocessorModelingCreateElementsAuto NumberedThru Nodes穿越结点命令。 创建结点对话框 3.7 ANSYS桁架结构计算示例 第三章 通过结点