Code for a concentrated load on square plate bending?

I got this image from my abaqus analysis .The arrow in the image is showing the direction of applied concentrated load on a plate which is at the middle of 1/7th length of the plate near applied boundary condition. Can some one help me in how to mention this concentrated load in my matlab FE code?

The FE code is attached with this question.

[dpb]

Looks to be a very problem-specific solution with many (essentially all) BCs hardcoded as "magic numbers"; probably going to be somewhat difficult to generalize.

Only forces mentioned are initialized to zero; not clear, but I'm guessing that is solution, not input.

In the section

% converting uniform pressure into equivalent nodal forces. For this
% integration we are using one point Guass quadrature.
fe = Force(nnel,shape,P) ;             % Force vector

there's a computation of a force vector from (one presumes a) function Force but there is no definition for Force in the posted code.

This appears as is to be written only for a uniform pressure load, P; you could approximate the solution by instead of using a uniform pressure over the full plate, turn P into a 2D array and set it to zero excepting at the element containing the desired point and setting P=F_applied/A_element for that element, zero elsewhere.

Where would be based on the fixed nodes as are defined at the beginning as to there only being 20x20 nodes. Of course, it wouldn't actually be a point load, but the approximation of a uniform pressure of same total force spread over the area bounding the desired location.

To actually include a point load the formulation would have to be modified significantly.

[zainab malik]

yes Force is a function that i was using. This code was for uniform pressure but now i want to replace this part of code with a concentrated load.

[dpb]

Well, you'll have to rewrite the code entirely to do that it would seem; it doesn't have anything in it for point loads at all.

[zainab malik]

Can you please check this attached code if its fine for point load.

[dpb]

"Fine"? I've no idea.

It looks like it does have a point load added at a specific node by the comments, however, yes.