• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.6
• /
• pp.414-421
• /
• 2019

In this paper, thermomechanical coupled analysis of carbon/phenolic composites structures in reentry environment was performed. The interface of thermomechanical coupled analysis was constructed using commercial software. The governing equations of temperature and displacement fields were considered to simulate change of physical behavior due to pyrolysis and ablation effects. The results of thermomechanical coupled analysis were compared with the results of ablation test using arc-heated wind tunnel. Also, the structural stability of reentry capsule was analyzed using the analysis interface. The excellent ablation characteristics and thermal protection effects of the carbon/phenolic composites were confirmed and the constructed analysis interface can be effectively used to perform thermal protection system design.

• Coupled systems mechanics
• /
• v.3 no.1
• /
• pp.1-25
• /
• 2014

In this work, we present the theoretical formulation, operator split solution procedure and partitioned software development for the coupled thermomechanical systems. We consider the general case with nonlinear evolution for each sub-system (either mechanical or thermal) with dedicated time integration scheme for each sub-system. We provide the condition that guarantees the stability of such an operator split solution procedure for fully nonlinear evolution of coupled thermomechanical system. We show that the proposed solution procedure can accommodate different evolution time-scale for different sub-systems, and allow for different time steps for the corresponding integration scheme. We also show that such an approach is perfectly suitable for parallel computations. Several numerical simulations are presented in order to illustrate very satisfying performance of the proposed solution procedure and confirm the theoretical speed-up of parallel computations, which follow from the adequate choice of the time step for each sub-problem. This work confirms that one can make the most appropriate selection of the time step with respect to the characteristic time-scale, carry out the separate computations for each sub-system, and then enforce the coupling to preserve the stability of the operator split computations. The software development strategy of direct linking the (existing) codes for each sub-system via Component Template Library (CTL) is shown to be perfectly suitable for the proposed approach.

• Tunnel and Underground Space
• /
• v.4 no.2
• /
• pp.102-118
• /
• 1994

Thermomechanical analysis is conducted on the radioactive repository in deep rock mass considering the in-situ stress, excavation and thermal loading of a radioactive waste. Thermomechanical properties of a discontinuous rock mass are estimated by a theoretical method so called sequential analysis. Using the estimated properties as input for finite element analysis, the influence on temperature distribution and thermal stress is analyzed within the scope of 2-dimensional steady state and transient heat transfer and coupled thermal elastic plastic behaviour. Granitic rock mass is taken for this analysis. The analysis is done for two different rock mass conditions, i.e. continuous-homogeneous and highly jointed conditions, for the purpose of comparison. In the case of steady state, the extent of disturbed zone around the storage tunnel due to the heat production of the spent-fuel canister varies depending on the thermomechanical properties of the rock mass. In the case of transient analyses, the response of the jointed rock mass to the thermal loading after radioactive waste disposal varies significantly with time, resulting in dramatic changes in the both size and location of disturbed zone.

• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.45-65
• /
• 2002

A constitutive model on oorthotropic thermo-elasto-viscoplasticity for fiber-reinforced composite materials Is illustrated, and their thermomechanical responses are predicted with the fully-coupled finite element formulation. The unmixing-mixing scheme can be adopted with the multipartite matrix method as the constitutive model. Basic assumptions based upon the composite micromechanics are postulated, and the strain components of thermal expansion due to temperature change are included In the formulation. Also. more than two sets of mechanical variables, which represent the deformation states of multipartite matrix can be introduced arbitrarily. In particular, the unmixing-mixing scheme can be used with any well-known isotropic viscoplastic theory of the matrix material. The scheme unnecessitates the complex processes for developing an orthotropic viscoplastic theory. The governing equations based on fully-coupled thermomechanics are derived with constitutive arrangement by the unmixing-mixing concept. By considering some auxiliary conditions, the Initial-boundary value problem Is completely set up. As a tool of numerical analyses, the finite element method Is used with isoparametric Interpolation fer the displacement and the temperature fields. The equation of mutton and the energy conservation equation are spatially discretized, and then the time marching techniques such as the Newmark method and the Crank-Nicolson technique are applied. To solve the ultimate nonlinear simultaneous equations, a successive iteration algorithm is constructed with subincrementing technique. As a numerical study, a series of analyses are performed with the main focus on the thermomechanical coupling effect in composite materials. The progress of viscoplastic deformation, the stress-strain relation, and the temperature History are careful1y examined when composite laminates are subjected to repeated cyclic loading.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.6 s.249
• /
• pp.705-712
• /
• 2006

A finite element method based on the penalized subdomain-interface framework is proposed for fully-coupled, nonlinear thermomechanical analyses with thermal contact anuor radiation boundaries. In the variational formulation, a well-known penalty functional scheme is adopted for connecting subdomains and interfaces that satisfy various continuity requirements. As a logical consequence, the whole domain can be arbitrarily divided into independently-modeled subdomains without considering the conformity of meshes along their interfaces. Since the nonlinearities due to the contact and radiation boundaries can be localized within a few subdomains, the computational efficiency of the present method is greatly increased with appropriate solution algorithms. By solving some numerical problems, these advantageous features are confirmed carefully.

• Structural Engineering and Mechanics
• /
• v.13 no.6
• /
• pp.591-614
• /
• 2002

This paper presents experimental and numerical analyses of the thermomechanical behaviour that takes place in SAE1020 mild steel cylindrical specimens during the conventional tensile test. A set of experiments has been carried out in order to obtain the stress-strain curve and the diameter evolution at the neck which allow, in turn, to derive the elastic and hardening parameters characterizing the material response. Temperature evolutions have also been measured for a high strain rate situation. Moreover, a finite element large strain thermoelastoplasticity-based formulation is proposed and used to simulate the deformation process during the whole test. Some important aspects of this formulation are discussed. Finally, the results provided by the simulation are experimentally validated.

• Steel and Composite Structures
• /
• v.10 no.2
• /
• pp.129-149
• /
• 2010

The objective of this study is to formulate a general 3D material-structural analysis framework for the thermomechanical behavior of steel-concrete structures in a fire environment. The proposed analysis framework consists of three sequential modeling parts: fire dynamics simulation, heat transfer analysis, and a thermomechanical stress analysis of the structure. The first modeling part consists of applying the NIST (National Institute of Standards and Technology) Fire Dynamics Simulator (FDS) where coupled CFD (Computational Fluid Dynamics) with thermodynamics are combined to realistically model the fire progression within the steel-concrete structure. The goal is to generate the spatial-temporal (ST) solution variables (temperature, heat flux) on the surfaces of the structure. The FDS-ST solutions are generated in a discrete form. Continuous FDS-ST approximations are then developed to represent the temperature or heat-flux at any given time or point within the structure. An extensive numerical study is carried out to examine the best ST approximation functions that strike a balance between accuracy and simplicity. The second modeling part consists of a finite-element (FE) transient heat analysis of the structure using the continuous FDS-ST surface variables as prescribed thermal boundary conditions. The third modeling part is a thermomechanical FE structural analysis using both nonlinear material and geometry. The temperature history from the second modeling part is used at all nodal points. The ABAQUS (2003) FE code is used with external user subroutines for the second and third simulation parts in order to describe the specific heat temperature nonlinear dependency that drastically affects the transient thermal solution especially for concrete materials. User subroutines are also developed to apply the continuous FDS-ST surface nodal boundary conditions in the transient heat FE analysis. The proposed modeling framework is applied to predict the temperature and deflection of the well-documented third Cardington fire test.

• Journal of Powder Materials
• /
• v.4 no.2
• /
• pp.83-89
• /
• 1997

A finite element analysis to solve the coupled thermomechanical problem in the plane strain upsetting of the porous metals was performed. The analysis was formulated using the yield function advanced by Lee and kim and developed using the thermo-elasto-plastic time integration procedure. The density and temperature dependent thermal and mechanical properties of porous metals were considered. The internal heat generation by the plastic deformation and the changing thermal boundary conditions corresponding to the geometry were incorporated in the program. The distributions of the stress, strain, pressure, density and temperature were predicted during the free resting period, deformation period and dwelling period of the forging process.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.11
• /
• pp.7-14
• /
• 2005

A subdomain-interface finite element method is suggested to solve a class of fully- coupled thermomechanical problems with contact boundaries. The penalty method is used for connecting subdomains that satisfy interface compatibility conditions. As a result, effective stiffness matrices are always positive definite, and computational efficiency can be improved to a considerable degree. Moreover, any complex-shaped domain can be divided into independently modeled subdomains without considering the conformity of meshes on interfaces. Using a computer code based on the present method, these advantageous features are shown through a set of numerical studies.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1995.10a
• /
• pp.16-29
• /
• 1995

A full three-dimensional thermo-coupled rigid-viscoplastic finite element method and the currently developed microstructural evolution system which includes semi-empirical mathematical equations suggested by different research groups were used together to form an integrated system of process and microstructure simulation of hot rolling. The distribution and time history of thermomechanical variables such as temperature, strain, strain rate, and time during pass and between passes were obtained FEM analysis of multipass hot rolling processes. Then distribution of metallurgical variables were calculated successfully on the basis of instantaneous thermomechanical data. For the verification of this method the evolution of microstructure in plate rolling and shape rolling was simulated and their results were compared with the data available in literature. Consequently, this approach makes it passible to describe the realistic evolution of microstructure by avoiding the use of erroneous average value and can be used in CAE of multipass hot rolling.