找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
楼主: 明月山河

解一个阀门的数学模型

[复制链接]
发表于 2015-5-24 21:17:22 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:08
+ Z6 }1 H7 h) l6 x我来简化一下:$ v# f3 N# U1 Q+ ^* Q" a
看似是非线性方程组,实际上不过是一元二次方程求解,初中生足以。
& P6 \2 z3 \$ b5 ~% r最后一个方程:S6,α ...

% L( N; F! }' D8 s5 n; g. ]1 d$ O>> syms a b c x5 Z# r# X, v0 _
>> f=sym('a*x^2+b*x+c')1 ^: c9 O3 t+ j, z- q1 h8 U
( {# b0 n: E* k! q9 h  y- p7 L
f =
5 M( l9 J$ h) ]4 `6 h& f, Z- }
! z4 L( ?2 k+ i9 ua*x^2 + b*x + c
7 s- d7 L! c* O
5 d7 V; w* g. m8 y>> finverse(f,x)
# c+ Q6 D& I- M$ M6 F, a% ^! b! G. D4 I- D
ans =
2 G# A5 K# U% o. N
: |* v! n, v5 L7 {  f( k; |-(b + (b^2 - 4*a*c + 4*a*x)^(1/2))/(2*a)% M9 S& i7 ~: t( `! {8 ]

9 d, J5 i9 Q  S. c4 @我用matlab 推导 一元二次方程求根公式      
0 i6 y, Y2 ~% i2 ^% F! ^            
9 r, K7 \$ A! d9 s4 c1 P/ m1 {               2
  v/ i2 j7 D/ q' x) t  b + sqrt(b  - 4 a c + 4 a x)$ R+ Z  Q; X2 `9 X, _% |7 l
- ----------------------------  w# s, d; v/ I. t
               2 a' D! L% p  X  F+ h8 a

9 }" {4 i) c! W$ F% X
% I9 q. O0 h$ F5 P" z0 w8 }7 T这个也含有X  估计  可以删除4 K( T) n) r3 P( D

2 u. h% U+ e: r% H' |% r
发表于 2015-5-24 21:46:45 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
发表于 2015-5-24 22:48:17 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:46 ' h- Q  C2 k" C5 j
糊涂了,既然如此,那就是迭代吧。
# ?- d/ M  K+ k! ^& s1 k7 m7 m2 R不管含有多少个ri,哪怕是离散的映射关系,反正有边界条件可以控制。
3 a9 ]1 K' U  j具 ...
7 b7 W8 y, Y& g0 V- U1 T
不想要开平方根的S   可以用权函数  之类的表达    可以吧' T' |- k% f4 r; q5 R- h; ^  ]

点评

试试看?  发表于 2015-5-25 00:27
 楼主| 发表于 2015-5-25 00:24:06 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:46 # F! f( l: j+ U* M5 x" }" R
糊涂了,既然如此,那就是迭代吧。1 ~/ V2 A5 o- E$ C
不管含有多少个ri,哪怕是离散的映射关系,反正有边界条件可以控制。
$ Z; _8 f4 V: y1 U+ _' k具 ...

; o9 R5 W' O4 P! j流量系数只有40,但是这个数字没有什么意义,这个阀是用来清洗物体的换能阀,流体高频振动的,惯性力起主要作用,而且与频率有关, 流动阻力相比之下很小,对阻抗的关注更大一些。收缩系数暂定取0.9以下。修正系数?好像用不到。
# h# R' F: J/ O- X5 _( @另外图中的通道也不是主要水流流道,而是脉冲波的传输通道。六个格子是分隔不同相位的脉冲,防止其能量混合相消。所以格子的位置关系非常严格,才有前面的模型。由于脉冲波的扩散方式还没有选定,所以Si的函数也没有确定。所以提取了一个子问题出来单独解决。其他的参数还在设计中。* n. z. M5 j* _) l
发表于 2015-6-24 14:37:12 | 显示全部楼层
看不懂,阀门好复杂的样子呢~

点评

是比较复杂,相依参数很多。  发表于 2015-6-24 23:40
发表于 2015-7-8 18:18:44 | 显示全部楼层
你给的公式真复杂。
4 Y5 ^1 q" _. V: Y* U- `把 R6-δ 定义为R0
5 p2 F1 R& s" A2 I, qS6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ7 @" I) N# o+ Y+ [  w& z5 E' r1 b
还是3个未知,再简化下,a角把δ/2去掉,
 楼主| 发表于 2016-1-19 14:29:25 | 显示全部楼层
zms9439 发表于 2015-7-8 18:18 $ z6 z7 J$ b6 |$ s/ ?
你给的公式真复杂。
; e7 `( q- N. _- ~+ _$ C  @把 R6-δ 定义为R0
5 w. T: c1 w; U  l" |S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ

; j6 n$ K: ?  ]' H- d谢谢关注,我已经解决了,用级数展开到4次项。
! I# j0 t: C: [+ N; O0 Q7 _
发表于 2016-1-19 15:08:26 | 显示全部楼层
明月山河 发表于 2016-1-19 14:29
7 I2 S# P' U. I- c$ o谢谢关注,我已经解决了,用级数展开到4次项。
/ Q  `2 j3 d! G" T/ Z
呵呵,时隔半年,终于解决了。恭喜- _# ?4 b5 [+ Y. j9 k
计算的结果可以发上来看看呀?

点评

哦,看来楼主还是很用了一番心思呀,不错,支持一下  发表于 2016-1-19 15:22
可以呀,不过我更改个好几个设计版本,这种已经是比较老的了。 倒不用半年,用了大概几天吧。  发表于 2016-1-19 15:12
发表于 2016-1-19 18:02:47 | 显示全部楼层
好高端的样子,真心不懂。希望以后能懂

点评

会懂的。 好好学习技术,但不要学那些缺德的东西。  发表于 2016-1-19 18:19
发表于 2016-12-14 16:21:10 | 显示全部楼层
牛人就是多
您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

Archiver|手机版|小黑屋|机械社区 ( 京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号 )

GMT+8, 2025-6-17 09:23 , Processed in 0.085213 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

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