fluidpipelinesimulinkmodels.mdl

These Simulink blocks contain transfer functions that model the pressure and flow transients for axi

Model {
  Name			  "pipelines"
  Version		  2.20
  SimParamPage		  Solver
  SampleTimeColors	  off
  InvariantConstants	  off
  WideVectorLines	  off
  ShowLineWidths	  off
  StartTime		  "0.0"
  StopTime		  "10.0"
  Solver		  ode45
  RelTol		  "1e-3"
  AbsTol		  "1e-6"
  Refine		  "1"
  MaxStep		  "auto"
  InitialStep		  "auto"
  FixedStep		  "auto"
  MaxOrder		  5
  OutputOption		  RefineOutputTimes
  OutputTimes		  "[]"
  LoadExternalInput	  off
  ExternalInput		  "[t, u]"
  SaveTime		  on
  TimeSaveName		  "tout"
  SaveState		  off
  StateSaveName		  "xout"
  SaveOutput		  on
  OutputSaveName	  "yout"
  LoadInitialState	  off
  InitialState		  "xInitial"
  SaveFinalState	  off
  FinalStateName	  "xFinal"
  LimitMaxRows		  off
  MaxRows		  "1000"
  Decimation		  "1"
  AlgebraicLoopMsg	  warning
  MinStepSizeMsg	  warning
  UnconnectedInputMsg	  warning
  UnconnectedOutputMsg	  warning
  UnconnectedLineMsg	  warning
  ConsistencyChecking	  off
  ZeroCross		  on
  SimulationMode	  normal
  RTWSystemTargetFile	  "grt.tlc"
  RTWInlineParameters	  off
  RTWPlaceOutputsASAP	  off
  RTWRetainRTWFile	  off
  RTWTemplateMakefile	  "grt_default_tmf"
  RTWMakeCommand	  "make_rtw"
  RTWGenerateCodeOnly	  off
  ExtModeMexFile	  "ext_comm"
  ExtModeBatchMode	  off
  BlockDefaults {
    Orientation		    right
    ForegroundColor	    black
    BackgroundColor	    white
    DropShadow		    off
    NamePlacement	    normal
    FontName		    "Helvetica"
    FontSize		    10
    FontWeight		    normal
    FontAngle		    normal
    ShowName		    on
  }
  AnnotationDefaults {
    HorizontalAlignment	    center
    VerticalAlignment	    middle
    ForegroundColor	    black
    BackgroundColor	    white
    DropShadow		    off
    FontName		    "Helvetica"
    FontSize		    10
    FontWeight		    normal
    FontAngle		    normal
  }
  LineDefaults {
    FontName		    "Helvetica"
    FontSize		    9
    FontWeight		    normal
    FontAngle		    normal
  }
  System {
    Name		    "pipelines"
    Location		    [438, 330, 879, 564]
    Open		    on
    ToolBar		    on
    StatusBar		    on
    ScreenColor		    white
    PaperOrientation	    landscape
    PaperPositionMode	    auto
    PaperType		    usletter
    PaperUnits		    inches
    Block {
      BlockType		      SubSystem
      Name		      "P-pipeline"
      Ports		      [2, 3, 0, 0, 0]
      Position		      [195, 41, 230, 89]
      ShowPortLabels	      off
      MaskType		      "P-pipeline"
      MaskDescription	      "Linear dynamic model of a rigid uniform fluid t"
"ransmission line with inputs = pressures P0 and P1 at the ends,  outputs = vo"
"lumetric flow rate at the ends (inflow positive) and pressure mode shapes p. "
"The pressure at x is P0*(1-x/length) + P1*x/length + p*sin(pi*(0:8)'*x/length"
")"
      MaskHelp		      "The model is documented in J. Makinen, R. Piche"
" and A. Ellman (2000), Fluid Transmission Line Modeling Using a Variational M"
"ethod, ASME Journal of Dynamic Systems, Measurement, and Control,122 (1), 153"
"-162.\n"
      MaskPromptString	      "fluid density|kinematic viscosity|length|radius"
"|speed of sound|initial pressure|filtering"
      MaskStyleString	      "edit,edit,edit,edit,edit,edit,popup(None|Hann|H"
"amming|Riemann|Blackman-Harris)"
      MaskTunableValueString  "on,on,on,on,on,on,on"
      MaskVariables	      "p1=@1;p2=@2;p3=@3;p4=@4;p5=@5;p6=@6;p7=@7;"
      MaskInitialization      "T=p3/p5;ee=8*p2*T/p4^2;Z0=p1*p5/(pi*p4^2);jj=1:"
"8;e=sqrt(jj*pi*ee)/2+7*ee/16;w=jj*pi-sqrt(jj*pi*ee)/4+ee/16;Pinit=pi*p4^2*p6/"
"p1/p5; if p7==1, k=ones(1,8);elseif p7==2, k=cos(jj*pi/9/2).^2;elseif p7==3, "
"k=0.54+.46*cos(jj*pi/9);elseif p7==4, k= 9 ./(jj*pi).*sin(jj*pi/9);elseif p7="
"=5, k= 0.42323+0.49755*cos(jj*pi/9)+0.07922*cos(2*jj*pi/9);else,  k=ones(1,8)"
";end; b=1 ./(8*sum(k(1:2:7) ./w(1:2:7).^2));"
      MaskDisplay	      "plot(-4,-10,24,16,[0 0 20 20],[0 2 2 0],[0 0 20"
" 20],[6 4 4 6])"
      MaskIconFrame	      on
      MaskIconOpaque	      on
      MaskIconRotate	      none
      MaskIconUnits	      autoscale
      MaskValueString	      "870|8e-5|19.74|6.17e-3|1400|0|Blackman-Harris"
      System {
	Name			"P-pipeline"
	Location		[5, 81, 931, 731]
	Open			off
	ToolBar			off
	StatusBar		off
	ScreenColor		white
	PaperOrientation	landscape
	PaperPositionMode	auto
	PaperType		usletter
	PaperUnits		inches
	Block {
	  BlockType		  Inport
	  Name			  "P0"
	  Position		  [15, 20, 35, 40]
	  Port			  "1"
	  PortWidth		  "-1"
	  SampleTime		  "-1"
	}
	Block {
	  BlockType		  Inport
	  Name			  "P1"
	  Position		  [15, 365, 35, 385]
	  Port			  "2"
	  PortWidth		  "-1"
	  SampleTime		  "-1"
	}
	Block {
	  BlockType		  Constant
	  Name			  "Constant1"
	  Position		  [15, 59, 55, 81]
	  ShowName		  off
	  Value			  "Pinit"
	}
	Block {
	  BlockType		  Constant
	  Name			  "Constant2"
	  Position		  [15, 334, 55, 356]
	  ShowName		  off
	  Value			  "Pinit"
	}
	Block {
	  BlockType		  Gain
	  Name			  "Gain"
	  Position		  [150, 352, 175, 378]
	  ShowName		  off
	  Gain			  "1/2"
	}
	Block {
	  BlockType		  Gain
	  Name			  "Gain1"
	  Position		  [145, 37, 170, 63]
	  ShowName		  off
	  Gain			  "1/2"
	}
	Block {
	  BlockType		  Gain
	  Name			  "Gain3"
	  Position		  [635, 156, 670, 184]
	  ShowName		  off
	  Gain			  "1/Z0"
	}
	Block {
	  BlockType		  Gain
	  Name			  "Gain4"
	  Position		  [635, 438, 670, 472]
	  ShowName		  off
	  Gain			  "1/Z0"
	}
	Block {
	  BlockType		  Mux
	  Name			  "Mux"
	  Ports			  [9, 1, 0, 0, 0]
	  Position		  [575, 594, 605, 726]
	  Inputs		  "9"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [115, 349, 130, 376]
	  ShowName		  off
	  Inputs		  "++"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum1"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [115, 33, 130, 62]
	  ShowName		  off
	  Inputs		  "+-"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum2"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [720, 162, 745, 193]
	  ShowName		  off
	  Inputs		  "++"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum3"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [725, 433, 750, 462]
	  ShowName		  off
	  Inputs		  "-+"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum4"
	  Ports			  [5, 1, 0, 0, 0]
	  Position		  [590, 13, 605, 327]
	  ShowName		  off
	  Inputs		  "+++++"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum5"
	  Ports			  [4, 1, 0, 0, 0]
	  Position		  [590, 330, 605, 580]
	  ShowName		  off
	  Inputs		  "++++"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum6"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [75, 354, 90, 381]
	  ShowName		  off
	  Inputs		  "-+"
	}
	Block {
	  BlockType		  Sum
	  Name			  "Sum7"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [80, 24, 95, 51]
	  ShowName		  off
	  Inputs		  "+-"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn"
	  Position		  [215, 338, 420, 392]
	  ShowName		  off
	  Numerator		  "[4*k(1)*T 0]"
	  Denominator		  "[T^2 e(1)*T w(1)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn1"
	  Position		  [215, 397, 420, 453]
	  ShowName		  off
	  Numerator		  "[4*k(3)*T 0]"
	  Denominator		  "[T^2 e(3)*T w(3)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn2"
	  Position		  [215, 81, 420, 139]
	  ShowName		  off
	  Numerator		  "[4*k(2)*T 0]"
	  Denominator		  "[T^2 e(2)*T w(2)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn3"
	  Position		  [215, 143, 420, 197]
	  ShowName		  off
	  Numerator		  "[4*k(4)*T 0]"
	  Denominator		  "[T^2 e(4)*T w(4)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn4"
	  Position		  [215, 201, 420, 259]
	  ShowName		  off
	  Numerator		  "[4*k(6)*T 0]"
	  Denominator		  "[T^2 e(6)*T w(6)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn5"
	  Position		  [215, 262, 420, 318]
	  ShowName		  off
	  Numerator		  "[4*k(8)*T 0]"
	  Denominator		  "[T^2 e(8)*T w(8)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn6"
	  Position		  [215, 457, 420, 513]
	  ShowName		  off
	  Numerator		  "[4*k(5)*T 0]"
	  Denominator		  "[T^2 e(5)*T w(5)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn7"
	  Position		  [215, 518, 420, 572]
	  ShowName		  off
	  Numerator		  "[4*k(7)*T 0]"
	  Denominator		  "[T^2 e(7)*T w(7)^2 ]"
	}
	Block {
	  BlockType		  TransferFcn
	  Name			  "Transfer Fcn8"
	  Position		  [285, 32, 370, 68]
	  ShowName		  off
	  Numerator		  "[2]"
	  Denominator		  "[T ee ]"
	}
	Block {
	  BlockType		  Outport
	  Name			  "Q0"
	  Position		  [785, 164, 825, 196]
	  Port			  "1"
	  OutputWhenDisabled	  held
	  InitialOutput		  "0"
	}
	Block {
	  BlockType		  Outport
	  Name			  "Q1"
	  Position		  [790, 435, 835, 465]
	  Port			  "2"
	  OutputWhenDisabled	  held
	  InitialOutput		  "0"
	}
	Block {
	  BlockType		  Outport
	  Name			  "p"
	  Position		  [725, 645, 770, 675]
	  Port			  "3"
	  OutputWhenDisabled	  held
	  InitialOutput		  "0"
	}
	Line {
	  SrcBlock		  "Transfer Fcn8"
	  SrcPort		  1
	  Points		  [130, 0]
	  Branch {
	    Points		    [0, 550]
	    DstBlock		    "Mux"
	    DstPort		    1
	  }
	  Branch {
	    DstBlock		    "Sum4"
	    DstPort		    1
	  }
	}
	Line {
	  SrcBlock		  "Sum5"
	  SrcPort		  1
	  DstBlock		  "Gain4"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Sum4"
	  SrcPort		  1
	  DstBlock		  "Gain3"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Sum"
	  SrcPort		  1
	  DstBlock		  "Gain"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Gain3"
	  SrcPort		  1
	  Points		  [20, 0]
	  Branch {
	    DstBlock		    "Sum2"
	    DstPort		    1
	  }
	  Branch {
	    Points		    [0, 270]
	    DstBlock		    "Sum3"
	    DstPort		    1
	  }
	}
	Line {
	  SrcBlock		  "Gain4"
	  SrcPort		  1
	  Points		  [0, 0]
	  Branch {
	    DstBlock		    "Sum3"
	    DstPort		    2
	  }
	  Branch {
	    Points		    [0, -270]
	    DstBlock		    "Sum2"
	    DstPort		    2
	  }
	}
	Line {
	  SrcBlock		  "Gain1"
	  SrcPort		  1
	  Points		  [0, 0]
	  Branch {
	    Points		    [5, 0; 0, 60]
	    Branch {
	      Points		      [0, 60]
	      Branch {
		DstBlock		"Transfer Fcn3"
		DstPort			1
	      }
	      Branch {
		Points			[0, 60]
		Branch {
		  DstBlock		  "Transfer Fcn4"
		  DstPort		  1
		}
		Branch {
		  Points		  [0, 60]
		  DstBlock		  "Transfer Fcn5"
		  DstPort		  1
		}
	      }
	    }
	    Branch {
	      DstBlock		      "Transfer Fcn2"
	      DstPort		      1
	    }
	  }
	  Branch {
	    DstBlock		    "Transfer Fcn8"
	    DstPort		    1
	  }
	}
	Line {
	  SrcBlock		  "Gain"
	  SrcPort		  1
	  Points		  [0, 0]
	  Branch {
	    Points		    [2, 60]
	    Branch {
	      Points		      [2, 60]
	      Branch {
		DstBlock		"Transfer Fcn6"
		DstPort			1
	      }
	      Branch {
		Points			[1, 60]
		DstBlock		"Transfer Fcn7"
		DstPort			1
	      }
	    }
	    Branch {
	      DstBlock		      "Transfer Fcn1"
	      DstPort		      1
	    }
	  }
	  Branch {
	    DstBlock		    "Transfer Fcn"
	    DstPort		    1
	  }
	}
	Line {
	  SrcBlock		  "Sum3"
	  SrcPort		  1
	  DstBlock		  "Q1"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Sum2"
	  SrcPort		  1
	  DstBlock		  "Q0"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Sum7"
	  SrcPort		  1
	  Points		  [0, 0]
	  Branch {
	    DstBlock		    "Sum1"
	    DstPort		    1
	  }
	  Branch {
	    DstBlock		    "Sum"
	    DstPort		    1
	  }
	}
	Line {
	  SrcBlock		  "P0"
	  SrcPort		  1
	  DstBlock		  "Sum7"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Constant1"
	  SrcPort		  1
	  DstBlock		  "Sum7"
	  DstPort		  2
	}
	Line {
	  SrcBlock		  "Constant2"
	  SrcPort		  1
	  DstBlock		  "Sum6"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "P1"
	  SrcPort		  1
	  DstBlock		  "Sum6"
	  DstPort		  2
	}
	Line {
	  SrcBlock		  "Sum6"
	  SrcPort		  1
	  Points		  [0, 0]
	  Branch {
	    DstBlock		    "Sum"
	    DstPort		    2
	  }
	  Branch {
	    Points		    [0, -315]
	    DstBlock		    "Sum1"
	    DstPort		    2
	  }
	}
	Line {
	  SrcBlock		  "Mux"
	  SrcPort		  1
	  DstBlock		  "p"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Transfer Fcn5"
	  SrcPort		  1
	  Points		  [0, 0]
	  Branch {
	    DstBlock		    "Sum4"
	    DstPort		    5
	  }
	  Branch {
	    Points		    [0, 430]
	    DstBlock		    "Mux"
	    DstPort		    9
	  }
	}
	Line {
	  SrcBlock		  "Transfer Fcn3"
	  SrcPort		  1
	  Points		  [40, 0]
	  Branch {
	    DstBlock		    "Sum4"
	    DstPort		    3
	  }
	  Branch {
	    Points		    [0, 490]
	    DstBlock		    "Mux"
	    DstPort		    5
	  }
	}
	Line {
	  SrcBlock		  "Transfer Fcn2"
	  SrcPort		  1
	  Points		  [60, 0]
	  Branch {
	    DstBlock		    "Sum4"
	    DstPort		    2
	  }
	  Branch {
	    Points		    [0, 520]
	    DstBlock		    "Mux"
	    DstPort		    3
	  }
	}
	Line {
	  SrcBlock		  "Transfer Fcn4"
	  SrcPort		  1
	  Points		  [20, 0]
	  Branch {
	    DstBlock		    "Sum4"
	    DstPort		    4
	  }
	  Branch {
	    Points		    [0, 460]
	    DstBlock		    "Mux"
	    DstPort		    7
	  }
	}
	Line {