2016-05-15 13 views
-4

直径dの変数dを解く際に実際に複雑な方程式があります。以前はMatlabを使用していませんが、これは私が利用できる唯一のリソースです。私が前にVBAでプログラミングを行っているMatlabを使用して代数方程式を解く

eq = 40; 
eq = 1+((2*d)/3); 
answer = solve(d) 

...基本たとえばだから、私は助けのためのMatlab

+0

「d」の式を解析的に見つけるのを止めるのは......? –

答えて

-2

感謝の正しい構文はよく分からないが、私はそれを考え出しました。

syms d eq; 
eq = 40 
eqn = 1+((2*d)/3)-eq==0; 
solve(eqn,d) 

完全な問題。

syms d eq Ktaxial Ktbending Kttorsion p1 p2 T sigxpoint1 sigxpoint2 Transhear torxypoint1 torxypoint2 torshear pie M Length Bending Area Q I; 
%Loading level 1 
eq = 41904.7617; %Constant 
Ktaxial = 1.95; %Constant 
Ktbending = 1.7; %Constant 
Kttorsion = 1.45; %Constant 
P1 = 9453; %Change 
P2 = 493; %Change 
T = 298; %Change 
Length = 22; %Constant 

%Constants, do not change 
pie = 3.141592654; 
Area = ((pie/4)*(d)^2); 
Ybar = ((2*d)/(3*pie)); 
Q = (Area/2)*(Ybar); 
I = ((pie/64)*(d)^4); 
M = P2*Length; 
Bending = Ktbending*((32*M)/(pie*(d)^3)); 
sigaxial = ((P1/Area))*Ktaxial; 
torshear = ((16*T)/(pie*(d)^3))*Kttorsion; 
Transhear = ((P2*Q)/(I*d)); 
%Calculate Stresses 
sigxpoint1 = (sigaxial+Bending); 
torxypoint1 = (torshear); 
sigxpoint2 = (sigaxial); 
torxypoint2 = (Transhear+torshear); 
%Final Equations 
point1 = (sqrt(((sigxpoint1)^2)+3*(torxypoint1)^2))-eq == 0; 
point2 = (sqrt(((sigxpoint2)^2)+3*(torxypoint2)^2))-eq == 0; 
%Solve 
double(solve(point1,d)) 
double(solve(point2,d)) 
関連する問題