湖南大学课程考试试卷.docVIP

诚信应考,考试作弊将带来严重后果！ 湖南大学课程考试试卷 课程名称：实用有限元方法；课程编码 04615 试卷编号： 1 ；考试时间：120分钟 题 号 总分 应得分 100 实得分 评卷人 Galerkin法推导该微分方程弱形式（10分）； （2）若将整个问题域离散为个节点，场变量。试给出此二 维问题的有限元离散平衡方程，并写出刚度矩阵的具体表达式（10分）。 2.（30分）下图所示为一平面四边形单元，各节点在物理坐标系及自然坐标系中的坐标值已标示完毕。试就下列问题给出相应解答： (1) 写出自然坐标系下，该等参元的形函数；并基于等参变换，给出雅克比矩阵J的表达式（6分）； (2) 积分点I，II，III和IV在自然坐标系中的坐标分别为（-0.5,0），（0，-0.5），（0.5,0）和（0,0.5），试给出第一个积分点（I）和第二个积分点（II）位置处的应变矩阵B（14分）。 (3)若已求得四个积分点处的应变分别为，，，，求四个节点处的应变（10分）。 （20分）考虑如下图所示的一个平面应力问题，假设厚度，材料各向同性，杨氏模量，泊松比，相关力和位移边界条件已在图中标示（左端固支约束，右端施加大小为1竖直向下的集中载荷）。试用3个三角形单元（已划分完毕）分析此问题，给出本问题组装后的总体刚度矩阵和总体力向量。 4.（30分）下图所示为一由两根杆组成的桁架系统，各杆件的横截面积，弹性模量。相关力和位移边界条件已在图中标示。 试写出全局刚度矩阵及载荷向量，并求节点2竖直位移（10分）。 下面是计算节点2竖直位移的部分matlab代码，请补充完整其中缺失的部分（20分）。 相关变量定义如下： emodule：弹性模量 edof：单元自由度数 ff：系统载荷向量 nd：单元节点 area：截面面积 sdof：系统自由度数 kk：系统刚度矩阵 xcoord：单元x坐标 P：外载荷幅值 gcoord：节点坐标数组 len：单元长度 ycoord：单元y坐标 nel：单元总数 nodes：单元节点数组 Tr：单元坐标转换矩阵 cosa：单元方向余弦 nnode：节点总数 bcdof：约束的自由度 ke：单元刚度矩阵 sina：单元方向正弦 nnel：单元节点数目 bcval：自由度约束值 Bmat：单元应变矩阵 %------------------------------------------------------ % 桁架系统，计算节点2的竖直位移 %------------------------------------------------------ clear format long %------------------------------------------------------ % 输入相关控制参数 %------------------------------------------------------ emodule=1; area=1; P=1; nel=2; nnode=3; nnel=2; ndof=2; edof= ; sdof= ; %------------------------------------------------------ % 输入节点坐标信息 %------------------------------------------------------ gcoord=[0,0;1,1;2,0]; %------------------------------------------------------ % 输入单元节点信息 %------------------------------------------------------ nodes=[1,2;2,3]; %------------------------------------------------------ % 输入边界信息 %------------------------------------------------------ bcdof=[