The application of polynomial reproducing methods has been explored in the context of linear and non linear problems. Of specific interest is the application of a recently developed reproducing scheme, referred to as the error reproducing kernel method (ERKM), which uses non-uniform rational B-splines (NURBS) to construct the basis functions, an aspect that potentially helps bring in locall support, convex approximation and variation diminishing properties in the functional approximation. Polynomial reproducing methods have been applied to solve problems coming under the class of a simplified theory called Cosserat theory. Structures such as a rod which have special geometric properties can be modeled with the aid of such simplified theories. It has been observed that the application of mesh-free methods to solve the aforementioned problems has the advantage that large deformations and exact cross-sectional deformations in a rod could be captured exactly by modeling the rod just in one dimension without the problem of distortion of elements or element locking which would have had some effect if the problem were to be solved using mesh based methods. Polynomial reproducing methods have been applied to problems in fracture mechanics to study the propagation of crack in a structure. As it is often desirable to limit the use of the polynomial reproducing methods to some parts of the domain where their unique advantages such as fast convergence, good accuracy, smooth derivatives, and trivial adaptivity are beneficial, a coupling procedure has been adopted with the objective of using the advantages of both FEM and polynomial reproducing methods. Exploration of SMW (Sherman-Morrison-Woodbury) in the context of polynomial reproducing methods has been done which would assist in calculating the inverse of a perturbed matrix (stiffness matrix in our case). This would to a great extent reduce the cost of computation. In this thesis, as a first step attempts have been made to apply Mesh free cosserat theory to one dimensional problems. The idea was to bring out the advantages and limitations of mesh free cosserat theory and then extend it to 2D problems.