量纲分析与无量纲化的实例及matlab求解.pptVIP

数学模型 量纲分析与无量纲化实例及 MATLAB求解 数学模型 29量纲分析与无量纲化 291量纲齐次原则 物长度的量纲记L= 动力学中 理质量m的量纲记M=m 口基本量纲 量时间x的量纲记r 1,M,T 的速度v的量纲=LT1 量加速度a的量纲[a=T 口导出量纲 纲力∫的量纲=LM2 引力常数k的量纲[k]=/[m]2=L3M172 对无量纲量a,[a=1(=LMT)f=kmm2 数学模型 量纲齐次原则等式两端的量纲一致 量纲分析-利用量纲齐次原则寻求物理量之间的关系 例:单摆运动求摆动周期t的表达式 设物理量4m8t=mn12g(1) 之间有关系式 a,a2a3为待定系数,4为无量纲量 (1)的量纲表达式[t]=[m][][g 口T=Ma+aT2 ng 对比 a2+a1=0→{a2=1/2 2 1/2 数学模型 量纲齐次原则等式两端的量纲一致 量纲分析-利用量纲齐次原则寻求物理量之间的关系 例:单摆运动求摆动周期t的表达式 设物理量4m8t=mn12g(1) 之间有关系式 a,a2a3为待定系数,4为无量纲量 (1)的量纲表达式[t]=[m][][g 口T=Ma+aT2 ng 对比 a2+a1=0→{a2=1/2 2 1/2 数学模型 t=Mm1g为什么假设这种形式 设p=f(x,y,x 对xyz的两组测量值x1ypz1和x2y22 ∫(x1ypx1),P2=f(x2y xy,z的量纲单 位缩小abc倍B=f(ax,by,C3),p2=f(ax2,by2c2 f(,yu, z f(ax,by, cz,) f(ax, by2, cz,) 日p=/x形式为f(x,y,z)=xyx 数学模型 单摆运动中tm,lg的一般表达式f(t,1,,g)=0 口t”m2lgx=xy1y4为待定常数,z为无量纲量 It=LM T (MT)(LMT°)2(LMT”) [m]=LMT° Z=LM To (LMT )4=LMTO [g]=LMT-2 MT2=LMoTo y3+y=0基本解 lg=丌F(x)=0 y2 0 y-2y=0=(2,0,-1,1 (t=4V/g) 数学模型 Pi定理( Buckingham)设/(q1,q2…,qm)=0 是与量纲单位无关的物理定律,X12,…,Xn是基本量 纲,n≤m,q1,q2,…,qm的量纲可表为 [q,]=∏x“,j=1,2,…,m i=1 量纲矩阵记作A={ ai1,若 rankA 线性齐次方程组Ay=0有mr个基本解,记作 y=(ys,y2…,m),s=1,2,…,m-r ∏G 为m-r个相互独立的无量纲量,且 F(x1,m2…,xm)=0与f(q1q2,…,qn)=0等价,F未定 数学模型 量纲分析示例:波浪对航船的阻力 航船阻力f 航船速度v,船体尺寸L,浸没面积s 海水密度a重力加速度g f(q1,q2,…qn)=0q(g,l,p,v,S,f)=0 [q]=X l=LT2,团=L,[p]=L3M, Iv]=LT-, [s]=L,f]= LMT 11-3121(L A={c2} A=001001(M) i1n×1m 200-10-2(7) m=6,n=3 (g)(D)(p)(v)(s)(f) 数学模 f(q1,q2,…,qn)=0q(g,l,p,v,S,f)=0 ranka=r ranka=3 Ay=0有m个基本解y=0有m=3个基本解 y=(y31,y32,…,ysm)7 y=(-12,-1/2.0,1,0,0 1.2 y2=(0.-2,0,0,1,0) y3=( 3,-1,0,0,1) mr个无量纲量 8 丌,=∏ π tpf 数学模 F(,,)=0与 a(g,p,s/=0等价 F(丌1,2,,xn)=0与 丌1=g27 ∫(q1,q2,…,qn)=0等价 兀 S f 为得到阻力∫的显式表达式F=0dz3=v(z1,z2) 日f=1g9(x,丌2),兀1=,丌2 12W未定