|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clc;%总润滑方程求解
8 V6 [1 {- V. L6 iclear;
3 ?% ~5 J, Y! ^7 c! Vclose all& r( g* Z) x0 G) x9 M3 f1 y( y
syms a b c d e x L O P U t o p/ h; @! n' I* \9 s! W/ [. ^% W
K=5;%常量9 u+ e6 Z( [% ?( ~
h=2;%常量 多孔储液介质厚度
5 `! j+ F# k+ U9 [& LD=0.012;%半径* J7 B, }; C) {5 T$ q% w8 j
H=0.2;%压缩后的高度1 ^) N6 F9 N8 O- X. D
L=@(x)0.2.*x;- F/ X. [$ r7 }$ V
% % % % 计算与孔隙率相关的brinkman数
) V/ i {, b/ dc=@(t)((t.^3.*D.^2)./(36.*(1-t).^2.*K)).^2;%渗透率求解6 g% L7 ]* E" R, z; V
t1=0.0033;%孔隙率) O* J. D& {- u
a=L./c(t1);%a为brinkman数
! W" w3 |! p/ jb=H./c(t1);
0 H& s+ e( p2 s; [' \0 B! t% % % % 计算方程中的f,控制此处f不与x一起变化
, r) f) @, }0 q8 ~9 Z9 If=(cos(h.*b.*H)-1)/(b.*sin(h.*b.*H));0 z$ W% U Q5 ]/ u0 p8 v7 ?& S& C/ h
% % % % 计算方程中常数余项C,此处x只与a有关' I$ m3 T/ A' g- @3 t
o=@(x)(a(x).^2.*(H-f)./(2.*f-H));5 t" H/ U4 M. ]: G. o
p=@(x)(a(x).^2./(2.*f-H));0 s: {0 K+ d# c3 o6 {
O=int(o(x),0,1);
' A8 @2 L* u% t# B0 M: VP=int(p,0,1);
9 \4 c. ?5 U7 x( @c=O./P;0 V4 C% B- E8 b5 ~
%%%%计算润滑升力 ?/ A" m3 ]- A, F# t6 E# Z. f% X% `
d=@(x)(1-x).*L./c(t1).^2.*(H-f+c)./(2.*f-H);. n% G% n: j% P
U=int(d,0,1);
" ]; ]8 s. a' ]4 z3 I/ P! L0 u我的代码一直在报错,请求高手帮忙看看哪里出了问题。; z- m' J1 o6 C0 k. G
错误提示:未定义与 'function_handle' 类型的输入参数相对应的函数 'rdivide'。% q2 A. D& n3 K2 N/ P1 t
3 Q. G/ }0 E# E0 i. i出错 runhuashengli (line 13)
: g/ Q( y! A! c& B% U* \5 v1 y4 ]a=L./c(t1);%a为brinkman数* y2 @+ g( r2 u. }7 o) u/ f
|
|