A new nonlinear fracture analysis framework is developed for the mode-I and II fracture response of thick-section fiber reinforced polymeric (FRP) composites. This framework employs 3D micromechanical constitutive models for the nonlinear material behavior along with cohesive elements for crack growth. Fracture tests on various cracked geometries are used to verify the prediction of the failure loads and the crack growth behavior. A commercially available pultruded E-glass/polyester and vinylester thick-section FRP composite material was used to demonstrate the proposed fracture approach along with the nonlinear constitutive modeling. A new Infra-red thermography technique is derived to measure the surface strain field near the crack tip in the linear response range. Mode I and II fracture toughness tests for pultruded composites are also examined using the eccentrically loaded, single-edge-notch tension, ESE(T), single-edge-notch tension, SEN(T), and a butterfly specimen with an Arcan-type fixture. Material nonlinearity and crack growth effects were observed during the tests and investigated using the proposed analysis framework. The effect of material orthotropy on the stress intensity factor solutions was addressed using the virtual crack closure technique. The analytic and experimental results support the use of the ESE(T) specimen for the measuring the mode-I fracture toughness and the butterfly shaped specimen for measuring the mode-II toughness. The calibrated cohesive models were able to predict the measured crack growth in both modes I and II for various crack geometries. A mixed mode failure criterion is proposed and verified with test results. Examples are presented for using this criterion and crack growth analyses. The experimental and analytical results of this study can form a foundation for using fracture-based methods for the design of structures using these materials.

Advances in the design and manufacturing technologies have greatly enhanced the utility of fiber reinforced composite materials in aircraft, helicopter and space- craft structural components. The special characteristics of composites such as high strength and stiffness, light-weight corrosion resistance make them suitable sub- stitute for metals/metallic alloys. However, composites are very sensitive to the anomalies induced during their fabrication and service life. Also, they are suscepti- ble to the impact and high frequency loading conditions because the epoxy matrix is at-least an order of magnitude weaker than the embedded reinforced carbon fibers. On the other hand, the carbon based matrix posses high electrical conductivity which is often undesirable. Subsequently, the metal matrix produces high brittleness. Var- ious forms of damage in composite laminates can be identified as indentation, fiber breakage, matrix cracking, fiber-matrix debonding and interply disbonding (delam- ination). Among all the damage modes mentioned above, delamination has been found to be serious for all cases of loading. They are caused by excessive interlaminar shear and normal stresses. The interlaminar stresses that arise in the case of composite materials due to the mismatch in the elastic constants across the plies. Delamination in composites reduce it’s tensile and compressive strengths by consid- erable margins. Hence the knowledge of these stresses is the most important aspect to be looked into. Basic theories like the Euler-Bernoulli’s theory and Timoshenko beam theory are based on many assumptions which poses limitation to determine these stresses accurately. Hence the determination of these interlaminar stresses accurately requires higher order theories to be considered. Most of the conventional methods of determination of the stresses are through the solutions, involving the trigonometric series, which are available only to small and simple problems. The most common method of solution is by Finite Element (FE) Method. There are only few elements existing in the literature and very few in the commercially available finite element software to determine the interlaminar stresses accurately in the composite laminates. Accuracy of finite element solution depends on the choice of functions to be used as interpolating polynomials for the field variable. In-appropriate choice will manifest in the form of delayed convergence. This delayed convergence and accuracy in predicting these stresses necessiates a formulation of elements with a completely new concept. The delayed convergence is sometimes attributed to the shear locking phenomena, which exist in most finite element formulation based on shear deformation theories. The present work aims in developing finite elements based on higher order theories, that alleviates the slow convergence and achieves the solutions at a faster rate without compromising on the accuracy. The accuracy primarily depends on the theory used to model the problem. Thus the basic theories (such as Elementary Beam theory and Timoshenko Beam theory) does not suffice the condition to accuratley determine the interlaminar stresses through the thickness, which is the primary cause for delamination in composites. Two different elements developed on the principle of super-convergence has been presented in this work. These elements are subjected to several numerical experiments and their performance is assessed by comparing the solutions with those available in literature. Spacecraft and aircraft structures are light in weight and are also lightly damped because of low internal damping of the material of construction. This increased exibility may allow large amplitude vibration, which might cause structural instability. In addition, they are susceptible to impact loads of very short duration, which excites many structural modes. Hence, structural dynamics and wave propagation study becomes a necessity. The wave based techniques have found appreciation in many real world problems such as in Structural Health Monitoring (SHM). Wave propagation problems are characterized by high frequency loads, that sets up stress waves to propagate through the medium. At high frequency, the wave lengths are small and from the finite element point of view, the element sizes should be of the same order as the wave lengths to prevent free edges of the element to act as a free boundary and start reflecting the stress waves. Also longer element size makes the mass distribution approximate. Hence for wave propagation problems, very large finite element mesh is an absolute necessity. However, the finite element problems size can be drastically reduced if we characterize the stiffness of the structure accurately. This can accelerate the convergence of the dynamic solution significantly. This can be acheived by the super-convergent formulation. Numerical results are presented to illustrate the efficiency of the new approach in both the cases of dynamic studies viz., the free vibration study and the wave propagation study. The thesis is organised into five chapters. A brief organization of the thesis is presented below, Chapter-1 gives the introduction on composite material and its constitutive law. The details of shear locking phenomena and the interlaminar stress distribution across the thickness is brought out and the present methods to avoid shear locking has been presented. Chapter-2 presents the different displacement based higher order shear deformation theories existing in the literature their advantages and limitations. Chapter-3 presents the formulation of a super-convergent finite element formulation, where the effect of lateral contraction is neglected. For this element static and free vibration studies are performed and the results are validated with the solution available in the open literature. Chapter-4 presents yet another super-convergent finite element formulation, wherein the higher order effects due to lateral contraction is included in the model. In addition to static and free vibration studies, wave propagation problems are solved to demonstrate its effectiveness. In all numerical examples, the super-convergent property is emphasized. Chapter-5 gives a brief summary of the total research work performed and presents further scope of research based on the current research.

Carbon nanotubes, graphene and nano sized core shell rubber particles have all been extensively researched for their capability to improve mechanical properties of thermoset resins. However, there has been a lack of research on their evaluation for energy absorption in high velocity impact scenarios, and the fundamental mechanics of their failure mechanisms during highly dynamic stress transfer through the matrix. This fundamental research is essential for laying the foundation for improvement in ballistic performance in composite armor. In hard armor applications, energy absorption is largely accomplished through delamination between plies of the composite laminate. This energy absorption is accomplished through two mechanisms. The first being the elongation of the fiber reinforcement contained in the resin matrix, and the second is the propagation of the crack in between the discreet fabric plies. This research aims to fundamentally study the energy absorption characteristics of various nano-particles as reinforcements in thermoset resin for high velocity impact applications. Multiple morphologies will be evaluated through use of platelet, tubular and spherical shaped nano-particles. Evaluations of the effect on stress transfer through the matrix due to the combination of nano sized and micro scale particles of milled fiber is conducted. Three different nano-particles are utilized, specifically, multi-walled carbon nanotubes, graphene, and core shell rubber particles. The difference in surface area, aspect ratio and molecular structure between the tube, platelet and spherical nano-particles causes energy absorption through different failure mechanisms. This changes the impact performance of composite panels enhanced with the nanoparticle fillers. Composite panels made through the use of dispersing the various nano-particles iv in a non-contact planetary mixer, are evaluated through various dynamic and static testing, including unnotched cantilever beam impact, mixed mode fracture toughness, split-Hopkinson bar, and ballistic V50 testing. The unnotched cantilever beam testing showed that the addition of milled fiber degraded the impact resistance of the samples. Addition of graphene nano platelets unilaterally degraded impact resistance through the unnotched cantilever beam testing. 1.5% loading of MWCNT showed the greatest increase in impact resistance, with a 43% increase over baseline. Determining the critical load for mixed mode interlaminar shear testing can be difficult for composite panels that bend without breaking. An iterative technique of optimizing the coefficient of determination, R2 , in linear regression is developed for objectively determining the point of non-linearity for critical load. This allows for a mathematical method of determination; thereby eliminating any subjective decision of choosing where the data becomes non-linear. The core shell rubber nano particles showed the greatest strain energy release rate with an exponential improvement over the baseline results. Synergistic effects between nano and micro sized particles in the resin matrix during transfer of the stress wave were created and evaluated. Loadings of 1% milled carbon fiber enhanced the V50 ballistic performance of both carbon nanotube and core shell rubber particles in the resin matrix. However, the addition of milled carbon fiber degrades the impact resistance of all nano-particle enhanced resin matrices. Therefore, benefits gained from the addition of microsized particles in combination with nano-sized particles, are only seen in high energy impact scenarios with micro second durations. v Loadings of 1% core shell rubber particles and 1% milled carbon fiber have an improvement of 8% in V50 ballistic performance over the baseline epoxy sample for 44 mag single wad cutter gas check projectiles. Loadings of 1% multi-walled carbon nanotubes with 1% milled carbon fiber have an improvement of 7.3% in V50 ballistic performance over the baseline epoxy sample. The failure mechanism of the various nano-particle enhanced resin matrices during the ballistic event is discussed through the use of scanning electron microscope images and Raman spectroscopy of the panels after failure. The Raman spectroscopy data shows a Raman shift for the fibers that had an enhancement in the V50 performance through the use of nano-particles. The Raman band for Kevlar® centered at 1,649 cm-1 stemming from the stretching of the C==O bond of the fiber shows to be more sensitive to the residual axial strain, while the Raman band centered at 1,611 cm-1 stemming from the C-C phenyl ring is minimally affected for the CSR enhanced panels due to the failure mechanism of the CSR particles during crack propagation.