代码下载链接
图示结构,单元划分如图,已知,μ、E, a,承受均布压力q (μ、E, a, q自己给数字,每人不同)
试用有限元法(平面三角形单元)划分四个单元如图,解此平面应力问题。求出应力,应变及支座反力。
(可用MATLAB, maple等编程求解,提交程序及结果)
% %
clc
clear;
%—鐗╃悊鍙傛暟------------------
E=7e10; %寮规�фā閲�,鍗曚綅Pa
t=0.01; %鍗曞厓鍘氬害,鍗曚綅m
a=1; %灏哄a锛屽崟浣峬
q=10e6; %鍧囧竷鍔涳紝鍗曚綅N/m2
miu=0.3; %娉婃澗姣�
%-------------------------------
node=[1 0 0 0;
2 1.5a 0 0;
3 3a 0 0;
4 0 a 0;
5 1.5a a 0;
6 0 2a 0;]; %鑺傜偣淇℃伅锛岀涓�鍒椾负鑺傜偣缂栧彿锛�2~4鍒楀垎鍒负x,y,z鏂瑰悜鍧愭爣
ele=[1 6 4 5;
2 1 2 4;
3 2 4 5;
4 2 3 5;]; %鍗曞厓淇℃伅锛岀涓�鍒椾负鍗曞厓缂栧彿锛屽悗闈㈠悇鍒椾负鍗曞厓涓婄殑鑺傜偣鍙风爜
n_ele=length(ele(:,1)); %鍗曞厓鏁�
代码截图
CAE有限元仿真服务