没有合适的资源?快使用搜索试试~ 我知道了~
三角形三节点有限元分析,三角形三节点,任意编号均可,整刚存储方式为旋转存储,采用四维数组
资源推荐
资源详情
资源评论
之前上传过一篇几乎一模一样的文章,
但本人经后期使用时发现其中还有错误
之处,本想删掉以前的重新发一遍,奈
何在我的账号上显示已删,但是在文库
里仍能搜出,故本人只能在此篇修改之
后的文档开头声明此问题,前一篇文档
给大家带来的不便望大家见谅,并诚恳
地希望大家提出文档中的不足之处
惯例---先附上本人调试结果:
调试源数据如下:
结构结点总数 NJ 10
结构单元总数 NE 9
问题类型代码 L*M 0/1
弹性模量 EO 1
泊松比 AMU 0.25
结构容重 ALOU 0/1
单元厚度 TE 1
结点坐标数组 0,0/2,0/4,0/6,0/1,2/3,2/5,2/2,4/4,4/3,6
单元结点码数组 1,2,5/2,3,6/3,4,7/3,7,6/2,6,5/5,6,8/6,7,9/6,9,8/8,9,10
支杆(约束)位置数组 1,2,7,8
结点载荷位置及大小数组 10,20
四次计算所得结果如下:
第一次
"结点个数 NJ 为 10"
"单元个数 NE 为 9"
"问题类型代码 L * M(平面应力问题:L*M=0,平面应变问题 L*M=1)为 1"
"弹性模量 EO 为 1"
"泊松比 AMU 为 0.25"
"结构容重 ALOU 为 0"
"单元厚度为 1"
"结点坐标数组为为"
"结点 1","(0,0)"
"结点 2","(2,0)"
"结点 3","(4,0)"
"结点 4","(6,0)"
"结点 5","(1,2)"
"结点 6","(3,2)"
"结点 7","(5,2)"
"结点 8","(2,4)"
"结点 9","(4,4)"
"结点 10","(3,6)"
"单元结点码数组(局部 i,j,k)为"
"单元 1","(1,2,5)"
"单元 2","(2,3,6)"
"单元 3","(3,4,7)"
"单元 4","(3,7,6)"
"单元 5","(2,6,5)"
"单元 6","(5,6,8)"
"单元 7","(6,7,9)"
"单元 8","(6,9,8)"
"单元 9","(8,9,10)"
"支杆位置数组为"
"1,2,7,8"
"结点载荷大小及位置数组"
"10,20"
"结点位移数组为(x 方向,y 方向)"
"结点 1","(0,0)"
"结点 2","(.174064491200528,6.09891965373474)"
"结点 3","(4.62569442396194E-15,-8.4802828571676E-17)"
"结点 4","(0,0)"
"结点 5","(-1.15496193389917,5.22835193717885)"
"结点 6","(-.49299322794274,4.88555124978286)"
"结点 7","(2.14314760774482E-16,9.62753833846531E-16)"
"结点 8","(.296586527465448,9.23523333937133)"
"结点 9","(2.46163910513264,8.76708689404178)"
"结点 10","(1.84725926162858,24.4048794464007)"
"单元应力大小为"
"单元 1","σx=4.23870978951537E+15,σy=4.23870978951537E+15,τxy=.971385094847063"
"单元 2","σx=-.734436569166194,σy=6.96257964802094E-02,τxy=-1.33578902545555"
" 单元 3"," σ x=-1.12616993441886E-15, σ y=1.12616993441886E-15, τ
xy=-4.02745924526963E-16"
"单元 4","σx=-.19423855586746,σy=.585674402912016,τxy=-1.02640957275085"
"单元 5","σx=.302787417812317,σy=-.378787288671809,τxy=-.268168551903495"
"单元 6","σx=-.122432553765182,σy=1.28409266220844,τxy=.155552684198082"
"单元 7","σx=-.257289675295777,σy=2.07523754118727,τxy=-.43548310613577"
"单元 8","σx=.507126091375968,σy=1.28734860707857,τxy=.280791919782448"
"单元 9","σx=-1.24175643011917E-02,σy=5.28304913650959,τxy=-1.7409576384686E-16"
第二次
"结点个数 NJ 为 10"
"单元个数 NE 为 9"
"问题类型代码 L * M(平面应力问题:L*M=0,平面应变问题 L*M=1)为 1"
"弹性模量 EO 为 1"
"泊松比 AMU 为 0.25"
"结构容重 ALOU 为 1"
"单元厚度为 1"
"结点坐标数组为为"
"结点 1","(0,0)"
"结点 2","(2,0)"
"结点 3","(4,0)"
"结点 4","(6,0)"
"结点 5","(1,2)"
"结点 6","(3,2)"
"结点 7","(5,2)"
剩余13页未读,继续阅读
资源评论
ruxin1991
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功