Advanced Search
All Posts : Product Tips & Techniques

Implicit equation plotting

cabs360 12 Maple
May 01 2006
0
Can anyone help me please?! Cash Reward! My program solves my implicit equation but i am having trouble getting more than one result, ie, varying one of the variables. The best i can do is make a loop which gives a table of results which requires individually copying the results across into Excel for anysort of graph. Ideally i want to be able to plot the results in Maple, for any variable. This is for my 3rd year engineering project, any help much appreciated. > restart; > with(LinearAlgebra): > # IMPACT OF A FLUID FILLED SPHERICAL SHELL ("NON-HERTZIAN") WITH A "HERTZIAN" BALL- > # All SI units metric > # Input Parameters fluid filled spherical shell > > Rsh:=0.0725: # radius brain > h:=0.00561: # thickness skull > Esh:=1.6e9: # Young's modulus skull > mush:=0.3: # Poisson's ratio skull > vsh:=0: # velocity of head > msh:=1.96: # mass of head > ksh:=evalf(2.3*Esh*h^2/(Rsh*(1-mush^2)^(1/2))): # ovalstiff > # Input Parameters for hertzian ball > mub:=0.25: # Poisson's ratio ball > Eb:=207e9: # Young's modulus ball > Rb:=0.009: # Radius curvature ball > vb:=1.98: # Velocity of ball > # Variable > for i from 0 to 2 do > mb:=0.045*i; # Mass ball (make it infinite for ground) > # Combined derived parameters > m:=1/(1/msh+1/mb); # Effective mass m* > R:=1/(1/(Rsh+h)+1/Rb); # Effective curvature R* > Ep:=Esh*Eb/((1-mush^2)*Esh+(1-mub^2)*Eb):# Effective Modulus E* > v:=abs(vsh-vb); # Delta V > kb:=4/3*R^(1/2)*Ep: #non-linear hertzian contact spring > # > alpha:=4/5*ksh/kb^(2/3): > beta:=ksh*m*v^2: > eq:=PP+alpha*PP^(5/6)-beta: > PoP:=fsolve(eq=0,PP): > > # General Solution > P[i]:=PoP^(1/2); # Peak Load > > # > od:; > P(i) > ; TABLE([0 = 0., 1 = 437.1559281, 2 = 623.2774814])(i)
engineering plot graph linear

Please Wait...