BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sundararajan Natarajan

Sundararajan Natarajan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencecond-mat.mtrl-sciphysics.class-ph

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2023arXiv

On the fractional transversely isotropic functionally graded nature of soft biological tissues

This paper focuses on the origin of the poroelastic anisotropic behaviour of the meniscal tissue and its spatially varying properties. We present confined compression creep test results on samples extracted from three parts of the tissue (Central body, Anterior horn and Posterior horn) in three orientations (Circumferential, Radial and Vertical). We show that a poroelastic model in which the fluid flow evolution is ruled by non-integer order operators (fractional Darcy's law) provides accurate agreement with the experimental creep data. The model is validated against two additional sets of experimental data: stress relaxation and fluid loss during the consolidation process measured as weight reduction. Results show that the meniscus can be considered as a transversely isotropic poroelastic material. This behaviour is due to the fluid flow rate being about three times higher in the circumferential direction than in the radial and vertical directions in the body region of the meniscus. In the anterior horn, the elastic properties are transversely isotropic, with the aggregate modulus higher in the radial direction than in the circumferential and vertical directions. The 3D fractional poroelastic model is implemented in finite element software and quantities such as flux of interstitial fluid during the consolidation process, a non-trivial experimental measure, are determined.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive phase field modelling of crack propagation in orthotropic functionally graded materials

In this work, we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials (FGMs). A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement. The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency. The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations. The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle. The study is then extended to the analysis of orthotropic FGMs. It is observed that, if the gradation in fracture properties is neglected, the material gradient plays a secondary role, with the fracture behaviour being dominated by the orthotropy of the material. However, when the toughness increases along the crack propagation path, a substantial gain in fracture resistance is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Continuum modelling of stress diffusion interactions in an elastoplastic medium in the presence of geometric discontinuity

Chemo-mechanical coupled systems have been a subject of interest for many decades now. Previous attempts to solve such models have mainly focused on elastic materials without taking into account the plastic deformation beyond yield, thus causing inaccuracies in failure calculations. This paper aims to study the effect of stress-diffusion interactions in an elastoplastic material using a coupled chemo-mechanical system. The induced stress is dependent on the local concentration in a one way coupled system, and vice versa in a two way coupled system. The time-dependent transient coupled system is solved using a finite element formulation in an open-source finite element solver FEniCS. This paper attempts to computationally study the interaction of deformation and diffusion and its effect on the localization of plastic strain. We investigate the role of geometric discontinuities in scenarios involving diffusing species, namely, a plate with a notch/hole/void and particle with a void/hole/core. We also study the effect of stress concentrations and plastic yielding on the diffusion-deformation. The developed code can be from https://github.com/mrupeshkumar/Elastoplastic-stress-diffusion-coupling

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A hybrid T-Trefftz polygonal finite element for linear elasticity

In this paper, we construct hybrid T-Trefftz polygonal finite elements. The displacement field within the polygon is repre- sented by the homogeneous solution to the governing differential equation, also called as the T-complete set. On the boundary of the polygon, a conforming displacement field is independently defined to enforce continuity of the displacements across the element boundary. An optimal number of T-complete functions are chosen based on the number of nodes of the polygon and degrees of freedom per node. The stiffness matrix is computed by the hybrid formulation with auxiliary displacement frame. Results from the numerical studies presented for a few benchmark problems in the context of linear elasticity shows that the proposed method yield highly accurate results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Application of higher-order structural theory to bending and free vibration analysis of sandwich plates with CNT reinforced composite facesheets

In this paper, the bending and free flexural vibration behaviour of sandwich plates with carbon nanotube (CNT) reinforced facesheets are investigated using QUAD-8 shear flexible element developed based on higher-order structural theory. This theory accounts for the realistic variation of the displacements through the thickness, and the possible discontinuity in slope at the interface, and the thickness stretch affecting the transverse deflection. The in-plane and rotary inertia terms are considered in the formulation. The governing equations obtained using Lagrange's equation of motions are solved for static and dynamic analyses considering a sandwich plate with homogeneous core and CNT reinforced face sheets. The accuracy of the present formulation is tested considering the problems for which solutions are available. A detailed numerical study is carried out based on various higher-order models deduced from the present theory to examine the influence of the volume fraction of the CNT, core-to-face sheet thickness and the plate thickness ratio on the global/local response of different sandwich plates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Assessment of certain higher-order structural models based on global approach for bending analysis of curvilinear composite laminates

In this paper, the performance of different structural models based on global approach in evaluating the static response of curvilinear fibre composite laminates is analyzed. A Co shear flexible Quad-8 element developed based on higher-order structural theory is employed for the present study. The structural theory accounting for the realistic variation of displacements through the thickness and the possible discontinuity in the slope at the interface is considered. Four alternate discrete structural models, deduced from the generic structural model by retaining various terms in the displacement functions are examined for their applicability. The accuracy of the present formulation is demonstrated considering the problems for which analytical solutions are available. A systematic numerical study, assuming different ply-angle and lay-up, is conducted in bringing out the influence of various structural models on the static response of composite laminates with curvilinear fibres.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical evaluation of stress intensity factors and T-stress for interfacial cracks and cracks terminating at the interface without asymptotic enrichment

In this paper, we extend the recently proposed extended scaled boundary finite element method (xSBFEM)~\cite{natarajansong2013} to study fracture parameters of interfacial cracks and cracks terminating at the interface. The approach is also applied to crack growth along the interface and crack deflecting into the material within the context of linear elastic fracture mechanics. Apart from the stress intensity factors, the T-stress can be computed directly from the definitions, without any requirement of path independent integrals. The method aims at improving the capability of the extended finite element method in treating crack tip singularities of cracks at interfaces. An optimum size of the scaled boundary region is presented for multimaterial junctions. The proposed method: (1) does not require special numerical integration technique; (2) does not require a priori knowledge of the asymptotic fields and (3) the stiffness of the region containing the crack tip is computed directly. The robustness of the proposed approach is demonstrated with a few examples in the context of linear elastic fracture mechanics. A discussion on the crack growth along the interface and crack deflecting into the material is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the equivalence between the cell-based smoothed finite element method and the virtual element method

We revisit the cell-based smoothed finite element method (SFEM) for quadrilateral elements and extend it to arbitrary polygons and polyhedrons in 2D and 3D, respectively. We highlight the similarity between the SFEM and the virtual element method (VEM). Based on the VEM, we propose a new stabilization approach to the SFEM when applied to arbitrary polygons and polyhedrons. The accuracy and the convergence properties of the SFEM are studied with a few benchmark problems in 2D and 3D linear elasticity. Later, the SFEM is combined with the scaled boundary finite element method to problems involving singularity within the framework of the linear elastic fracture mechanics in 2D.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Displacement based finite element formulations over polygons: a comparison between Laplace interpolants, strain smoothing and scaled boundary polygon formulation

Three different displacement based finite element formulations over arbitrary polygons are studied in this paper. The formulations considered are: the conventional polygonal finite element method (FEM) with Laplace interpolants, the cell-based smoothed polygonal FEM with simple averaging technique and the scaled boundary polygon formulation. For the purpose of numerical integration, we employ the sub-traingulation for the polygonal FEM and classical Gaussian quadrature for the smoothed FEM and for the scaled boundary polygon formulation. The accuracy and the convergence properties of these formulations are studied with a few benchmark problems in the context of linear elasticity and the linear elastic fracture mechanics. The extension of scaled boundary polygon to higher order polygons is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite element computations on quadtree meshes: strain smoothing and semi-analytical formulation

This short communication discusses two alternate techniques to treat hanging nodes in a quadtree mesh. Both the techniques share similarities, in that, they require only boundary information. Moreover, they do not require an explicit form of the shape functions, unlike the conventional approaches, for example, as in the work of Gupta \cite{gupta1978} or Tabarraei and Sukumar \cite{tabarraeisukumar2005}. Hence, no special numerical integration technique is required. One of the techniques relies on the strain projection procedure, whilst the other is based on the scaled boundary finite element method. Numerical examples are presented to demonstrate the accuracy and the convergence properties of the two techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hygrothermal effects on free vibration and buckling of laminated composites with cutouts

The effect of moisture concentration and the thermal gradient on the free flexural vibration and buckling of laminated composite plates are investigated. The effect of a centrally located cutout on the global response is also studied. The analysis is carried out within the framework of the extended finite element method. A Heaviside function is used to capture the jump in the displacement and an enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The formulation takes into account the transverse shear deformation and accounts for the lamina material properties at elevated moisture concentrations and temperature. The influence of the plate geometry, the geometry of the cutout, the moisture concentration, the thermal gradient and the boundary conditions on the free flexural vibration is numerically studied.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear flutter analysis of functionally graded panels using cell based smoothed finite element method and discrete shear gap technique

In this paper, a cell-based smoothed finite element method with discrete shear gap technique for triangular ele- ments is employed to study the linear flutter characteristics of functionally graded material (FGM) flat panels. The influence of thermal environment, the presence of a centrally located circular cutout and the aerodynamic damping on the supersonic flutter characteristics of flat FGM panels is also investigated. The structural for- mulation is based on the first-order shear deformation theory and the material properties are assumed to be temperature dependent and graded only in the thickness direction according to power law distribution in terms of the volume fraction of its constituent materials. The aerodynamic force is evaluated by considering the first order high mach number approximation to linear potential flow theory. The formulation includes transverse shear deformation and in-plane and rotary inertia effects. The influence of the plate thickness, aspect ratio, boundary conditions, material gradient index, temperature dependent material properties, damping, cutout size, skewness of the plate and boundary conditions on the critical aerodynamic pressure is numerically studied.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Smoothed Finite Element and Genetic Algorithm based optimization for Shape Adaptive Composite Marine Propellers

An optimization scheme using the Cell-based Smoothed Finite Element Method (CS-FEM) combined with a Genetic Algorithm (GA) framework is proposed in this paper to design shape adaptive laminated composite marine propellers. The proposed scheme utilise the bend-twist coupling characteristics of the composites to achieve the required performance. An iterative procedure to evaluate the unloaded shape of the propeller blade is proposed, confirming the manufacturing requirements at the initial stage. The optimization algorithm and codes developed in this work were implemented under a variety of parameter settings and compared against the requirement to achieve an ideally passive pitch varying propeller. Recommendations for the required thickness of the propeller blade to achieve optimum bend-twist coupling performance without resulting large rake deformations are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vibration of functionally graded material plates with cutouts & cracks in thermal environment

In this paper, the effect of a centrally located cutout (circular and elliptical) and cracks emanating from the cutout on the free flexural vibration behaviour of functionally graded material plates in thermal environment is studied. The discontinuity surface is represented independent of the mesh by exploiting the partition of unity method framework. A Heaviside function is used to capture the jump in the displacement across the discontinuity surface and asymptotic branch functions are used to capture the singularity around the crack tip. An enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction. The effective properties of the functionally graded material are estimated using the Mori- Tanaka homogenization scheme and the plate kinematics is based on the first order shear deformation theory. The influence of the plate geometry, the geometry of the cutout, the crack length, the thermal gradient and the boundary conditions on the free flexural vibration is numerically studied.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient recovery-based error estimation for the smoothed finite element method for smooth and singular linear elasticity

An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown to provide significant advantages compared to conventional finite element approximations. In particular, a widely cited strength of such methods is improved accuracy for the same computational cost. Yet, few attempts have been made to directly assess the quality of the results obtained during the simulation by evaluating an estimate of the discretization error. Here we propose a recovery type error estimator based on an enhanced recovery technique. The salient features of the recovery are: enforcement of local equilibrium and, for singular problems a "smooth+singular" decomposition of the recovered stress. We evaluate the proposed estimator on a number of test cases from linear elastic structural mechanics and obtain precise error estimations whose effectivities, both at local and global levels, are improved compared to recovery procedures not implementing these features.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

NURBS-based finite element analysis of functionally graded plates: static bending, vibration, buckling and flutter

In this paper, a non-uniform rational B-spline based iso-geometric finite element method is used to study the static and dynamic characteristics of functionally graded material (FGM) plates. The material properties are assumed to be graded only in the thickness direction and the effective properties are computed either using the rule of mixtures or by Mori-Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation plate theory (FSDT). The shear correction factors are evaluated employing the energy equivalence principle and a simple modification to the shear correction factor is presented to alleviate shear locking. Static bending, mechanical and thermal buckling, linear free flexural vibration and supersonic flutter analysis of FGM plates are numerically studied. The accuracy of the present formulation is validated against available three-dimensional solutions. A detailed numerical study is carried out to examine the influence of the gradient index, the plate aspect ratio and the plate thickness on the global response of functionally graded material plates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Representation of singular fields without asymptotic enrichment in the extended finite element method

In this paper, we replace the asymptotic enrichments around the crack tip in the extended finite element method (XFEM) with the semi-analytical solution obtained by the scaled boundary finite element method (SBFEM). The proposed method does not require special numerical integration technique to compute the stiffness matrix and it improves the capability of the XFEM to model cracks in homogeneous and/or heterogeneous materials without a priori knowledge of the asymptotic solutions. A heaviside enrichment is used to represent the jump across the discontinuity surface. We call the method as the extended scaled boundary finite element method (xSBFEM). Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics show that the proposed method yields accurate results with improved condition number. A simple MATLAB code is annexed to compute the terms in the stiffness matrix, which can easily be integrated in any existing FEM/XFEM code.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method (XFEM) allow discontinuities to be simulated independently of the mesh [1]. This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains [2] to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem, show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the approximation in the smoothed finite element method (SFEM)

This letter aims at resolving the issues raised in the recent short communication [1] and answered by [2] by proposing a systematic approximation scheme based on non-mapped shape functions, which both allows to fully exploit the unique advantages of the smoothed finite element method (SFEM) [3, 4, 5, 6, 7, 8, 9] and resolve the existence, linearity and positivity deficiencies pointed out in [1]. We show that Wachspress interpolants [10] computed in the physical coordinate system are very well suited to the SFEM, especially when elements are heavily distorted (obtuse interior angles). The proposed approximation leads to results which are almost identical to those of the SFEM initially proposed in [3]. These results that the proposed approximation scheme forms a strong and rigorous basis for construction of smoothed finite element methods.

0 citations0 reviewsNo code linkedOpen access