[wave modeling 2] three dimensional wave modeling and wave generator modeling matlab simulation

1. Software version

matlab2017b

2. Core source code

function func_power_gen_machine(HHup,R_buoy,seepart,seeall);


if seepart == 0 & seeall == 1
    % First generate the fixed part 
    % First of all ： Cylinder on the seafloor 
    t  = 0:pi/20:2*pi;
    RR =  100;
    HH = -6.5;
    x= 600+RR*sin(t);
    y= 300+RR*cos(t)/2;
    z=linspace(-8,HH,length(t));
    X=meshgrid(x);
    Y=meshgrid(y);
    Z=[meshgrid(z)]';
    surf(X,Y,Z);
    xlabel('x'),ylabel('y'),zlabel('z')
    shading interp;
    colormap([255/255,100/255,103/255]);
    hold on
    clear X Y Z


    % second ： The upper half of the column 
    t  = 0:pi/20:2*pi;
    RR2 = 60;
    HH2 = -6;
    x= 600+RR2*sin(t);
    y= 300+RR2*cos(t)/2;
    z=linspace(-8,HH2,length(t));
    X=meshgrid(x);
    Y=meshgrid(y);
    Z=[meshgrid(z)]';
    surf(X,Y,Z);
    xlabel('x'),ylabel('y'),zlabel('z')
    shading interp;
    colormap([255/255,100/255,103/255]);
    hold on
    clear X Y Z

    % Third ： coil 
    HH3 = -8:0.5:-3;
    for i = 1:floor(length(HH3)/2)
        t  = 0:pi/20:2*pi;
        RR3 = 10;

        x= 600+RR3*sin(t);
        y= 300+RR3*cos(t)/2;
        z=linspace(HH3(2*i-1),HH3(2*i),length(t));
        X=meshgrid(x);
        Y=meshgrid(y);
        Z=[meshgrid(z)]';
        plot3(X,Y,Z,'b');
        xlabel('x'),ylabel('y'),zlabel('z')
        shading interp;
        colormap([255/255,100/255,103/255]);
        hold on;
    end

    alpha(0.5);
    clear X Y Z


    HHup = HHup + 0.1;
    % Fourth ： Slide the piston 
    t  = 0:pi/20:2*pi;
    RR2  = 15;
    x= 600+RR2*sin(t);
    y= 300+RR2*cos(t)/2;

    z=linspace(-6+HHup,-5+HHup,length(t));
    X=meshgrid(x);
    Y=meshgrid(y);
    Z=[meshgrid(z)]';
    surf(X,Y,Z);
    xlabel('x'),ylabel('y'),zlabel('z')
    shading interp;
    colormap([255/255,100/255,103/255]);
    hold on
    alpha(0.5);
    clear X Y Z

    % The fifth ： Connect leads 
    t  = 0:pi/20:2*pi;
    RR2  = 0.01;
    x= 600+RR2*sin(t);
    y= 300+RR2*cos(t)/2;

    z=linspace(-5+HHup,0+HHup,length(t));
    X=meshgrid(x);
    Y=meshgrid(y);
    Z=[meshgrid(z)]';
    plot3(X,Y,Z,'k');
    xlabel('x'),ylabel('y'),zlabel('z')
    shading interp;
    colormap([255/255,100/255,103/255]);
    hold on
    alpha(0.5);
    clear X Y Z

    % The sixth ： Surface buoy , spherical 
    t=linspace(0,pi,25);
    p=linspace(0,2*pi,25);
    [theta,phi]=meshgrid(t,p);
    x= 600+R_buoy*sin(theta).*sin(phi);
    y= 300+R_buoy*sin(theta).*cos(phi)/2;
    z=R_buoy*cos(theta)/28 + HHup;
    surf(x,y,z);
    shading interp;
    colormap([255/255,100/255,103/255]);
    hold on
    alpha(0.5);

    alpha(0.5);
end 


if seepart == 1 & seeall == 0
    % The sixth ： Surface buoy , spherical 
    t=linspace(0,pi,25);
    p=linspace(0,2*pi,25);
    [theta,phi]=meshgrid(t,p);
    x= 600+R_buoy*sin(theta).*sin(phi);
    y= 300+R_buoy*sin(theta).*cos(phi)/2;
    z=R_buoy*cos(theta)/28 + HHup;
    surf(x,y,z);
    shading interp;
    colormap([255/255,100/255,103/255]);
    hold on
    alpha(0.5);

    alpha(0.5);
end 

3. test result

On the basis of the last achievement 【 Wave modeling 1】, Add the dynamic effect of wave generator , And simply calculate the power generation . According to the design requirements , After the whole generator , The simulation results are as follows ：

        You can see it here , The control interface has been adjusted , Three buttons are added , among ,part_View It is the dynamic effect of the sea buoy ,All_View It is the dynamic effect of the whole generator ,Close the power mach It is the dynamic effect of turning off the generator , Only show the effect of waves ,electricity It shows the effect of power generation . The simulation results of the whole system are as follows ：

The initial wave modeling ：

Click on Part_View Key ：

Adjust the view slider , The following effect ：

Click on All_View Key ：

Adjust the view slider , The following effect ：

Click on Close the power Mach Key ：

A19-07