Model Drivetrain Noise

This example shows how to inject a fault into a drivetrain using aTorque Noise Source块。将噪声注入你的模型允许你to predict how your actual physical system responds when it experiences environmental or internal disturbances. It also allows you to test the robustness and responsiveness of your control system.

In the example, the noise is injected into one of the sprockets on a metal sheet feeder that is moving a piece of stock. Although you can perform most of the steps in this example using tools that the Simulink^®andSimscape™ Driveline™user interfaces provide, programmatic commands are supplied. You can combine the programmatic commands to create a script for parameter sweeps.

Open the model. At the MATLAB^®command prompt, enter:
```
model ='sdl_sheet_metal_feeder'; open_system(model)
```

Simulate the model and plot the results.

Script for Generating and Plotting Simulation Results

Add, configure, connect, and arrange these blocks as shown:

Stepblock — Specify7.5for theStep timeparameter and300for theFinal valueparameter.
Simulink-PS Converterblock
Mechanical Rotational Referenceblock
Torque Noise Sourceblock — For theSample time, specify1e-1, forRepeatability, selectSpecify seed, and forSeedspecify0.

Script for Adding, Configuring, Connecting, and Arranging Blocks

%% Expand the model window to make room for new blocksset_param(model,'Location',[108 73 881 682])%% Add and configure Step block% Define blockstepPath = [model,'/Step']; stepLib ='simulink/Sources/Step'; stepPosition = [-195 205 -165 235]; stepTime ='7.5'; stepValue ='300';% Add blockadd_block(stepLib,stepPath,'Position',stepPosition)% Configure Step blockset_param(stepPath,'Time',stepTime,...'After',stepValue)% Check block configurationopen_system(stepPath) pause(3); close_system(stepPath)% Get output portstepPortHandle = get_param(stepPath,'PortHandles');stepOutport = stepPortHandle.Outport;%% Add and configure S-PS Converter block% Define blocksPSConvPath = [model,'/Simulink-PS Converter']; sPSConvLib =...'nesl_utility/Simulink-PS Converter'; sPSConvPathPosition = [-115 205 -85 235];% Add blockadd_block(sPSConvLib,sPSConvPath,'Position',sPSConvPathPosition)% Get output portsPSConvPortHandle = get_param(sPSConvPath,'PortHandles');sPSConvInport = sPSConvPortHandle.Inport; sPSConvConPort = sPSConvPortHandle.RConn(1,1);%% Add Mechanical Rotational Reference block% Define blockmechRotRefPath = [model,'/Mechanical Rotational Reference']; mechRotRefLib =...'fl_lib/Mechanical/Rotational Elements/Mechanical Rotational Reference'; mechRotRefPosition = [-40 255 -20 275];% Add blockadd_block(mechRotRefLib,mechRotRefPath,'Position',mechRotRefPosition)%% Get output portmechRotRefPortHandle = get_param(mechRotRefPath,'PortHandles');mechRotRefConPort = mechRotRefPortHandle.LConn(1,1);%% Add and configure Torque Noise Source block% Define blocktorqueNoisePath = [model,'/Torque Noise Source']; torqueNoiseLib ='sdl_lib/Sources/Torque Noise Source'; torqueNoisePosition = [25 210 65 250];%% Add blockadd_block(torqueNoiseLib,torqueNoisePath,'Position',torqueNoisePosition)%% Get output porttorqueNoisePortHandle = get_param(torqueNoisePath,'PortHandles');torqueNoiseConPort1 = torqueNoisePortHandle.LConn(1,1); torqueNoiseConPort2 = torqueNoisePortHandle.LConn(1,2); torqueNoiseConPort3 = torqueNoisePortHandle.RConn(1,1);%% Configure blockset_param(torqueNoisePath,'sample_time','1e-1',...'repeatability','3')%% Check block configurationopen_system(torqueNoisePath) pause(3); close_system(torqueNoisePath)%% Get block points for connecting as a branched linetorqueNoisePortPoints = get_param(torqueNoisePath,'PortConnectivity');[torqueNoiseLConnPoints1,torqueNoiseLConnPoints2,torqueNoiseRConnPoints]...= torqueNoisePortPoints.Position; sprocketPath = [model,'/Sprocket to Belt R']; sprocketPortHandle = get_param(sprocketPath,'PortHandles');sprocketConPort1 = sprocketPortHandle.RConn(1,1); sproketPortPoints = get_param(sprocketPath,'PortConnectivity');[sproketLConnPoints,sproketRConnPoints] = sproketPortPoints.Position;%% Connect blocksadd_line(model,stepOutport,sPSConvInport) add_line(model,sPSConvConPort,torqueNoiseConPort1) add_line(model,mechRotRefConPort,torqueNoiseConPort2) add_line(model,[sproketLConnPoints;torqueNoiseRConnPoints])%% Hold to check connectionspause(5);

Simulate the model and plot the results.

Script for Generating and Plotting Simulation Results

The noise injected at simulation time,t= 7.5 seconds introduces torque disturbances.

Zoom in to see the effects of the disturbance.
Script for Zooming In
%% Zoom to examine data at peak noisexlim(axes3,[12.0 14.50]); ylim(axes3,[-120 120]);
When the noise source torque exceeds ±50 N*m, it most significantly effects the torque applied by the sprocket and, therefore, the position of the stock.

Simscape Driveline Documentation

金宝app

MATLAB、Simulink、その他の製品をお試しください

無料評価版のダウンロード