How can I simulate motion of buoy due ocean waves in 2D or 3D ?

hi, I'm new in matlab. can any one help me please? ocean waves is sine waves, I'm not understand how to simulate the motion of bouy due ocean waves in matlab, maybe any code or toolbox can represent about that.

[Andrew Newell]

Start with TUTORIAL: how to ask a question (on Answers) and get a fast answer.

Answer 1

[Chad Greene]

Run this:

x = 1:.1:100; % horizontal distance
lambda = 20;  % wavelength
c = 5;        % wave speed
A = 1.5;      % wave amplitude 

% Initialize plot: 
hold on
xlabel('distance (m)')
ylabel('wave height (m)') 
box off
axis([0 100 -5 5])

% Make an animation: 
for t= 0:.1:10

      % Calculate wave height along x at time t: 
      y = A*sin((2*pi/lambda)*(x-c*t)); 

      % Plot wave: 
      hwave = fill([x,100,0],[y,-5,-5],[.01 .44 .61],'edgecolor','b'); 

      % Get buoy height if buoy is at x = 75 m:  
      ybuoy = y(x==75)+.1; 
      hbuoy = plot(75,ybuoy,'ko','markersize',16,'markerfacecolor','r');

      % Draw a frame, wait 1/10 second, delete objects and repeat loop: 
      drawnow
      pause(.1)
      delete([hwave hbuoy])
  end

[Muhammad Luqman Amiruddin]

thank you for this your answer Mr. Chad Greene. can I ask again ? how to make line in the center of buoy to seabed ?

[Chad Greene]

Is this what you mean?

x = 1:.1:100; % horizontal distance
lambda = 20;  % wavelength
c = 5;        % wave speed
A = 1.5;      % wave amplitude 

% Initialize plot: 
hold on
xlabel('distance (m)')
ylabel('wave height (m)') 
box off
axis([0 100 -5 5])

% Make an animation: 
for t= 0:.1:10

        % Calculate wave height along x at time t: 
        y = A*sin((2*pi/lambda)*(x-c*t)); 

        % Plot wave: 
        hwave = fill([x,100,0],[y,-5,-5],[.01 .44 .61],'edgecolor','b'); 

        % Get buoy height if buoy is at x = 75 m:  
        ybuoy = y(x==75)+.1; 
        h_line = plot([75 75],[-5 ybuoy],'k-'); 
        hbuoy = plot(75,ybuoy,'ko','markersize',16,'markerfacecolor','r');

        % Draw a frame, wait 1/10 second, delete objects and repeat loop: 
        drawnow
        pause(.1)
        delete([hwave hbuoy h_line])
    end

[Muhammad Luqman Amiruddin]

yes. thank you very much Mr. Chad.

[Muhammad Luqman Amiruddin]

hello Mr. Chad Greene, can i ask you some question again? do you understand about simulink? especially SimHydraulic in Simscape. i have problem to make velocity hydraulic motor to constant. i will attach my system printscreen. thank you.

Answer 2

[Toby]

For anyone that stumbles on this question, there is now a MATLAB/Simulink simulation package for bodies in waves, distributed by the US Dept. of Energy:

http://wec-sim.github.io/WEC-Sim/

If you'd rather code up the basic equations yourself, you can do something like this:

https://www.researchgate.net/publication/272092752_Dynamics_of_a_Floating_Platform_Mounting_a_Hydrokinetic_Turbine

[rees adah]

Thank you for this resource