2-2 图示悬臂板,属于平面应力问题,其网格图及单元、节点编号见图2-1,E=2、1×1011,u=0、28,演算其单刚阵到总刚阵得组集过程,并用MATLAB 软件计算总刚阵。

图2-1答:根据图2-1所示列出单元节点列表:i j k 1 3 5 4 2 2 5 3 3 2 6 5 4162(1)计算单元刚度阵 单元1得刚度矩阵: ,; 单元2得刚度矩阵:,; 单元3得刚度矩阵:,; 单元4得刚度矩阵:,; 总刚度矩阵:[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡++++++++++++=46,636,635,642,632,641,636,535,525,515,514,523,513,532,522,515,414,413,425,315,314,323,313,322,346,236,235,225,223,242,232,222,241,246,142,141,100000000000k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K 节点单元Matlab 程序语言得编写:function Idexglobal gNode gElement gMaterialgNode=[0、0 0、010、5 0、011、0 0、011、0 0、00、5 0、00、0 0、0]%gNode 同样就是一个矩阵,每一行表示一个结点,第1 列就是结点得x 坐标,第2 列就是结点得y坐标gElement=[3 4 52 3 52 5 61 2 6 ];%gElement 就是一个矩阵,每一行表示一个单元,第1 行就是单元得第1 个结点号,第2 行就是单元得第2个结点号。

Returnfunction k=StiffnessMatrix(ie)%计算单元刚度矩阵函数global gNode gElementk=zeros(6,6); %6x6单元刚阵E=2、1*10^11; %材料特性u=0、28 ; %材料特性t=0、01; %材料特性xi=gNode(gElement(ie,1),1);yi=gNode(gElement(ie,1),2);xj=gNode(gElement(ie,2),1);yj=gNode(gElement(ie,2),2);xm=gNode(gElement(ie,3),1);ym=gNode(gElement(ie,3),2); %计算节点坐标分量ai=xj*ym-xm*yj;aj=xm*yi-xi*ym;am=xi*yj-xj*yi;bi=yj-ym;bj=ym-yi;bm=yi-yj;ci=-(xj-xm);cj=-(xm-xi);cm=-(xi-xj);d=[1,xi,yi;1,xj,yj;1,xm,ym];area=det(d); %计算单元面积B=[bi 0 bj 0 bm 0 ;0 ci 0 cj 0 cm;ci bi cj bj cm bm];B=B/2/area;D=[1 u 0;u 1 0;0 0 (1-u)/2];D=D*E/(1-u^2);k=transpose(B)*D*B*t*abs(area); %计算单元刚度矩阵Returnfunction gK=AssembleStiffnessMatrix% 计算总刚阵global gElement gK iegK=zeros(12,12);for ie =1:1:4 %单元循环k=StiffnessMatrix(ie);for i=1:1:3 %节点循环for j=1:1:3 %节点循环for p=1:1:2 %自由度循环for q=1:1:2 %自由度循环m=(i-1)*2+p; %每个节点有2个自由度,i节点得第p个自由度为(i-1)*2+pn=(j-1)*2+q; %每个节点有2个自由度,i节点得第p个自由度为(i-1)*2+pM=(gElement(ie,i)-1)*2+p;N=(gElement(ie,j)-1)*2+q;gK(M,N)=gK(M,N)+k(m,n);endendendendendReturn则单元1得刚度矩阵为>> StiffnessMatrix(1)ans =1、0e+010 *2、0508 0 -2、0508 0、0410 0 -0、04100 5、6966 0、0319 -5、6966 -0、0319 0-2、0508 0、0319 2、0531 -0、0729 -0、0023 0、04100、0410 -5、6966 -0、0729 5、6974 0、0319 -0、00080 -0、0319 -0、0023 0、0319 0、0023 0-0、0410 0 0、0410 -0、0008 0 0、0008 单元2得刚度矩阵>> StiffnessMatrix(2)ans =1、0e+010 *2、0531 -0、0729 -2、0508 0、0319 -0、0023 0、0410-0、0729 5、6974 0、0410 -5、6966 0、0319 -0、0008 -2、0508 0、0410 2、0508 0 0 -0、04100、0319 -5、6966 0 5、6966 -0、0319 0-0、0023 0、0319 0 -0、0319 0、0023 00、0410 -0、0008 -0、0410 0 0 0、0008 单元3得刚度矩阵为>> StiffnessMatrix(3)ans =1、0e+010 *0、0023 0 -0、0023 0、0319 0 -0、03190 0、0008 0、0410 -0、0008 -0、0410 0-0、0023 0、0410 2、0531 -0、0729 -2、0508 0、03190、0319 -0、0008 -0、0729 5、6974 0、0410 -5、69660 -0、0410 -2、0508 0、0410 2、0508 0-0、0319 0 0、0319 -5、6966 0 5、6966单元4得刚度矩阵>> StiffnessMatrix(4)ans =1、0e+010 *2、0531 -0、0729 -2、0508 0、0319 -0、0023 0、0410-0、0729 5、6974 0、0410 -5、6966 0、0319 -0、0008-2、0508 0、0410 2、0508 0 0 -0、04100、0319 -5、6966 0 5、6966 -0、0319 0-0、0023 0、0319 0 -0、0319 0、0023 00、0410 -0、0008 -0、0410 0 0 0、0008总刚度矩阵为ans =1、0e+011 *Columns 1 through 80、2053 -0、0073 -0、0002 0、0041 0 00 0-0、0073 0、5697 0、0032 -0、0001 0 00 0-0、0002 0、0032 0、4106 -0、0073 -0、0002 0、0041 0 00、0041 -0、0001 -0、0073 1、1395 0、0032 -0、00010 00 0 -0、0002 0、0032 0、2053 0 -0、2051 0、00410 0 0、0041 -0、0001 0 0、5697 0、0032 -0、56970 0 0 0 -0、2051 0、0032 0、2053 -0、00730 0 0 0 0、0041 -0、5697 -0、0073 0、56970 0 -0、4102 0、0073 0 -0、0073 -0、0002 0、00320 0 0、0073 -1、1393 -0、0073 0 0、0041 -0、0001-0、2051 0、0041 0 -0、0073 0 00 00、0032 -0、5697 -0、0073 0 0 00 0Columns 9 through 120 0 -0、2051 0、00320 0 0、0041 -0、5697-0、4102 0、0073 0 -0、00730、0073 -1、1393 -0、0073 00 -0、0073 0 0-0、0073 0 0 0-0、0002 0、0041 0 00、0032 -0、0001 0 00、4106 -0、0073 -0、0002 0、0041-0、0073 1、1395 0、0032 -0、0001-0、0002 0、0032 0、2053 00、0041 -0、0001 0 0、56972-3 在平面问题有限元分析中,(1)用到了哪些弹性力学中得基本方程？答:平衡微分方程、几何方程、相容方程(形变协调方程)。

(2)力得平衡条件就是如何满足得？答:根据能量守恒原理,有外力所作虚功应该等于内力虚功。

也就就是结构在外载荷作用下处于平衡状态则在结构上得力在任意虚功位移上所作得虚功之与等于零。

以下就是用到得方程:(3)变形协调条件就是如何满足得？答:对材料进行线弹性与各向同性得假设,用弹性力学中应力-应变之间得关系得到变形协调条件。

下面就是形变协调方程。

2-4 在平面三角形单元中得位移、应变、应力具有什么特征？位移特征:(1)必须包含单元得刚体位移;(2)必须包含单元得常应变状态;(3)必须保证不偏惠各坐标轴;(4)必须保证单元内位移连续。

应力特征:(1)三角形单元其应力仅与单元材料与几何尺寸有关,与节点位移有关,而与单元内位置坐标无关,也即这类单元内得应力就是常量。

(2)三角形单元内应力连续,但在公共边界上应力有突变,密布网格可以减少这种冲突得不合理性。

应变特征:由于简单三角形单元取线性位移模式,其应变矩阵为常数矩阵,即在这样得位移模式下,三角形单元内得应变为某一常量。

2-5 在平面三角形单元中,当尺寸逐步缩小,单元中得位移、应变、应力具有什么特征？当单元尺寸逐步减小时,单元各点得应变趋于相等,这时常量应变成为主要成分,因此,位移应能反应这种常应变状态,由于应力矩阵也就是常数矩阵,单元应力也就是常量。

但就是相邻单元一般具有不同得力,在单元得公共边上会有应力突变,随着单元尺寸得逐步减小,这种突变会急剧降低,从而不会妨碍有线单元法得简答收敛于精确解。