首页 | 家园 | 百科 | 下载 | 书库 | 仪器 | 展会 |人才 | 公司 | 网址 | 问答 | 论坛 | 注册 | 通行证登录 | English

生物行

生物行   进展 | 摘要 | 人物 | 医药 | 疾病 | 技术 | 健康 | 能源   生物航   供应 | 求购 | 公司 | 展会 | 要发布
论坛   神经科学 | 神经系统疾病 | 实验技术 | 求职招聘 | 器材试剂 | 文献检索 | 读书笔记 | 考试招生 | 软件使用

[图文] 8节点等参单元matlab编程

时间:2006-08-12 20:17来源:Internet 作者:admin 阅读:

function u=u()                    % 本程序采用8结点四边形单元,分析固定位移边界条件的弹性力学平面问题   

 

fprintf('请输入弹性模量E的值(单位均采用国际单位制,下同) \n ');

E=input('');

fprintf('请输入泊松比NU的值\n ');

NU=input('');

fprintf('请输入结构的厚度hou\n ');

hou=input('');

fprintf('请输入分析类型p,平面应力问题请输入1,平面应变问题请输入2\n ');

p=input('');

while p~=1&p~=2

    fprintf('输入有误。请重新输入分析类型p,平面应力问题请输入1,平面应变问题请输入2\n ');

    p=1;

end

 

if p==1                           % 根据问题类型确定弹性矩阵D

   D=(E/(1-NU*NU))*[1 NU 0; NU 1 0;0 0 (1-NU)/2];

elseif p==2

   D=(E/(1+NU)/(1-2*NU))*[1-NU NU 0;NU 1-NU 0;0 0 (1-2*NU)/2];

end

 

fprintf('请输入结点个数nn(nn为不小于8的正整数) \n ');

nn=input('');

while nn<8

    fprintf('输入有误,请重新输入结点个数nn(nn为不小于8的正整数) \n ');

    nn=input('');

end

 

fprintf('请输入单元个数ne(ne为不小于1的正整数) \n ');

ne=input('');

while ne<1

(责任编辑:泉水)
------分隔线----------------------------
发表评论
请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
评价:
表情:
用户名: 验证码:点击我更换图片
推荐内容
  • MATLAB命令大全

    点击浏览该文件 MATLAB 命令大全 管理命令和函数 help 在线帮助文件 doc 装入超文本说...