[数学]第三章 MATLAB有限元分析与应用
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
步骤2:形成单元刚度矩阵
弹簧元
调用 function y = SpringElementStiffness(k)函数
k1=SpringElementStiffness(100);
k1 =
100 -100 -100 100
k2=SpringElementStiffness(200);
k2 = 200 -200 2018/10/4 -200 200
已知边界条件:
U1 0, F2 0, F3 15
0 0 F1 100 100 100 300 200 U 0 2 U 15 200 200 3 0
2018/10/4
10
§3-1 5、实例计算分析应用
-100 100 0
K = SpringAssemble(K,k2,2,3)
2018/10/4
Βιβλιοθήκη Baidu
-100 300 -200
0 -200 200
9
§3-1 4、实例计算分析应用
步骤4:引入边界条件
弹簧元
0 U1 F1 100 100 100 300 200 U F 2 2 U F 200 200 3 0 3
function y = SpringElementForces(k,u)
%SpringElementForces This function returns the element nodal force % vector given the element stiffness matrix k % and the element nodal displacement vector u.
步骤5:解方程
弹簧元
U=zeros(2,1); F=[0;15]; K = K(2:3,2:3);
K(1,:)=[]; K(:,1)=[];
300 200 U 2 0 200 200 U 15 3
U=K\F
2018/10/4
U=inv(K)*F
弹簧元
弹簧元是总体和局部坐标一致的一维有限单元 每个弹簧元有两个节点(node) 单刚矩阵为:
k k k
k k
2 2
总刚矩阵: 结构方程:
n n
K U F
单元节点力: f k u
2018/10/4 3
§3-1 3、MATLAB函数编写
8
§3-1 4、实例计算分析应用
步骤3:集成整体刚度矩阵
弹簧元
调用 function y = SpringAssemble(K,k,i,j)函数
n=3; K = zeros(n,n);
K= 0 0 0
K= 100 -100 0 K= 100 -100 0
0 0 0
0 0 0 0 0 0
K = SpringAssemble(K,k1,1,2)
y = [k -k ; -k k];
2018/10/4
4
§3-1 3、MATLAB函数编写
3.2 整体刚度矩阵的形成
弹簧元
function y = SpringAssemble(K,k,i,j)
%SpringAssemble This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % This function returns the global stiffness matrix K % after the element stiffness matrix k is assembled.
U= 0.1500 0.2250
11
§2-1 5、实例计算分析应用
步骤6:后处理
U=[0;U] U=
0 0.1500 0.2250 F=K*U F= -15.0000 0.0000 15.0000
2018/10/4
弹簧元
u1=U(1:2); f1=SpringElementForces(k1,u1); f1 = -15.0000 15.0000
u2=U(2:3); f2=SpringElementForces(k2,u2); f2 = -15.0000 15.0000
12
§3-1 5、实例计算分析应用
第三章 MATLAB有限元分析与应用
§3-1
§3-2 §3-3 §3-4 §3-5 §3-6
弹簧元
线性杆元 二次杆元 平面桁架元 空间桁架元 梁元
2018/10/4
1
§3-1 1、有限元方法的步骤:
弹簧元
离散化域
形成单刚矩阵 集成整体刚度矩阵 引入边界条件
求解方程
后处理
2018/10/4 2
§3-1 2、基本方程
3.1 单元刚度矩阵的形成
弹簧元
function y = SpringElementStiffness(k)
%SpringElementStiffness This function returns the element stiffness %matrix for a spring with stiffness k. %The size of the element stiffness matrix is 2 x 2.
K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K;
2018/10/4 5
§3-1 3、MATLAB函数编写
3.3 节点载荷计算
弹簧元
y = k * u;
2018/10/4
6
§3-1 4、实例计算分析应用
弹簧元
如图所示二弹簧元结构,假定k1=100kN/m,k2=200kN/m,P=15kN。 求:系统的整体刚度矩阵; 节点2、3的位移; 节点1的支反力; 每个弹簧的内力
解:
步骤1:离散化域
2018/10/4
7
§3-1 4、实例计算分析应用
弹簧元
调用 function y = SpringElementStiffness(k)函数
k1=SpringElementStiffness(100);
k1 =
100 -100 -100 100
k2=SpringElementStiffness(200);
k2 = 200 -200 2018/10/4 -200 200
已知边界条件:
U1 0, F2 0, F3 15
0 0 F1 100 100 100 300 200 U 0 2 U 15 200 200 3 0
2018/10/4
10
§3-1 5、实例计算分析应用
-100 100 0
K = SpringAssemble(K,k2,2,3)
2018/10/4
Βιβλιοθήκη Baidu
-100 300 -200
0 -200 200
9
§3-1 4、实例计算分析应用
步骤4:引入边界条件
弹簧元
0 U1 F1 100 100 100 300 200 U F 2 2 U F 200 200 3 0 3
function y = SpringElementForces(k,u)
%SpringElementForces This function returns the element nodal force % vector given the element stiffness matrix k % and the element nodal displacement vector u.
步骤5:解方程
弹簧元
U=zeros(2,1); F=[0;15]; K = K(2:3,2:3);
K(1,:)=[]; K(:,1)=[];
300 200 U 2 0 200 200 U 15 3
U=K\F
2018/10/4
U=inv(K)*F
弹簧元
弹簧元是总体和局部坐标一致的一维有限单元 每个弹簧元有两个节点(node) 单刚矩阵为:
k k k
k k
2 2
总刚矩阵: 结构方程:
n n
K U F
单元节点力: f k u
2018/10/4 3
§3-1 3、MATLAB函数编写
8
§3-1 4、实例计算分析应用
步骤3:集成整体刚度矩阵
弹簧元
调用 function y = SpringAssemble(K,k,i,j)函数
n=3; K = zeros(n,n);
K= 0 0 0
K= 100 -100 0 K= 100 -100 0
0 0 0
0 0 0 0 0 0
K = SpringAssemble(K,k1,1,2)
y = [k -k ; -k k];
2018/10/4
4
§3-1 3、MATLAB函数编写
3.2 整体刚度矩阵的形成
弹簧元
function y = SpringAssemble(K,k,i,j)
%SpringAssemble This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % This function returns the global stiffness matrix K % after the element stiffness matrix k is assembled.
U= 0.1500 0.2250
11
§2-1 5、实例计算分析应用
步骤6:后处理
U=[0;U] U=
0 0.1500 0.2250 F=K*U F= -15.0000 0.0000 15.0000
2018/10/4
弹簧元
u1=U(1:2); f1=SpringElementForces(k1,u1); f1 = -15.0000 15.0000
u2=U(2:3); f2=SpringElementForces(k2,u2); f2 = -15.0000 15.0000
12
§3-1 5、实例计算分析应用
第三章 MATLAB有限元分析与应用
§3-1
§3-2 §3-3 §3-4 §3-5 §3-6
弹簧元
线性杆元 二次杆元 平面桁架元 空间桁架元 梁元
2018/10/4
1
§3-1 1、有限元方法的步骤:
弹簧元
离散化域
形成单刚矩阵 集成整体刚度矩阵 引入边界条件
求解方程
后处理
2018/10/4 2
§3-1 2、基本方程
3.1 单元刚度矩阵的形成
弹簧元
function y = SpringElementStiffness(k)
%SpringElementStiffness This function returns the element stiffness %matrix for a spring with stiffness k. %The size of the element stiffness matrix is 2 x 2.
K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K;
2018/10/4 5
§3-1 3、MATLAB函数编写
3.3 节点载荷计算
弹簧元
y = k * u;
2018/10/4
6
§3-1 4、实例计算分析应用
弹簧元
如图所示二弹簧元结构,假定k1=100kN/m,k2=200kN/m,P=15kN。 求:系统的整体刚度矩阵; 节点2、3的位移; 节点1的支反力; 每个弹簧的内力
解:
步骤1:离散化域
2018/10/4
7
§3-1 4、实例计算分析应用