有限单元法习题.pdfVIP

第一题 用matlab 编程解决该问题。 程序如下： 程序 1：计算应变矩阵 function B = FEM_findB(p) %%% find the strain matrix according to certain three %%% points given %%% p1,p2,p3three points coordinates. must be vector. p1 = p(1,:); p2 = p(2,:); p3 = p(3,:); b1 = p2(2) - p3(2); c1 = p3(1) - p2(1); b2 = p3(2) - p1(2); c2 = p1(1) - p3(1); b3 = p1(2) - p2(2); c3 = p2(1) - p1(1); B11 = [b1,0;0,c1;c1,b1]; B12 = [b2,0;0,c2;c2,b2]; B13 = [b3,0;0,c3;c3,b3]; B = [B11,B12,B13]; 程序2 ：计算弹性矩阵 function D = FEM_findD(E,v) %%% find the elasticity matrix according to certain E and v %%% D elasticity matrix %%% E modulus of elasticity %%% v Poisson ratio if nargin 2 v = 0; if nargin == 0 E = 1; end end D = E/(1-v^2)*[1,v,0;v,1,0;0,0,(1-v)/2]; 程序3 ：计算总刚，单刚 function [K,Km,BE] = FEM_findK(P,lamda,mo,t,A) %%% find the total stiffness matrix of the structure %%% firstly find the element stiffness matrix %%% secondly find the structures total stiffness matrix %%% K the structures total stiffness matrix %%% Km modified total stiffness matrix %%% BE the matrix of B of all elements %%% P points coordinates matrix %%% lamda nodes number matrix of elements %%% mo modification vector indicating the removed ranks number %%% t thickness of the elements %%% A area of the elements %%% dp nodes displacement matrix of elements if nargin 5 A = 0.5; if nargin 4 t = 1; if nargin 3 error(please input P and lamda and mo!); end end end Ke = zeros(6); BE = []; K = zeros(12); dp = [2*lamda(:,1)-1,2*lamda(:,1),2*lamda(:,2)-1,2*lamda(:,2),... 2*lamda(:,3)-1,2*l