计算力学双语编程大作业-YongliangWang.pdf

计算力学(双语)编程大作业 Programming Homework 班级：分析力学15-1 姓名： 徐逍遥 学号： 1510610119 时间： 2018 年6 月 目录 1. 问题说明4 2. 前处理过程4 3. 求解问题5 3.1 控制方程5 3.2 离散方程5 3.4 弹性矩阵6 3.5 应变位移矩阵6 3.6 单元刚度矩阵7 3.7 数值积分9 3.8 边界荷载 10 3.9 位移边界条件 13 3.10 节点位移函数 13 3.11 PHW-2018 问题的计算程序 14 4. 算例 14 4.1 HW21 检验 14 4.2 手算结果 15 4.3 MATLAB 程序运行结果 15 4.4 结果对比 17 5. PHW-2018 题目计算 17 5.1 输入程序 17 5.2 位移结果分析 18 5.3 应变结果分析（以各单元中心点为例分析）20 5.4 应力结果分析（以各单元中心点为例分析）20 6. 感想21 7. 完整的代码附录22 1. 问题说明 计算如图所示二维平面结构，可以简化为平面应力问题求解。该模型上部施加分布集度 为F 的载荷，右部施加最大集度为F 的倒三角形载荷，左侧在图示位置施加两个大小为Q 的集中力。板件的尺寸如图示。模型尺寸如下表所示。 Table.1 模型材料参数表 宽度a 高度b 集中力Q 荷载F 弹性模量E 泊松比 （m ） （m ） （kN ） （kN/m ） （GPa ） μ 7 4 40 7 21.92 0.28 Fig.1 待分析模型 采用有限单元法，编程采用 MATLAB ，算法依据课本上平面问题的弹性力学解法。思 路是将《THE FINITE ELEMENT METHOD》中涉及到的求解问题的公式全部转化为矩阵的 运算，编制程序，求解问题。 2. 前处理过程 将上图所示结构进行有限元网格划分，结点和单元编号如图所示。 Fig.2 单元划分后的模型示意图 如图所示的单元划分与结点编号不满足带宽最小的条件，因此不是最优划分，但是本问 题计算量小，这种划分计算简便，而带来的效率问题可以忽略，故采用此种划分方法。 每个单元的定位向量如下： cn {1}=[7 8 9 10 3 4 1 2]; cn {2}=[9 10 11 12 5 6 3 4]; cn {3}=[13 14 15 16 9 10 7 8]; cn {4}=[15 16 17 18 11 12 9 10]; cn {5}=[19 20 21 22 15 16 13 14]; cn {6}=[21 22 23 24 17 18 15 16]. 3. 求解问题 3.1 控制方程 图示平面问题的控制方程为： ′′ − − = 对于静力问题，不考虑惯性项，上述方程化简为： + = 3.2 离散方程 经推导得单元离散代数方程