找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
楼主: 明月山河

解一个阀门的数学模型

  [复制链接]
发表于 2015-5-24 21:17:22 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:08
5 |$ h/ D' [1 t5 {; O我来简化一下:- {, A4 e; o2 `' e" u; l7 n  a, e. x
看似是非线性方程组,实际上不过是一元二次方程求解,初中生足以。
$ R) T+ d3 T, m1 R8 W* Y8 p  {最后一个方程:S6,α ...
3 t) }# G0 c( p! b
>> syms a b c x! q6 Q3 E2 D, l5 e' G$ _1 g
>> f=sym('a*x^2+b*x+c')& g4 [- |, d% U2 v; w4 G/ W, r

4 F8 ]4 }, |: B5 ~' wf =
# s& p% |0 r6 k. h9 `& T) K$ _* W
$ T9 l6 G2 N# s. `a*x^2 + b*x + c
& k, j3 i" a! J1 @3 p
6 @7 q# b0 j" x6 A4 ]>> finverse(f,x)
0 R: k/ `6 l6 n" K( d0 d* ]( |0 w% m6 K. [7 n( N" J
ans =
2 S( o+ r1 {/ S% ~5 k5 F( }9 B$ u5 r. ^. W4 R
-(b + (b^2 - 4*a*c + 4*a*x)^(1/2))/(2*a)
5 g7 U, |8 h( J! v( n" o& a$ O; d" C' v4 q. v  O8 {
我用matlab 推导 一元二次方程求根公式      . P$ B  z; p- Y% B( }
            
3 _* I1 m1 e6 W* h4 c               2
# h& c& T4 R5 R; J1 k  b + sqrt(b  - 4 a c + 4 a x)
4 f! P0 r* y  P0 @" i8 {- ----------------------------  }) Y4 T+ z% Y& P
               2 a8 {: f# T# P, C, X8 c  n: M5 ]

% {! w* ?: d9 V9 a/ y3 A; |: Q+ n
+ Y& x+ W" G& ?2 r/ O' e0 |. k这个也含有X  估计  可以删除; h& d8 a* K5 Y3 Z  C) {. z; S

, _# T2 z  D1 H) V$ v' L; Z4 P
发表于 2015-5-24 21:46:45 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
发表于 2015-5-24 22:48:17 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:46 + ^. t. g2 K7 I) O
糊涂了,既然如此,那就是迭代吧。1 B. Y9 W. [8 R+ K! x/ i
不管含有多少个ri,哪怕是离散的映射关系,反正有边界条件可以控制。
" V# q  B7 n8 u% ]具 ...

- p# S8 P& K! S* f+ B不想要开平方根的S   可以用权函数  之类的表达    可以吧' p, j$ l# A. P2 y

点评

试试看?  发表于 2015-5-25 00:27
 楼主| 发表于 2015-5-25 00:24:06 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:46 : C) U9 Y; V/ W4 z
糊涂了,既然如此,那就是迭代吧。
( D$ Q2 Z3 f3 E不管含有多少个ri,哪怕是离散的映射关系,反正有边界条件可以控制。- n# Q  w: ~/ c& O" p% d8 B& u
具 ...

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

点评

是比较复杂,相依参数很多。  发表于 2015-6-24 23:40
发表于 2015-7-8 18:18:44 | 显示全部楼层
你给的公式真复杂。
7 c$ {7 o" P& Z2 i7 ^$ u( i把 R6-δ 定义为R0
" O% k* [; _: |# g# p+ {S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ# d* \" X4 [; y" N8 E( m
还是3个未知,再简化下,a角把δ/2去掉,
 楼主| 发表于 2016-1-19 14:29:25 | 显示全部楼层
zms9439 发表于 2015-7-8 18:18 0 S$ f/ v, y. Z7 Q
你给的公式真复杂。
/ K$ k4 N7 v8 X) N  w( Q1 w把 R6-δ 定义为R0: H. f" Z; X+ i* k0 I* v. ^
S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ
/ u: U; l* V$ k+ A
谢谢关注,我已经解决了,用级数展开到4次项。) U, h; l4 G$ D8 n
发表于 2016-1-19 15:08:26 | 显示全部楼层
明月山河 发表于 2016-1-19 14:29
* T* x: o+ J: P+ d; z) w8 f% \谢谢关注,我已经解决了,用级数展开到4次项。

) h, |- S  I+ {/ q8 v, `/ w, @1 ?呵呵,时隔半年,终于解决了。恭喜
' L1 n* `' K% K* k( W. p$ J计算的结果可以发上来看看呀?

点评

哦,看来楼主还是很用了一番心思呀,不错,支持一下  发表于 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-10-20 00:26 , Processed in 0.064394 second(s), 14 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

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