医工互联

 找回密码
 注册[Register]

手机动态码快速登录

手机号快速登录

微信登录

微信扫一扫,快速登录

QQ登录

只需一步,快速开始

查看: 293|回复: 0
收起左侧

外科手术器械设计 超声刀设计思路和原理 - 04

[复制链接]

  离线 

发表于 2022-10-24 10:26:53 | 显示全部楼层 |阅读模式 <
113014hwrly18lexykpywp.png

根据上文牛顿第二定律的谐振频率模型可以得到:
113014eleglrmmnr1mnlna.png

又有,在无阻尼和无外力的情况下,可简化为:
113015tlubgaagnutvf9zg.png

根据简谐振动的情况,里面的t可以用x代替;根据二阶齐次微分方程的通解可以求得振速,和应力的关系:
  1. syms k1 k2 k3 k4 k5 k6;
  2. syms l1 l2 l3 l4 l5 l6;
  3. syms A1 A2 A3 A4 A5 A6;
  4. syms B1 B2 B3 B4 B5 B6;
  5. syms Z1 Z2 Z3 Z4 Z5 Z6;
  6. syms x1 x2 x3 x4 x5 x6; % variables
  7. syms vb vf; %input velocity and output velocity
  8. v1 = A1*sin(k1*x1)+B1*cos(k1*x1); %velocity equation
  9. v2 = A2*sin(k2*x2)+B2*cos(k2*x2);
  10. v3 = A3*sin(k3*x3)+B3*cos(k3*x3);
  11. v4 = A4*sin(k4*x4)+B4*cos(k4*x4);
  12. v5 = A5*sin(k5*x5)+B5*cos(k5*x5);
  13. v6 = A6*sin(k6*x6)+B6*cos(k6*x6);
  14. f1 = -1i*Z1*(A1*cos(k1*x1)-B1*sin(k1*x1)); %force equation
  15. f2 = -1i*Z2*(A2*cos(k2*x2)-B2*sin(k2*x2));
  16. f3 = -1i*Z3*(A3*cos(k3*x3)-B3*sin(k3*x3));
  17. f4 = -1i*Z4*(A4*cos(k4*x4)-B4*sin(k4*x4));
  18. f5 = -1i*Z5*(A5*cos(k5*x5)-B5*sin(k5*x5));
  19. f6 = -1i*Z6*(A6*cos(k6*x6)-B6*sin(k6*x6));
复制代码
因为力和速度是连续的,所以可以有下面边界条件:
  1. tempA3 = solve(subs(f3,x3,0),A3);
  2. tempB3 = solve(subs(v3,x3,0)==vb, B3);
  3. tempB2 = solve(subs(v2,x2,0)==subs(subs(subs(v3,x3,l3),A3,tempA3),B3,tempB3),B2);
  4. tempA2 = solve(subs(f2,x2,0)==subs(subs(subs(f3,x3,l3),A3,tempA3),B3,tempB3),A2);
  5. tempB1 = solve(subs(v1,x1,0)==subs(subs(subs(v2,x2,l2),A2,tempA2),B2,tempB2),B1);
  6. tempA1 = solve(subs(subs(subs(v1,x1,l1),B1,tempB1),B2,tempB2)==0,A1);
  7. tempB4 = 0;
  8. tempB5 = vf*cos(k5*l5);
  9. tempA4 = vf*cos(k5*l5)/sin(k4*l4);
  10. tempA5 = vf*sin(k5*l5);
  11. %
  12. coefficient(1,1) = -(cos(k1*l1)*(vb*cos(k2*l2)*cos(k3*l3) - (Z3*vb*sin(k2*l2)*sin(k3*l3))/Z2))/sin(k1*l1);
  13. coefficient(1,2) = -(Z3*vb*sin(k3*l3))/Z2;
  14. coefficient(1,3) = 0;
  15. coefficient(1,4) = (vf*cos(k5*l5))/sin(k4*l4);
  16. coefficient(1,5) = vf*sin(k5*l5);
  17. coefficient(2,1) = vb*cos(k2*l2)*cos(k3*l3) - (Z3*vb*sin(k2*l2)*sin(k3*l3))/Z2;
  18. coefficient(2,2) = vb*cos(k3*l3);
  19. coefficient(2,3) = vb;
  20. coefficient(2,4) = 0;
  21. coefficient(2,5) = vf*cos(k5*l5);
  22. %
复制代码
再根据边界条件中多余的两个等式可求得节点平面左边和右边的频率方程 
  1. FreqEquationLeft = Z2*(vb*cos(k3*l3)*sin(k2*l2) + (Z3*vb*cos(k2*l2)*sin(k3*l3))/Z2) - (Z1*cos(k1*l1)*(vb*cos(k2*l2)*cos(k3*l3) - (Z3*vb*sin(k2*l2)*sin(k3*l3))/Z2))/sin(k1*l1);
  2. %if l1=0 the equation can simplify:
  3. FreqEquationLeft = (Z3*tan(k2*l2)*tan(k3*l3))/Z2 == 1
  4. %
  5. FreqEquationRight = (Z4*vf*cos(k4*l4)*cos(k5*l5))/sin(k4*l4) - Z5*vf*sin(k5*l5);
复制代码
 至此我们可以进行简单的计算了:
  1. %input parameters%
  2. tempZ3 = 7900*sqrt(19.5E+10/7900)*1E-6; %SS
  3. tempZ2 = 7750*sqrt(8E+10/7750)*1E-6; %PZT
  4. tempk2 = 2*pi*55E+3/sqrt(8E+10/7750);
  5. templ2 = 0.0025*4;
  6. tempk3 = 2*pi*55E+3/sqrt(19.5E+10/7900);
  7. %=>
  8. templ3 = atan(tempZ2/(tempZ3*tan(tempk2*templ2)))/tempk3
  9. %0.047%
  10. %
  11. %input parameters%
  12. ratio = vf/vb == -(tempZ2/tempZ5)*(cos(k4l4)/sin(k2l2))*sqrt((1+(tempZ5/tempZ4)^2*tan(k4*l4)^2)/(1+(tempZ2/tempZ3)^2*atan(k2*l2)^2));
  13. templ4 = 0.015;
  14. tempZ4 = 2790*(sqrt(7.15E+10/2790))*1E-6;
  15. tempZ5 = 2790*(sqrt(7.15E+10/2790))*1E-6
  16. R4/R5 = 2.4; %the radius reduction
  17. %=>
  18. templ5 = atan(tempZ4/(tempZ5*tan(tempk4*templ4)))/tempk5
  19. %0.009%
  20. vf/vb
  21. %10.7%
  22. %
复制代码
113015jlp1ty377uiu7jsm.png



来源:https://blog.csdn.net/weixin_43904073/article/details/122180126
免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!
回复

使用道具 举报

提醒:禁止复制他人回复等『恶意灌水』行为,违者重罚!
您需要登录后才可以回帖 登录 | 注册[Register] 手机动态码快速登录 微信登录

本版积分规则

发布主题 快速回复 收藏帖子 返回列表 客服中心 搜索
简体中文 繁體中文 English 한국 사람 日本語 Deutsch русский بالعربية TÜRKÇE português คนไทย french

QQ|RSS订阅|小黑屋|处罚记录|手机版|联系我们|Archiver|医工互联 |粤ICP备2021178090号 |网站地图

GMT+8, 2024-11-22 14:45 , Processed in 0.270885 second(s), 66 queries .

Powered by Discuz!

Copyright © 2001-2023, Discuz! Team.

快速回复 返回顶部 返回列表