| 查看: 958 | 回复: 2 | ||
[求助]
悬臂梁的模态分析【传递矩阵法】 已有1人参与
|
|
悬臂梁受定常轴向力的作用，将悬臂梁分成50段，每一段悬臂梁的传递矩阵方程：Z(i+1)=[tr1]i*[tr2]i* Z(i); 其中[tr1]i=[sin(a1*l) cos(a1*l) sinh(a2*l) cosh(a2*l); a1*cos(a1*l) -a1*sin(a1*l) a2*cosh(a2*l) a2*sinh(a2*l); -E*I*a1^2*sin(a1*l) -E*I*a1^2*cos(a1*l) E*I*a2^2*sinh(a2*l) E*I*a2^2*cosh(a2*l); (-E*I*a1^3)*cos(a1*l) (E*I*a1^3)*sin(a1*l) (E*I*a2^3)*cosh(a2*l) (E*I*a2^3)*sinh(a2*l)]; [tr2]i=(a1^2+a2^2).*[0 a2^2/a1 0 -1/(E*I*a1); a2^2 0 -1/(E*I) 0; 0 a1^2/a2 0 1/(E*I*a2); a1^2 0 1/(E*I) 0]; 那么总传递矩阵是不是 trr=(tr1.*tr2)^50; 不知理解正确否，请各位大虾指点。。 如果正确，我按照这种思路编了相应的程序： clc clear syms E I a1 a2 m l s omiga f pi tr1 tr2 trr trr1 ; %文献悬臂梁的参数 E=2.1856e11; I=5e-8; l=0.44; m=3/l; s=1e4; omiga=2*pi*f; %传递矩阵中的各个参数 a1=((sqrt(s^2+4*m*omiga^2*E*I)-s)/(2*E*I))^(1/2); a2=((sqrt(s^2+4*m*omiga^2*E*I)+s)/(2*E*I))^(1/2); %由边界条件得到的总传递矩阵（这里只是trr=tr1.*tr2，个人感觉trr的50次与一次得到的特征方程的解相差不大） trr=[0, a2*cosh(a2*l)*(a1^16 + a2^16); -(cosh(a2*l)*(- E*I*a2^3 + s*a2)*(a1^16 + a2^16))/(E*I), 0]; trr1=det(trr); solve(vpa(trr1,2)) 得到的结果为 84.071，与文献中116.267相差很大，不知道是怎么回事，，望有能力的大虾们指点迷津，非常感激！望回复。。 |
回复此楼» 猜你喜欢
英国巴斯大学 University of Bath 博士后及博士招聘
已经有5人回复
-
求调剂收留
已经有1人回复
-
金属材料论文润色/翻译怎么收费?
已经有193人回复
-
罕见病药小批量生产，高腐蚀工况下钽衬里是否最优解？
已经有0人回复
-
085500机械专硕 331分求调剂
已经有1人回复
-
机械专硕311求调剂
已经有0人回复
-
求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！！ANSYS带质量的模态分析。
已经有4人回复
- 【已完结】anasys瞬态分析 （悬臂梁） 命令流如下：
已经有4人回复
- 有限元分析基础教程(ANSYS算例)(曾攀)
已经有96人回复
- 底座钢结构分析 死活算不出来 求助啊
已经有10人回复
- 岩土工程ANSYS经典72例（含命令代码）
已经有284人回复
|
本帖内容被屏蔽 |
2楼2022-10-26 20:11:48
|
本帖内容被屏蔽 |
3楼2022-10-26 20:13:20
