Journal Papers : Scipedia

azrr tt
. Amir

Journal Papers (2022). 2

Abstract
aaaaaaaaaaaaaaa
- 16
- read
CHILDREN’S RIGHTS TO EDUCATION IN SOUTH AFRICA TWENTY YEARS AFTER DEMOCRACY: A REFLECTION ON ACHIEVEMENTS, PROBLEMS AND AREAS FOR FUTURE ACTION
B. Chitsamatanga, B. CHITSAMATANGA

(2020). B BANDA

Abstract
Education is the most powerful weapon that can be used to change the world. These noble words by Nelson Mandela provide a yardstick against which South Africa can measure whether the rights of the child, key of which is the right to education, are progressively being realised. Towards this end, education is seen as the ultimate foundation and solution for the attainment of freedom, democracy and human rights; social justice and equality for South African children. Thus, the question remains: after more than twenty years of democracy, to what extent can South Africa pride herself of having achieved this noble goal? This article reflects on the national effort to realise the rights of children to education against the frameworks provided for in international and regional instruments. The latter require governments’ commitment to ensure that the rights of children to education in particular, are implemented and respected. An analysis will be made on how the new democratic government has succeeded to reverse the discriminatory policies, law and practice of the apartheid era. It further discusses how the existing education system translates the attainment of children rights to education in terms of the standard and quality of teaching and learning; existing school facilities; racial segregation; human and financial resources; and safe school learning environment. The article concludes by identifying measures that need to be taken by all stakeholders in order to promote the realisation of children’s right to education.
Abstract
Education is the most powerful weapon that can be used to change the world. These noble words by Nelson Mandela provide a yardstick against which South Africa can measure [...]
- 5
- read
CLOUD TECHNOLOGY: A CATALYST FOR EFFECTIVE ENTREPRENEURSHIP EDUCATION IN NIGERIA
. AKPANOBONG, U. AKPANOBONG

Nigerian Journal of Business Education (2017). Alghali, M., Alwi, N. H. M., & Ismail, R. (2013). Privacy in cloud based e-Learning. The Second International Conference on Informatics Engineering & Information Science (ICIEIS2013). pp. 355-362. The Society of Digital Informationand Wireless Communication. Alghali, M., Alwi, N. H. M., and Ismail, R. (2014). Challenges and benefits of implementing cloud-based e-learning in developing countries. Proceeding of the Social Sciences Research (ICSSR-2014). pp. 9-10. Amue, G. J., Igwe, S. R. & Abiye, H. (2014). ICT entrepreneurship and small business innovation: A mechanism for sustainability. European Journal of Business and Social Sciences, 3 (6), pp. 103 -112. Retrieved May 20, 2016 from http:// www.ejbss.com/recent.aspx Armbrust, M., Fox, A., Griffith, R., Joseph, A. D., Katz, R. H., Konwinski, A. Zaharia, M. (2009, February 10). Above the clouds: A Berkeley view of cloud computing. Retrieved from http://www.eecs.berkeley.edu/Pubs/TechRpts /2009/EECS-2009- 28.html Eken, H. (2014). Cloud computing security issues and recommendations. International Journal of Intelligent Computing Research (IJICR), 5 (1/2). Ercan, T. (2010). Effective use of cloud computing in educational institutions. Procedia Social and Behavioral Sciences 2, pp. 938-942. Tadapaneni, N. R. (2017). Different Types of Cloud Service Models. Available at SSRN 3614630. Gkikas, D. (2014). The impact of cloud computing on entrepreneurship and Start-ups: case of Greece. Sweden: KTH Industrial Engineering and Management. Larry, C. K. & Best, Z. Y. (2011). The use of ICT in entrepreneurial development. Journal of Research in Entrepreneurship, 16 (6), 339-352. Hamburg, I. & Brien, E. O. (2014). Cloud computing to support e-learning and business in SMEs. International Journal for e-learning Security (IJeLS),4 (3,4). Hamburg, I., Bucksch, S. & Brien, E. O. (2015). ICT–based approaches for Entrepreneurship education. The Tenth International Conference on Internet and Web Application and Services ICIW. Marston, S., Li, Z., Bandyopadhyay, S., Zhang, J. & Ghalsasi, A. (2011). “Cloud computing – The business perspective”. Decision Support Systems, 51, pp. 176- 189. Mathew, S. (2012). Implementation of cloud computing in education – A revolution. International Journal of Computer Theory and Engineering, 4 (3). Mell, P. & Grance, T. (2011). The NIST definition of cloud computing, Recommendation of the National Institute of Standards and Technology (draft). NIST Special Publication, 800 (145), 7. Odunaike, S. A., Olugbara, O. O., & Ojo, S. O. (2012). Using cloud computing to mitigate rural e-learning sustainability and challenges. Proceedings of the World Congress on Engineering and Computer Science, 1. San Francisco, USA. Seleman, Y. E. A. (2015). Modern methods in the development of education using cloud computing applications. Global Journal of Computer Science and Technology: B Cloud and Distributed, 15 (2, 1) Staten, J., Schadler, T., Rymer, J. R. & Wang, C. (2009). Q & A : By 2011, CIOs must answer the question, “Why not run in the cloud ?” (Leaver, S. and Viglianti, A. Interviewers). Sclater, N. (2010). Cloud computing in education. Institute for Information Technologies in Education (IITE). Policy Brief. UNESCO. Retrieved August, 26, 2016 from http://iite.unesco.org/pics/publications/en/files/3214674.pdf Yadav, K. (2014). Role of cloud computing in education. International Journal of Innovative Research in Computer and Communication Engineering, 2 (2).

Abstract
Information technology plays a major role in organising and operating a business by any entrepreneur. Entrepreneurship involves the ability to be creative, innovative and taking risk as well as the capability to plan and manage projects to achieve the set objectives. Entrepreneurship education contributes to economic growth in developed nations. The financial crisis and the economic recession being experienced now in the country has really affected the entrepreneurs who operate Small and Medium Scale Enterprises (SMEs) in Nigeria due to inability to explore modern technology in their businesses. Therefore, there is need for cloud technology to be adopted as a powerful catalyst that would enhance job creation through entrepreneurship education. The potential for cloud computing to boost entrepreneurs as a platform for entrepreneurship education is significant. This technology can be adapted to support entrepreneurship education as it can be scaled over time to offer the required resources. Research has shown that cloud computing has played a catalytic role in boosting education in the world. Thus, the paper examined the potential benefits and issues involved in adopting cloud computing to support entrepreneurship education. It was concluded that the adoption of cloud computing would help to generate thousands of new jobs in Nigeria. It was also recommended among others that the National Information Technology Development Agency should organise workshops, seminars and training for existing entrepreneurs on the benefits of using cloud technology. Cloud technology should also be introduced as one of the courses in business education and allied fields.
Abstract
Information technology plays a major role in organising and operating a business by any entrepreneur. Entrepreneurship involves the ability to be creative, innovative and [...]
- 14
- read
Effect of the tool tilt angle on the material flow and heat generation in Friction Stir Welding
N. Dialami, M. Cervera, M. Chiumenti

Metals (2019). Vol. 9 (1), pp. 28-45

Abstract
This work studies the effect of the tool tilt angle on the generated heat and the material flow in the work pieces joint by Friction Stir Welding (FSW). An apropos kinematic framework together with a two-stage speed-up strategy is adopted to simulate the FSW problem. The effect of tilt angle on the FSWelds is modeled through the contact condition by modifying an enhanced friction model. A rotated friction shear stress is proposed, the angle of rotation depending on the process parameters and the tilt angle. The proposed rotation angle is calibrated by the experimental data provided for a tilt angle 2.5◦ . The differences of generated heat and material flow for the cases of tool with tilt angle of 0◦ and 2.5◦ are discussed. It is concluded that due to the higher temperature, softer material and greater frictional force in the trailing side of the tool, the material flow in the rear side of the FSW tool with the title angle is considerably enhanced, which assists to prevent the generation of defect.
Abstract
This work studies the effect of the tool tilt angle on the generated heat and the material flow in the work pieces joint by Friction Stir Welding (FSW). An apropos kinematic [...]
- 66
- read
Numerical modelling and experimental validation in Selective Laser Melting
M. Chiumenti, E. Neiva, E. Salsi, M. Cervera, E. Badia, J. Moya, Z. Chen, C. Lee, C. Davies

Additive Manufacturing (2017). (preprint) Vol. 18, pp. 171-185

Abstract
In this work a finite-element framework for the numerical simulation of the heat transfer analysis of additive manufacturing processes by powder-bed technologies, such as Selective Laser Melting, is presented. These kind of technologies allow for a layer-by-layer metal deposition process to cost-effectively create, directly from a CAD model, complex functional parts such as turbine blades, fuel injectors, heat exchangers, medical implants, among others. The numerical model proposed accounts for different heat dissipation mechanisms through the surrounding environment and is supplemented by a finite-element activation strategy, based on the born-dead elements technique, to follow the growth of the geometry driven by the metal deposition process, in such a way that the same scanning pattern sent to the numerical control system of the AM machine is used. An experimental campaign has been carried out at the Monash Centre for Additive Manufacturing using an EOSINT-M280 machine where it was possible to fabricate different benchmark geometries, as well as to record the temperature measurements at different thermocouple locations. The experiment consisted in the simultaneous printing of two walls with a total deposition volume of 107 cm3 in 992 layers and about 33,500 s build time. A large number of numerical simulations have been carried out to calibrate the thermal FE framework in terms of the thermophysical properties of both solid and powder materials and suitable boundary conditions. Furthermore, the large size of the experiment motivated the investigation of two different model reduction strategies: exclusion of the powder-bed from the computational domain and simplified scanning strategies. All these methods are analysed in terms of accuracy, computational effort and suitable applications
Abstract
In this work a finite-element framework for the numerical simulation of the heat transfer analysis of additive manufacturing processes by powder-bed technologies, [...]
- 49
- read
Numerical modelling of concrete curing, regarding hydration and temperature phenomena
M. Cervera, R. Faria, T. Prato, J. Oliver

Computers and Structures (2002). Vol. 80 (18–19), pp. 1511-1521

Abstract
A numerical model that accounts for the hydration and aging phenomena during the early ages of concrete curing is presented in a format suitable for a finite element implementation. Assuming the percolation of water through the hydrates already formed as the dominant mechanism of cement hydration, the model adopts an internal variable called hydration degree, whose evolution law is easily calibrated and allows an accurate prediction of the hydration heat production. Compressive strength evolution is related to the aging degree, a concept that accounts for the influences of the hydration and curing temperature on the final mechanical properties of concrete. The model capabilities are illustrated by means of a wide set of experimental tests involving ordinary and high performance concretes, and through the simulation of the concrete curing on a viaduct deck of the Öresund Link.
Abstract
A numerical model that accounts for the hydration and aging phenomena during the early ages of concrete curing is presented in a format suitable for a finite element implementation. [...]
- 48
- read
A strain-based plastic viscous-damage model for massive concrete structures
R. Faria, M. Cervera, J. Oliver

Int. J. Solids and Structures (1998). Vol. 35 (14), pp. 1533-1558

Abstract
Within the framework of continuum damage mechanics, a new constitutive damage model for massive concrete is presented, mainly intended for the seismic analysis of gravity and arch dams. Consistent with thermodynamic requirements, a straindriven formalism is adopted, improving the algorithmic efficiency as much as required for the analysis of large scale problems to become feasible. Two scalar damage variables are introduced as internal variables, as well as a plastic-strain tensor. An extension to account for the concrete strain-rate dependency, suitable for seismic analysis, is presented at the end. The efficiency of numerical predictions from the constitutive model is illustrated through numerical applications and algorithmic implementation is also detailed.
Abstract
Within the framework of continuum damage mechanics, a new constitutive damage model for massive concrete is presented, mainly intended for the seismic analysis of gravity [...]
- 46
- read
Thermo-chemo-mechanical model for concrete. I: Hydration and aging
M. Cervera, T. Prato, J. Oliver

Journal of Engineering Mechanics (1999). Vol. 125 (9), pp. 1018-1027

Abstract
In this work a coupled thermo-chemo-mechanical model for the behavior of concrete at early ages is proposed. The model allows simulation of the observed phenomena of hydration, aging, damage, and creep. It is formulated within an appropriate thermodynamic framework, from which the state equations are derived. In this first part, the formulation and assessment of the thermochemical aspects of the model are presented. It is based on the reactive porous media theory, and it can accurately predict the evolution in time of the hydration degree and the hydration heat production. The evolution of the compressive and tensile strengths and elastic moduli is related to the aging degree, a concept introduced to account for the effect of the curing temperature in the evolution of the mechanical properties. The short- and long-term mechanical behavior is modeled by means of a viscoelastic damage model that accounts for the aging effects. The formulation and assessment of the mechanical part of the model are relegated to a companion paper.
Abstract
In this work a coupled thermo-chemo-mechanical model for the behavior of concrete at early ages is proposed. The model allows simulation of the observed phenomena of hydration, [...]
- 35
- read
Numerical simulation and experimental calibration of Additive Manufacturing by blown powder technology. Part I: thermal analysis
M. Chiumenti, X. Lin, M. Cervera, W. Lei, Y. Zheng, W. Huang

Rapid Prototyping Journal (2017). (preprint) Vol. 23 (2), pp. 448-463

Abstract
In this work the numerical simulation of Additive Manufacturing (AM) processes is addressed. The numerical results are compared with the experimental campaign carried out at SKLSP laboratories, where a Laser Solid Forming (LSF) machine, also referred to as Laser Engineered Net Shaping (LENS), is used to fabricate metal parts directly from CAD models. Ti-6Al-4V metal powder is injected into the molten pool created by a focused, high-energy laser beam and a layer of added material is sinterized according to the laser scanning pattern speciÖed by the user. The objectives of this paper are twofold: Örstly, this work is intended to calibrate the software for the numerical simulation of the AM process, to achieve high accuracy. Secondly, the sensitivity of the numerical model to the process parameters and modelling data is analysed. The numerical model adopts an apropos FE activation technology, which reproduces the same scanning pattern set for the numerical control system of the AM machine. This consists of a complex sequence of polylines, used to deÖne the contour of the component, and hatches patterns to Öll the inner section. The full sequence is given through the Common Layer Interface (CLI) format, a standard format for di§erent manufacturing processes such as Rapid Prototyping (RP), Shape Metal Deposition (SMD) or machining processes, among others. The result is a layer-by-layer metal deposition which can be used to build-up complex structures for components such as turbine blades, aircraft sti§eners, cooling systems, or medical implants, among others.
Abstract
In this work the numerical simulation of Additive Manufacturing (AM) processes is addressed. The numerical results are compared with the experimental campaign carried out [...]
- 32
- read
Derivation of thin plate bending elements with one degree of freedom per node: a simple three node triangle
E. Oñate, M. Cervera

Engineering Computations (1993). Vol. 10 (6), pp. 543-561

Abstract
A general methodology for deriving thin plate bending elements with a single degree of freedom per node is presented. The formulation is based on the combination of a standard C0 finite element interpolation for the deflection field with an independent approximation of the curvatures which are expressed in terms of the deflection gradient along the sides using a finite volume‐like approach. The formulation is particularized for the simplest element of the family, i.e. the three node triangle with three degrees of freedom. The potential of the new element is shown through different examples of application.
Abstract
A general methodology for deriving thin plate bending elements with a single degree of freedom per node is presented. The formulation is based on the combination of a standard C0 [...]
- 31
- read
Numerical analysis of stereolithography processes using the finite element method
G. Bugeda, M. Cervera, G. Lombera, E. Oñate

Rapid Prototyping Journal (1995). Vol. 1 (2), pp. 13-23

Abstract
Stereolithography (SLA) is one of the most important techniques used in rapid prototyping processes. It has a great industrial interest because it allows for dramatic time savings with respect to traditional manufacturing processes. One of the main sources of error in the final dimensions of the prototype is the curl distortion effect owing to the shrinkage of the resin during the SLA process. Presents a study of the influence of different constructive and numerical parameters in the curl distortion, an analysis which was made using the computer code stereolithography analysis program, developed to model SLA processes using the finite element method. Also briefly presents this code.
Abstract
Stereolithography (SLA) is one of the most important techniques used in rapid prototyping processes. It has a great industrial interest because it allows for dramatic time [...]
- 28
- read
The intrinsic time for the streamline upwind Petrov-Galerkin formulation using quadratic elements
R. Codina, E. Oñate, M. Cervera

Comput. Methods Appl. Mech. Engrg., (1992). Vol. 94, pp. 239-262

Abstract
In this paper the functions of the Péclet number that appear in the intrinsic time of the streamline upwind/Petrov-Galerkin (SUPG) formulation are analyzed for quadratic elements. Some related issues such as the computation of the characteristic element length and the introduction of source terms in the one-dimensional model problem are also addressed.
Abstract
In this paper the functions of the Péclet number that appear in the intrinsic time of the streamline upwind/Petrov-Galerkin (SUPG) formulation are analyzed for quadratic [...]
- 24
- read
A computational model for progressive cracking in large dams due to the swelling of concrete
M. Cervera, J. Oliver, E. Herrero, E. Oñate

Engineering Fracture Mechanics (1990). Vol. 35 (1-3), pp. 573-585

Abstract
The paper presents a computational model for the analysis of large concrete gravity dams subjected to severe damage due to internal actions. It has been observed in several operating concrete dams that the combined effects of water intrusion and concrete expansion produce a time-advancing deteriorating process that may endanger the global stability of the construction. The present study was undertaken to simulate numerically the observed phenomena and to enable to predict future developments. Tensile cracking of the concrete is modelled using an elasticfracturing constitutive model. The model is able to simulate in a realistic manner the phenomena of primary and secondary crack initiation, elastic degradation, crack closing and reopening. The triggering of the volumetric expansion due to water intrusion is linked to the onset of cracking at each point, assuming that water enters the dam mostly through the cracks. Temperature and pore-water pressure effects are included using assumed distributions based on available field data. Construction joints are modelled using a frictional joint element, although a “constitutive alternative” is outlined. Close surveillance of the behaviour of a dam that presented this sort of problem was used to calibrate the numerical model and to confirm the obtained results.
Abstract
The paper presents a computational model for the analysis of large concrete gravity dams subjected to severe damage due to internal actions. It has been observed in several [...]
- 21
- read
Structural Analysis of Masonry Historical Constructions. Classical and Advanced Approaches
P. Roca, M. Cervera, G. Gariup, L. Pelà

Archives of Comp. Meths. Engng. (2010). 17 (3), pp. 299-325

Abstract
A review of methods applicable to the study of masonry historical construction, encompassing both classical and advanced ones, is presented. Firstly, the paper offers a discussion on the main challenges posed by historical structures and the desirable conditions that approaches oriented to the modeling and analysis of this type of structures should accomplish. Secondly, the main available methods which are actually used for study masonry historical structures are referred to and discussed. The main available strategies, including limit analysis, simplified methods, FEM macro- or micro-modeling and discrete element methods (DEM) are considered with regard to their realism, computer efficiency, data availability and real applicability to large structures. A set of final considerations are offered on the real possibility of carrying out realistic analysis of complex historic masonry structures. In spite of the modern developments, the study of historical buildings is still facing significant difficulties linked to computational effort, possibility of input data acquisition and limited realism of methods.
Abstract
A review of methods applicable to the study of masonry historical construction, encompassing both classical and advanced ones, is presented. Firstly, the paper offers a discussion [...]
- 21
- read
A penalty finite element method for non‐Newtonian creeping flows
R. Codina, M. Cervera, E. Oñate

Int. J. Numer. Meth. Engng. (1993). Vol. 36 (8), pp. 1395-1412

Abstract
In this paper we present an iterative penalty finite element method for viscous non‐Newtonian creeping flows. The basic idea is solving the equations for the difference between the exact solution and the solution obtained in the last iteration by the penalty method. For the case of Newtonian flows, one can show that for sufficiently small penalty parameters the iterates converge to the incompressible solution. The objective of the present work is to show that the iterative penalization can be coupled with the iterative scheme used to deal with the non‐linearity arising from the constitutive law of non‐Newtonian fluids. Some numerical experiments are conducted in order to assess the performance of the approach for fluids whose viscosity obeys the power law.
Abstract
In this paper we present an iterative penalty finite element method for viscous non‐Newtonian creeping flows. The basic idea is solving the equations for the difference [...]
- 18
- read
Prediction of joint line remnant defect in friction stir welding
N. Dialami, M. Cervera, M. Chiumenti, A. Segatori

Int. J. Mechanical Sciences (2019). Vol. 151, pp. 61-69

Abstract
This work studies computationally and experimentally the joint line remnant defects emerging in friction stir welding (FSW) due to oxide layer propagation into the weld. A finite element-based model is used for the computational simulation. To follow the evolution of the oxide layer originally placed between the butted surfaces, a material tracing technique is incorporated in the numerical model. This approach allows tracking the position of the tracers representing oxide layer particles knowing the nodal velocities. A robust and fast two-stage numerical strategy is adopted for the analysis of FSW process to solve the underlying thermo-mechanical problem. The first stage is a speed-up stage solved on a fixed mesh that allows to quickly obtain the steady state. Oxide layer evolution is traced in the second stage where the rotation of the tool is modelled. Experimentally, to produce a clearly dispersed oxide line in the weld, one of the workpieces is anodized while the other one is taken as extruded. The computationally obtained oxide layer patterns are compared to those obtained experimentally using macrograph analysis of the joint cross section. The effect of the pin features and the process parameters on the final result is studied. The results show that with appropriate modelling of the material tracers in FSW, significant agreement can be attained between the computed and measured post-FSW oxide layer evolved results.
Abstract
This work studies computationally and experimentally the joint line remnant defects emerging in friction stir welding (FSW) due to oxide layer propagation into the weld. A [...]
- 17
- read
Numerical Analysis of Dams with Extensive Cracking due to Concrete Hydration: Simulation of a Real Case
M. Cervera, M. Galindo, J. Oliver

Dam Engng. (1992). Vol. 3 (1), pp. 1-22.

Abstract
The paper presents the computational model developed to study an existing large gravity dam suffering from observed "pathological behavior" resulting concrete hydration. The combined effect of water intrusion and chemical and numerical studies. The present study was undertaken to simulate numerically the observed phenomenon and to predict future developments. Concrete is modelled using an isotropic continuum damage model including visco-elastic effects. The triggering of the volumetric expansion resulting from water intrusion is linked to the onset of damage at each a point, assuming that water enters dam mostly through the damaged areas. Transversal construction joints are modelled using smeared joints. Temperature and porewater pressure effects are included using assumed distributions based on available field data. Results are presented which match reasonably well with the available filed measurements.
Abstract
The paper presents the computational model developed to study an existing large gravity dam suffering from observed "pathological behavior" resulting concrete hydration. [...]
- 16
- read
Simulation of Construction of RCC Dams. I: Temperature and Aging
M. Cervera, J. Oliver, T. Prato

Journal of Structural Engineering (2000). 126 (9), pp. 1053-1060

Abstract
The increasing number of roller compacted concrete (RCC) dams being built around the world demands accurate methodologies for the realistic short- and long-term evaluations of the risk of thermally induced cracking in these constructions. In this work a numerical procedure for the simulation of the construction process of RCC dams is presented. It takes into account the more relevant features of the behavior of concrete at early ages, such as hydration, aging, creep, and damage. A 2D model of the Urugua-í RCC Dam, built in Argentina, is used to perform the corresponding analyses. In this first part only the thermochemical aspects of the simulation of the construction process are presented. The temperature distribution and evolution inside the dam are obtained before and after the completion of the dam. The evolution of the compressive and tensile strengths and elastic moduli and their final distribution inside the dam can also be predicted. Results from 2D and simplified vertical 1D models are compared to assess the validity of the latter, and several parametric studies are carried out. The simulation and discussion of the mechanical aspects of the construction process are relegated to a companion paper that follows.
Abstract
The increasing number of roller compacted concrete (RCC) dams being built around the world demands accurate methodologies for the realistic short- and long-term evaluations [...]
- 16
- read
Numerical modelling of microstructure evolution in Friction Stir Welding (FSW)"
N. Dialami, M. Cervera, M. Chiumenti

Metals (2018). (preprint) Vol. 8 (3), pp. 183

Abstract
This work studies the metallurgical and microstructural aspects of Friction Stir Welding (FSW) in terms of grain size and microhardness. The modelling is based on the combination of an apropos kinematic framework for the local simulation of FSW processes and a material particle tracing technique for tracking the material flow during the weld. The resulting grain size and microhardness values are validated with experimental observations from an identical processed sample. A Sheppard-Wright constitutive relation is adopted to describe the mechanical behavior of AZ31 Mg alloy considered in this work. The strain rate and temperature histories obtained from the numerical model are stored on the tracers. The relationship among the grain size, microhardness, strain rate, and temperature is obtained using Zener-Hollomon parameter and Hall-Petch relationship. A linear description relates the logarithm of average grain size to the logarithm of Zener-Hollomon parameter. The relationship between microhardness and average grain size stands away from the linear trend.
Abstract
This work studies the metallurgical and microstructural aspects of Friction Stir Welding (FSW) in terms of grain size and microhardness. The modelling is based on the combination [...]
- 15
- read
Mixed stabilized finite element methods in nonlinear solid mechanics: Part II: Strain localization
M. Cervera, M. Chiumenti, R. Codina

Comput. Methods Appl. Mech. Engrg., (2010). Vol. 199 (37-40), pp. 2571-2589

Abstract
This paper deals with the question of strain localization associated with materials which exhibit softening due to tensile straining. A standard local isotropic Rankine damage model with strain-softening is used as exemplary constitutive model. Both the irreducible and mixed forms of the problem are examined and stability and solvability conditions are discussed. Lack of uniqueness and convergence difficulties related to the strong material nonlinearities involved are also treated. From this analysis, the issue of local discretization error in the pre-localization regime is deemed as the main difficulty to be overcome in the discrete problem. Focus is placed on low order finite elements with continuous strain and displacement fields (triangular P1P1 and quadrilateral Q1Q1), although the presented approach is very general. Numerical examples show that the resulting procedure is remarkably robust: it does not require the use of auxiliary tracking techniques and the results obtained do not suffer from spurious mesh-bias dependence.
Abstract
This paper deals with the question of strain localization associated with materials which exhibit softening due to tensile straining. A standard local isotropic Rankine damage [...]
- 14
- read
Finite element modeling of multi-pass welding and shaped metal deposition processes
M. Chiumenti, M. Cervera, A. Salmi, C. Agelet de Saracibar, N. Dialami, K. Matsui

Comput. Methods Appl. Mech. Engrg., (2010). Vol. 199 (37-40), pp. 2343-2359

Abstract
This paper describes the formulation adopted for the numerical simulation of the shaped metal deposition process (SMD) and the experimental work carried out at ITP Industry to calibrate and validate the proposed model. The SMD process is a novel manufacturing technology, similar to the multi-pass welding used for building features such as lugs and flanges on fabricated components (see Fig. 1a and b). A fully coupled thermo-mechanical solution is adopted including phase-change phenomena defined in terms of both latent heat release and shrinkage effects. Temperature evolution as well as residual stresses and distortions, due to the successive welding layers deposited, are accurately simulated coupling the heat transfer and the mechanical analysis. The material behavior is characterized by a thermo-elasto-viscoplastic constitutive model coupled with a metallurgical model. Nickel super-alloy 718 is the target material of this work. Both heat convection and heat radiation models are introduced to dissipate heat through the boundaries of the component. An in-house coupled FE software is used to deal with the numerical simulation and an ad-hoc activation methodology is formulated to simulate the deposition of the different layers of filler material. Difficulties and simplifying hypotheses are discussed. Thermo-mechanical results are presented in terms of both temperature evolution and distortions, and compared with the experimental data obtained at the SMD laboratory of ITP.
Abstract
This paper describes the formulation adopted for the numerical simulation of the shaped metal deposition process (SMD) and the experimental work carried out at ITP Industry [...]
- 13
- read
Smeared crack approach: back to the original track
M. Cervera, M. Chiumenti

Int. J. Numer. Anal. Meth. Geomech. (2006). Vol. 30 (12), pp. 1173-1199.

Abstract
This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches.
Abstract
This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack [...]
- 13
- read
Continuum damage model for orthotropic materials: Application to masonry
L. Pelà, M. Cervera, P. Roca

Comput. Methods Appl. Mech. Engrg., (2011). Vol. 200 (9-12), pp. 917-930

Abstract
This paper contributes to the formulation of continuum damage models for orthotropic materials under plane stress conditions. Two stress transformation tensors, related to tensile and compressive stress states, respectively, are used to establish a one-to-one mapping relationship between the orthotropic behaviour and an auxiliary model. This allows the consideration of two individual damage criteria, according to different failure mechanisms, i.e. cracking and crushing. The constitutive model adopted in the mapped space makes use of two scalar variables which monitor the local damage under tension and compression, respectively. The model affords the simulation of orthotropic induced damage, while also accounting for unilateral effects, thanks to a stress tensor split into tensile and compressive contributions. The fundamentals of the method are presented together with the procedure utilized to adjust the model in order to study the mechanical behaviour of masonry material. The validation of the model is carried out by means of comparisons with experimental results on different types of orthotropic masonry at the material level.
Abstract
This paper contributes to the formulation of continuum damage models for orthotropic materials under plane stress conditions. Two stress transformation tensors, related to [...]
- 13
- read
A stabilized formulation for incompressible plasticity using linear triangles and tetrahedra
M. Chiumenti, Q. Valverde, C. Agelet de Saracibar, M. Cervera

International Journal of Plasticity (2004). 20 (8-9), pp. 1487-1504

Abstract
In this paper, a stabilized finite element method to deal with incompressibility in solid mechanics is presented. Both elastic and J2-plastic constitutive behavior have been considered. A mixed formulation involving pressure and displacement fields is used and a continuous linear interpolation is considered for both fields. To circumvent the Babuška–Brezzi condition a stabilization technique based on the orthogonal sub-scale method is introduced. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate for real industrial applications. Results are compared with standard Galerkin and Q1P0 mixed formulations in either elastic or elasto-plastic incompressible problems.
Abstract
In this paper, a stabilized finite element method to deal with incompressibility in solid mechanics is presented. Both elastic and J2-plastic constitutive behavior have been [...]
- 13
- read
Simulation of Construction of RCC Dams. II: Stress and Damage
M. Cervera, T. Prato, J. Oliver

Journal of Structural Engineering (2000). 126 (9), pp. 1062-1068

Abstract
The increasing number of roller compacted concrete (RCC) dams being built around the world demands accurate methodologies for the realistic short- and long-term evaluations of the risk of thermally induced cracking in these constructions. In this work a numerical procedure for the simulation of the construction process of RCC dams is presented. It takes into account the more relevant features of the behavior of concrete at early ages, such as hydration, aging, creep, and damage. A 2D model of the Urugua-;aaı RCC Dam, built in Argentina, is used to perform the corresponding analyses. In this second part of the paper, the mechanical aspects of the simulation are presented; long-term effects are included by incorporating a creep model that naturally accounts for the aging effects, and the risk of tensile damage is also considered. The methodology determines the stress field inside the dam at any time during the construction and in the following years. Results for the reference case assess the suitability of the adopted design. This is compared to alternative studies considering different construction schedules to conclude that for these cases changes should be introduced in the dam design.
Abstract
The increasing number of roller compacted concrete (RCC) dams being built around the world demands accurate methodologies for the realistic short- and long-term evaluations [...]
- 13
- read
Modeling of microstructure evolution of Ti6Al4V for Additive Manufacturing
E. Salsi, M. Chiumenti, M. Cervera

Metals (2018). Vol. 8 (8), pp. 633-657

Abstract
AM processes are characterized by complex thermal cycles that have a deep influence on the microstructural transformations of the deposited alloy. In this work, a general model for the prediction of microstructure evolution during solid state transformations of Ti6Al4V is presented. Several formulations have been developed and employed for modeling phase transformations in other manufacturing processes and, particularly, in casting. The proposed model is mainly based on the combination and modification of some of these existing formulations, leading to a new overall model specifically dedicated to AM. The accuracy and suitability of the integrated model is enhanced, introducing new dedicated features. In fact the model is designed to deal with fast cooling and re-heating cycles typical of AM processes because: (a) it is able to consider incomplete transformations and varying initial content of phases and (b) it can take into account simultaneous transformations.The model is implemented in COMET, an in-house Finite Element (FE)-based framework for the solution of thermo-mechanical engineering problems. The validation of the microstructural model is performed by comparing the simulation results with the data available in the literature. The sensitivity of the model to the variation of material parameters is also discussed.
Abstract
AM processes are characterized by complex thermal cycles that have a deep influence on the microstructural transformations of the deposited alloy. In this work, a general [...]
- 13
- read
In situ measurements and thermo-mechanical simulation of Ti–6Al–4V laser solid forming processes
X. Lu, X. Lin, M. Chiumenti, M. Cervera, Y. Hu, X. Ji, L. Ma, W. Huang

Int. J. Mechanical Sciences (2019). (preprint) Vol. 153-154, pp. 119-130

Abstract
Residual stresses and distortions are two technical obstacles for popularizing the additive manufacturing (AM) technology. The evolution of the stresses in AM components during the thermal cycles of the metal depositing process is not yet clear, and more accurate in situ measurements are necessary to calibrate and validate the numerical tools developed for its simulation. In this work a fully coupled thermo-mechanical analysis to simulate the laser solid forming (LSF) process is carried out. At the same time, an exhaustive experimental campaign is launched to measure the temperature evolution at different locations, as well as the distortions and both the stress and strain fields. The thermal and mechanical responses of single-wall coupons under different process parameters are recorded and compared with the numerical models. Good agreement between the numerical results and the experimental measurements is obtained. Sensitivity analysis demonstrates that the AM process is significantly affected by the laser power and the feeding rate, while poorly influenced by the scanning speed.
Abstract
Residual stresses and distortions are two technical obstacles for popularizing the additive manufacturing (AM) technology. The evolution of the stresses in AM components during [...]
- 13
- read
Residual stress and distortion of rectangular and S-shaped Ti-6Al-4V parts by Directed Energy Deposition: Modelling and experimental calibration
X. Lu, X. Lin, M. Chiumenti, M. Cervera, Y. Hu, X. Ji, L. Ma, H. Yang, W. Huang

Additive Manufacturing (2019). (preprint) Vol. 26, pp. 166-179

Abstract
Residual stresses and distortion in Additive Manufactured (AM) parts are two key obstacles which seriously hinder the wide application of this technology. Nowadays, understanding the thermomechanical behavior induced by the AM process is still a complex task which must take into account the effects of both the process and the material parameters, the microstructure evolution as well as the pre-heating strategy. One of the challenges of this work is to increase the complexity of the geometries used to study the thermomechanical behavior induced by the AM process. The reference geometries are a rectangular and a S-shaped structures made of 44-layers each. The samples have been fabricated by Directed Energy Deposition (DED). In-situ thermal and distortion histories of the substrate are measured in order to calibrate the 3D coupled thermo-mechanical model. Once the numerical results showed a good agreement with the temperature measurements, the validated model has been used to predict the residual stresses and distortions. Different process parameters have been analyzed to study their sensitivity to the process assessment. Different preheating strategies have been also analyzed to check their effectiveness on the mitigation of both distortions and residual stresses. Finally, some simplifications of the actual scanning sequence are proposed to reduce the computational cost without loss of the accuracy of the simulation framework.
Abstract
Residual stresses and distortion in Additive Manufactured (AM) parts are two key obstacles which seriously hinder the wide application of this technology. Nowadays, understanding [...]
- 13
- read
Thermo‐mechanical analysis of industrial solidification processes
M. Cervera, C. Agelet de Saracibar, M. Chiumenti

Int. J. Numer. Meth. Engng. (1999). Vol. 46 (9), pp. 1575-1591

Abstract
The paper presents an up‐to‐date finite element numerical model for fully coupled thermo‐mechanical problems, focussing in the simulation of solidification processes of industrial metal parts. The proposed constitutive model is defined by a thermo‐visco‐elasto‐(visco)plastic free energy function which includes a contribution for thermal multiphase changes. Mechanical and thermal properties are assumed to be temperature‐dependent, and viscous‐like strains are introduced to account for the variation of the elastic moduli during the cooling process. The continuous transition between the initial fluid‐like and the final solid‐like behaviour of the part is modelled by considering separate viscous and elasto‐plastic responses as a function of the solid fraction. Thermo‐mechanical contact conditions between the mould and the part are specifically considered, assuming that the heat flux is a function of the normal pressure and the thermal and mechanical gaps. A fractional step method arising from an operator split of the governing equations is used to solve the non‐linear coupled system of equations, leading to a staggered product formula solution algorithm suitable for large‐scale computations. Representative simulations of industrial solidification processes are shown, and comparison of computed results using the proposed model with available experimental data is given.
Abstract
The paper presents an up‐to‐date finite element numerical model for fully coupled thermo‐mechanical problems, focussing in the simulation of solidification processes [...]
- 12
- read
An orthotropic damage model for the analysis of masonry structures
L. Pelà, M. Cervera, P. Roca

Construction and Building Materials (2013). (preprint) Vol. 41, pp. 957-967.

Abstract
This paper presents a numerical model for nonlinear analysis of masonry structural elements based on Continuum Damage Mechanics. The material is described at the macro-level, i.e. it is modeled as a homogeneous orthotropic continuum. The orthotropic behavior is simulated by means of an original methodology, resulting from the concept of mapped tensors from the anisotropic field to an auxiliary workspace. The application of this idea to strain-based Continuum Damage Models is innovative and leads to several computational benefits. The suitability of the model for representing the behavior of different types of brickwork masonry is shown via the simulation of experimental tests.
Abstract
This paper presents a numerical model for nonlinear analysis of masonry structural elements based on Continuum Damage Mechanics. The material is described at the macro-level, [...]
- 12
- read
Numerical modeling of friction stir welding processes
M. Chiumenti, M. Cervera, C. Agelet de Saracibar, N. Dialami

Comput. Methods Appl. Mech. Engrg., (2013). (preprint) Vol. 254, pp. 353-369

Abstract
This work describes the formulation adopted for the numerical simulation of the friction stir welding (FSW) process. FSW is a solid-state joining process (the metal is not melted during the process) devised for applications where the original metallurgical characteristics must be retained. This process is primarily used on aluminum alloys, and most often on large pieces which cannot be easily heat treated to recover temper characteristics. Heat is either induced by the friction between the tool shoulder and the work pieces or generated by the mechanical mixing (stirring and forging) process without reaching the melting point (solid-state process). To simulate this kind of welding process, a fully coupled thermo-mechanical solution is adopted. A sliding mesh, rotating together with the pin (ALE formulation), is used to avoid the extremely large distortions of the mesh around the tool in the so called stirring zone while the rest of the mesh of the sheet is fixed (Eulerian formulation). The orthogonal subgrid scale (OSS) technique is used to stabilize the mixed velocity–pressure formulation adopted to solve the Stokes problem. This stabilized formulation can deal with the incompressible behavior of the material allowing for equal linear interpolation for both the velocity and the pressure fields. The material behavior is characterized either by Norton–Hoff or Sheppard–Wright rigid thermo-visco-plastic constitutive models. Both the frictional heating due to the contact interaction between the surface of the tool and the sheet, and the heat induced by the visco-plastic dissipation of the stirring material have been taken into account. Heat convection and heat radiation models are used to dissipate the heat through the boundaries. Both the streamline-upwind/Petrov–Galerkin (SUPG) formulation and the OSS stabilization technique have been implemented to stabilize the convective term in the balance of energy equation. The numerical simulations presented are intended to show the accuracy of the proposed methodology and its capability to study real FSW processes where a non-circular pin is often used.
Abstract
This work describes the formulation adopted for the numerical simulation of the friction stir welding (FSW) process. FSW is a solid-state joining process (the metal is not [...]
- 12
- read
Benchmarking on bifurcation and localization in J2 plasticity for plane stress and plane strain conditions
M. Cervera, M. Chiumenti, D. Di Capua

Comput. Methods Appl. Mech. Engrg., (2012). Vol. 241-244, pp. 206-224

Abstract
This paper studies the phenomenon of strain bifurcation and localization in J2plasticity under plane stress and plane strain conditions. Necessary conditions for the outcome of bifurcation, localization, stress boundedness and decohesion are analytically established. It is shown that the explicit consideration of these conditions allows for the determination of localization angles in certain situations of interest that can be used to conduct benchmark tests on finite element formulations. The relative merits of irreducible, (stabilized) mixed and (displacement and/or strain) enhanced formulations are discussed. Numerical examples show that the mixed displacement/pressure formulation is to be preferred to the standard irreducible schemes in order to predict correct failure mechanisms with localized patterns of plastic deformation. Mixed elements are shown to be practically free from mesh directional bias dependence.
Abstract
This paper studies the phenomenon of strain bifurcation and localization in J2plasticity under plane stress and plane strain conditions. Necessary conditions for the outcome [...]
- 12
- read
Seismic evaluation of concrete dams via continuum damage models
M. Cervera, R. Faria, J. Oliver

Earthquake Engineering Structural Dynamics (1995). Vol. 24 (9), pp. 1125-1245

Abstract
In this paper a general methodology for the analysis of large concrete dams subjected to seismic excitation is outlined. It is valid both for gravity dams (2D representation) and arch dams (3D representation). The method allows for non‐linear material behaviour of the dam, ‘transparent fictitious boundaries’ for dealing properly with in‐coming and out‐going seismic waves, and an efficient procedure to deal with dam‐soil‐fluid interaction. The mechanical behaviour of concrete is modelled using an isotropic damage model which allows for tension and compression damage, and exhibits stiffness recovery upon load reversals. Emphasis is placed in the treatment of fluid‐structure interaction, regarding both formulation and efficiency aspects. A gravity dam and an arch dam are analysed subjected to artificially generated earthquakes of different intensities, and the results are used to study the degree of (un)safety of the dams.
Abstract
In this paper a general methodology for the analysis of large concrete dams subjected to seismic excitation is outlined. It is valid both for gravity dams (2D representation) [...]
- 12
- read
On the constitutive modeling of coupled thermomechanical phase-change problems
C. Agelet de Saracibar, M. Cervera, M. Chiumenti

International Journal of Plasticity (2001). 17 (12), pp. 1565-1622.

Abstract
This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient and robust numerical model, able for the numerical simulation of industrial solidification processes. Some of the current issues addressed in the paper are the following. A fractional step method arising from an operator split of the governing differential equations has been used to solve the nonlinear coupled system of equations, leading to a staggered product formula solution algorithm. Nonlinear stability issues are discussed and isentropic and isothermal operator splits are formulated. Within the isentropic split, a strong operator split design constraint is introduced, by requiring that the elastic and plastic entropy, as well as the phase-change induced elastic entropy due to the latent heat, remain fixed in the mechanical problem. The formulation of the model has been consistently derived within a thermodynamic framework. All the material properties have been considered to be temperature dependent. The constitutive behavior has been defined by a thermoviscous/elastoplastic free energy function, including a thermal multiphase change contribution. Plastic response has been modeled by a J2 temperature dependent model, including plastic hardening and thermal softening. The constitutive model proposed accounts for a continuous transition between the initial liquid state, the intermediate mushy state and the final solid state taking place in a solidification process. In particular, a pure viscous deviatoric model has been used at the initial fluid-like state. A thermomecanical contact model, including a frictional hardening and temperature dependent coupled potential, is derived within a fully consistent thermodinamical theory. The numerical model has been implemented into the computational finite element code COMET developed by the authors. Numerical simulations of solidification processes show the good performance of the computational model developed.
Abstract
This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The [...]
- 12
- read
A rate-dependent isotropic damage model for the seismic analysis of concrete dams
M. Cervera, O. Manzoli, J. Oliver

Engineering Computations (1996). Vol. 25 (9), pp. 987-1010.

Abstract
In this paper a rate‐dependent isotropic damage model developed for the numerical analysis of concrete dams subjected to seismic excitation is presented. The model is shown to incorporate two features essential for seismic analysis: stiffness degradation and stiffness recovery upon load reversals and strain‐rate sensitivity. The issue of mesh objectivity is addressed using the concept of the ‘characteristic length’ of the fracture zone, to show that both the softening modulus and the fluidity parameter must depend on it to provide consistent results as the computational mesh is refined. Some aspects of the numerical implementation of the model are also treated, to show that the model can be easily incorporated in any standard non‐linear finite element code. The application of the proposed model to the seismic analysis of a large gravity concrete dam shows that the structural response may vary significantly in terms of the development of damage. The inclusion of rate sensitivity is able to reproduce the experimental observation that the tensile peak strength of concrete can be increased up to 50 percent for the range of strain rates that appear in a structural safety analysis of a dam subjected to severe seismic actions.
Abstract
In this paper a rate‐dependent isotropic damage model developed for the numerical analysis of concrete dams subjected to seismic excitation is presented. The model is shown [...]
- 12
- read
A novel positive/negative projection in energy norm for the damage modeling of quasi-brittle solids
J. Wu, M. Cervera

Int. J. Solids and Structures (2018). (preprint) Vol. 139–140, pp. 250–269

Abstract
The asymmetric tensile/compressive material behavior and microcracks closure-reopening (MCR) effects exhibited by quasi-brittle solids are of significant importance to the nonlinear responses of engineering structures under cyclic loading, e.g., earthquake excitations. Based on our previous work (Cervera et al., 1995; Faria et al., 1998; Wu et al., 2006) this work addresses a novel thermodynamically consistent unilateral damage model for concrete. In particular, the positive/negative projection (PNP) of the effective stress tensor and the additive bi-scalar damage constitutive relation are maintained owing to the conceptual simplicity and computational efficiency. It is found that the classical PNP widely adopted in the literature is not optimal for this damage model, since the resulting stiffness is not always of major symmetry. Consequently, a well-defined free energy potential does not exist in general cases and the model cannot be cast into the framework of thermodynamics with internal variables. Furthermore, the damage induced anisotropy cannot be captured, exhibiting excessive lateral deformationsunder uniaxial tension. To overcome the above issues, a novel PNP, variationally interpreted as the closest point projection of the effective stress in energy norm, is proposed with closed-form solution. With the novel PNP, the secant stiffness tensorof the proposed unilateral damage model always possesses major symmetry and exhibits orthotropic behavior under uniaxial tension and mixed tension/compression. The corresponding thermodynamics framework is then given, resulting in an energy release rate based rounded-Rankine type damage criterion appropriate for tensile failure in quasi-brittle solids. Several numerical examples of single-point verifications and benchmark tests are presented. It is demonstrated that the proposed model is capable of characterizing localized failureof concrete under proportional and non-proportional static loading, as well as the MCR effects under seismic cyclic loading.
Abstract
The asymmetric tensile/compressive material behavior and microcracks closure-reopening (MCR) effects exhibited by quasi-brittle solids are of significant [...]
- 12
- read
A finite volume format for structural mechanics
E. Oñate, M. Cervera, O. Zienkiewicz

Int. J. Numer. Meth. Engng. (1994). Vol. 37 (2), pp. 181-201

Abstract
A general Finite Volume Method (FVM) for the analysis of structural problems is presented. It is shown that the FVM can be considered to be a particular case of finite elements with a non‐Galerkin weighting. For structural analysis this can readily be interpreted as equivalent to the unit displacement method which involves mainly surface integrals. Both displacement and mixed FV formulations are presented for static and dynamic problems.
Abstract
A general Finite Volume Method (FVM) for the analysis of structural problems is presented. It is shown that the FVM can be considered to be a particular case of finite elements [...]
- 11
- read
Softening, localization and stabilization: capture of discontinuous solutions in J2 plasticity
M. Cervera, M. Chiumenti, C. Agelet de Saracibar

Int. J. Numer. Anal. Meth. Geomech. (2004). 28 (5), pp. 373-393.

Abstract
This paper exploits the concept of stabilization techniques to improve the behaviour of mixed linear/linear simplicial elements (triangles and tetrahedra) in incompressible or nearly incompressible situations. Elasto‐J2‐plastic constitutive behaviour has been considered with linear and exponential softening. Two different stabilization methods are used to attain global stability of the corresponding discrete finite element formulation. Implementation and computational aspects are also discussed, showing that a robust application of the proposed formulation is feasible. Numerical examples show that the formulation derived is free of volumetric locking and spurious oscillations of the pressure. The results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with the standard, non‐stabilized, approaches.
Abstract
This paper exploits the concept of stabilization techniques to improve the behaviour of mixed linear/linear simplicial elements (triangles and tetrahedra) in incompressible [...]
- 11
- read
Local–global strategy for the prediction of residual stresses in FSW processes
N. Dialami, M. Cervera, M. Chiumenti, C. Agelet de Saracibar

International Journal of Advanced Manufacturing Technology (2017). (preprint) Vol. 88 (9), pp. 3099–3111

Abstract
This work describes the local–global strategy proposed for the computation of residual stresses in friction stir welding (FSW) processes. A coupling strategy between the analysis of the process zone nearby the pin tool (local level analysis) and the simulation carried out for the entire structure to be welded (global level analysis) is implemented to accurately predict the temperature histories and, thereby, the residual stresses in FSW. As a first step, the local problem solves the material stirring as well as the heat generation induced by the pin and shoulder rotation at the heat affected zone. The Arbitrary Lagrangian Eulerian (ALE) formulation is adopted to deal with the rotation of complex pin shapes. A thermo-rigid-viscoplastic constitutive law is employed to characterize the viscous flow of the material, driven by the high-strain rates induced by the FSW process. A mixed temperature–velocity–pressure finite element technology is used to deal with the isochoric nature of the strains. The output of this local analysis is the heat generated either by plastic dissipation or by friction, and it is used as the power input for the welding analysis at structural (global) level. The global problem is tackled within the Lagrangian framework together with a thermo-elasto-viscoplastic constitutive model. In addition, in this case, the mixed temperature–displacement–pressure format is introduced to deal with the deviatoric nature of the plastic strains. The outcomes of this analysis are the distortions and the residual stresses after welding. The material used in this work is stainless steel 304 L; however, the methodology presented is applicable to a wide range of materials. The proposed numerical strategy is validated by the experimental evidence.
Abstract
This work describes the local–global strategy proposed for the computation of residual stresses in friction stir welding (FSW) processes. A coupling strategy between [...]
- 11
- read
Numerical modeling of the electron beam welding and its experimental validation
M. Chiumenti, M. Cervera, N. Dialami, B. Wu, L. Jinwei, C. Agelet de Saracibar

Finite Elements in Analysis and Design (2016). (preprint) Vol. 121, pp. 118–133

Abstract
Electron Beam Welding (EBW) is a highly efficient and precise welding method increasingly used within the manufacturing chain and of growing importance in different industrial environments such as the aeronautical and aerospace sectors. This is because, compared to other welding processes, EBW induces lower distortions and residual stresses due to the lower and more focused heat input along the welding line. This work describes the formulation adopted for the numerical simulation of the EBW process as well as the experimental work carried out to calibrate and validate it. The numerical simulation of EBW involves the interaction of thermal, mechanical and metallurgical phenomena. For this reason, in this work the numerical framework couples the heat transfer process to the stress analysis to maximize accuracy. An in-house multi-physics FE software is used to deal with the numerical simulation. The definition of an ad hoc moving heat source is proposed to simulate the EB power surface distribution and the corresponding absorption within the work-piece thickness. Both heat conduction and heat radiation models are considered to dissipate the heat through the boundaries of the component. The material behavior is characterized by an apropos thermo-elasto-viscoplastic constitutive model. Titanium-alloy Ti6A14V is the target material of this work. From the experimental side, the EB welding machine, the vacuum chamber characteristics and the corresponding operative setting are detailed. Finally, the available facilities to record the temperature evolution at different thermo-couple locations as well as to measure both distortions and residual stresses are described. Numerical results are compared with the experimental evidence.
Abstract
Electron Beam Welding (EBW) is a highly efficient and precise welding method increasingly used within the manufacturing chain and of growing importance in different industrial [...]
- 11
- read
Computational Modeling and Sub-Grid Scale Stabilization of Incompressibility and Convection in the Numerical Simulation of Friction Stir Welding Processes
C. Agelet de Saracibar, M. Chiumenti, M. Cervera, N. Dialami, A. Seret

Archives of Comp. Meths. Engng. (2014). (preprint) Vol. 21 (1), pp. 3-37

Abstract
This paper deals with the computational modeling and sub-grid scale stabilization of incompressibility and convection in the numerical simulation of the material flow around the probe tool in a friction stir welding (FSW) process. Within the paradigmatic framework of the multiscale stabilization methods, suitable pressure and convective derivative of the temperature sub-grid scale stabilized coupled thermomechanical formulations have been developed using an Eulerian description. Norton-Hoff and Sheppard-Wright thermo-rigid-viscoplastic constitutive material models have been considered. Constitutive equations for the sub-grid scale models have been proposed and an approximation of the sub-grid scale variables has been given. In particular, algebraic sub-grid scale (ASGS) and orthogonal sub-grid scale (OSGS) methods for mixed velocity, pressure and temperature P1/P1/P1 linear elements have been considered. Furthermore, it has been shown that well known classical stabilized formulations, such as the Galerkin least-squares (GLS) for incompressible (or quasi-incompressible) problems or the Streamline Upwind/Petrov-Galerkin (SUPG) method for convection dominant problems, can be recovered as particular cases of the multiscale stabilization framework considered. Using a product formula algorithm for the solution of the coupled thermomechanical problem, the resulting algebraic system of equations has been solved using a staggered procedure in which a mechanical problem, defined by the linear momentum balance equation, under quasi-static conditions, and the incompressibility equation, is solved first at constant temperature. Then a thermal problem, defined by the energy balance equation, is solved keeping constant the mechanical variables, i.e. velocity and pressure. The computational model has been implemented in an enhanced version of the finite element software COMET, developed by the authors at the International Center for Numerical Methods in Engineering (CIMNE). Two numerical examples have been considered. The first one deals with the numerical simulation of a coupled thermomechanical flow in a 2D rectangular domain. Steady-state and transient conditions have been considered. The goal of this numerical example has been the comparison between different sub-grid scale stabilization methods for the velocity and temperature equations. In particular, using a GLS stabilization method for the pressure equation, a comparison between SUPG and OSGS convective stabilization methods has been performed. Additionally, using a SUPG stabilization method for the temperature equation, a comparison between GLS and OSGS pressure stabilization methods has been done. The second example deals with the 3D numerical simulation of a representative FSW process. Numerical results obtained have been compared with experimental results available in the literature. A good agreement on the temperature distribution has been obtained and predicted peak temperatures compare well, both in value and position, with the experimental results available.
Abstract
This paper deals with the computational modeling and sub-grid scale stabilization of incompressibility and convection in the numerical simulation of the material flow around [...]
- 10
- read
Size effect and localization in J2 plasticity
M. Cervera, M. Chiumenti

Int. J. Solids and Structures (2009). Vol. 46 (17), pp. 3301-3312.

Abstract
This paper studies the phenomenon of structural size effect and strain localization in J2 plasticity. Size effect is here understood as the change in the response of a given structure when the spatial dimensions are scaled up or down while the geometry and other relevant properties of the structure are preserved. The work exploits the advantages of the mixed displacement–pressure formulation in incompressible or quasi-incompressible situations. Elasto-J2-plastic constitutive behaviour with regularized softening is considered. Stability issues are discussed to ensure existence and uniqueness of the solution of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is able to solve a wide range of structural scales, including real life engineering applications. The results obtained do not suffer from spurious mesh dependence. Furthermore, the formulation includes the classical theories of perfect plasticity and linear fracture mechanics as limit cases.
Abstract
This paper studies the phenomenon of structural size effect and strain localization in J2 plasticity. Size effect is here understood as the change in the response of a given [...]
- 10
- read
Numerical prediction of temperature and density distributions in selective laser sintering processes
G. Bugeda, M. Cervera, G. Lombera

Rapid Prototyping Journal (1999). Vol. 5 (1), pp. 21-26

Abstract
A finite element model has been developed for the 3D simulation of the sintering of a single track during a selective laser sintering process (SLS). The model takes into account both the thermal and the sintering phenomena involved in the process. Owing to the continuous movement of the laser beam the model takes also into account the transient nature of the problem. This is transformed into a pseudo‐static one through a transformation of the coordinates system of the equations. Nevertheless, this transformation introduces a convective term into the heat equations that produces instabilities in the solution. These instabilities have been solved by using a stream upwind Petrov Galerkin (SUPG) strategy together with a shock capturing scheme. Finally, a fixed point strategy is used for the solution of the analysis. The model has been tested through the solution of some examples.
Abstract
A finite element model has been developed for the 3D simulation of the sintering of a single track during a selective laser sintering process (SLS). The model takes into account [...]
- 10
- read
Thermo-Chemo-Mechanical Model for Concrete. II: Damage and Creep
M. Cervera, J. Oliver, T. Prato

Journal of Engineering Mechanics (1999). Vol. 125 (9), pp. 1028-1038

Abstract
In this work a coupled thermo-chemo-mechanical model for the behavior of concrete at early ages is proposed. This paper presents the formulation and assessment of the mechanical aspects of the model. Short- and long-term mechanical behaviors are modeled via a viscoelastic damage model that accounts for the aging effects. The short-term model is based on the framework of the continuum damage mechanics theory. A novel normalized format of the damage model is proposed, so that the phenomenon of aging is accounted for in a natural fashion. Long-term effects are included by incorporating a creep model inspired in the microprestress-solidification theory.
Abstract
In this work a coupled thermo-chemo-mechanical model for the behavior of concrete at early ages is proposed. This paper presents the formulation and assessment of the mechanical [...]
- 10
- read
Challenges in thermo-mechanical analysis of Friction Stir Welding processes
N. Dialami, M. Chiumenti, M. Cervera, C. Agelet de Saracibar

Archives of Comp. Meths. Engng. (2017). (preprint) Vol. 24 (1), pp. 198-225

Abstract
This paper deals with the numerical simulation of friction stir welding (FSW) processes. FSW techniques are used in many industrial applications and particularly in the aeronautic and aerospace industries, where the quality of the joining is of essential importance. The analysis is focused either at global level, considering the full component to be jointed, or locally, studying more in detail the heat affected zone (HAZ). The analysis at global (structural component) level is performed defining the problem in the Lagrangian setting while, at local level, an apropos kinematic framework which makes use of an efficient combination of Lagrangian (pin), Eulerian (metal sheet) and ALE (stirring zone) descriptions for the different computational sub-domains is introduced for the numerical modeling. As a result, the analysis can deal with complex (non-cylindrical) pin-shapes and the extremely large deformation of the material at the HAZ without requiring any remeshing or remapping tools. A fully coupled thermo-mechanical framework is proposed for the computational modeling of the FSW processes proposed both at local and global level. A staggered algorithm based on an isothermal fractional step method is introduced. To account for the isochoric behavior of the material when the temperature range is close to the melting point or due to the predominant deviatoric deformations induced by the visco-plastic response, a mixed finite element technology is introduced. The Variational Multi Scale method is used to circumvent the LBB stability condition allowing the use of linear/linear P1/P1 interpolations for displacement (or velocity, ALE/Eulerian formulation) and pressure fields, respectively. The same stabilization strategy is adopted to tackle the instabilities of the temperature field, inherent characteristic of convective dominated problems (thermal analysis in ALE/Eulerian kinematic framework). At global level, the material behavior is characterized by a thermo–elasto–viscoplastic constitutive model. The analysis at local level is characterized by a rigid thermo–visco-plastic constitutive model. Different thermally coupled (non-Newtonian) fluid-like models as Norton–Hoff, Carreau or Sheppard–Wright, among others are tested. To better understand the material flow pattern in the stirring zone, a (Lagrangian based) particle tracing is carried out while post-processing FSW results. A coupling strategy between the analysis of the process zone nearby the pin-tool (local level analysis) and the simulation carried out for the entire structure to be welded (global level analysis) is implemented to accurately predict the temperature histories and, thereby, the residual stresses in FSW.
Abstract
This paper deals with the numerical simulation of friction stir welding (FSW) processes. FSW techniques are used in many industrial applications and particularly in the aeronautic [...]
- 10
- read
Finite element analysis and experimental validation of the thermomechanical behavior in laser solid forming of Ti-6Al-4V
X. Lu, X. Lin, M. Chiumenti, M. Cervera, J. Li, L. Ma, L. Wei, Y. Hu, W. Huang

Additive Manufacturing (2018). (preprint) Vol. 21, pp. 30-40

Abstract
A three-dimensional (3D) thermomechanical coupled model for Laser Solid Forming (LSF) of Ti-6Al-4V alloy has been calibrated through experiments of 40-layers metal deposition using different scanning strategies. The sensitivity analysis of the mechanical parameters shows that the thermal expansion coefficient as well as the elastic limit of Ti-6Al-4V have a great impact on the mechanical behavior. Using the validated model and optimal mechanical parameters, the evolution of thermo-mechanical fields in LSF has been analyzed. It has been found that the stresses and distortions develop in two stages, after the deposition of the first layer and during the cooling phase after the manufacturing of the component. The cooling phase is the responsible of 70% of the residual stresses and 60% of the total distortions. The analyses indicate that by controlling the initial substrate temperature (pre-heating phase) and the final cooling phase it is possible to mitigate both distortion and residual stresses. Hence, the influence of different pre-heating procedures on the mechanical fields has been analyzed. The results show that increasing the pre-heating temperature of the substrate is the most effective way to reduce the distortions and residual stresses in Additive Manufacturing.
Abstract
A three-dimensional (3D) thermomechanical coupled model for Laser Solid Forming (LSF) of Ti-6Al-4V alloy has been calibrated through experiments of 40-layers metal [...]
- 10
- read
An orthotropic mesh corrected crack model
M. Cervera

Comput. Methods Appl. Mech. Engrg., (2008). Vol. 197 (17-18), 1603-1619

Abstract
This paper recovers the original spirit of the continuous crack approaches, where displacements jumps across the crack are smeared over the affected elements and the behaviour is established through a softening stress–(total) strain law, using standard finite element displacement interpolations and orthotropic localconstitutive models. The paper focuses on the problem of shear locking observed in the discrete problem when orthotropic models are used. The solution for this drawback is found in the form of a mesh corrected crack model where the structure of the inelastic strain tensor is linked to the geometry of the cracked element. The discrete model is formulated as a non-symmetric orthotropic local damage constitutive model, in which the softening modulus is regularized according to the material fracture energy and the element size. The resulting formulation is easily implemented in standard non-linear FE codes and suitable for engineering applications. Numerical examples show that the results obtained using this crack model do not suffer from dependence on the mesh directional alignment, comparing very favourably with those obtained using related standard isotropic or orthotropic damage models.
Abstract
This paper recovers the original spirit of the continuous crack approaches, where displacements jumps across the crack are smeared over the affected elements and the behaviour [...]
- 9
- read
Mixed stabilized finite element methods in nonlinear solid mechanics: Part I: Formulation
M. Cervera, M. Chiumenti, R. Codina

Comput. Methods Appl. Mech. Engrg., (2010). Vol. 199 (37-40), pp. 2559-2570

Abstract
This paper exploits the concept of stabilized finite element methods to formulate stable mixed stress/displacement and strain/displacement finite elements for the solution of nonlinear solid mechanics problems. The different assumptions and approximations used to derive the methods are exposed. The proposed procedure is very general, applicable to 2D and 3D problems. Implementation and computational aspects are also discussed, showing that a robust application of the proposed formulation is feasible. Numerical examples show that the results obtained compare favorably with those obtained with the corresponding irreducible formulation.
Abstract
This paper exploits the concept of stabilized finite element methods to formulate stable mixed stress/displacement and strain/displacement finite elements for the solution [...]
- 9
- read
A localized mapped damage model for orthotropic materials
L. Pelà, M. Cervera, S. Oller, M. Chiumenti

Engineering Fracture Mechanics (2014). (preprint) Vol. 124-125, pp. 196-216

Abstract
This paper presents an implicit orthotropic model based on the Continuum Damage Mechanics isotropic models. A mapping relationship is established between the behaviour of the anisotropic material and that of an isotropic one. The proposed model is used to simulate the failure loci of common orthotropic materials, such as masonry, fibre-reinforced composites and wood. The damage model is combined with a crack-tracking technique to reproduce the propagation of localized cracks in the discrete FE problem. The proposed numerical model is used to simulate the mixed mode fracture in masonry members with different orientations of the brick layers.
Abstract
This paper presents an implicit orthotropic model based on the Continuum Damage Mechanics isotropic models. A mapping relationship is established between the behaviour of [...]
- 9
- read
Challenges, tools and applications of tracking algorithms in the numerical modelling of cracks in concrete and masonry structures
S. Salaustros, M. Cervera, L. Pelà

Archives of Comp. Meths. Engng. (2018). (preprint) Published online, pp. 1-45

Abstract
The importance of crack propagation in the structural behaviour of concrete and masonry structures has led to the development of a wide range of finite element methods for crack simulation. A common standpoint in many of them is the use of tracking algorithms, which identify and designate the location of cracks within the analysed structure. In this way, the crack modelling techniques, smeared or discrete, are applied only to a restricted part of the discretized domain. This paper presents a review of finite element approaches to cracking focusing on the development and use of tracking algorithms. These are presented in four categories according to the information necessary for the definition and storage of the crack-path. In addition to that, the most utilised criteria for the selection of the crack propagation direction are summarized. The various algorithmic issues involved in the development of a tracking algorithm are discussed through the presentation of a local tracking algorithm based on the smeared crack approach. Challenges such as the modelling of arbitrary and multiple cracks propagating towards more than one direction, as well as multi-directional and intersecting cracking, are detailed. The presented numerical model is applied to the analysis of small- and large-scale masonry and concrete structures under monotonic and cyclic loading.
Abstract
The importance of crack propagation in the structural behaviour of concrete and masonry structures has led to the development of a wide range of finite element methods for [...]
- 9
- read
Mixed linear/linear simplicial elements for incompressible elasticity and plasticity
M. Cervera, M. Chiumenti, Q. Valverde, C. Agelet de Saracibar

Comput. Methods Appl. Mech. Engrg., (2003). Vol. 192 (49-50), pp. 5249-5263

Abstract
This paper exploits the concept of orthogonal sub-grid scales to stabilize the behaviour of mixed linear/linear simplicial elements (triangles and tetrahedra) in incompressible or nearly incompressible situations. Both incompressible elastic and J2-plastic constitutive behaviours have been considered. The different assumptions and approximations used to derive the method are exposed. Implementation and computational aspects are also discussed, showing that a robust application of the proposed formulation is feasible. Numerical examples show that the elements derived are free of volumetric locking and spurious oscillations of the pressure, and that the results obtained compare favourably with those obtained with the Q1P0 quadrilateral.
Abstract
This paper exploits the concept of orthogonal sub-grid scales to stabilize the behaviour of mixed linear/linear simplicial elements (triangles and tetrahedra) in incompressible [...]
- 8
- read
Nonlinear analysis of reinforced concrete plate and shell structures using 20-noded isoparametric brick elements
M. Cervera, E. Hinton, O. Hassan

Computers and Structures (1987). Vol. 25 (6), pp. 845-869

Abstract
This paper describes an efficient and accurate, 3-D finite element model which may be adopted in the nonlinear analysis of reinforced concerete plate and shell structures. Benchmark tests are undertaken to check the computational model.
Abstract
This paper describes an efficient and accurate, 3-D finite element model which may be adopted in the nonlinear analysis of reinforced concerete plate and shell structures. [...]
- 8
- read
Mixed stabilized finite element methods in nonlinear solid mechanics. Part III: Compressible and incompressible plasticity
M. Cervera, M. Chiumenti, L. Benedetti, R. Codina

Comput. Methods Appl. Mech. Engrg., (2015). (preprint) Vol. 285, pp. 752-775

Abstract
This paper presents the application of a stabilized mixed strain/displacement finite element formulation for the solution of nonlinear solid mechanics problems involving compressible and incompressible plasticity. The variational multiscale stabilization introduced allows the use of equal order interpolations in a consistent way. Such formulation presents two advantages when compared to the standard, displacement based, irreducible formulation: (a) it provides enhanced rate of convergence for the strain (and stress) field and (b) it is able to deal with incompressible situations. The first advantage also applies to the comparison with the mixed pressure/displacement formulation. The paper investigates the effect of the improved strain and stress fields in problems involving strain softening and localization leading to failure, using low order finite elements with continuous strain and displacement fields (P1P1 triangles or tetrahedra and Q1Q1quadrilaterals, hexahedra, and triangular prisms) in conjunction with an associative frictional Drucker–Prager plastic model. The performance of the strain/displacement formulation under compressible and nearly incompressible deformation patterns is assessed and compared to a previously proposed pressure/displacement formulation. Benchmark numerical examples show the capacity of the mixed formulation to predict correctly failure mechanisms with localized patterns of strain, virtually free from any dependence of the mesh directional bias. No auxiliary crack tracking technique is necessary.
Abstract
This paper presents the application of a stabilized mixed strain/displacement finite element formulation for the solution of nonlinear solid mechanics problems involving compressible [...]
- 8
- read
Mesh objective tensile cracking via a local continuum damage model and a crack tracking technique
M. Cervera, M. Chiumenti

Comput. Methods Appl. Mech. Engrg., (2006). Vol. 196 (1-3), pp. 304-320

Abstract
This paper describes a procedure for the solution of problems involving tensile cracking using the so-called smeared crack approach, that is, standard finite elements with continuous displacement fields and a standard local constitutive model with strain-softening. An isotropic Rankine damage model is considered. The softening modulus is adjusted according to the material fracture energy and the element size. The resulting continuum and discrete mechanical problems are analyzed and the question of predicting correctly the direction of crack propagation is deemed as the main difficulty to be overcome in the discrete problem. It is proposed to use a crack tracking technique to attain the desired stability and convergence properties of the corresponding formulation. Numerical examples show that the resulting procedure is well-posed, stable and remarkably robust; the results obtained do not seem to suffer from spurious mesh-size or mesh-bias dependence.
Abstract
This paper describes a procedure for the solution of problems involving tensile cracking using the so-called smeared crack approach, that is, standard finite elements [...]
- 8
- read
Modeling material failure in concrete structures under cyclic actions
R. Faria, M. Cervera, J. Oliver

Journal of Structural Engineering (2004). 130 (12), pp. 1997-2005.

Abstract
A constitutive model devised for the analysis of concrete structures, and suitable for generic two- or three-dimensional applications, is presented and validated. For plain concrete a tension-compression distinguishing stress split is performed, and two scalar damage variables account for the degradation induced by the tensile and compressive stress components. As outcomes the model reproduces the stiffness recovery upon load reversal, and it captures the strength enhancement under multiaxial compression. Besides, the simple formulation as well as the extremely reduced number of parameters involved in the concrete model makes it quite suitable for the analysis of real structures, and constitutes a useful design tool. As regards to the nonlinear performance of the steel reinforcement, the explicit Giuffrè–Menegotto–Pinto model is adopted. Efficiency of the global model is illustrated via two seismic applications: one concerning an arch dam, and the other a six-floor reinforced concrete wall. The latter application is presented for validation purposes.
Abstract
A constitutive model devised for the analysis of concrete structures, and suitable for generic two- or three-dimensional applications, is presented and validated. For plain [...]
- 8
- read
A crack-tracking technique for localized damage in quasi-brittle materials
M. Cervera, L. Pelà, R. Clemente, P. Roca

Engineering Fracture Mechanics (2010). Vol. 77 (13), pp. 2431-2450.

Abstract
This work presents a procedure to simulate the growth and propagation of localized tensile cracks on quasi-brittle materials. The so-called smeared damage approach, which consists in standard finite elements and local nonlinear constitutive laws, is recovered and improved in order to represent crack localization and avoid spurious mesh-bias dependence in the discrete problem. This is achieved by means of the implementation of a local crack-tracking algorithm which can reproduce individual (discrete) cracks and ensure objectivity of the finite element problem solution. The performance of the localized damage model is stressed by means of the analyses of structural case-studies. Compared to the Smeared Crack Approach in its original form, the presented procedure shows clearly a better capacity to predict realistic collapse mechanisms. The proposed tracking technique is relatively inexpensive.
Abstract
This work presents a procedure to simulate the growth and propagation of localized tensile cracks on quasi-brittle materials. The so-called smeared damage approach, which [...]
- 8
- read
3D numerical models of FSW processes with non-cylindrical pin
P. Bussetta, N. Dialami, M. Chiumenti, C. Agelet de Saracibar, M. Cervera, R. Boman, J. Ponthot

Advanced Modeling and Simulation in Engineering Sciences (2015). (preprint) Vol. 2 (27), pp. 1-19

Abstract
Friction stir welding process is a relatively recent welding process (patented in 1991). FSW is a solid-state joining process during which materials to be joined are not melted. During the FSW process, the behaviour of the material is at the interface between solid mechanics and fluid mechanics. In this paper, a 3D numerical model of the FSW process with a non-cylindrical tool based on a solid formulation is compared to another one based on a fluid formulation. Both models use advanced numerical techniques such as the Arbitrary Lagrangian Eulerian formulation, remeshing or the Orthogonal Sub-Grid Scale method. It is shown that these two formulations essentially deliver the same results.
Abstract
Friction stir welding process is a relatively recent welding process (patented in 1991). FSW is a solid-state joining process during which materials to be joined are not melted. [...]
- 8
- read
Strain localization analysis of elastic-damaging frictional-cohesive materials: Analytical results and numerical verification
J. Wu, M. Cervera

Materials (2017). (preprint) Vol. 10 (4), pp. 434

Abstract
Damage-induced strain softening is of vital importance for the modeling of localized failure in frictional-cohesive materials. This paper addresses strain localization of damaging solids and the resulting consistent frictional-cohesive crack models. As a supplement to the framework recently established for stress-based continuum material models in rate form (Wu and Cervera 2015, 2016), several classical strain-based damage models, expressed usually in total and secant format, are considered. Upon strain localization of such damaging solids, Maxwell’s kinematics of a strong (or regularized) discontinuity has to be reproduced by the inelastic damage strains, which are defined by a bounded characteristic tensor and an unbounded scalar related to the damage variable. This kinematic constraint yields a set of nonlinear equations from which the discontinuity orientation and damage-type localized cohesive relations can be derived. It is found that for the “Simó and Ju 1987” isotropic damage model, the localization angles and the resulting cohesive model heavily depend on lateral deformations usually ignored in classical crack models for quasi-brittle solids. To remedy this inconsistency, a modified damage model is proposed. Its strain localization analysis naturally results in a consistent frictional-cohesive crack model of damage type, which can be regularized as a classical smeared crack model. The analytical results are numerically verified by the recently-proposed mixed stabilized finite element method, regarding a singly-perforated plate under uniaxial tension. Remarkably, for all of the damage models discussed in this work, the numerically-obtained localization angles agree almost exactly with the closed-form results. This agreement, on the one hand, consolidates the strain localization analysis based on Maxwell’s kinematics and, on the other hand, illustrates versatility of the mixed stabilized finite element method.
Abstract
Damage-induced strain softening is of vital importance for the modeling of localized failure in frictional-cohesive materials. This paper addresses strain localization of [...]
- 8
- read
A stabilized formulation for incompressible elasticity using linear displacement and pressure interpolations
M. Chiumenti, Q. Valverde, C. Agelet de Saracibar, M. Cervera

Comput. Methods Appl. Mech. Engrg., (2002). Vol. 191 (46), pp. 5253-5264

Abstract
In this paper a stabilized finite element method to deal with incompressibility in solid mechanics is presented. A mixed formulation involving pressure and displacement fields is used and a continuous linear interpolation is considered for both fields. To overcome the Babuška–Brezzi condition, a stabilization technique based on the orthogonal sub-scale method is introduced. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate for real industrial applications. Results are compared with standard Galerkin and Q1P0 mixed formulations for nearly incompressible problems in the context of linear elasticity.
Abstract
In this paper a stabilized finite element method to deal with incompressibility in solid mechanics is presented. A mixed formulation involving pressure and displacement fields [...]
- 7
- read
On the formulation of coupled thermoplastic problems with phase-change
C. Agelet de Saracibar, M. Cervera, M. Chiumenti

International Journal of Plasticity (1999). Vol. 15 (1), pp. 1-34

Abstract
This paper deals with a numerical formulation for coupled thermoplastic problems including phase-change phenomena. The final goal is to get an accurate, efficient and robust numerical model, allowing the numerical simulation of solidification processes in the metal casting industry. Some of the current issues addressed in the paper are the following. A fractional step method arising from an operator split of the governing differential equations has been used to solve the nonlinear coupled system of equations, leading to a staggered product formula solution algorithm. Nonlinear stability issues are discussed and isentropic and isothermal operator splits are formulated. Within the isentropic split, a strong operator split design constraint is introduced, by requiring that the elastic and plastic entropy, as well as the phase-change induced elastic entropy due to the latent heat, remain fixed in the mechanical problem. The formulation of the model has been consistently derived within a thermodynamic framework. The constitutive behavior has been defined by a thermoelastoplastic free energy function, including a thermal multiphase change contribution. Plastic response has been modeled by a J2 temperature dependent model, including plastic hardening and thermal softening. A brief summary of the thermomechanical frictional contact model is included. The numerical model has been implemented into the computational Finite Element code COMET developed by the authors. A numerical assessment of the isentropic and isothermal operator splits, regarding the nonlinear stability behavior, has been performed for weakly and strongly coupled thermomechanical problems. Numerical simulations of solidification processes show the performance of the computational model developed.
Abstract
This paper deals with a numerical formulation for coupled thermoplastic problems including phase-change phenomena. The final goal is to get an accurate, efficient and robust [...]
- 7
- read
A thermodynamically consistent plastic-damage framework for localized failure in quasi-brittle solids: Material model and strain localization analysis
J. Wu, M. Cervera

Int. J. Solids and Structures (2016). (preprint) Vol. 88–89, pp. 227–247

Abstract
Aiming for the modeling of localized failure in quasi-brittle solids, this paper addresses a thermodynamically consistent plastic-damage framework and the corresponding strain localization analysis. A unified elastoplastic damage model is first presented based on two alternative kinematic decompositions assuming infinitesimal deformations, with the evolution laws of involved internal variables characterized by a dissipative flow tensor. For the strong (or regularized) discontinuity to form in such inelastic quasi-brittle solids and to evolve eventually into a fully softened one, a novel strain localization analysis is then suggested. A kinematic constraint more demanding than the classical discontinuous bifurcation condition is derived by accounting for the traction continuity and the loading/unloading states consistent with the kinematics of a strong (or regularized) discontinuity. More specifically, the strain jumps characterized by Maxwell’s kinematic condition have to be completely inelastic (energy dissipative). Reproduction of this kinematics implies vanishing of the aforesaid dissipative flow tensorial components in the directions orthogonal to the discontinuity orientation. This property allows naturally developing a localized plastic-damage model for the discontinuity (band), with its orientation and the traction-based failure criterion consistently determined a posteriori from the given stress-based counterpart. The general results are then particularized to the 2D conditions of plane stress and plane strain. It is found that in the case of plane stress, strain localization into a strong (or regularized) discontinuity can occur at the onset of strain softening. Contrariwise, owing to an extra kinematic constraint, in the condition of plane strain some continuous inelastic deformations and substantial re-orientation of principal strain directions in general have to take place in the softening regime prior to strain localization. The classical Rankine, Mohr–Coulomb, von Mises (J2) and Drucker–Prager criteria are analyzed as illustrative examples. In particular, both the closed-form solutions for the discontinuity angles validated by numerical simulations and the corresponding traction-based failure criteria are obtained.
Abstract
Aiming for the modeling of localized failure in quasi-brittle solids, this paper addresses a thermodynamically consistent plastic-damage framework and the corresponding strain [...]
- 7
- read
Tracking multi-directional intersecting cracks in numerical modelling of masonry shear walls under cyclic loading
S. Salaustros, M. Cervera, L. Pelà

Meccanica (2017). (preprint) Vol. 53 (7), pp. 1757–1776

Abstract
In-plane cyclic loading of masonry walls induces a complex failure pattern composed of multiple diagonal shear cracks, as well as flexural cracks. The realistic modelling of such induced localized cracking necessitates the use of costly direct numerical simulations with detailed information on both the properties and geometry of masonry components. On the contrary, computationally efficient macro-models using standard smeared-crack approaches often result in a poor representation of fracture in the simulated material, not properly localized and biased by the finite element mesh orientation. This work proposes a possible remedy to these drawbacks of macro-models through the use of a crack-tracking algorithm. The macro-modelling approach results in an affordable computational cost, while the tracking algorithm aids the mesh-bias independent and localized representation of cracking. A novel methodology is presented that allows the simulation of intersecting and multi-directional cracks using tracking algorithms. This development extends the use of localized crack approaches using tracking algorithms to a wider field of applications exhibiting multiple, arbitrary and interacting cracking. The paper presents also a novel formulation including into an orthotropic damage model the description of irreversible deformations under shear loading. The proposed approach is calibrated through the comparison with an experimental test on a masonry shear wall against in-plane cyclic loading.
Abstract
In-plane cyclic loading of masonry walls induces a complex failure pattern composed of multiple diagonal shear cracks, as well as flexural cracks. The realistic modelling [...]
- 7
- read
Numerical analysis of the manufacturing processes of a mock-up of the ITER NHF First Wall Panel
J. González, M. Chiumenti, M. Cervera, C. Agelet de Saracibar, F. Samaniego, I. Cobo

Fusion Engineering and Design (2018). (preprint) Vol. 135, Part A, pp. 65-73

Abstract
The objective of ITER is to build a new Tokamak, with the goal of demonstrating the scientific and technical feasibility of fusion power. The First Wall Panels are the inner component of the reactor, built with different materials that must support high heat flux levels inside the vacuum vessel. The manufacturing processes of the First Wall are a complex procedure including bending, hipping and cutting procedures which, in general, lead to residual stresses and distortions of the fabricated component. In this work, the analysis of the thermo-mechanical response of a simplified prototype of the ITER NHF First Wall Panel is presented from the numerical point of view. The experimental procedure within each phase of the whole manufacturing process is described. Residual stresses and distortions have been measured and analyzed. The numerical simulation of the manufacturing process includes the description of the main hypothesis, the applied loads and the boundary conditions assumed at every stage of the process. Special attention is paid to the simulation of machining and cutting by means of an ad-hoc element deactivation strategy. The numerical results are compared with the experimental evidence to show the prediction capability and the limitations of the proposed numerical model.
Abstract
The objective of ITER is to build a new Tokamak, with the goal of demonstrating the scientific and technical feasibility of fusion power. The First Wall Panels are the inner [...]
- 7
- read
Experimental validation of a FSW model with an enhanced friction model: application to a threaded cylindrical pin tool
N. Dialami, M. Cervera, M. Chiumenti, A. Segatori, W. Osikowicz

Metals (2017). (preprint) Vol. 7 (11), pp. 491

Abstract
This work adopts a fast and accurate two-stage computational strategy for the analysis of FSW (Friction stir welding) processes using threaded cylindrical pin tools. The coupled thermo-mechanical problem is equipped with an enhanced friction model to include the effect of non-uniform pressure distribution under the pin shoulder. The overall numerical strategy is successfully validated by the experimental measurements provided by the industrial partner (Sapa). The verification of the numerical model using the experimental evidence is not only accomplished in terms of temperature evolution but also in terms of torque, longitudinal, transversal and vertical forces.
Abstract
This work adopts a fast and accurate two-stage computational strategy for the analysis of FSW (Friction stir welding) processes using threaded cylindrical pin tools. The coupled [...]
- 7
- read
A mixed three-field FE formulation for stress accurate analysis including the incompressible limit
M. Chiumenti, M. Cervera, R. Codina

Comput. Methods Appl. Mech. Engrg., (2015). (preprint) Vol. 238, pp. 1095-1116

Abstract
In previous works, the authors have presented the stabilized mixed displacement/pressure formulation to deal with the incompressibility constraint. More recently, the authors have derived stable mixed stress/displacement formulations using linear/linear interpolations to enhance stress accuracy in both linear and non-linear problems. In both cases, the Variational Multi Scale (VMS) stabilization technique and, in particular, the Orthogonal Subgrid Scale (OSS) method allows the use of linear/linear interpolations for triangular and tetrahedral elements bypassing the strictness of the inf–sup condition on the choice of the interpolation spaces. These stabilization procedures lead to discrete problems which are fully stable, free of volumetric locking or stress oscillations. This work exploits the concept of mixed finite element methods to formulate stable displacement/stress/pressure finite elements aimed for the solution of nonlinear problems for both solid and fluid finite element (FE) analyses. The final goal is to design a finite element technology able to tackle simultaneously problems which may involve isochoric behavior (preserve the original volume) of the strain field together with high degree of accuracy of the stress field. These two features are crucial in nonlinear solid and fluid mechanics, as used in most numerical simulations of industrial manufacturing processes. Numerical benchmarks show that the results obtained compare very favorably with those obtained with the corresponding mixed displacement/pressure formulation.
Abstract
In previous works, the authors have presented the stabilized mixed displacement/pressure formulation to deal with the incompressibility constraint. More recently, [...]
- 6
- read
Viscoelasticity and damage model for creep behaviour of historical masonry structures
P. Roca, M. Cervera, L. Pelà, R. Clemente, M. Chiumenti

The Open Civil Engineering Journal (2012). Vol. 6 (Suppl. 1-M7), pp. 188-199.

Abstract
This paper presents a continuum model for the simulation of the viscous effects and the long-term damage accumulation in masonry structures. The rheological model is based on a generalized Maxwell chain representation with a constitutive law utilizing a limited number of internal variables. Thanks to its computational efficiency, this approach is suitable for the analysis of large and complex structures. In the paper, the viscous and damage models are presented and their coupling is discussed. The FE simulation of the construction process of the representative bay of Mallorca Cathedral is presented, together with the analysis of the long-term effects. The parameters of the model are tentatively calibrated on the basis of the time-dependent viscous deformations detected during the cathedral monitoring.
Abstract
This paper presents a continuum model for the simulation of the viscous effects and the long-term damage accumulation in masonry structures. The rheological model is based [...]
- 6
- read
Current Developments on the Coupled Thermomechanical Computational Modeling of Metal Casting Processes
C. Agelet de Saracibar, M. Chiumenti, M. Cervera

Computer Methods in Material Science (2006). Vol. 6 (1-2), pp. 15-25

Abstract
In this paper, current developments on the coupled thermomechanical computational simulation of metal casting processes are presented A thermodynamically consistent constitutive material model is derived from a thermoviscoplastic free energy function. A continuous transition between the initial fluid-like and the final solid-like is modeled by considering a J2 thermoviscoplastic model. Thus, an thermoelastoviscoplastic model, suitable for the solid-like phase, degenerates into a pure thermoviscous model, suitable for the liquid-like phase, according to the evolution of the solid fraction function. A thermomechanical contact model, taking into account the insulated effects of the air-gap due to thermal shrinkage of the part during solidification and cooling, is introduced. A fractional step method, arising from an operator split of the governing differential equations, is considered to solve the coupled problem using a staggered scheme. Within a finite element setting, using low-order interpolation elements, a multiscale stabilization technique is introduced as a convenient framework to overcome the Babuska-Brezzi condition and avoid volumetric locking and pressure instabilities arising in incompressible or quasi-incompressible problems. Computational simulation of industrial castings show the good performance of the model.
Abstract
In this paper, current developments on the coupled thermomechanical computational simulation of metal casting processes are presented A thermodynamically consistent constitutive [...]
- 6
- read
Material flow visualization in Friction Stir Welding via particle tracing
N. Dialami, M. Chiumenti, M. Cervera, C. Agelet de Saracibar, J. Ponthot

Int. J. of Material Forming (2015). (preprint) Vol. 8 (2), pp. 167-181

Abstract
This work deals with the modeling of the material flow in Friction Stir Welding (FSW) processes using particle tracing method. For the computation of particle trajectories, three accurate and computationally efficient integration methods are implemented within a FE model for FSW process: the Backward Euler with Sub-stepping (BES), the 4-th order Runge–Kutta (RK4) and the Back and Forth Error Compensation and Correction (BFECC) methods. Firstly, their performance is compared by solving the Zalesak’s disk benchmark. Later, the developed methodology is applied to some FSW problems providing a quantitative 2D and 3D view of the material transport in the process area. The material flow pattern is compared to the experimental evidence.
Abstract
This work deals with the modeling of the material flow in Friction Stir Welding (FSW) processes using particle tracing method. For the computation of particle trajectories, [...]
- 6
- read
3D numerical models of FSW processes with non-cylindrical pin
P. Bussetta, N. Dialami, M. Chiumenti, C. Agelet de Saracibar, M. Cervera, J. Ponthot

Advances in Materials and Processing Technologies (2016). (preprint) Vol. 1 (3-4), pp. 275-287

Abstract
Friction stir welding (FSW) process is a relatively recent welding process (patented in 1991). FSW is a solid-state joining process during which materials to be joined are not melted. During the FSW process, the behaviour of the material is at the interface between solid mechanics and fluid mechanics. In this paper, a 3D numerical model of the FSW process with a non-cylindrical tool based on a solid formulation is compared to another one based on a fluid formulation. Both models use advanced numerical techniques such as the arbitrary Lagrangian–Eulerian formulation, remeshing or the orthogonal sub-grid scale method. It is shown that these two formulations essentially deliver the same results.
Abstract
Friction stir welding (FSW) process is a relatively recent welding process (patented in 1991). FSW is a solid-state joining process during which materials to be joined are [...]
- 6
- read
Analysis of the effect of provisional ties in the construction and current deformation of Mallorca cathedral
L. Pelà, J. Bourgeois, P. Roca, M. Cervera, M. Chiumenti

International Journal of Architectural Heritage (2016). (preprint) Vol. 10 (4), pp. 418–437

Abstract
This article presents the analysis of the structure of Mallorca Cathedral taking into account the influence on the structural behavior of auxiliary iron ties used during the construction process. Recent studies (Roca et al. 2012, 2013) presented some hypotheses about the construction process of the cathedral. This article complements the previous results by considering the use of auxiliary ties as temporary stabilizing device during the construction. Evidence of the use of ties during the construction has been recognized after a comprehensive survey. The study of the role of such ties and the effect of their later removal are studied by a FE analysis carried out on a representative bay of the structure. The study includes a time-dependent FE analysis after the removal of the ties to assess the long-term structural behavior. The results of the numerical analysis are compared with the deformation trends identified by means of a recent monitoring campaign.
Abstract
This article presents the analysis of the structure of Mallorca Cathedral taking into account the influence on the structural behavior of auxiliary iron ties used during the [...]
- 6
- read
Modelling of Bingham and Herschel–Bulkley flows with mixed P1/P1 finite elements stabilized with orthogonal subgrid scale
E. Moreno, A. LARESE, M. Cervera

Journal of Non-Newtonian Fluid Mechanics (2016). (preprint) Vol. 228, pp. 1-16

Abstract
This paper presents the application of a stabilized mixed pressure/velocity finite element formulation to the solution of viscoplastic non-Newtonian flows. Both Bingham and Herschel–Bulkley models are considered. The detail of the discretization procedure is presented and the Orthogonal Subgrid Scale (OSS) stabilization technique is introduced to allow for the use of equal order interpolations in a consistent way. The matrix form of the problem is given. A series of examples is presented to assess the accuracy of the method by comparison with the results obtained by other authors. The extrusion in a Bingham fluid and the movement of a moving and rotating cylinder are analyzed in detail. The evolution of the streamlines, the yielded and unyielded regions, the drag and lift forces are presented. These benchmark examples show the capacity of the mixed OSS formulation to reproduce the behavior of a Bingham and Herschel–Bulkley flows with the required accuracy.
Abstract
This paper presents the application of a stabilized mixed pressure/velocity finite element formulation to the solution of viscoplastic non-Newtonian flows. Both Bingham and [...]
- 6
- read
An enhanced finite element macro-model for the realistic simulation of localized cracks in masonry structures: A large-scale application
S. Salaustros, L. Pelà, M. Cervera, P. Roca

International Journal of Architectural Heritage (2017). (preprint) Vol. 12 (3), pp. 432-447

Abstract
Finite element macro-modeling approaches are widely used for the analysis of large-scale masonry structures. Despite their efficiency, they still face two important challenges: the realistic representation of damage and a reasonable independency of the numerical results to the used discretization. In this work, the classical smeared crack approach is enhanced with a crack-tracking algorithm, originating from the analysis of localized cracking in quasi-brittle materials. The proposed algorithm is for the first time applied to a large-scale wall exhibiting multiple shear and flexural cracking. Discussion covers structural aspects, as the response of the structure under different assumptions regarding the floor rigidity, but also numerical issues, commonly overlooked in the simulation of large structures, such the mesh-dependency of the numerical results.
Abstract
Finite element macro-modeling approaches are widely used for the analysis of large-scale masonry structures. Despite their efficiency, they still face two important challenges: [...]
- 6
- read
Enhanced friction model for Friction Stir Welding (FSW) analysis: Simulation and experimental validation
N. Dialami, M. Chiumenti, M. Cervera, A. Segatori, W. Osikowicz

Int. J. Mechanical Sciences (2017). (preprint) Vol. 133, pp. 555-567

Abstract
Friction is one of the main heat generation mechanisms in Friction Stir Welding (FSW). This phenomenon occurs between the pin and the workpiece as the rotating tool moves along the weld line. An accurate friction model is essential for obtaining realistic results in a FSW simulation in particular temperature, forces and torque. In this work, a modified Norton's friction law is developed. The suggested enhanced friction model aims at providing not only the realistic temperature field but also the forces and torque. This model does not exclusively relate the frictional shear stress to the sliding velocity; conversely it takes into account the effect of non-uniform pressure distribution under the shoulder, as this latter has an important role in the process of heat generation. Longitudinal, transversal and vertical forces and torque are numerically calculated. The effect of the enhanced friction model is reflected in these forces. In particular, it leads to a more realistic estimation of the transversal and longitudinal forces in comparison with the results obtained using former models. The friction model is successfully validated by the experimental measurements provided by the industrial partner (Sapa). The experimental analysis is performed for the material characterization, the calibration of the friction model and, more generally, the assessment of the overall numerical strategy proposed for the FSW simulation.
Abstract
Friction is one of the main heat generation mechanisms in Friction Stir Welding (FSW). This phenomenon occurs between the pin and the workpiece as the rotating tool moves [...]
- 6
- read
Shear band localization via local J2 continuum damage mechanics
M. Cervera, M. Chiumenti, C. Agelet de Saracibar

Comput. Methods Appl. Mech. Engrg., (2004). Vol. 193 (9-11), pp. 849-880

Abstract
This paper describes a novel formulation for the solution of problems involving shear band localization using a local isotropic J2 continuum damage model and mixed linear simplex (triangles and tetrahedra). Stabilization methods are used to ensure existence and uniqueness of the solution, attaining global and local stability of the corresponding discrete finite element formulation. Consistent residual viscosity is used to enhance robustness and convergence of the formulation. Implementation and computational aspects are also discussed. A simple isotropic local J2 damage constitutive model is considered, either with linear or exponential softening. The softening modulus is regularized according to the material mode II fracture energy and the element size. Numerical examples show that the formulation derived is fully stable and remarkably robust, totally free of volumetric locking and spurious oscillations of the pressure. As a consequence, the results obtained do not suffer from spurious mesh-size or mesh-bias dependence, comparing very favourably with those obtained with the ill-posed standard approaches.
Abstract
This paper describes a novel formulation for the solution of problems involving shear band localization using a local isotropic J2 continuum damage model and mixed linear [...]
- 6
- read
On the computational efficiency and implementation of block-iterative algorithms for nonlinear coupled Problems
M. Cervera, R. Codina, M. Galindo

Engineering Computations (1996). Vol. 13 (6), pp. 4-30

Abstract
Outlines a general methodology for the solution of the system of algebraic equations arising from the discretization of the field equations governing coupled problems. Considers that this discrete problem is obtained from the finite element discretization in space and the finite difference discretization in time. Aims to preserve software modularity, to be able to use existing single field codes to solve more complex problems, and to exploit computer resources optimally, emulating parallel processing. To this end, deals with two well-known coupled problems of computational mechanics – the fluid-structure interaction problem and thermally-driven flows of incompressible fluids. Demonstrates the possibility of coupling the block-iterative loop with the nonlinearity of the problems through numerical experiments which suggest that even a mild nonlinearity drives the convergence rate of the complete iterative scheme, at least for the two problems considered here. Discusses the implementation of this alternative to the direct coupled solution, stating advantages and disadvantages. Explains also the need for online synchronized communication between the different codes used as is the description of the master code which will control the overall algorithm.
Abstract
Outlines a general methodology for the solution of the system of algebraic equations arising from the discretization of the field equations governing coupled problems. Considers [...]
- 5
- read
On the Numerical modeling of the thermomechanical contact for metal casting analysis
M. Chiumenti, C. Agelet de Saracibar, M. Cervera

J. Heat Transfer (2008). Vol. 130 (6), 10 pp.

Abstract
The paper shows the intrinsic difficulties found in the numerical simulation of industrial casting processes using finite element (FE) analysis. Up until now, uncoupled pure thermal simulations have been mostly considered to model solidification and cooling phenomena. However, a fully coupled thermomechanical analysis provides a more complete insight of the casting process and the final outcome regarding the quality of the part. In this type of analysis, the thermomechanical model used plays a role of paramount importance, as the problem is coupled both ways through contact between part and mould. The paper presents the full statement of the problem regarding contact, and it considers the difficulties associated with FE mesh generation and time integration strategy. It also reviews soft and hard algorithms for mechanical contact presenting some new alternatives. Evaluation of coefficients used for thermal contact is also discussed, and a new proposal is presented. Finally, some numerical applications are presented to assess the performance of the proposed strategies both in benchmark and industrial problems.
Abstract
The paper shows the intrinsic difficulties found in the numerical simulation of industrial casting processes using finite element (FE) analysis. Up until now, uncoupled pure [...]
- 5
- read
Continuum FE models for the analysis of a representative bay in Mallorca Cathedral
P. Roca, M. Cervera, L. Pelà

Engineering Structures (2013). (preprint) Vol. 46, pp. 653-670

Abstract
From the theoretical point of view, systems composed by masonry arches or vaults would require, during construction, the simultaneous activation of all structural elements in order to reach the optimum balance of thrusts. This is not obviously the case of complex ancient masonry constructions, whose long and gradual building process may have contributed to their deformed condition and even to damage.In this paper, the possible influence of the construction process as well as that of later long-term deformation on the final condition of the building is investigated in the case of a complex and large historical structure, namely Mallorca Cathedral. A FE code has been specifically developed for the present study. The code is able to account for construction processes through sequential-evolutionary analyses, with the description of masonry mechanical damage and long-term deformation. The representative bay of the cathedral is analyzed taking into account different construction phases, as emerged from historical research. The response of such substructure to transverse earthquake equivalent forces is then investigated. In this case, the damage model is improved with a local crack-tracking algorithm. This numerical strategy models the tensile damage as distinct cracks, leading to a better prediction of realistic collapsing mechanisms.
Abstract
From the theoretical point of view, systems composed by masonry arches or vaults would require, during construction, the simultaneous activation of all structural elements [...]
- 5
- read
An apropos kinematic framework for the numerical modelling of Friction Stir Welding
N. Dialami, M. Chiumenti, M. Cervera, C. Agelet de Saracibar

Computers and Structures (2013). (prepint) Vol. 117, pp. 48-57.

Abstract
This paper describes features of a fully coupled thermo-mechanical model for Friction Stir Welding (FSW) simulation. An apropos kinematic setting for different zones of the computational domain is introduced and an efficient coupling strategy is proposed. Heat generation via viscous dissipation as well as frictional heating is considered. The results of the simulation using the proposed model are compared with the experimental evidence. The effect of slip and stick condition on non-circular pin shapes is analyzed. Simulation of material stirring is also carried out via particle tracing, providing insight of the material flow pattern in the vicinity of the pin.
Abstract
This paper describes features of a fully coupled thermo-mechanical model for Friction Stir Welding (FSW) simulation. An apropos kinematic setting for different zones of the [...]
- 5
- read
Comparison of a Fluid and a Solid Approach for the Numerical Simulation of Friction Stir Welding with a Non‐Cylindrical Pin
P. Bussetta, N. Dialami, R. Boman, M. Chiumenti, C. Agelet de Saracibar, M. Cervera, J. Ponthot

Steel Research International (2014). (preprint) Vol. 85 (6), pp. 968-979

Abstract
Friction stir welding (FSW) process is a solid‐state joining process during which materials to be joined are not melted. As a consequence, the heat‐affected zone is smaller and the quality of the weld is better with respect to more classical welding processes. Because of extremely high strains in the neighborhood of the tool, classical numerical simulation techniques have to be extended in order to track the correct material deformations. The Arbitrary Lagrangian–Eulerian (ALE) formulation is used to preserve a good mesh quality throughout the computation. With this formulation, the mesh displacement is independent from the material displacement. Moreover, some advanced numerical techniques such as remeshing or a special computation of transition interface is needed to take into account non‐cylindrical tools. During the FSW process, the behavior of the material in the neighborhood of the tool is at the interface between solid mechanics and fluid mechanics. Consequently, a numerical model of the FSW process based on a solid formulation is compared to another one based on a fluid formulation. It is shown that these two formulations essentially deliver the same results in terms of pressures and temperatures.
Abstract
Friction stir welding (FSW) process is a solid‐state joining process during which materials to be joined are not melted. As a consequence, the heat‐affected zone is smaller [...]
- 5
- read
High-fidelity prediction of crack formation in 2D and 3D pullout tests
L. Benedetti, M. Cervera, M. Chiumenti

Computers and Structures (2016). (preprint) Vol. 172, pp. 93–109

Abstract
This paper presents the 2D and 3D numerical analysis of pullout tests on steel anchorages in concrete blocks using standard and mixed finite elements. A novel (stabilized) mixed formulation in the variables of total strain ε and displacements u is introduced to overcome the intrinsic deficiencies of the standard displacement-based one in the context of localization of strains, such as mesh dependency. The quasi-brittle behavior of concrete is described through an elastoplastic constitutive law with a local Rankine yielding criterion. The proposed formulation is shown to be a reliable and accurate tool, sensitive to the physical parameters of the pullout tests, but objective with respect to the adopted FE mesh. Furthermore, the mixed ε/u finite element is able to capture the correct failure mechanism with relatively coarse discretizations. At the same time, the spurious behavior of the standard formulation is not alleviated by mesh-refinement.
Abstract
This paper presents the 2D and 3D numerical analysis of pullout tests on steel anchorages in concrete blocks using standard and mixed finite elements. A novel (stabilized) [...]
- 5
- read
An energy-equivalent d+/d− damage model with enhanced microcrack closure-reopening capabilities for cohesive-frictional materials
M. Cervera, C. Tesei

Materials (2017). (preprint) Vol. 10 (4), pp. 433

Abstract
In this paper, an energy-equivalent orthotropic d+/d− damage model for cohesive-frictional materials is formulated. Two essential mechanical features are addressed, the damage-induced anisotropy and the microcrack closure-reopening (MCR) effects, in order to provide an enhancement of the original d+/d−model proposed by Faria et al. 1998, while keeping its high algorithmic efficiency unaltered. First, in order to ensure the symmetry and positive definiteness of the secant operator, the new formulation is developed in an energy-equivalence framework. This proves thermodynamic consistency and allows one to describe a fundamental feature of the orthotropic damage models, i.e., the reduction of the Poisson’s ratio throughout the damage process. Secondly, a “multidirectional” damage procedure is presented to extend the MCR capabilities of the original model. The fundamental aspects of this approach, devised for generic cyclic conditions, lie in maintaining only two scalar damage variables in the constitutive law, while preserving memory of the degradation directionality. The enhanced unilateral capabilities are explored with reference to the problem of a panel subjected to in-plane cyclic shear, with or without vertical pre-compression; depending on the ratio between shear and pre-compression, an absent, a partial or a complete stiffness recovery is simulated with the new multidirectional procedure.
Abstract
In this paper, an energy-equivalent orthotropic d+/d− damage model for cohesive-frictional materials is formulated. Two essential mechanical features are addressed, [...]
- 5
- read
Appraisement of planar, bending and twisting cracks in 3D with isotropic and orthotropic damage models
G. Barbat, M. Cervera, M. Chiumenti

International Journal of Fracture (2018). (preprint) Vol. 210 (1-2), pp. 45–79

Abstract
This paper discusses the modeling of cracking in quasi-brittle materials using isotropic and orthotropic damage constitutive laws. A mixed strain/displacement finite element formulation is used, taking advantage of its enhanced precision and its enforced interelemental strain continuity. On the one hand, this formulation avoids the spurious mesh dependency of the computed solution associated to standard elements and does not require the use of tracking techniques. On the other hand, it greatly alleviates the spurious stress locking associated to the use of orthotropic models on standard finite elements. The performance of several isotropic and orthotropic damage constitutive laws is assessed through an extensive comparison with analytical solutions, numerical tests and experimental evidence reported in the literature. The behavior of the different damage models in terms of crack surface, collapse mechanism and force displacement curves is investigated performing 3D analyses in several conditions including Mode I, Mixed Mode and Mode III fracture. When performing the appraisement of planar, bending and twisting cracks, the enhanced accuracy of the mixed formulation allows for a distinct assessment of the several damage models considered. Aspects related to the behavior of damage models, such as the influence of Poisson’s ratio, the shape of the damage surface and the adoption of isotropic and orthotropic models are investigated and noteworthy conclusions are drawn
Abstract
This paper discusses the modeling of cracking in quasi-brittle materials using isotropic and orthotropic damage constitutive laws. A mixed strain/displacement finite element [...]
- 5
- read
A Phenomenological Model for the Solidification of Eutectic and Hypoeutectic Alloys Including Recalescence and Undercooling
M. Chiumenti, M. Cervera, E. Salsi, A. Zonato

J. of Heat Transfer (2018). (preprint) Vol. 140 (8), pp. 12

Abstract
In this work, a novel phenomenological model is proposed to study the liquid-to-solid phase change of eutectic and hypoeutectic alloy compositions. The objective is to enhance the prediction capabilities of the solidification models based on a-priori definition of the solid fraction as a function of the temperature field. However, the use of models defined at the metallurgical level is avoided to minimize the number of material parameters required. This is of great industrial interest because, on the one hand, the classical models are not able to predict recalescence and undercooling phenomena, and, on the other hand, the complexity as well as the experimental campaign necessary to feed most of the microstructure models available in the literature make their calibration difficult and very dependent on the chemical composition and the treatment of the melt. Contrarily, the proposed model allows for an easy calibration by means of few parameters. These parameters can be easily extracted from the temperature curves recorded at the hot spot of the quick cup test, typically used in the differential thermal analysis (DTA) for the quality control of the melt just before pouring. The accuracy of the numerical results is assessed by matching the temperature curves obtained via DTA of eutectic and hypoeutectic alloys. Moreover, the model is validated in more complex casting experiments where the temperature is measured at different thermocouple locations and the metallurgical features such as grain size and nucleation density are obtained from an exhaustive micrography campaign. The remarkable agreement with the experimental evidence validates the predicting capabilities of the proposed model.
Abstract
In this work, a novel phenomenological model is proposed to study the liquid-to-solid phase change of eutectic and hypoeutectic alloy compositions. The objective is to enhance [...]
- 5
- read
Out-of-plane seismic response and failure mechanism of masonry structures using finite elements with enhanced strain accuracy
G. Vlachakis, M. Cervera, G. Barbat, S. Salaustros

Engineering Failure Analysis (2019). (preprint) Vol. 97, pp. 534-555

Abstract
The out-of-plane response is a complex and at the same time key aspect of the seismic vulnerability of masonry structures. It depends on several factors, some of which are the material properties, the quality of the walls, the geometry of the structure, the connections between structural elements and the stiffness of the diaphragms. During the last years, a wide variety of numerical methods has been employed to assess the out-ofplane behaviour of unreinforced masonry structures. Finite element macro-modelling approaches are among the most famous as they allow modelling large structures at a reasonable computational cost. However, macro-modelling approaches may result in a non-realistic representation of localized cracks and a dependency of the numerical solution on the finite element mesh. Mixed strain/displacement finite elements have been recently proposed as a remedy to the above numerical pathologies. Due to the independent interpolation of strains and displacements these finite element formulations are characterized by an enhanced accuracy in strain localization and crack propagation problems, being at the same time practically mesh independent. For these reasons, mixed finite elements are employed in this work for the out-of-plane assessment of unreinforced masonry structures, being at the same time their first real-scale application. A full-scale experimental campaign of two masonry structures, a stone and a brick one, subjected to shaking-table tests is chosen as reference benchmark. Their structural response under seismic actions is numerically assessed through nonlinear static analysis. The proposed approach is validated through the comparison of the numerical results with the experimental ones, as well as with the results obtained using standard irreducible finite elements.
Abstract
The out-of-plane response is a complex and at the same time key aspect of the seismic vulnerability of masonry structures. It depends on several factors, some of which are [...]
- 5
- read
On the orthogonal subgrid scale pressure stabilization of finite deformation J2 plasticity
C. Agelet de Saracibar, M. Chiumenti, Q. Valverde, M. Cervera

Comput. Methods Appl. Mech. Engrg., (2006). Vol. 195 (9-12), pp. 1224-1251

Abstract
The use of stabilization methods is becoming an increasingly well-accepted technique due to their success in dealing with numerous numerical pathologies that arise in a variety of applications in computational mechanics. In this paper a multiscale finite element method technique to deal with pressure stabilization of nearly incompressibility problems in nonlinear solid mechanics at finite deformations is presented. A J2-flow theory plasticity model at finite deformations is considered. A mixed formulation involving pressure and displacement fields is used as starting point. Within the finite element discretization setting, continuous linear interpolation for both fields is considered. To overcome the Babuška–Brezzi stability condition, a multiscale stabilization method based on the orthogonal subgrid scale (OSGS) technique is introduced. A suitable nonlinear expression of the stabilization parameter is proposed. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate and, therefore, very convenient for practical industrial applications. Numerical results obtained using the OSGS stabilization technique are compared with results provided by the P1 standard Galerkin displacements linear triangular/tetrahedral element, P1/P1 standard mixed linear displacements/linear pressure triangular/tetrahedral element and Q1/P0 mixed bilinear/trilinear displacements/constant pressure quadrilateral/hexahedral element for 2D/3D nearly incompressible problems in the context of a nonlinear finite deformation J2 plasticity model.
Abstract
The use of stabilization methods is becoming an increasingly well-accepted technique due to their success in dealing with numerous numerical pathologies that arise in a variety [...]
- 4
- read
Stress-accurate Mixed FEM for soil failure under shallow foundations involving strain localization in plasticity
L. Benedetti, M. Cervera, M. Chiumenti

Computers and geotechnics (2015). (preprint) Vol. 64, pp. 32-47

Abstract
The development of slip lines, due to strain localization, is a common cause for failure of soil in many circumstances investigated in geotechnical engineering. Through the use of numerical methods – like finite elements – many practitioners are able to take into account complex geometrical and physical conditions in their analyses. However, when dealing with shear bands, standard finite elements display lack of precision, mesh dependency and locking. This paper introduces a (stabilized) mixed finite element formulation with continuous linear strain and displacement interpolations. Von Mises and Drucker–Prager local plasticity models with strain softening are considered as constitutive law. This innovative formulation succeeds in overcoming the limitations of the standard formulation and provides accurate results within the vicinity of the shear bands, specifically without suffering from mesh dependency. Finally, 2D and 3D numerical examples demonstrate the accuracy and robustness in the computation of localization bands, without the introduction of additional tracking techniques as usually required by other methods.
Abstract
The development of slip lines, due to strain localization, is a common cause for failure of soil in many circumstances investigated in geotechnical engineering. Through the [...]
- 4
- read
Explicit mixed strain-displacement finite element for dynamic geometrically non-linear solid mechanics
N. Lafontaine, R. Rossi, M. Cervera, M. Chiumenti

Computational Mechanics (2015). (preprint) Vol. 55 (3), pp. 543-559

Abstract
Low-order finite elements face inherent limitations related to their poor convergence properties. Such difficulties typically manifest as mesh-dependent or excessively stiff behaviour when dealing with complex problems. A recent proposal to address such limitations is the adoption of mixed displacement-strain technologies which were shown to satisfactorily address both problems. Unfortunately, although appealing, the use of such element technology puts a large burden on the linear algebra, as the solution of larger linear systems is needed. In this paper, the use of an explicit time integration scheme for the solution of the mixed strain-displacement problem is explored as an alternative. An algorithm is devised to allow the effective time integration of the mixed problem. The developed method retains second order accuracy in time and is competitive in terms of computational cost with the standard irreducible formulation.
Abstract
Low-order finite elements face inherent limitations related to their poor convergence properties. Such difficulties typically manifest as mesh-dependent or excessively stiff [...]
- 4
- read
On the conformity of strong, regularized, embedded and smeared discontinuity approaches for the modeling of localized failure in solids
M. Cervera, J. Wu

Int. J. Solids and Structures (2015). (preprint) Vol. 71 (1), pp. 19-38

Abstract
Once strain localization occurs in softening solids, inelastic loading behavior is restricted within a narrow band while the bulk unloads elastically. Accordingly, localized failure in solids can be approached by embedding or smearing a traction-based inelastic discontinuity (band) within an (equivalent) elastic matrix along a specific orientation. In this context, the conformity of the strong/regularized and embedded/smeared discontinuity approaches are investigated, regarding the strategies dealing with the kinematics and statics. On one hand, the traction continuity condition imposed in weak form results in the strong and regularized discontinuity approaches, with respect to the approximation of displacement and strain discontinuities. In addition to the elastic bulk, consistent plastic-damage cohesive models for the discontinuities are established. The conformity between the strong discontinuity approach and its regularized counterpart is shown through the fracture energy analysis. On the other hand, the traction continuity condition can also be enforced point-wisely in strong form so that the standard principle of virtual work applies. In this case, the static constraint resulting from traction continuity can be used to eliminate the kinematic variable associated with the discontinuity (band) at the material level. This strategy leads to embedded and smeared discontinuity models for the overall weakened solid which can also be cast into the elastoplastic degradation framework with a different kinematic decomposition. Being equivalent to the kinematic constraint guaranteeing stress continuity upon strain localization, Mohr’s maximization postulate is adopted for the determination of the discontinuity orientation. Closed-form results are presented in plane stress conditions, with the classical Rankine, Mohr–Coulomb, von Mises and Drucker–Prager criteria as illustrative examples. The orientation of the discontinuity (band) and the stress-based failure criteria consistent with the given traction-based counterparts are derived. Finally, a generic failure criterion of either elliptic, parabolic or hyperbolic type, appropriate for the modeling of mixed-mode failure, is analyzed in a unified manner. Furthermore, a novel method is proposed to calibrate the involved mesoscopic parameters from available macroscopic test data, which is then validated against Willam’s numerical test.
Abstract
Once strain localization occurs in softening solids, inelastic loading behavior is restricted within a narrow band while the bulk unloads elastically. Accordingly, localized [...]
- 4
- read
Explicit mixed strain-displacement finite elements for compressible and quasi-incompressible elasticity and plasticity
M. Cervera, N. Lafontaine, R. Rossi, M. Chiumenti

Computational Mechanics (2016). (Preprint) Vol. 58 (3), pp. 511–532

Abstract
This paper presents an explicit mixed finite element formulation to address compressible and quasi-incompressible problems in elasticity and plasticity. This implies that the numerical solution only involves diagonal systems of equations. The formulation uses independent and equal interpolation of displacements and strains, stabilized by variational subscales. A displacement sub-scale is introduced in order to stabilize the mean-stress field. Compared to the standard irreducible formulation, the proposed mixed formulation yields improved strain and stress fields. The paper investigates the effect of this enhancement on the accuracy in problems involving strain softening and localization leading to failure, using low order finite elements with linear continuous strain and displacement fields (P1P1 triangles in 2D and tetrahedra in 3D) in conjunction with associative frictional Mohr–Coulomb and Drucker–Prager plastic models. The performance of the strain/displacement formulation under compressible and nearly incompressible deformation patterns is assessed and compared to analytical solutions for plane stress and plane strain situations. Benchmark numerical examples show the capacity of the mixed formulation to predict correctly failure mechanisms with localized patterns of strain, virtually free from any dependence of the mesh directional bias. No auxiliary crack tracking technique is necessary.
Abstract
This paper presents an explicit mixed finite element formulation to address compressible and quasi-incompressible problems in elasticity and plasticity. This implies that [...]
- 4
- read
A fast and accurate two-stage strategy to evaluate the effect of the pin tool profile on metal flow, torque and forces during friction stir welding
N. Dialami, M. Cervera, M. Chiumenti, C. Agelet de Saracibar

Int. J. Mechanical Sciences (2017). (preprint) Vol. 122, pp. 215–227

Abstract
Pin geometry is a fundamental consideration in friction stir welding (FSW). It influences the thermal behaviour, material flow and forces during the weld and reflects on the joint quality. This work studies four pin tools with circular, triflute, trivex, and triangular profiles adopting a validated model of FSW process developed by the authors. The effect of the rotating tool geometry on the flow behaviour and process outcomes is analysed. Additionally, longitudinal and transversal forces and torque are numerically calculated and compared for the different pin shapes. The study is carried out for slip and stick limiting friction cases between pin and workpiece. The main novelties of the paper are a “speed-up” two-stage simulation methodology and a piecewise linear version of the constitutive model, both of them conceived for the use in real case industrial applications, where the achievement of accuracy with affordable simulation times is of importance. The Norton-Hoff constitutive model is adopted to characterize the material behaviour during the weld. The piecewise linear version of the model developed by the authors greatly facilitates the convergence of the numerical solution ensuring both computational efficiency and accuracy. A two-stage computational procedure is applied. In the first stage, a forced transient is carried out; in the second one, the magnitudes of interest are computed. The study shows that the proposed modelling approach can be used to predict and interpret the FSW behaviour for a specific pin geometry. Moreover, the reduction of the simulation time using the two-stage strategy can be up to 90%, compared to a standard single stage strategy.
Abstract
Pin geometry is a fundamental consideration in friction stir welding (FSW). It influences the thermal behaviour, material flow and forces during the weld and [...]
- 4
- read
3D numerical modelling of twisting cracks under bending and torsion of skew notched beams
L. Benedetti, M. Cervera, M. Chiumenti

Engineering Fracture Mechanics (2017). (preprint) Vol. 176, pp. 235–256

Abstract
The testing of mode III and mixed mode failure is every so often encountered in the dedicated literature of mechanical characterization of brittle and quasi-brittle materials. In this work, the application of the mixed strain displacement ε-u finite element formulation to three examples involving skew notched beams is presented. The use of this FE technology is effective in problems involving localization of strains in softening materials. The objectives of the paper are: (i) to test the mixed formulation in mode III and mixed mode failure and (ii) to present an enhancement in terms of computational time given by the kinematic compatibility between irreducible displacement-based and the mixed strain-displacement elements. Three tests of skew-notched beams are presented: firstly, a three point bending test of a PolyMethyl MethaAcrylate beam; secondly, a torsion test of a plain concrete prismatic beam with square base; finally, a torsion test of a cylindrical beam made of plain concrete as well. To describe the mechanical behavior of the material in the inelastic range, Rankine and Drucker-Prager failure criteria are used in both plasticity and isotropic continuum damage formats. The proposed mixed formulation is capable of yielding results close to the experimental ones in terms of fracture surface, peak load and global loss of carrying capability. In addition, the symmetric secant formulation and the compatibility condition between the standard irreducible method and the strain-displacement one is exploited, resulting in a significant speedup of the computational procedure.
Abstract
The testing of mode III and mixed mode failure is every so often encountered in the dedicated literature of mechanical characterization of brittle and quasi-brittle [...]
- 4
- read
Mesh objective modeling of cracks using continuous linear strain and displacement interpolations
M. Cervera, M. Chiumenti, R. Codina

Int. J. Numer. Meth. Engng. (2011). Vol. 87 (10), pp. 962-987

Abstract
The paper addresses the problem of tensile and mixed‐mode cracking within the so‐called smeared crack approach. Because lack of point‐wise convergence on stresses is deemed as the main difficulty to be overcome in the discrete problem, a (stabilized) mixed formulation with continuous linear strain and displacement interpolations is used. The necessary convergence rate can be proved for such a formulation, at least in the linear problem. Two standard local isotropic Rankine damage models with strain‐softening, differing in the definition of the damage criteria, are used as discrete constitutive model. Numerical examples demonstrate the application of the proposed formulation using linear triangular P1P1 and bilinear quadrilateral Q1Q1 mixed elements. The results obtained do not suffer from spurious mesh‐bias dependence without the use of auxiliary tracking techniques.
Abstract
The paper addresses the problem of tensile and mixed‐mode cracking within the so‐called smeared crack approach. Because lack of point‐wise convergence on stresses is [...]
- 3
- read
On the equivalence between traction- and stress-based approaches for the modeling of localized failure in solids
J. Wu, M. Cervera

Journal of the Mechanics and Physics of Solids (2015). (preprint) Vol. 82, pp. 137-163

Abstract
This work investigates systematically traction- and stress-based approaches for the modeling of strong and regularized discontinuities induced by localized failure in solids. Two complementary methodologies, i.e., discontinuities localized in an elastic solid and strain localization of an inelastic softening solid, are addressed. In the former it is assumed a priori that the discontinuity forms with a continuous stress field and along the known orientation. A traction-based failure criterion is introduced to characterize the discontinuity and the orientation is determined from Mohr's maximization postulate. If the displacement jumps are retained as independent variables, the strong/regularized discontinuity approaches follow, requiring constitutive models for both the bulk and discontinuity. Elimination of the displacement jumps at the material point level results in the embedded/smeared discontinuity approaches in which an overall inelastic constitutive model fulfilling the static constraint suffices. The second methodology is then adopted to check whether the assumed strain localization can occur and identify its consequences on the resulting approaches. The kinematic constraint guaranteeing stress boundedness and continuity upon strain localization is established for general inelastic softening solids. Application to a unified stress-based elastoplastic damage model naturally yields all the ingredients of a localized model for the discontinuity (band), justifying the first methodology. Two dual but not necessarily equivalent approaches, i.e., the traction-based elastoplastic damage model and the stress-based projected discontinuity model, are identified. The former is equivalent to the embedded and smeared discontinuity approaches, whereas in the later the discontinuity orientation and associated failure criterion are determined consistently from the kinematic constraint rather than given a priori. The bi-directional connections and equivalence conditions between the traction- and stress-based approaches are classified. Closed-form results under plane stress condition are also given. A generic failure criterion of either elliptic, parabolic or hyperbolic type is analyzed in a unified manner, with the classical von Mises (J2), Drucker–Prager, Mohr–Coulomb and many other frequently employed criteria recovered as its particular cases.
Abstract
This work investigates systematically traction- and stress-based approaches for the modeling of strong and regularized discontinuities induced by localized failure in [...]
- 3
- read
Finite element modelling of internal and multiple localized cracks
S. Salaustros, L. Pelà, M. Cervera, P. Roca

Computational Mechanics (2017). (preprint) Vol. 59 (2), pp. 299–316

Abstract
Tracking algorithms constitute an efficient numerical technique for modelling fracture in quasi-brittle materials. They succeed in representing localized cracks in the numerical model without mesh-induced directional bias. Currently available tracking algorithms have an important limitation: cracking originates either from the boundary of the discretized domain or from predefined “crack-root” elements and then propagates along one orientation. This paper aims to circumvent this drawback by proposing a novel tracking algorithm that can simulate cracking starting at any point of the mesh and propagating along one or two orientations. This enhancement allows the simulation of structural case-studies experiencing multiple cracking. The proposed approach is validated through the simulation of a benchmark example and an experimentally tested structural frame under in-plane loading. Mesh-bias independency of the numerical solution, computational cost and predicted collapse mechanisms with and without the tracking algorithm are discussed.
Abstract
Tracking algorithms constitute an efficient numerical technique for modelling fracture in quasi-brittle materials. They succeed in representing localized cracks in the numerical [...]
- 3
- read
Cracking of quasi-brittle structures under monotonic and cyclic loadings: A d+/d− damage model with stiffness recovery in shear
M. Cervera, C. Tesei, G. Ventura

Int. J. Solids and Structures (2017). (preprint) Vol. 135, pp. 148–171

Abstract
In the present paper, a new d+/d− damage model apt for quasi-brittle materials is described and its validation is carried out considering unreinforced concrete, reinforced concrete and masonry elements. Two independent scalar damage variables, d+ and d−, in combination with the split of the reversible strain tensor into its positive and negative counterparts, are adopted in order to simulate the pronounced dissimilar response under tension and compression, typical of these materials. An energy-equivalent framework is considered for representing the orthotropy induced in the material by the degradation process, with the consequence that a thermodynamically consistent constitutive operator, positive definite, symmetric and strain-driven, is derived. In addition to the degradation parameters, the permanent strain tensor is also contemplated by the model and a modification of the exponential softening modulus is proposed in order to treat the evolution of the two causes of dissipation, damage and irreversible deformations, in a coupled way. The formulation is integrated with a multidirectional damage procedure, addressed to extend the microcrack closure-reopening (MCR) capabilities of the model to shear cyclic conditions, characterized by orthogonal (or however intersecting) sets of cracks. Maintaining unaltered the dependence of the constitutive law from two scalar indeces, d+ and d−, this approach activates or deactivates a tensile(compressive) damage value on the base of the current maximum (minimum) principal strain direction. In correspondence with damage activation (crack opening) or deactivation (crack closure), a smooth transition is introduced, in order to avoid abrupt changes in stiffness and enhance the numerical performance and robustness of the multidirectional procedure. The adequacy of the proposed constitutive model in reproducing experimental results has been proven for both monotonic and cyclic loading conditions. The two examples of application involving cyclic loads, dominated by shear, constitute a validation of the multidirectional damage approach, showing how the suitable representation of unilateral effects and permanent deformations is essential to model the observed structural response in terms of maximum resistance, evolution of stiffness degradation and dissipation capacity.
Abstract
In the present paper, a new d+/d− damage model apt for quasi-brittle materials is described and its validation is carried out considering unreinforced [...]
- 3
- read
Finite element modeling of quasi-brittle cracks in 2D and 3D with enhanced strain accuracy
M. Cervera, G. Barbat, M. Chiumenti

Computational Mechanics (2017). (preprint) Vol. 60 (5), pp. 767–796

Abstract
This paper discusses the finite element modeling of cracking in quasi-brittle materials. The problem is addressed via a mixed strain/displacement finite element formulation and an isotropic damage constitutive model. The proposed mixed formulation is fully general and is applied in 2D and 3D. Also, it is independent of the specific finite element discretization considered; it can be equally used with triangles/tetrahedra, quadrilaterals/hexahedra and prisms. The feasibility and accuracy of the method is assessed through extensive comparison with experimental evidence. The correlation with the experimental tests shows the capacity of the mixed formulation to reproduce the experimental crack path and the force–displacement curves with remarkable accuracy. Both 2D and 3D examples produce results consistent with the documented data. Aspects related to the discrete solution, such as convergence regarding mesh resolution and mesh bias, as well as other related to the physical model, like structural size effect and the influence of Poisson’s ratio, are also investigated. The enhanced accuracy of the computed strain field leads to accurate results in terms of crack paths, failure mechanisms and force displacement curves. Spurious mesh dependency suffered by both continuous and discontinuous irreducible formulations is avoided by the mixed FE, without the need of auxiliary tracking techniques or other computational schemes that alter the continuum mechanical problem.
Abstract
This paper discusses the finite element modeling of cracking in quasi-brittle materials. The problem is addressed via a mixed strain/displacement finite element formulation [...]
- 3
- read
A smeared‐embedded mesh‐corrected damage model for tensile cracking
M. Cervera

Int. J. Numer. Meth. Engng. (2008). Vol. 76 (12), pp. 1930-1954

Abstract
Traditional smeared orthotropic models display an unacceptable dependence of the solution on the alignment of the mesh, which usually manifests as stress locking. A solution for this drawback is proposed in this paper by adopting the concept of embedded inelastic strains, rather than displacement jumps, and by linking the structure of the inelastic strain to the geometry of the cracked element. The resulting model, applicable to linear 3‐noded triangles, is formulated as a non‐symmetric orthotropic local damage constitutive model, with the softening modulus regularized according to the material fracture energy and the element size. Analytical and numerical results show that this approach is effective in removing the locking problem as well as efficient from the computational point of view.
Abstract
Traditional smeared orthotropic models display an unacceptable dependence of the solution on the alignment of the mesh, which usually manifests as stress locking. A solution [...]
- 2
- read
Preconditioned conjugate gradient method for the non‐linear finite element analysis with particular reference to 3D reinforced concrete structures
M. Cervera, Y. Liu, E. Hinton

Engineering Computations (1986). Vol. 3 (3), pp.235-242

Abstract
A hierarchically preconditioned conjugate gradient (PCG) method for finite element analysis is presented. Its use is demonstrated for the difficult problem of the non‐linear analysis of 3D reinforced concrete structures. Examples highlight the dramatic savings in computer storage and more modest savings in solution times obtained using PCG especially for large problems.
Abstract
A hierarchically preconditioned conjugate gradient (PCG) method for finite element analysis is presented. Its use is demonstrated for the difficult problem of the non‐linear [...]
- 2
- read
A crack-tracking technique for localized cohesive-frictional damage
S. Saloustros, L. Pelà, M. Cervera

Engineering Fracture Mechanics (2015). (preprint) Vol. 150, pp. 96-114

Abstract
This paper extends the use of crack-tracking techniques within the smeared crack approach for the numerical simulation of cohesive–frictional damage on quasi-brittle materials. The mechanical behaviour is described by an isotropic damage model with a Mohr–Coulomb failure surface. The correct crack propagation among the two alternative fracture planes proposed by the Mohr–Coulomb theory is selected with the use of an energy criterion based on the total elastic strain energy. The simulation of three benchmark problems of mixed-mode fracture in concrete demonstrates that the proposed methodology can reproduce the material’s frictional characteristics, showing robustness, as well as mesh-size and mesh-bias independence.
Abstract
This paper extends the use of crack-tracking techniques within the smeared crack approach for the numerical simulation of cohesive–frictional damage on quasi-brittle [...]
- 2
- read
Strong discontinuities and continuum plasticity models: the strong discontinuity approach
M. Cervera, O. Manzoli, J. Oliver

International Journal of Plasticity (1999). Vol. 15 (3), pp. 319-351

Abstract
The paper presents the Strong Discontinuity Approach for the analysis and simulation of strong discontinuities in solids using continuum plasticity models. Kinematics of weak and strong discontinuities are discussed, and a regularized kinematic state of discontinuity is proposed as a mean to model the formation of a strong discontinuity as the collapsed state of a weak discontinuity (with a characteristic bandwidth) induced by a bifurcation of the stress–strain field, which propagates in the solid domain. The analysis of the conditions to induce the bifurcation provides a critical value for the bandwidth at the onset of the weak discontinuity and the direction of propagation. Then a variable bandwidth model is proposed to characterize the transition between the weak and strong discontinuity regimes. Several aspects related to the continuum and, their associated, discrete constitutive equations, the expended power in the formation of the discontinuity and relevant computational details related to the finite element simulations are also discussed. Finally, some representative numerical simulations are shown to illustrate the proposed approach.
Abstract
The paper presents the Strong Discontinuity Approach for the analysis and simulation of strong discontinuities in solids using continuum plasticity models. Kinematics [...]
- 0
- read

editors

Cervera, Miguel 543
Universitat Politècnica de...

information

How to submit

Open access

Contact

editors

Cervera, Miguel 543

information

Categories