7结构动力分解.ppt

7.6.2 ANSYS算例2 （2）模态分析 分析类型：modal; 分析选项：子空间法subspace，求解8阶，扩展8阶。 模态分析：当前分析 模态显示：读入结果/by pick选择（set,1,n)）需显示的模态；绘制变形图（PLDISP,2） 7.6.2 ANSYS算例2 （3）谱分析 分析类型：specturm; 分析选项：单点，模态数8个 模态阻尼：各阶0.02； 地震反应谱定义： 1）激励设置setting：地震加速度，方向X：sedx=1 2）频率表 3）反应谱值 7 度多遇，地震影响系数 0.08 ，III 场地， 第一组 ：TG＝0.45 FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1.00,1.25,1.667 SV,0.05,0.1544,0.1916,0.2102,0.2242,0.2507,0.329,0.5149,0.4844,0.5987,0.7869, FREQ,2.222,10,100000 SV,0.05,1.0342,1.0342,0.3528 模态组合：SRSS 7.6.2 ANSYS算例2 反应谱分析:当前分析 结果显示： 1）读入结果：File菜单/Read input from: *.mcom 2）绘制变形图：PLDSP，2 （4）时程分析 分析类型：transient；方法：full 分析选项：输入瑞雷阻尼α、δ=2ξf1： 2*dampr*freq1*6.2832 荷载步默认1000修改：FINISH; /CONFIG,NRES,10000 地震波定义(parameters)： 1）参数定义：点数NT=1001；间隔DT=0.02； 数组AC,行数NT 2)读入地震波数据：数组参数中指定文件名。格式填入(F8.4) 7.6.2 ANSYS算例2 时程分析： *DO,I,1,1001 ACEL,AC(I),0,0 ！X加速度激励 TIME,I*0.02 ！时间不长0.02s OUTRES,ALL,ALL ！输出所有结果 SOLVE *ENDDO 结果显示(TimeHist Postpro) 1）定义时程变量（+）：X向位移 2）绘制时程图： 选择节点20，时程变量定义窗口（3th按钮）绘制 第7章 结构动力问题的有限元法 7.1 结构的动力方程 单元内任意点位移： 单元内任意点速度： 单元内任意点加速度： 单位体积惯性力： 单位体积阻尼力： 等价结点惯性力： 等价结点阻尼力： 将式（7-1a）、（7-1b）代入上式，有 单元质量矩阵 单元阻尼矩阵 结构动力方程 将惯性力、阻力作为载荷，得到用结点位移表示的单元平衡方程 按有限元的集合方法，得到结构动力平衡方程 7.2 单元质量矩阵 两种计算方法 1）集中质量矩阵：一种是用静力等效的方法将单元质量团聚在各个结点上，形成一个对角质量矩阵，计算简单，但只适用于低阶单元 2）一致质量矩阵：采用式（7-2）进行计算，得到的质量矩阵称，计算较为复杂，但可适用于各种不同单元 。 7.2.1 杆单元 1．集中质量矩阵 按静力等效方法将质量平均分配到结点 i, j i j ρ，A，l 2．一致质量矩阵 杆单元的形函数为 杆单元的一致质量矩阵 7.2.2 梁单元 1．集中质量矩阵 按静力等效方法将质量平均分配到结点，并略去转动项，则梁单元的集中质量矩阵为 2．一致质量矩阵 梁单元的形函数为 7.2.3 平面三结点三角形单元 1．集中质量矩阵 三角形单元面积为A，厚度为h，密度为ρ,按照静力等效方法，将单元质量平均分配给每一个结点，得到单元的集中质量矩阵为 2．一致质量矩阵 单元的形函数[N] 7.3 阻尼矩阵 （1）当阻尼力正比于质点运动速度时，单位体积上作用的阻尼力为 单元结点上的阻尼力为 单元阻尼矩阵正比于单元质量矩阵 （2）当阻尼力正比于应变速度时，阻尼应力可表示为 单元中有初应力引起单元中的势能为 初应力引起单元中的结点力为 由阻尼应力引起的单元结点上的阻尼力为 单元阻尼矩阵正比于单元刚度矩阵 （3）一般情形 单元阻尼矩阵可看作是质量矩阵和刚度矩阵的线性组合 （瑞雷阻尼） 取两阶频率ωi， ωj ，带入上式 7.4 特征值问题 结构中各点的振幅不能全为零，则上面齐次线性方程组中系数矩阵的行列式必须等于零，即特征值方程 算例：图示的阶梯型悬臂杆，试求其固有频率和特征向量。 将杆划分为2个单元3个结点，单元刚度矩阵为： 2A A L/2 L/2