Mathematical Modeling of Displacements and Thermal Stresses in Anisotropic Materials (Sapphire) in Cooling
|
|
- Kory Phelps
- 6 years ago
- Views:
Transcription
1 Mathematical Modeling of Displacements and Thermal Stresses in Anisotropic Materials (Sapphire) in Cooling Timo Tiihonen & Tero Tuovinen September 11, 2015 European Study Group with Industry, ESGI 112, LUT, Sept 7-11, 2015 Problem setup by Jari Järvinen Silicom Ltd. Abstract In this study we are mathematically modeling a sapphire crystal growth by heat exchanger method (HEM). In the system, solid material is initially at temperature close to its melting point temperature, and the cooling takes place until the solid is at room temperature. Cooling is organized with decreasing heater power (heat source) in the system. We have modeled the system that is isolated thermally in every boundary except at the bottom and a cone, and there we have defined constant heat out-flow. Sapphire belongs to the hexagonal system and rhombohedra single crystals. In our analysis, the aim is to study displacements and thermal stresses for isotropic/anisotropic 3D and axi-symmetric solid material (sapphire) during the cooling? In particularly, we will answer the reliability issues related to how valid the axi-symmetric approximation compared to 3D computations. We will consider the validity of the isotropic assumption compared to the thermally and structurally anisotropic case. The numerical simulations are implemented by using Comsol Multiphysics software. 1 Introduction Sapphire crystal growth with heat exchanger method is schematized in Figure 1. The process takes place in an isolated chamber filled with inert gas. Due to high temperatures involved (over 2000C) radiativeheattransferispredominant in the gas. As our focus here is primarily in the prediction of thermally induced stresses in the crystal we will focus on modeling the crystal and assume its thermal environment known. By simple estimation it can be concluded that the role of gravitation (self weight and related forces against supporting surfaces) is 1
2 Figure 1: Sapphire Crystal Growth by Heat Exchanger Method (HEM). quite minor compared to expected level of thermal stresses. This, together with the fact that sapphire shrinks during the cooling faster than the supporting materials, enables us to ignore also the surrounding vessel from mechanical considerations. Hence the modeling can be restricted to sapphire alone. The modeled sapphire crystal is sketched in Figure 2. It has three regions: the seed crystal, the expanding cone and the cylindrical part that will be the actual produced crystal to be processed further. In what follows we briefly summarize the governing equations with focus on the material characteristic properties, sketch the simulation scenarios and summarize the main findings. Sample numerical data/graphs are included as appendix. Additional information related to the topic can be found in [Buchanan, 2000], [Chen et al., 2014], [Dobrovinskaya et al., 2009] and [Fang et al., 2013]. 2
3 2 Governing equations 2.1 Heat transfer in solids Modeling the heat transfer in solids, we have used the following representation = r (krt ) (1) where C p denotes the heat capacity (constant in our study), k the heat conductivity (likewise treated as constant and isotropic in our study) and T is temperature. We have assumed that system can be controlled outside using heaters so that boundary conditions in producing chamber can be considered thermally insulated. The equation for thermal insulation is defined as n ( krt )=0 (2) The computations have been executed for two separate boundary conditions for cooling, heat flux and constant out-flow. However, we have reported only the results from out-flow condition. The equation for heat flux boundary condition is described as follows n ( krt )=h (T ext T ) (3) and boundary conditions for constant out-flow can be represented as n ( krt )=q 0. (4) T ext is external temperature and we have assumed that to be 273K. 2.2 Solid mechanics Modeling of mechanical deformation and stress analysis, we have considered the following equation r + F v =0 (5) in the domain. F v stands for volumetric force (gravitation) that is ignored in our study. One-side coupling with thermal effects is done by using following condition = C :( T ) (6) where = 1 ru +(ru) T and T = (T T ref ). (7) 2 and denotes the (matrix of) heat expansion coefficients that is slightly anisotropic in our case. Our aim was to disturb the stresses of the system as little as possible using prescribed displacement and other constraints that are necessary for numerical convergence. In particular, we have used spring foundation approach where we 3
4 have defined very weak springs in the upper corners of the domain. conditions for spring foundation can be expressed as Line n = k A (u u 0 ) (8) where k A is spring constant per unit area. The anisotropic elastic properties are expressed in matrix notation, also called Voigt notation. To do this we take advantage of the symmetry of the stress and strain tensors and express them as six-dimensional vectors in an orthonormal coordinate system (e1, e2, e3) as [ ]= 11, 22 33, 23, 13, 12 1, 2, 3, 4, 5, 6 (9) and [ ] = 11, 22, 33, 2 23, 2 13, , 2 3, 4, 5, 6 (10) Then the stiffness tensor Ccan be expressed as 2 3 c 1111 c 1122 c 1133 c 1123 c 1131 c 1112 c 2211 c 2222 c 2233 c 2223 c 2231 c 2212 C = c 3311 c 3322 c 3333 c 3323 c 3331 c c 2311 c 2322 c 2333 c 2323 c 2331 c (11) 7 4 c 3111 c 3122 c 3133 c 3123 c 3131 c c 1211 c 1222 c 1233 c 1223 c 1231 c 1212 Sapphire is transversely isotropic material. That is, it is symmetric with respect to a rotation about an axis of symmetry. For such a material, if e 3 is the axis of symmetry, Hooke s law can be expressed as = C 11 C 12 C C 12 C 11 C C 13 C 13 C C C (C 11 C 12 ) (12) For our analysis we split C to isotropic part C iso and anisotropic C an for which C an 33 =C 33 -C 11 C an 13 =C 13 -C 12 and (13) C an 44 =C 44 -C 66. This allows convenient parametrization of the level of anisotropy as C( ) =C iso + C an. (14) 4
5 Similar anisotropy is present also in the thermal expansion that differs in e 3 -direction from the others. For Sapphire the main anisotropy relates to C13 an that is about 25% of the corresponding isotropic value. For other components of C the anisotropic effects are smaller. Transversely isotropic material allows axisymmetric treatment only if the axis of symmetry and the axis of anisotropy coincide. Otherwise a non axisymmetric effect appears even if all the loading is axisymmetric. To estimate the significance of this effect we have to consider the change in material properties expressed in the axisymmetric frame when the plane of isotropy is tilted by angle. The components of the stiffness tensor in the rotated configuration are related to the components in the reference configuration by the relation c pqrs = l pi l qj l rk l sl c ijkl (15) where l ab are the components of an orthogonal rotation matrix [L( )]. Tilting of the material can be modeled by rotation in, say yz-plane. By definition the rotation leaves isotropic properties intact and for the considered anisotropic part the leading order of changes is sin( )sin( ) for small values of. Hence for relatively small rotations the axisymmetric approximation can still be used quite safely. 2.3 Simulated scenarios We first simulate the system in axisymmetric setting. The material is initially at constant temperature (2000C) and will be cooled with constant rate either from the seed crystal boundaries only or from seed crystal and the conical surface with controlled distribution of heat fluxes simulation parameter. Other boundaries are taken as thermally insulated. For the mechanical part the system is fixed with very weak springs from the boundaries of the uppermost surface to avoid artificial fixture induced stresses in the thermally loaded zone. It should be noted that any kind of support of the seed crystal (fixing either the z-displacements or all displacements at the lowermost surface) leads to stress peaks at the fixed boundary. The following scenarios were simulated: equal cooling power to seed crystal and cone isotropic material anisotropic material 75% of cooling to seed crystal, 25 % to cone surface isotropic material anisotropic material 5
6 Figure 2: Geometry and mesh used in computations. The temperature histories were the same for all comparable cases as the temperature evolution does not depend on mechanical deformations (in our simplified model) and hence presented only once. Between the scenarios we compared the following quantities: peak stress stress distribution in the symmetry plane stress distribution on the surface distribution of total displacement All these were evaluated after 3600 s of cooling. Moreover the evolution of temperature and stress on the symmetry axis was monitored for longer time. For 3D simulations the above scenarios were repeated and in addition we monitored stress distribution in a cross section and minima and maxima of stresses on the circumference as indicator of the numerics induced loss of symmetry. Finally for both cooling scenarios the 3D simulation was done for anisotropic material that was rotated 90 degrees (changing y- and z- coordinates). 3 Numerical considerations Figure 2 presents the used geometries and defined mesh for the computations. Numerical results have been simulated using Comsol Multiphysics software and many of the predefined configurations have been used. Meshes have been build using definition finer and solver configuration have been accepted as it is by default. Computational time used in axi-symmetrical setup was around 20 s and full 3D computing takes time approximately 10 minutes. Some of the used parameters are represented in the table 1. 6
7 Property Name Value Unit Density 3980 kg/m 3 Thermal conductivity k 35 W/(m K) Coefficient of the thermal expansion E 6 1/K Heat capacity at constant pressure C p 800 J/(kg/K) Table 1: Material parameters. 4 Analysis of the results and conclusions We have presented the main results in the Table 2. In our experiments we found up to 8% difference in the stress peak between isotropic and anisotropic material models. As a summary we can conclude that anisotropic material model should be used. If the material is aligned with the symmetry axis, the simulation can be easily done with axisymmetric model. In the full 3D-model we noticed that the grid effects cause artificial asymmetry and hide the anisotropic effects partially. As a whole 3D simulation would require quite fine grid to be able to give reliable results. The results depend, naturally, quite strongly on the boundary conditions. More information of actual cooling strategies would be needed. The simulated scenarios are quite extreme in that respect concentrating the heat flux to a small part of the crystal boundary. More even distribution of heat flux is likely to produce lower peak stresses. Due to strong shrinking of sapphire even minor blocking of deformations induces significant stresses. Hence for realistic modeling the supports and contact to them must be modeled carefully (contact surfaces allowing detachment during shrinking). Situations where the material anisotropy is not aligned with symmetry axis were studied in two ways. Numerically we studied the case where crystal is turned 90 degrees. In this case clear symmetry loss was anticipated. We noticed over20% asymmetry in peak stresses and higher Von Mises stress levels. For small deviations between the material and symmetry axes we concluded with sketchy perturbation analysis that for small rotations loss of symmetry is proportional to sin 2. Hence, given the challenges in 3D simulations, it may be more reliable to use anisotropic axisymmetric model even if the positioning of the seed crystal and its orientation are not perfect. As a summary of conclusions is that material anisotropy should be taken into the model. It has 10% effect to z-displacements and may effect simulations with realistic supporting mechanisms. Moreover, it effects to peak stresses. Material anisotropy can be implemented in axisymmetric code and it is not overly sensitive to perfect alignment to symmetry axis. However, there is no immediate need for 3D model. First of all, it is costly compared to axisymmetric model is same level of reliability is aimed to. Secondly, there are no dominant 3D effects present in the system apart from growing in wrong direction or from Grasshoff convection in melt phase. 7
8 Simulation Cooling Isotropy Peak stress Displacement Axisymmetric 50 Isotropic 3,07E+07 1,16E-04 Axisymmetric 50 Anisotropic 3,05E+07 1,26E-04 Axisymmetric 75 Isotropic 4,72E+07 1,41E-04 Axisymmetric 75 Anisotropic 4,71E+07 1,52E-04 3D 50 Isotropic 3,51E+07 1,21E-04 3D 50 Anisotropic 3,25E+07 1,32E-04 3D 75 Isotropic 5,36E+07 1,47E-04 3D 75 Anisotropic 4,97E+07 1,57E-04 5 Discussion Table 2: Peak stresses and displacements. Because this was small workshop research related to event of mathematical modeling of industrial problems, it is convenient to express some notes about the our modeling process. At the first, we focused on the question of comparison between axisymmetric and full 3D computations. The starting point was 3D simulation of thermal expansion induced stresses. Using extreme scenarios, we estimate the role of the (asymmetric) anisotropy. The lack of practical knowledge of the dimensions of real system drives us to working with world record scale dimensions. Next, we build up axisymmetric model with smaller dimensions. After some analytical approximations, we succeeded to estimate semi-realistic parameters to heat flows etc. Practically the estimation relies on the information that the product took one week time to cool down. In the workshop, we had limited time resources, and therefore some details are excluded from this report. 8
9 References [Buchanan, 2000] Buchanan, G. R. (2000). Vibration of truncated conical cylinders of crystal class 6/mmm. Journal of Vibration and Control, 6(7), [Chen et al., 2014] Chen, C.-H., Chen, J.-C., Chiue, Y.-S., Chang, C.-H., Liu, C.-M., & Chen, C.-Y. (2014). Thermal and stress distributions in larger sapphire crystals during the cooling process in a kyropoulos furnace. Journal of Crystal Growth, 385, [Dobrovinskaya et al., 2009] Dobrovinskaya, E. R., Lytvynov, L. A., & Pishchik, V. (2009). Sapphire: material, manufacturing, applications. Springer Science & Business Media. [Fang et al., 2013] Fang, H., Pan, Y., Zheng, L., Zhang, Q., Wang, S., & Jin, Z. (2013). To investigate interface shape and thermal stress during sapphire single crystal growth by the cz method. Journal of Crystal Growth, 363,
10 APPENDIX Temperature Figure 3: Temperature fields using full 3D computations. Surface and cut-plane. 50/50 cooling. 10
11 Figure 4: Temperature evolution versus time. The behavior will be similar until the end. 50/50 and 75/25 cooling. 11
12 Figure 5: Temperature fields in surface and cut-off plane using full 3D computations. 25/75 cooling. 12
13 Von Mises Stress Figure 6: Isotropic and anisotropic Von Mises stress fields using axial symmetric approach. 50/50 cooling. 13
14 Figure 7: Isotropic and anisotropic Von Mises stress fields using axial symmetric approach, 3D-view. 50/50 cooling. 14
15 Figure 8: Isotropic and anisotropic Von Mises stress fields using full 3D computations. 50/50 cooling. 15
16 Figure 9: Isotropic and anisotropic Von Mises stress fields in cut-off plane using full 3D computations. 50/50 cooling. 16
17 Figure 10: Isotropic and anisotropic Von Mises stress fields using axial symmetric approach. 25/75 cooling. 17
18 Figure 11: Isotropic and anisotropic Von Mises stress fields using axial symmetric approach, 3D-view. 25/75 cooling. 18
19 Figure 12: Isotropic and anisotropic Von Mises stress fields using full 3D computations. 25/75 cooling. 19
20 Figure 13: Isotropic and anisotropic Von Mises stress fields in cut-off plane using full 3D computations. 25/75 cooling. 20
21 Figure 14: Asymmetric and symmetric anisotropic Von Mises stress fields using full 3D computations. 25/75 cooling. 21
22 Figure 15: Asymmetric and symmetric anisotropic Von Mises stress fields in cut-off plane using full 3D computations. 25/75 cooling. 22
23 Figure 16: Asymmetric anisotropicvon Mises stress fields using full 3D computations. 75/25 and 50/50 cooling. 23
24 Figure 17: Asymmetric anisotropic Von Mises stress fields in cut-off plane using full 3D computations. 75/25 and 50/50 cooling. 24
25 Figure 18: asymmetric and symmetric anisotropic Von Mises stress cut using full 3D computations. 25
26 Displacement Figure 19: Isotropic and anisotropic displacement fields using axial symmetric approach. 50/50 cooling. 26
27 Figure 20: Isotropic and anisotropic displacements using full 3D computations. 50/50 cooling. 27
28 Figure 21: Isotropic and anisotropic displacement fields using axial symmetric approach. 25/75 cooling. 28
29 Figure 22: Isotropic and anisotropic displacements using full 3D computations. 25/75 cooling. 29
30 Figure 23: Asymmetric and symmetric anisotropic displacements using full 3D computations. 25/75 cooling. 30
31 Figure 24: Asymmetric anisotropic displacements using full 3D computations. 75/25 and 50/50 cooling. 31
32 Evolutionary behavior Figure 25: Isotropic and anisotropic Von Mises stress evolution in the middle line of the system using full 3D computations. 50/50 cooling. 32
33 Figure 26: Isotropic and anisotropic Von Mises stress evolution in the middle line of the system using full 3D computations. 25/75 cooling. 33
34 Figure 27: Asymmetric and symmetric anisotropic Von Mises stress evolution in the middle line of the system using full 3D computations. 25/75 cooling. 34
35 Figure 28: Asymmetric anisotropic Von Mises stress evolution in the middle line of the system using full 3D computations. 75/25 and 50/50 cooling. 35
Introduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationUsing Thermal Boundary Conditions in SOLIDWORKS Simulation to Simulate a Press Fit Connection
Using Thermal Boundary Conditions in SOLIDWORKS Simulation to Simulate a Press Fit Connection Simulating a press fit condition in SOLIDWORKS Simulation can be very challenging when there is a large amount
More informationGeneralized Plane Piezoelectric Problem: Application to Heterostructure Nanowires
Generalized Plane Piezoelectric Problem: Application to Heterostructure Nanowires H. T. Mengistu and A. Garcia-Cristobal Institute of Material Science, University of Valencia, P.O.Box 22085, E-46071, Valencia,
More informationPeristaltic Pump. Introduction. Model Definition
Peristaltic Pump Introduction In a peristaltic pump, rotating rollers are squeezing a flexible tube. As the pushed down rollers move along the tube, the fluid in the tube follows the motion. The main advantage
More informationFinal Design Project: Biodiesel Settling Tank Analysis
MESSIAH COLLEGE ENGR 495 Finite Element Methods Tuesday, December 16, 2003 : Biodiesel Settling Tank Analysis Brandon Apple Jon Bitterman Becky Gast Kyle McNamara ENGR 495 Finite Element Methods 1 Abstract
More informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
More information3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials
3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials Abstract: Talapady S. Bhat and T. A. Venkatesh Department of Material Science and Engineering Stony Brook University,
More informationLaminated Composite Plates and Shells
Jianqiao Ye Laminated Composite Plates and Shells 3D Modelling With 62 Figures Springer Table of Contents 1. Introduction to Composite Materials 1 1.1 Introduction 1 1.2 Classification of Composite Materials
More informationA Finite Element Model for Numerical Analysis of Sintering
A Finite Element Model for Numerical Analysis of Sintering DANIELA CÂRSTEA High-School Group of Railways, Craiova ION CÂRSTEA Department of Computer Engineering and Communication University of Craiova
More informationPredicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics
Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Neelam Mishra 1, and Kaushik Das *1 1 School of Minerals Metallurgical and Materials Engineering, Indian
More informationShape, Convection and Convergence
Presented at the COMSOL Conference 2008 Boston Shape, Convection and Convergence,VP Research, Inc. Who is? Who is? We are: A Technical Consultation firm in business since 1993 Who is? We are: A Technical
More informationConstitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
More informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationA Flexible Scheme for the Numerical Homogenisation in Linear Elasticity
Presented at the COMSOL Conference 2010 Paris A Flexible Scheme for the Numerical Homogenisation in Linear Elasticity A. Gerisch 1, S. Tiburtius 1, Q. Grimal 2, K. Raum 3 1 Technische Universität Darmstadt,
More informationStress analysis of deflection analysis flexure and obif Vertical Load orientation
Stress analysis of deflection analysis flexure and obif Vertical Load orientation Note: Do not base your design decisions solely on the data presented in this report. Use this information in conjunction
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationFEA Mechanical Modeling of Torque Transfer Components for Fully Superconducting Rotating Machines
FEA Mechanical Modeling of Torque Transfer Components for Fully Superconducting Rotating Machines Tingcheng.Wu 1, Guillaume.Escamez 1, Clement.Lorin 1, Philippe J. Masson 1 1 University of Houston *Mechanical
More information3.22 Mechanical Properties of Materials Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 3.22 Mechanical Properties of Materials Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Quiz #1 Example
More informationAEROELASTIC ANALYSIS OF COMBINED CONICAL - CYLINDRICAL SHELLS
Proceedings of the 7th International Conference on Mechanics and Materials in Design Albufeira/Portugal 11-15 June 2017. Editors J.F. Silva Gomes and S.A. Meguid. Publ. INEGI/FEUP (2017) PAPER REF: 6642
More informationTable of Contents. Preface...xvii. Part 1. Level
Preface...xvii Part 1. Level 1... 1 Chapter 1. The Basics of Linear Elastic Behavior... 3 1.1. Cohesion forces... 4 1.2. The notion of stress... 6 1.2.1. Definition... 6 1.2.2. Graphical representation...
More informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationDESIGN OF LAMINATES FOR IN-PLANE LOADING
DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily
More informationarxiv: v1 [physics.ins-det] 4 Jun 2018
Noname manuscript No. (will be inserted by the editor) Simulation of force-insensitive optical cavities in cubic spacers Eugen Wiens Stephan Schiller arxiv:186.2176v1 [physics.ins-det] 4 Jun 218 Received:
More informationFinite Element Modeling for Transient Thermal- Structural Coupled Field Analysis of a Pipe Joint
International Conference on Challenges and Opportunities in Mechanical Engineering, Industrial Engineering and Management Studies 88 Finite Element Modeling for Transient Thermal- Structural Coupled Field
More informationA Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core
Commun. Theor. Phys. 56 774 778 Vol. 56, No. 4, October 5, A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core JIANG
More informationEffect of Temperature on Mechanical Properties and Dynamic Characteristics of the Rotor in High Acceleration Mei-Hui JIA a,*, Cheng-Lin WANG b
2016 International Conference on Mechanics Design, Manufacturing and Automation (MDM 2016) ISBN: 978-1-609-34-0 Effect of Temperature on Mechanical Properties and Dynamic Characteristics of the Rotor in
More informationNonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su. ME 563: Nonlinear Finite Element Analysis.
ME 563: Nonlinear Finite Element Analysis Spring 2016 Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su Department of Mechanical and Nuclear Engineering,
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More information142. Determination of reduced mass and stiffness of flexural vibrating cantilever beam
142. Determination of reduced mass and stiffness of flexural vibrating cantilever beam Tamerlan Omarov 1, Kuralay Tulegenova 2, Yerulan Bekenov 3, Gulnara Abdraimova 4, Algazy Zhauyt 5, Muslimzhan Ibadullayev
More informationStresses Analysis of Petroleum Pipe Finite Element under Internal Pressure
ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT
More informationMATERIAL REDUCTION & SYMMETRY PLANES
MATERIAL REDUCTION & SYMMETRY PLANES ISSUES IN SIMULATING A CONNECTING ROD LIMITATIONS ON THE ANALYSES Recall from previous lectures that, in static stress analysis that we are subjecting the component
More informationDEVELOPMENT OF SOPHISTICATED MATERIALS USING THE COMPUTER SIMULATION
7 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE Baia Mare, Romania, May 17-18, 2007 ISSN-1224-3264 DEVELOPMENT OF SOPHISTICATED MATERIALS USING THE COMPUTER SIMULATION Zuzana MURČINKOVÁ Ing. Zuzana MURČINKOVÁ,
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam
MODELING OF CONCRETE MATERIALS AND STRUCTURES Class Meeting #1: Fundamentals Kaspar Willam University of Colorado at Boulder Notation: Direct and indicial tensor formulations Fundamentals: Stress and Strain
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More information1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine Spring 2004 LABORATORY ASSIGNMENT NUMBER 6
1.103 CIVIL ENGINEERING MATERIALS LABORATORY (1-2-3) Dr. J.T. Germaine MIT Spring 2004 LABORATORY ASSIGNMENT NUMBER 6 COMPRESSION TESTING AND ANISOTROPY OF WOOD Purpose: Reading: During this laboratory
More informationMultiphysics Modeling
11 Multiphysics Modeling This chapter covers the use of FEMLAB for multiphysics modeling and coupled-field analyses. It first describes the various ways of building multiphysics models. Then a step-by-step
More informationBuckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate
Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate Outline Introduction Representative Volume Element (RVE) Periodic Boundary Conditions on RVE Homogenization Method Analytical
More informationTwo-dimensional ternary locally resonant phononic crystals with a comblike coating
Two-dimensional ternary locally resonant phononic crystals with a comblike coating Yan-Feng Wang, Yue-Sheng Wang,*, and Litian Wang Institute of Engineering Mechanics, Beijing Jiaotong University, Beijing,
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationMECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
More informationPrediction of Transformer Core Noise
Prediction of Transformer Core Noise R. Haettel *,1, M. Kavasoglu 1, A. Daneryd 1 and C. Ploetner 2 1 ABB Corporate Research Sweden, 2 ABB Transformers Canada *Corresponding author: 721 78 Västerås Sweden,
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
More informationStrain Measurement Techniques for Composite Coupon Testing
Strain Measurement Techniques for Composite Coupon Testing Introduction Characterization of the properties of anisotropic and inhomogeneous composite materials for use in demanding structural applications
More informationFundamental Solution of 3D Isotropic Elastic Material
Page 6 Fundamental Solution of 3D Isotropic Elastic Material Md. Zahidul Islam, Md. Mazharul Islam *, Mahesh Kumar Sah Department of Mechanical Engineering, Khulna University of Engineering & Technology,
More informationSIMULATION AND TESTING OF STRUCTURAL COMPOSITE MATERIALS EXPOSED TO FIRE DEGRADATION
SIMULATION AND TESTING OF STRUCTURAL COMPOSITE MATERIALS EXPOSED TO FIRE DEGRADATION A. Lozano Martín a*, A. Fernández López a, A. Güemes a a Department of Aerospace Materials and Manufacturing, School
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationAE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF HOOP STRESSES FOR A THIN-WALL PRESSURE VESSEL
Objective AE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF OOP STRESSES FOR A TIN-WA PRESSURE VESSE This experiment will allow you to investigate hoop and axial stress/strain relations
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationSimulation and Experimental Validation of Induction Heating of MS Tube for Elevated Temperature NDT Application.
Simulation and Experimental Validation of Induction Heating of MS Tube for Elevated Temperature NDT Application. B. Patidar, M.M.Hussain, Sanjoy Das, D Mukherjee, A.P. Tiwari Bhabha Atomic Research Centre,
More informationExample-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium
Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic
More informationDesign Optimization of an Electronic Component with an Evolutionary Algorithm Using the COMSOL-MATLAB LiveLink
Design Optimization of an Electronic Component with an Evolutionary Algorithm Using the COMSOL-MATLAB LiveLink Eva Pelster 1,David Wenger,1 1 Wenger Engineering GmbH, Einsteinstr. 55, 8977 Ulm, mail@wenger-engineering.com
More informationPressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials
Pressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials Pressure Vessels: In the previous lectures we have discussed elements subjected
More informationModule III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina
Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the
More informationSupplementary Figures
Supplementary Figures 8 6 Energy (ev 4 2 2 4 Γ M K Γ Supplementary Figure : Energy bands of antimonene along a high-symmetry path in the Brillouin zone, including spin-orbit coupling effects. Empty circles
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationThe Use of COMSOL Multiphysics for Studying the Fracture Pressure of Rectangular Micro-Channels Embedded in Thin Silicon Substrates
The Use of COMSOL Multiphysics for Studying the Fracture Pressure of Rectangular Micro-Channels Embedded in Thin Silicon Substrates Kaitlin Howell 1, H. Georgiou 2, A. Mapelli 3, P. Petagna 3, G. Romagnoli
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationNUMERICAL ANALYSIS OF HEAT STORAGE AND HEAT CONDUCTIVITY IN THE CONCRETE HOLLOW CORE DECK ELEMENT
NUMERICAL ANALYSIS OF HEAT STORAGE AND HEAT CONDUCTIVITY IN THE CONCRETE HOLLOW CORE DECK ELEMENT Michal Pomianowski 1, Per Heiselberg 1, Rasmus Lund Jensen 1, and Hicham Johra 1 1 Aalborg University,
More informationConstitutive Relations
Constitutive Relations Andri Andriyana, Ph.D. Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationMultiphysics Simulation and Optimization for Thermal Management of Electronics Systems
Multiphysics Simulation and Optimization for Thermal Management of Electronics Systems Ercan M. Dede, Jaewook Lee, & Tsuyoshi Nomura Toyota Research Institute of North America Ann Arbor, MI APEC 2012 Industry
More informationA Semi-Analytical Thermal Elastic Model for Directional Crystal Growth with Weakly Anisotropy
A Semi-Analytical Thermal Elastic Model for Directional Crystal Growth with Weakly Anisotropy School of Mathematical Sciences Peking University with Sean Bohun, Penn State University Huaxiong Huang, York
More informationModelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases
Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of
More informationReceived 23 January 2016; received in revised form 20 May 2016; accepted 21 June 2016 DOI
2084. A unified formulation for vibration analysis of open cylindrical shells coupled with annular sector plates under general boundary and coupling conditions Huimin Liu 1, Fanming Liu 2, Haoran Bai 3,
More informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationCoupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine
Coupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine Yasar Deger*, Burkhard Simperl*, Luis P. Jimenez** *Sulzer Innotec, Sulzer Markets and Technology Ltd, Winterthur, Switzerland **Guascor
More information1. Background. is usually significantly lower than it is in uniaxial tension
NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and
More informationAzimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen
Azimuthal AVO and Curvature Jesse Kolb* David Cho Kris Innanen Azimuthal AVO What is azimuthal AVO?: Analysis of incidence angle and azimuthal amplitude variations of reflection coefficients; Measures
More informationCOUPLED MAGNETO-MECHANICAL FINITE ELEMENT ANALYSIS OF A POWER TRANSFORMER IN SHORT CIRCUIT CONDITIONS
Journal of Energy VOLUME 63 2014 journal homepage: http://journalofenergy.com/ Bruno Bošnjak Siemens AG Power Transformers Nürnberg bruno.bosnjak@siemens.com Andreas Hauck SIMetris GmbH andreas.hauck@simetris.de
More informationSeismic Waves in Complex 3 D Structures, 26 (2016), (ISSN , online at
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationHeat Transfer Modeling using ANSYS FLUENT
Lecture 2 - Conduction Heat Transfer 14.5 Release Heat Transfer Modeling using ANSYS FLUENT 2013 ANSYS, Inc. March 28, 2013 1 Release 14.5 Agenda Introduction Energy equation in solids Equation solved
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationThermal and hydraulic modelling of road tunnel joints
Thermal and hydraulic modelling of road tunnel joints Cédric Hounyevou Klotoé 1, François Duhaime 1, Lotfi Guizani 1 1 Département de génie de la construction, École de technologie supérieure, Montréal,
More informationBack Calculation of Rock Mass Modulus using Finite Element Code (COMSOL)
Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL) Amirreza Ghasemi 1. Introduction Deformability is recognized as one of the most important parameters governing the behavior of rock
More informationMultiphysics Simulations for the design of a Superconducting magnet for proton therapy
WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN Paul Scherrer Institut Multiphysics Simulations for the design of a Superconducting magnet for proton therapy Ciro Calzolaio October 19, 2017 OUTLINE Superconducting
More informationAnalysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements
Acta Mech. Sin. (2016) 32(2):295 300 DOI 10.1007/s10409-015-0534-2 RESEARCH PAPER Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements V. H. Carneiro 1
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationModelling Anisotropic, Hyperelastic Materials in ABAQUS
Modelling Anisotropic, Hyperelastic Materials in ABAQUS Salvatore Federico and Walter Herzog Human Performance Laboratory, Faculty of Kinesiology, The University of Calgary 2500 University Drive NW, Calgary,
More informationChapter 6: Plastic Theory
OHP Mechanical Properties of Materials Chapter 6: Plastic Theory Prof. Wenjea J. Tseng 曾文甲 Department of Materials Engineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W. F.
More informationThermal Analysis. with SolidWorks Simulation 2013 SDC. Paul M. Kurowski. Better Textbooks. Lower Prices.
Thermal Analysis with SolidWorks Simulation 2013 Paul M. Kurowski SDC PUBLICATIONS Schroff Development Corporation Better Textbooks. Lower Prices. www.sdcpublications.com Visit the following websites to
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationExample Resistive Heating
Example Resistive Heating SOLVED WITH COMSOL MULTIPHYSICS 3.5a COPYRIGHT 2008. All right reserved. No part of this documentation may be photocopied or reproduced in any form without prior written consent
More informationMODELING THE EFFECTIVE ELASTIC MODULUS OF RC BEAMS EXPOSED TO FIRE
Journal of Marine Science and Technology, Vol., No., pp. -8 () MODELING THE EFFECTIVE ELASTIC MODULUS OF RC BEAMS EXPOSED TO FIRE Jui-Hsiang Hsu*, ***, Cherng-Shing Lin**, and Chang-Bin Huang*** Key words:
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationLoad Cell Design Using COMSOL Multiphysics
Load Cell Design Using COMSOL Multiphysics Andrei Marchidan, Tarah N. Sullivan and Joseph L. Palladino Department of Engineering, Trinity College, Hartford, CT 06106, USA joseph.palladino@trincoll.edu
More informationSolving the Inverse Problem of Resonant Ultrasound Spectroscopy on Dumbbell-shaped Compression Samples using COMSOL
Solving the Inverse Problem of Resonant Ultrasound Spectroscopy on Dumbbell-shaped Compression Samples using COMSOL Michael C. Golt U.S. Army Research Laboratory c/o Dynamic Science Inc., Aberdeen Proving
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationFrequently Asked Questions
Frequently Asked Questions Why do we have to make the assumption that plane sections plane? How about bars with non-axis symmetric cross section? The formulae derived look very similar to beam and axial
More informationCalculating Mechanical Transfer Functions with COMSOL Multiphysics. Yoichi Aso Department of Physics, University of Tokyo
Calculating Mechanical Transfer Functions with COMSOL Multiphysics Yoichi Aso Department of Physics, University of Tokyo Objective Suspension Point You have a pendulum like the one shown on the right.
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationMultiphysicss. 1 Indian Institute of. Machining. materials. In. for the anode tool. metal which leads to. form: made to study the.
Thermal Model for Single Discharge EDM Process Using COMSOL Multiphysicss K. Gajjar 1, U. Maradia 2, and K. Wegener 2 1 Indian Institute of Technology Bombay, Maharashtra, India 2 Institute for Process
More informationCoupled electromagnetic, thermal and stress analysis of large power electrical transformers
Coupled electromagnetic, thermal and stress analysis of large power electrical transformers DANIELA CÂRSTEA High-School Group of Railways, Craiova ROMANIA ALEXANDRU ADRIAN CÂRSTEA University of Craiova
More informationDesign and Analysis of Various Microcantilever Shapes for MEMS Based Sensing
ScieTech 014 Journal of Physics: Conference Series 495 (014) 01045 doi:10.1088/174-6596/495/1/01045 Design and Analysis of Various Microcantilever Shapes for MEMS Based Sensing H. F. Hawari, Y. Wahab,
More information