## Thin plate deflection formula

Hi. I have been doing some FEA modelling with solid works and am trying to calculate my deflection for a point load at the centre of a clamped periphery (not simply supported) circular thin plate. I need to calculate the theoretical values to make sure that my FEA is correct

The formula I have found is this :

w=(-W/16pieD)*(a^2-r^2*(1+2*ln(a/r)))
for r not = to 0

w= deflection
W= load in N
a= fixed max radius (m)
r = variable radius (m)

D=flexual rigidity = Eh^3/(1-v^2)
E=Young's modulus (Pa)
h=plate thickness (m)
v=poissons ratio

When I plug my relevent data into the formula I get stuck because i am using a point load at the centre....therefore my r=0...I cannot find the formula for when the load is at the centre. Can anyone please help?

I have been having no trouble doing this with a distruited load (pressure) but it's the point load that I have been having trouble with.

Recognitions:
You will get a problem using concentrated analytic loadings with erroneously high answers. Roark says in Chapter 11: Flat Plates, Section 1: Common Case:

 Concentrated Loading It will be noted that all formulas for maximum stress due to a load applied over a small area give very high values when the radius of the loaded area approaches zero. Analysis by a more precise method (Ref 12) shows that the actual maximum stress produced by a load concentrated on a very small area of radius $$r_0$$ can be found by replacing the actual $$r_0$$ by a so-called equivalent radius $$r'_0$$, which depends largely upon the thickness of the plate $$t$$ and to a lesser degree on its least transverse dimension. Holl (Ref. 13) shows how $$r'_0$$ varies with the width of a flat plate. Westergaard (Ref. 14) gives an approximate expression for this equivalent radius: $$r'_0 = \sqrt{ 1.6 r^2_0 + t^2} - 0.675t$$ This formula which applies to a plate of any form, may be used for all values of $$r_0$$ less than $$0.5t$$; for larger values the actual $$r_0$$ may be used. Use of the equivalent radius makes possible the calculation of the finite maximum stresses produced by a (nominal) point loading whereas the ordinary forumula would indicate that these stresses were infinite.
So, the application of a concentrated loading physically is erroneous. You can try to apply the loading in your FEA as a concentrated surface loading over a finite area. Then use the formula given to get an equivalent radius, thereby which you can get the stresses and deflections.

For uniform loading over a very small central circular area of radius r0, those are Roark cases 16 and 17 depending on the boundary conditions. I can supply those if you would like. I think you have case 17 though (edges fixed rather than simply supported). In that case, the maximum deflection at r=0 is:

$$y_{max} = \frac{ -W a^2}{16 \pi D}$$

Where:

$$W = q \pi r^2_0$$

q being the "pressure", and a being the radius of the flat plate.
 thanks for that minger. That solves my problem. :)

Recognitions: