Excitation force floating body

Discussion in 'Hydrodynamics and Aerodynamics' started by Troglo, Nov 30, 2012.

  1. Troglo
    Joined: Nov 2012
    Posts: 4
    Likes: 0, Points: 0, Legacy Rep: 10
    Location: New York

    Troglo New Member

    I am working on a heave model , I am trying to make a model in heave in the time domain. I am using wamit the identify the hydrodinamic parameters as added mass, radiation damping and excitation force. I have some problem with the result, indeed using the diffraction force and the haskind force coming from wamit I have a response of my model not reasonable, I have a amplitude of the motion 2 or 3 times the incident wave. I know that the parameters of the force are normalized by rho*g*H/2, but made a total excitation force due to the sum of Haskind+diffraction I could find any reasonable results. Can you explain how use this values?

    I write the equation that I am solving through the matlab function ODE45, where bvz is the viscous factor damping, brz is the radiation damping, Awp is the restoring coefficient, m and awz are mass and add mass and Fz0 is the total excitation force
    y(2) is the velocity in heave of my body and y(1) is the displacement.

    dydt=[y(2);
    (-(bvz*y(2)*abs(y(2))))-brz*y(2)-(((Awp))*y(1))+Fz0*cos((ome.*t))))/(m+awz)];

    I would be very grateful if someone could help me to understand where is the error!
     
  2. philSweet
    Joined: May 2008
    Posts: 2,682
    Likes: 451, Points: 83, Legacy Rep: 1082
    Location: Beaufort, SC and H'ville, NC

    philSweet Senior Member

    your parentheses are a mess. (check (your)) ttyping)));)
     
  3. Troglo
    Joined: Nov 2012
    Posts: 4
    Likes: 0, Points: 0, Legacy Rep: 10
    Location: New York

    Troglo New Member

    You are right, sorry!
    I fix it a bit

    dydt=[y(2);
    (-(bvz*(y(2)*|(y(2)|-brz*y(2)-Awp*y(1)+Fz0*cos(ome.*t))/(m+awz)];
     
  4. philSweet
    Joined: May 2008
    Posts: 2,682
    Likes: 451, Points: 83, Legacy Rep: 1082
    Location: Beaufort, SC and H'ville, NC

    philSweet Senior Member

    next item-

    The quantity on the left is an acceleration - y" or d^2y/dt^2. Are you integrating it twice to get heave? If I understand the code, it will return an acceleration, but I've never used Matlab.

    Is Fz0 a surface displacement, or a force? if the former, wouldn't it just get rolled into the y(1) term? ie Awp*(y(1)+Fz0*cos(ome.*t))

    when all else fails, write the whole thing out with units explicitly and make sure the units add up.
     
  5. Troglo
    Joined: Nov 2012
    Posts: 4
    Likes: 0, Points: 0, Legacy Rep: 10
    Location: New York

    Troglo New Member

    ok, I was not really clear, my fault. Fz0 is the amplitude of the excitation force.
    My question is: Should I consider the total amplitutde of the excitation force acting on my body as the sum of the two output coming from wamit, ie. HASKIND and DIFFRACTION excitation force, or these two parameters are two way to identify the same quantity?

    Thanks and sorry for the lack of clarity!
     
  6. CWTeebs
    Joined: Apr 2011
    Posts: 232
    Likes: 15, Points: 0, Legacy Rep: 171
    Location: Maine

    CWTeebs AnomalyGenerator

    Are you simulating in regular or irregular seas? Can you post your MatLab code? What do your variables mean? What coordinate system are you using?

    It's the "Haskind relations" not the "Haskind forces," meaning the pressure in the fluid domain at the location of the structure without the ship, plus the force required to scatter/diffract the water away from the structure. Together these constitute the excitation forces. Well, I guess you could call these the Haskind forces, I'm just used to seeing it referred to as the Haskind relations, they're just two distinct components that makeup the inputs to the system.

    Your problem may be related to how your undimensionalizing the WAMIT parameters. I've not personally used WAMIT. See the following thread:

    http://www.boatdesign.net/forums/de...tation-force-nondimensionalization-45057.html
     

  7. Troglo
    Joined: Nov 2012
    Posts: 4
    Likes: 0, Points: 0, Legacy Rep: 10
    Location: New York

    Troglo New Member

    I am simulating in regular waves, the coordinate system is refered to the center of gravity of the body.
    bvz=coeff. of viscous damping
    brz=coeff. of radiation damping
    Awp=waterplane area
    Fz0=Excitation force amplitude
    m=bodu mass
    awz=added mass;

    Mu question is, the output of wamit is:
    - HASKIND EXCITING FORCES AND MOMENTS;
    - DIFFRACTION EXCITING FORCES AND MOMENTS;

    In my model I use for Fz0, one of the two above forces,is it right?
    My personal e-mail is: t.rog.lo@hotmail.com, if you write to me I can send the code.
     
Loading...
Similar Threads
  1. JoshTruman
    Replies:
    8
    Views:
    1,119
  2. Ricky Larsen
    Replies:
    10
    Views:
    1,440
  3. captncoop
    Replies:
    5
    Views:
    1,454
  4. sun
    Replies:
    0
    Views:
    1,050
  5. Sailor Al
    Replies:
    40
    Views:
    6,104
  6. Sailor Al
    Replies:
    80
    Views:
    8,850
  7. Ittiandro
    Replies:
    2
    Views:
    6,749
  8. Toru
    Replies:
    7
    Views:
    2,536
  9. BlunderBus
    Replies:
    11
    Views:
    5,607
  10. rspawarr
    Replies:
    12
    Views:
    7,250
Forum posts represent the experience, opinion, and view of individual users. Boat Design Net does not necessarily endorse nor share the view of each individual post.
When making potentially dangerous or financial decisions, always employ and consult appropriate professionals. Your circumstances or experience may be different.