| 查看: 1014 | 回复: 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相差很大，不知道是怎么回事，，望有能力的大虾们指点迷津，非常感激！望回复。。 |
回复此楼» 猜你喜欢
广州双一流大学08代码学硕接受调剂
已经有7人回复
-
求助：能与1,3-二氯-2-丙醇（一氯甘油；CAS号: 96-24-2）共沸的有机溶剂？
已经有3人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有107人回复
-
河南科技大学动力工程及工程热物理学硕招收调剂生（急需）
已经有1人回复
-
硕士招生-能源动力、动力工程及工程热物理、新能源科学与工程、建环、材料等工科背景
已经有2人回复
-
硕士招生-能源动力、动力工程及工程热物理、新能源科学与工程、建环、材料等工科背景
已经有1人回复
-
【博士招生】英国拉夫堡大学博士招生 （LU-CSC scholarship, 暂免ATAS）
已经有0人回复
-
北京理工大学 空天学院-招收发动机相关专业（力学、燃烧传热等方向）博士生
已经有2人回复
-
已确定
已经有2人回复
-
北京航空航天大学宇航学院 张新华副教授课题组招收 2026年申请考核制博士生
已经有2人回复
-
上海交通大学智慧能源创新学院曾美容课题组招聘博士后
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！！ANSYS带质量的模态分析。
已经有4人回复
- 【已完结】anasys瞬态分析 （悬臂梁） 命令流如下：
已经有4人回复
- 有限元分析基础教程(ANSYS算例)(曾攀)
已经有96人回复
- 底座钢结构分析 死活算不出来 求助啊
已经有10人回复
- 岩土工程ANSYS经典72例（含命令代码）
已经有284人回复
|
本帖内容被屏蔽 |
2楼2022-10-26 20:11:48
|
本帖内容被屏蔽 |
3楼2022-10-26 20:13:20
10