湿/热/力耦合条件下复合材料结构渐进损伤仿真

发展了湿/热/静力耦合条件下复合材料结构渐进损伤仿真方法,考虑复合材料力学性能的随机分布特性建立有限元模型,改进了湿热条件下复合材料本构关系,以Hashin准则、最大应力准则作为失效判据,完整仿真了不同湿热条件、不同载荷类型作用下复合材料开孔层合板从损伤起始至最终失效的损伤演化全过程,极限强度预测结果与相应试验结果的误差在10%以内,验证了仿真分析方法的有效性.湿热效应对损伤起始和演化的影响表明:由于湿热条件作用,相对于室温干态而言,复合材料开孔层合板的损伤起始的时间提前、所需载荷降低,极限强度和极限应变减小.

热/力耦合作用下基于应力分析的冰破坏准则

根据飞机热除冰的物理过程,考虑外部空气动力和蒙皮表面加热的作用,建立了NACA 0012翼型前缘冰层应力计算模型。采用有限元方法和平面三角形单元对控制方程组进行了求解,获得了外部空气动力和蒙皮表面加热对冰层黏附界面应力的影响规律。研究表明:蒙皮表面不加热时,来流速度影响了黏附界面应力的强度,来流攻角影响了黏附界面应力的分布,冰-蒙皮间黏附界面切应力最大值随来流速度呈近似线性增大趋势,但外部空气动力很难造成冰层破坏。蒙皮表面加热时,冰-蒙皮间黏附界面的耦合应力和冰层内部的主应力随着热流密度的增大而增大,很容易超过剪切强度,这是造成冰破坏的关键因素。耦合冰-蒙皮剪切强度随界面温度的变化关系,初步建立了基于应力分析和热/力耦合作用的冰破坏判断准则。外部空气动力产生的界面应力和蒙皮表面加热产生的界面热应力之和,必须大于与蒙皮表面温度相关的剪切强度,则冰层发生破坏,破坏位置是耦合应力超过剪切强度的区域。

A thermo-mechanical damage constitutive model for deep rock considering brittleness-ductility transition characteristics

This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determined through the extremum features of stress−strain curve.Subsequently,the model predictions were compared with experimental results of marble samples.It is found that when the treatment temperature rises,the coupling damage evolution curve shows an S-shape and the slope of ascending branch gradually decreases during the coupling damage evolution process.At a constant temperature,confining pressure can suppress the expansion of micro-fractures.As the confining pressure increases the rock exhibits ductility characteristics,and the shape of coupling damage curve changes from an S-shape into a quasi-parabolic shape.This model can well characterize the influence of high temperature on the mechanical properties of deep rock and its brittleness-ductility transition characteristics under confining pressure.Also,it is suitable for sandstone and granite,especially in predicting the pre-peak stage and peak stress of stress−strain curve under the coupling action of confining pressure and high temperature.The relevant results can provide a reference for further research on the constitutive relationship of rock-like materials and their engineering applications.

Thermal-hydro-mechanical coupling stress intensity factor of brittle rock

A new calculation formula of THM coupling stress intensity factor was derived by the boundary collocation method, in which an additional constant stress function was successfully introduced for the cracked specimen with hydraulic pressure applied on its crack surface. Based on the newly derived formula, THM coupling fracture modes （including tensile, shear and mixed fracture mode） can be predicted by a new fracture criterion of stress intensity factor ratio, where the maximum axial load was measured by self-designed THM coupling fracture test. SEM analyses of THM coupling fractured surface indicate that the higher the temperature and hydraulic pressure are and the lower the confining pressure is, the more easily the intergranular （tension） fracture occurs. The transgranular （shear） fracture occurs in the opposite case while the mixed-mode fracture occurs in the middle case. The tested THM coupling fracture mechanisms are in good agreement with the predicted THM coupling fracture modes, which can verify correction of the newly-derived THM coupling stress intensity factor formula.

Numerical studies of effect of tool sizes and pin shapes on friction stir welding of AA2024-T3 alloy

Coupled thermo-mechanical model was used to investigate the effects of the pin diameter, the shoulder diameter and the in conical angle on the heat generations, the material deformations and the energy histories in friction stir welding（FSW） of AA2024-T3 alloy. Results indicate that the shoulder-plate contact area takes more important contribution to the heat generation than the pin-plate contact area. The increase of the shoulder diameter or the decrease of the pin diameter can lead to the increase of the welding temperature in FSW, but the change of shoulder size is more important. Compared to the cases in FSW of AA6061-T6, the input power is obviously increased in FSW of AA2024-T3 and the ratio of the plastic dissipation to the friction dissipation becomes decreased.

Phase field simulation for non-isothermal solidification of multicomponent alloys coupled with thermodynamics database

In order to quantitively model the real solidification process of industrial multicomponent alloys, a non-isothermal phase field model was studied for multicomponent alloy fully coupled with thermodynamic and diffusion mobility database, which can accurately predict the phase equilibrium, solute diffusion coefficients, specific heat capacity and latent heat release in the whole system. The results show that these parameters are not constants and their values depend on local concentration and temperature. Quantitative simulation of solidification in multicomponent alloys is almost impossible without such parameters available. In this model, the interfacial region is assumed to be a mixture of solid and liquid with the same chemical potentials, but with different composition. The anti-trapping current is also considered in the model. And this model was successfully applied to industrial A1-Cu-Mg alloy for the free equiaxed dendrite solidification process.

Thermo-mechanical coupled particle model for rock

A thermo-mechanical coupled particle model for simulation of thermally-induced rock damage based on the particle simulation method was proposed.The simulation results of three verification examples,for which the analytical solutions are available,demonstrate the correctness and usefulness of the thermo-mechanical coupled particle model.This model is applied to simulating an application example with two cases：one is temperature-independent elastic modulus and strength,while the other is temperature-dependent elastic modulus and strength.The related simulation results demonstrate that microscopic crack initiation and propagation process with consideration of temperature-independent and temperature-dependent elastic modulus and strength are different and therefore,the corresponding macroscopic failure patterns of rock are also different.On the contrary,considering the temperature-dependent elastic modulus and strength has no or little effect on the heating conduction behavior.Numerical results,which are obtained by using the proposed model with temperature-dependent elastic modulus and strength,agree well with the experimental results.This also reveals that the rock subjected to heating experiences much more cracking than the rock subjected to cooling.

Study on inhomogeneous cooling behavior of extruded profile with unequal and large thicknesses during quenching using thermo-mechanical coupling model

The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanical coupling model for simulating the on-line water quenching of extruded profile with unequal and large thicknesses was developed. The temperature field, residual stress field and distortion of profile during quenching were investigated systematically. The results show that heat transfer coefficient increases as water flow rate increases. The peak heat transfer coefficient with higher water flow rates appears at lower interface temperatures. The temperature distribution across the cross-section of profile during quenching is severe nonuniform and the maximum temperature difference is 300 ℃ at quenching time of 3.49 s. The temperature difference through the thickness of different parts of profile first increases sharply to a maximum value, and then gradually decreases. The temperature gradient increases obviously with the increase of thickness of parts. After quenching, there exist large residual stresses on the inner side of joints of profile and the two ends of part with thickness of 10 mm. The profile presents a twisting-type distortion across the cross-section under non-uniform cooling and the maximum twisting angle during quenching is 2.78°.

Coupled thermo-hydro-mechanical-migratory model for dual-porosity medium and numerical analysis

A coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturated-unsaturated ubiquitous-joint rockmass was established,in which the stress field and the temperature field were single,but the seepage field and the concentration field were double,and the influences of sets,spaces,angles,continuity ratios,stiffnesses of fractures on the constitutive relationship of the medium were considered.Also,the relative two-dimensional program of finite element method was developed.Taking a hypothetical nuclear waste repository as a calculation example,the case in which the rockmass was unsaturated dual-porosity medium and radioactive nuclide leak was simulated numerically,and the temperatures,negative pore pressures,saturations,flow velocities,nuclide concentrations and principal stresses in the rockmass were investigated.The results show that the negative pore pressures and nuclide concentrations in the porosity and fracture present different changes and distributions.Even though the saturation degree in porosity is only about 1/10 that in fracture,the flow velocity of underground water in fracture is about three times that in porosity because the permeability coefficient of fracture is almost four orders higher than that of porosity.The value of nuclide concentration in fracture is close to that in porosity.

Influence of radial forging process on strain inhomogeneity of hollow gear shaft using finite element method and orthogonal design

Due to the current trend towards lightweight design in automotive industry,hollow stepped gear shafts for automobile and its radial forging process are widely investigated.Utilizing coupled finite element thermo-mechanical model,radial forging process of a hollow stepped gear shaft for automobile was simulated.The optimal combination of three process parameters including initial temperature,rotation rate and radial reduction was also selected using orthogonal design method.To examine the strain inhomogeneity of the forging workpiece,the strain inhomogeneity factor was introduced.The results reveal that the maximum effective strain and the minimum effective strain appeared in the outermost and innermost zones of different cross sections for the hollow stepped gear shaft,respectively.Optimal forging parameters are determined as a combination of initial temperature of 780°C,rotation rate of 21°/stroke and radial reduction of 3 mm.

Railway subgrade thermal-hydro-mechanical behavior and track irregularity under the sunny-shady slopes effect in seasonal frozen regions

The sunny-shady slopes effect is a phenomenon that impacts the temperature distribution of high-speed railway subgrades,resulting in uneven frost heaving deformation on the subgrade surface,which in turn causes static irregularity in the slab track.Based on the hydraulics theory,a thermal-hydro-mechanical(THM)coupled model of frozen soil is established and verified.We explore the process and characteristics of the temperature field and deformation of soil during the freezing process of high-speed railway subgrades and analyze the track irregularity variation law of China Railway Track SystemⅢslab tracks under uneven frost heaving deformation.The results show that,because the left and right slopes of high-speed railway subgrade are exposed to different amounts of solar radiation,which is the key factor causing uneven frost heaving of subgrade.Different strike angles cause changes in temperature of the subgrade’s upper part and the frost heaving amount on the surface,leading to differences in the deformation of the slab track structure:Increased strike angle weakens the rail level irregularity of the down line and marginally increases the rail level irregularity of the up line,and these become consistent in north-south directions.Therefore,when selecting railway lines in seasonal frozen areas,the west-east direction should be avoided to prevent the extremes in sunny-shady slopes effect on subgrades.

Coupled thermo-hydro-mechanical process in buffer material and self-healing effects with joints

Within the multi-barrier system for high-level waste disposal,the technological gap formed by combined buffer material block becomes the weak part of buffer layer.In this paper,Gaomiaozi bentonite buffer material with technological gap was studied,the heat transfer induced by liquid water flow and water vapor was embedded into the energy conservation equation.Based on the Barcelona basic model,the coupled thermo-hydro-mechanical model of unsaturated bentonite was established by analyzing the swelling process of bentonite block and the compression process of joint material.The China-Mock-up test was adopted to compare the numerical calculation results with the test results so as to verify the rationality of the proposed model.On this basis,the effect of joint self-healing on dry density,thermal conductivity and permeability coefficient of buffer material was further analyzed.The results show that,with bentonite hydrating and swelling,the joint material gradually increases in dry density,and exhibits comparatively uniform hydraulic and thermal conductivity properties as compacted bentonite block.As a result,the buffer material gradually shifts to homogenization due to the coordinated deformation.

Crack initiation rate of brittle rock under thermal-hydro-mechanical coupling condition

A calculation formula of thermal-hydro-mechanical(THM)coupling crack initiation rate for brittle rock was derived based on the energy conservation law.The self-designed THM coupling fracture test with conductive adhesive electrical measurement method was applied to measuring the THM coupling crack propagation rate of brittle rock continuously.Research results show that both calculation and test results of crack initiation rate increased with increase of the temperature and the hydraulic pressure.They are almost in good agreement,which can prove validity of the calculation formula of THM coupling crack initiation rate.

Fracture prediction in non-isothermal viscous pressure bulging of aluminum alloy sheet using ductile fracture criterion

The failure of AA3003 aluminum alloy sheet metal was predicted for non-isothermal viscous pressure bulging (VPB). Utilizing the coupled thermo-mechanical finite element method combined with ductile fracture criterion, the calculations were carried out for non-isotherm VPB at various temperatures and the influences of the initial temperature of viscous medium on failure mode of bulge specimens were investigated. The results show that the failure modes are different for the non-isothermal VPB with different initial temperatures of viscous medium. For the non-isothermal VPB of AA3003 aluminum alloy sheet with initial temperature of 250 ℃, when the initial temperature of viscous medium ranges from 150 to 180 ℃, the formability of sheet metal can be improved to a full extent. The validity of the predictions is examined by comparing with experimental results.

A new measurement method of crack propagation rate for brittle rock under THMC coupling condition

A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.

2D FEM analysis for coupled thermo-hydro-mechanical-migratory processes in near field of hypothetical nuclear waste repository

In order to consider the influence of temperature and underground water movement, an elastoplastic model and a 2D FEM stress fields on the migration of radioactive nuclide with code for analysis of coupled thermo-hydro-mechanical （THM） processes in saturated and unsaturated porous media were extended and improved through introducing the percolation and migration equation, so that the code can be used for solving the temperature field, flow field, stress field and nuclide concentration field simultaneously. The states of temperatures, pore pressures and nuclide concentrations in the near field of a hypothetical nuclear waste repository were investigated. The influence of the half life of the radioactive nuclide on the temporal change of nuclide concentration was analyzed considering the thermo-hydro-mechanical-migratory coupling. The results show that, at the boundary of the vitrified waste, the concentration of radioactive nuclide with a half life of 10 a falls after a period of rising, with the maximum value of 0.182 mol/m3 and the minimum value of 0.181 mol/m^3 at the end of computation. For a half life of 1 000 a, the concentration of radioactive nuclide always increases with the increase of the time during the computation period; and the maximum value is 1.686 mol/m^3 at the end of the computation. Therefore, under the condition of THM coupling, the concentration of radioactive nuclide with a shorter half life will decrease more quickly with water flow; but for the radioactive nuclide with a longer half life, its concentration will keep at a higher level for a longer time in the migration process.

Cast-rolling force model in solid-liquid cast-rolling bonding(SLCRB) process for fabricating bimetal clad strips

Based on twin-roll casting, a cast-rolling force model was proposed to predict the rolling force in the bimetal solid-liquid cast-rolling bonding(SLCRB) process. The solid-liquid bonding zone was assumed to be below the kiss point(KP). The deformation resistance of the liquid zone was ignored. Then, the calculation model was derived. A 2D thermal-flow coupled simulation was established to provide a basis for the parameters in the model, and then the rolling forces of the Cu/Al clad strip at different rolling speeds were calculated. Meanwhile, through measurement experiments, the accuracy of the model was verified. The influence of the rolling speed, the substrate strip thickness, and the material on the rolling force was obtained. The results indicate that the rolling force decreases with the increase of the rolling speed and increases with the increase of the thickness and thermal conductivity of the substrate strip. The rolling force is closely related to the KP height. Therefore, the formulation of reasonable process parameters to control the KP height is of great significance to the stability of cast-rolling forming.

Rheological numerical simulation for thermo-hydro-mechanical coupling analysis for rock mass

Under the environment of seepage field, stress field and temperature field interaction and influence, the three fields will not only produce coupling effect, but also have deformation with time due to the rheological behavior of rock mass. In the paper, based on the fundamental theories of rock mass coupling theory and rheological mechanics, the rheological model for fully coupled thermo-hydro-mechanical analysis for rock mass was set up, and the corresponding constitutive relationship, the conservation equation of mass and the conservation equation of energy were given, and the finite element formulas were derived for coupling analysis of rock mass. During establishing governing equations, rock mass was assumed approximately as macro-equivalent continuum medium. The obtained rheological numerical model for fully coupled thermo-hydro-mechanical analysis can be used for analyzing and predicting the long-term stability of underground caverns and slope engineering under the condition of thermo-hydro-mechanical coupling with rheological deformation.

Research on numerical welding experiment of a thick spherical shell structure

In this paper, in order to predict the residual deformation of thick spherical structure, a welding program is compiled in APDL language based on Ansys and a numerical welding experiment of a welding example is carried out. The temperature field of welding was simulated firstly, then a thermal-structure coupling analysis was carried out, and at last the residual stress and deformation after welding were got. After that, the numerical experiment result was compared with physical experiment one. The comparative analysis shows that the numerical simulation fits well with physical experiment. On the basis of that, a three-dimensional numerical experiment of a thick spherical shell structure was carried out to get the changing rule of stress and deformation of a thick spherical shell structure during welding. The research is of great value to the prediction of residual deformation and high precision machining.

Aging Characteristics of Lithium-Ion Battery Under Fast Charging Based on Electrochemical-thermalmechanical Coupling Model

The aging characteristics of lithium-ion battery(LIB)under fast charging is investigated based on an electrochemical-thermal-mechanical(ETM)coupling model.Firstly,the ETM coupling model is established by COMSOL Multiphysics.Subsequently,a long cycle test was conducted to explore the aging characteristics of LIB.Specifically,the effects of charging(C)rate and cycle number on battery aging are analyzed in terms of nonuniform distribution of solid electrolyte interface(SEI),SEI formation,thermal stability and stress characteristics.The results indicate that the increases in C rate and cycling led to an increase in the degree of nonuniform distribution of SEI,and thus a consequent increase in the capacity loss due to the SEI formation.Meanwhile,the increases in C rate and cycle number also led to an increase in the heat generation and a decrease in the heat dissipation rate of the battery,respectively,which result in a decrease in the thermal stability of the electrode materials.In addition,the von Mises stress of the positive electrode material is higher than that of the negative electrode material as the cycling proceeds,with the positive electrode material exhibiting tensile deformation and the negative electrode material exhibiting compressive deformation.The available lithium ion concentration of the positive electrode is lower than that of the negative electrode,proving that the tensile-type fracture occurring in the positive material under long cycling dominated the capacity loss process.The aforementioned studies are helpful for researchers to further explore the aging behavior of LIB under fast charging and take corresponding preventive measures.