MSCNastran操作与实战培训教程概要.ppt

建立有限元模型的输入数据文件 按MSC/NASTRAN的输入数据格式，建立相应有限元模形的数据文件 第1步 确定分析类型，形成输入文件的执行控制语句集 ID MPM，CH 12 EXAMPLE SOL 101 TIME 100 CEND 第2步 离散物理模型以形成有限元模型 对承受集中力梁弯曲问题。因杆单元CROD不能承弯，故不能选用它。 CBAR是常用的等剖面直梁元，满足本问题的要求。 最后形成具有四个结点、三个CBAR单元、在结点3处作用集中力F的两端铰支梁模型 第3步 形成模型数据集 建立有限元模型的“几何”、“单元”、“材料”、“载荷”与“约束”五类数据 几何数据 GRID卡如下： 若采用自由域格式,则为: GRID , 1 , 0. , 0. , 0. GRID , 2 ,10. , 0. , 0. GRID , 3 , 20. , 0. , 0. GRID , 4 , 30. , 0. , 0. 单元数据 ??? CBAR单元卡中第6、7、8项为梁端点GA处的定位向量分量（X1，X2，X3）； 该向量确定梁元坐标系和主平面1的位置，据此计算梁剖面性质。尽管的选择有点任意，但其方向必须与梁剖面的惯性主平面相一致。 定义平面1和Y-单元轴 单元性质卡标识号为101，三个梁单元的CBAR卡数据 CBAR单元的性质PBAR卡： a)???? PBAR卡第5项和第6项。,两者不能颠倒。 ???? 应力恢复系数Ci、Di、Ei和Fi的选取，由用户在梁剖面（Y-Z平面）上任选的四点局部坐标值。 考虑横向剪切影响，需剪力面积系数K1和K2填入适当的值，对矩形剖面，K1=K2=5/6。若不填，其缺省值为无穷大，意味着剪切对挠度无影响。 对于无剪切影响，PBAR卡为 材料性质 梁的材料为钢材，弹性模量，泊桑比。采用MAT1卡来定义，材料标识号选为201 载荷数据 结点3处有100磅集中力，方向向下。用FORCE卡定义。选用载荷标识号选为10， 该FORCE卡如下 边界条件 （1）可用GRID卡第8项（PS）填入固定约束自由度分量。 （2）考虑到1和4结点约束相同，采用SPC1卡描述。SPC1卡 第4步 提出分析输出要求，建立情况控制指令集 位移和单元应力输出，如下的情况控制指令是必须的 DISP = ALL（输出全部结点位移） STRESS = ALL（打引所有单元的应力） 为了检查计算结果，采用如下指令是合适的： FORCE = ALL（打引所有单元力） SPCF = ALL（打引所有约束反力） 输入数据文件的输出采用如下指令： FCHO= BOTH （输入文件以分类和未分类两种形式输出） 输出结果页面的标题和子标题可用 TITLE和SUBTITLE指令： TITLE = HINGED BEAM SUBTITLE = WITH CONCENTRATED FORCE 对载荷和约束的选取如下： ? 标识号应与所选模型数据卡的标号相对应 最后，形成了如下情况控制指令集： CEND ECHO = BOTH DISP = ALL STRESS = ALL FORCE = ALL SPCF = ALL SPC = 100 LOAD = 10 TITLE = HINGED BEAM SUBTITLE = WITH CONCENTRATED FORCE 上述情况控知指令要求放在CEND之后，相互的顺序是任意的。 第5步 完成输入文件 运行MSC/NASTRAN并输出结果 MSC/NASTRAN的计算结果放在.F06文件中。 检验结果 错误信息 MSC/NASTRAN对有限元模型、运算过程错误信息，记录在文件.F06中。 残余力向量误差ε 残余力向量误差EPSILON是十分小的（～10-6），这表明稳定的数字特性，不存在大的舍入误差，所求解的刚度矩阵特性好（不存在奇异性、病态条件）。 位移量级判断 梁的挠度量级为英寸，这同梁的长度与剖面尺寸相比均是很小的量，符合小位移假设，结果是合理的。如果得到大于几英寸的挠度结果， 载荷太大，或是材料特性数据（弹性模量）不正确，或该问题根本不是线性问题。 平衡检查 为检验静力平衡，需计算支持点的反力 与经典理论结果的比较 第10章 模态分析 基本有限元方程 模态分析基本有限元方程 [M]和[K]分别为结构系统的质量矩阵和刚度矩阵,{u}和 分别为节点位移与加速度 解为如下的简谐运动 其中， 为模态形状， 为圆频率 等价为特征方程的非0解 有限元分析中，矩阵[K]和[M]实的对称矩阵，它们满