| 查看: 963 | 回复: 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相差很大，不知道是怎么回事，，望有能力的大虾们指点迷津，非常感激！望回复。。 |
回复此楼» 猜你喜欢
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有14人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有8人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有122人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有6人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有6人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有0人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有6人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有6人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有7人回复
-
【博导招生】上海大学张建华课题组-招收2026年秋季“申请-考核制”博士生一名
已经有0人回复
-
北京高校副校长团队招收机械类，环境类学硕和专硕
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！！ANSYS带质量的模态分析。
已经有4人回复
- 【已完结】anasys瞬态分析 （悬臂梁） 命令流如下：
已经有4人回复
- 有限元分析基础教程(ANSYS算例)(曾攀)
已经有96人回复
- 底座钢结构分析 死活算不出来 求助啊
已经有10人回复
- 岩土工程ANSYS经典72例（含命令代码）
已经有284人回复
|
本帖内容被屏蔽 |
3楼2022-10-26 20:13:20
|
本帖内容被屏蔽 |
2楼2022-10-26 20:11:48
