找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
楼主: 明月山河

解一个阀门的数学模型

[复制链接]
发表于 2015-5-24 21:17:22 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:08 ( ^' z5 l+ B4 d- X
我来简化一下:
1 a8 z/ w' R% @* N: g看似是非线性方程组,实际上不过是一元二次方程求解,初中生足以。' a/ |- A- `" Z
最后一个方程:S6,α ...
1 k) q" t3 V5 s  z4 }1 Y
>> syms a b c x( ^  S- f- [1 B% D
>> f=sym('a*x^2+b*x+c')
3 h9 R9 \4 t& `& q6 B& o
: i7 J! M6 q( l! T$ h) A* W8 Ff =  M, G8 u6 Y3 f) _1 Y4 x- q

* q6 g( w4 U+ n" M; s) va*x^2 + b*x + c
3 U) E  {/ L) Y0 C( y, O: @" D8 c3 x
6 g. w0 N4 r: W# @' R) `>> finverse(f,x)- ]$ I. n) w/ W4 T5 j) w

+ q+ Y% l. \7 Y: s1 B& `4 hans =8 Q: Z/ h$ c) [* U

: z' _8 f  Q# [5 F) l; G5 g-(b + (b^2 - 4*a*c + 4*a*x)^(1/2))/(2*a)' s5 n0 V; |% C3 c- q7 D
/ @% V- }" {% a, B9 j' V6 ]+ T
我用matlab 推导 一元二次方程求根公式      3 B% D' y! Z) k$ i+ E
            
- K% P3 y) Y8 ^" X& g$ z4 N               2; R9 \) j/ A2 i2 T7 @
  b + sqrt(b  - 4 a c + 4 a x)
0 k. F- K+ i5 j! v" s, \- ----------------------------7 y4 m# g3 H% I  L* ]* G
               2 a; V( f( Y6 }0 F' G1 m

2 c" f  r' t* C: q! B: N& k
( O, N9 h3 _5 R/ z4 p这个也含有X  估计  可以删除
; P$ _' i  k) y2 }4 y. Z2 k9 t1 e4 j% o4 h* z" O& E
发表于 2015-5-24 21:46:45 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
发表于 2015-5-24 22:48:17 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:46   P* x$ m0 b$ `% b# p( H
糊涂了,既然如此,那就是迭代吧。: Q, [, _* l  P7 _3 N
不管含有多少个ri,哪怕是离散的映射关系,反正有边界条件可以控制。0 [( K& v" ]/ U0 e- z
具 ...
* i1 o9 o7 c* b2 j7 P$ q
不想要开平方根的S   可以用权函数  之类的表达    可以吧
  o7 ^" [* V3 ^/ c, k

点评

试试看?  发表于 2015-5-25 00:27
 楼主| 发表于 2015-5-25 00:24:06 | 显示全部楼层
ngsxngtd 发表于 2015-5-24 21:46
, \$ @, w" s% Z3 r7 q9 Z糊涂了,既然如此,那就是迭代吧。
% L% ?+ T) Q8 P, b0 [/ h# Q不管含有多少个ri,哪怕是离散的映射关系,反正有边界条件可以控制。
6 m0 f; y9 H0 M) `& {具 ...
% a3 t) {& [5 U# J/ E0 j% @
流量系数只有40,但是这个数字没有什么意义,这个阀是用来清洗物体的换能阀,流体高频振动的,惯性力起主要作用,而且与频率有关, 流动阻力相比之下很小,对阻抗的关注更大一些。收缩系数暂定取0.9以下。修正系数?好像用不到。
- p/ ~; A1 x6 Y# A; v( T) b5 j另外图中的通道也不是主要水流流道,而是脉冲波的传输通道。六个格子是分隔不同相位的脉冲,防止其能量混合相消。所以格子的位置关系非常严格,才有前面的模型。由于脉冲波的扩散方式还没有选定,所以Si的函数也没有确定。所以提取了一个子问题出来单独解决。其他的参数还在设计中。
1 q, l- F9 v5 {1 o* T
发表于 2015-6-24 14:37:12 | 显示全部楼层
看不懂,阀门好复杂的样子呢~

点评

是比较复杂,相依参数很多。  发表于 2015-6-24 23:40
发表于 2015-7-8 18:18:44 | 显示全部楼层
你给的公式真复杂。! Q" n4 q1 S: Z5 u8 Q8 z; c% H: v3 ]) ^+ m
把 R6-δ 定义为R0
: Y7 G* J7 F& o2 f. ^" t$ o5 ]S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ
  b5 U7 {2 S. G0 j8 ~* W还是3个未知,再简化下,a角把δ/2去掉,
 楼主| 发表于 2016-1-19 14:29:25 | 显示全部楼层
zms9439 发表于 2015-7-8 18:18
& Y+ M+ b  f' v) f- E" n你给的公式真复杂。: X3 M9 d7 ?, z1 |! W  k) _
把 R6-δ 定义为R0. h8 _! f. L0 L# a
S6=1/2*(R0^2-R7^2)*a-(R0-R7)*δ

2 U' c- h1 p7 z* f) H8 j; Q谢谢关注,我已经解决了,用级数展开到4次项。* ]- X% R, d, T* N2 S. W6 K
发表于 2016-1-19 15:08:26 | 显示全部楼层
明月山河 发表于 2016-1-19 14:29
; }4 K( N6 ~# O5 G6 N谢谢关注,我已经解决了,用级数展开到4次项。

9 ?3 Q9 o  s& `9 W呵呵,时隔半年,终于解决了。恭喜8 ?, R* N; {+ C& B4 p! i. |8 T
计算的结果可以发上来看看呀?

点评

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

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

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