Based on the concept of constitutive relation error along with the residual of both origin and dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed in this paper. It lea...Based on the concept of constitutive relation error along with the residual of both origin and dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed in this paper. It leads to high quality local error bounds in the problem of fracture mechanics simulation with extended finite element method (XFEM), which involves enrichment to solve a stress singularity in the crack. Since goal-oriented error estimation with enriched degrees of freedom gives us a chance to evaluate the XFEM simulation, the stress intensity factor calculated by two kinds of XFEM programs developed by ourselves and by commercial code ABAQUS are compared in this work. By comparing the reliability of the stress intensity factor calculation, the accuracy of two programs in different cases is evaluated and the source of error is discussed. A 2-dimensional XFEM example is given to illustrate the computational procedure.展开更多
Using a recently established liquid crystal model for vesicles, we present a theoretical method to analyze the morphological stability of liquid crystal vesicles in an electric field. The coupled mechanical-electrical...Using a recently established liquid crystal model for vesicles, we present a theoretical method to analyze the morphological stability of liquid crystal vesicles in an electric field. The coupled mechanical-electrical effects associated with elastic bending, osmotic pressure, surface tension, Max- well pressure, as well as flexoelectric and dielectric proper- ties of the membrane are taken into account. The first and second variations of the free energy are derived in a com- pact form by virtue of the surface variational principle. The former leads to the shape equation of a vesicle embedded in an electric field, and the latter allows us to examine the stabil- ity of a given vesicle morphology. As an illustrative exam- ple, we analyze the stability of a spherical vesicle under a uniform electric field. This study is helpful for understanding and revealing the morphological evolution mechanisms of vesicles in electric fields and some associated phenomena of cells.展开更多
We develop a new hierarchical dislocation-grain boundary (GB) interaction model to predict the mechanical behavior of polycrystalline metals at micro and submicro scales by coupling 3D Discrete Dislocation Dynamics (D...We develop a new hierarchical dislocation-grain boundary (GB) interaction model to predict the mechanical behavior of polycrystalline metals at micro and submicro scales by coupling 3D Discrete Dislocation Dynamics (DDD) simulation with the Molecular Dynamics (MD) simulation.At the microscales,the DDD simulations are responsible for capturing the evolution of dislocation structures;at the nanoscales,the MD simulations are responsible for obtaining the GB energy and ISF energy which are then transferred hierarchically to the DDD level.In the present model,four kinds of dislocation-GB interactions,i.e.transmission,absorption,re-emission and reflection,are all considered.By this methodology,the compression of a Cu micro-sized bi-crystal pillar is studied.We investigate the characteristic mechanical behavior of the bi-crystal compared with that of the single-crystal.Moreover,the comparison between the present penetrable model of GB and the conventional impenetrable model also shows the accuracy and efficiency of the present model.展开更多
Hydraulic fracture(HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method(XFEM) coupling with Biot...Hydraulic fracture(HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method(XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures(NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk(KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.展开更多
In this study,we theoretically investigate the dynamic indentation for measuring the loss (damping) factor of a linear viscoelastic material from its indentation response.A rigid indenter with arbitrary tip profile is...In this study,we theoretically investigate the dynamic indentation for measuring the loss (damping) factor of a linear viscoelastic material from its indentation response.A rigid indenter with arbitrary tip profile is assumed to indent into a viscoelastic substrate with arbitrary shape.We perform a theoretical analysis and identify the conditions under which the loss factor of the material can be determined from the phase angle between the applied harmonic indentation load and the corresponding harmonic displacement,a directly measurable quantity in a dynamic indentation test.To validate the conclusion drawn from our theoretical analysis,a series of numerical experiments are performed,including the spherical indentation of a soft layer with irregular surface morphology bonded to a rigid substrate,a conical indenter with tip defects indenting into a half-spherical particle,and the indentation of porous materials.This study may facilitate the use of the dynamic indentation technique to evaluate the damping properties of linear viscoelastic materials,including some advanced polymers and biological soft tissues.展开更多
In order to analyze and simulate the impact collapse or seismic response of the reinforced concrete(RC)structures,a combined fiber beam model is proposed by dividing the cross section of RC beam into concrete fiber an...In order to analyze and simulate the impact collapse or seismic response of the reinforced concrete(RC)structures,a combined fiber beam model is proposed by dividing the cross section of RC beam into concrete fiber and steel fiber.The stress-strain relationship of concrete fiber is based on a model proposed by concrete codes for concrete structures.The stress-strain behavior of steel fiber is based on a model suggested by others.These constitutive models are implemented into a general finite element program ABAQUS through the user defined subroutines to provide effective computational tools for the inelastic analysis of RC frame structures.The fiber model proposed in this paper is validated by comparing with experiment data of the RC column under cyclical lateral loading.The damage evolution of a three-dimension frame subjected to impact loading is also investigated.展开更多
Most of the elements in structural design are subjected to the complex loading conditions that need to be predicted in advance for a reliable design. The failure prediction of these structural elements is based on the...Most of the elements in structural design are subjected to the complex loading conditions that need to be predicted in advance for a reliable design. The failure prediction of these structural elements is based on the failure theories that employ the experimental data of the uniaxial, the biaxial and the shear tests. Among these tests, the biaxial testing is used to model the behavior of the structure loaded in more than one direction. In this study, we discover an interesting fact for the biaxial tension specimen that the thickness reduction of the measuring region cannot increase the stress to a high level. This behavior of the biaxial testing is different from the uniaxial testing where the thickness reduction of the measuring region always increases its stress to a higher level. Therefore, we further investigate the biaxial testing specimen via related fundamental problem, i.e., an axisymmetric plate with varying thickness under axisymmetric tension. We perform the shape optimization based on our analytical solution and the FEM results. For an infinite axisymmetric disk, the optimized design should include a low slope diverging section connecting the inner to outer region. On the contrary, the larger inner measuring region for finite axisymmetric disk would result in higher corresponding stress. These conclusions can be used as a guideline for the design of a biaxial testing specimen.展开更多
In this paper,a crystal plasticity model considering the irradiation effect based on the thermal activation theory is established.The evolutions of screw dislocations,edge dislocations,and stacking fault tetrahedrals(...In this paper,a crystal plasticity model considering the irradiation effect based on the thermal activation theory is established.The evolutions of screw dislocations,edge dislocations,and stacking fault tetrahedrals(SFTs)(induced by irradiation)are included into the model.The interactions between dislocations and irradiation-induced SFTs are also considered.The constitutive model is numerically implemented on the ABAQUS platform through UMAT subroutine and applied to study the irradiation effect on the mechanical behavior of pure copper.The mechanical properties of single and polycrystalline copper are studied,and the simulation results show that the constitutive model can properly predict the mechanical behavior of irradiated pure copper.Especially for polycrystalline copper,the simulation results are in good agreement with the experimental data.展开更多
基金Project supported by the National Natural Science Foundation of China(No.10876100)
文摘Based on the concept of constitutive relation error along with the residual of both origin and dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed in this paper. It leads to high quality local error bounds in the problem of fracture mechanics simulation with extended finite element method (XFEM), which involves enrichment to solve a stress singularity in the crack. Since goal-oriented error estimation with enriched degrees of freedom gives us a chance to evaluate the XFEM simulation, the stress intensity factor calculated by two kinds of XFEM programs developed by ourselves and by commercial code ABAQUS are compared in this work. By comparing the reliability of the stress intensity factor calculation, the accuracy of two programs in different cases is evaluated and the source of error is discussed. A 2-dimensional XFEM example is given to illustrate the computational procedure.
基金supported by the National Natural Science Foundation of China(10972121,10732050 and 10525210)973 Program(2010CB631005)
文摘Using a recently established liquid crystal model for vesicles, we present a theoretical method to analyze the morphological stability of liquid crystal vesicles in an electric field. The coupled mechanical-electrical effects associated with elastic bending, osmotic pressure, surface tension, Max- well pressure, as well as flexoelectric and dielectric proper- ties of the membrane are taken into account. The first and second variations of the free energy are derived in a com- pact form by virtue of the surface variational principle. The former leads to the shape equation of a vesicle embedded in an electric field, and the latter allows us to examine the stabil- ity of a given vesicle morphology. As an illustrative exam- ple, we analyze the stability of a spherical vesicle under a uniform electric field. This study is helpful for understanding and revealing the morphological evolution mechanisms of vesicles in electric fields and some associated phenomena of cells.
基金the National Natural Science Foundation of China(Grant Nos.12172194,12222205,and 51971115)the National Key Research and Development Program of China(Grant No.2019YFE03130003).
基金supported by the National Natural Science Foundation of China(Grant No. 10772096)the National Basic Research Program of China(Grand No. 2010CB631005)
文摘We develop a new hierarchical dislocation-grain boundary (GB) interaction model to predict the mechanical behavior of polycrystalline metals at micro and submicro scales by coupling 3D Discrete Dislocation Dynamics (DDD) simulation with the Molecular Dynamics (MD) simulation.At the microscales,the DDD simulations are responsible for capturing the evolution of dislocation structures;at the nanoscales,the MD simulations are responsible for obtaining the GB energy and ISF energy which are then transferred hierarchically to the DDD level.In the present model,four kinds of dislocation-GB interactions,i.e.transmission,absorption,re-emission and reflection,are all considered.By this methodology,the compression of a Cu micro-sized bi-crystal pillar is studied.We investigate the characteristic mechanical behavior of the bi-crystal compared with that of the single-crystal.Moreover,the comparison between the present penetrable model of GB and the conventional impenetrable model also shows the accuracy and efficiency of the present model.
基金supported by the National Natural Science Foundation of China(Grant Nos.11532008,and 11372157)
文摘Hydraulic fracture(HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method(XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures(NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk(KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.
基金supported by the National Natural Science Foundation of China(Grant Nos. 10972112,10525210,and 10732050)the National Basic Research Program of China(Grant Nos. 2010CB631005)
文摘In this study,we theoretically investigate the dynamic indentation for measuring the loss (damping) factor of a linear viscoelastic material from its indentation response.A rigid indenter with arbitrary tip profile is assumed to indent into a viscoelastic substrate with arbitrary shape.We perform a theoretical analysis and identify the conditions under which the loss factor of the material can be determined from the phase angle between the applied harmonic indentation load and the corresponding harmonic displacement,a directly measurable quantity in a dynamic indentation test.To validate the conclusion drawn from our theoretical analysis,a series of numerical experiments are performed,including the spherical indentation of a soft layer with irregular surface morphology bonded to a rigid substrate,a conical indenter with tip defects indenting into a half-spherical particle,and the indentation of porous materials.This study may facilitate the use of the dynamic indentation technique to evaluate the damping properties of linear viscoelastic materials,including some advanced polymers and biological soft tissues.
基金supported by the National Natural Science Foundation of China(Grant No.90815026)
文摘In order to analyze and simulate the impact collapse or seismic response of the reinforced concrete(RC)structures,a combined fiber beam model is proposed by dividing the cross section of RC beam into concrete fiber and steel fiber.The stress-strain relationship of concrete fiber is based on a model proposed by concrete codes for concrete structures.The stress-strain behavior of steel fiber is based on a model suggested by others.These constitutive models are implemented into a general finite element program ABAQUS through the user defined subroutines to provide effective computational tools for the inelastic analysis of RC frame structures.The fiber model proposed in this paper is validated by comparing with experiment data of the RC column under cyclical lateral loading.The damage evolution of a three-dimension frame subjected to impact loading is also investigated.
基金supported by the National Natural Science Foundation of China(Grant Nos.11425208,11372158,11720101002,11890674)the Science Challenge Project(Grant No.TZ2018001)
文摘Most of the elements in structural design are subjected to the complex loading conditions that need to be predicted in advance for a reliable design. The failure prediction of these structural elements is based on the failure theories that employ the experimental data of the uniaxial, the biaxial and the shear tests. Among these tests, the biaxial testing is used to model the behavior of the structure loaded in more than one direction. In this study, we discover an interesting fact for the biaxial tension specimen that the thickness reduction of the measuring region cannot increase the stress to a high level. This behavior of the biaxial testing is different from the uniaxial testing where the thickness reduction of the measuring region always increases its stress to a higher level. Therefore, we further investigate the biaxial testing specimen via related fundamental problem, i.e., an axisymmetric plate with varying thickness under axisymmetric tension. We perform the shape optimization based on our analytical solution and the FEM results. For an infinite axisymmetric disk, the optimized design should include a low slope diverging section connecting the inner to outer region. On the contrary, the larger inner measuring region for finite axisymmetric disk would result in higher corresponding stress. These conclusions can be used as a guideline for the design of a biaxial testing specimen.
基金The support of the National Natural Science Foundation of China(NSFC)under Grant No.11202114Beijing Higher Education Young Elite Teacher Project under Grant No.YETP0156Tsinghua University Initiative Scientific Research Program under Grant No.2019Z08QCX06 are gratefully acknowledged。
文摘In this paper,a crystal plasticity model considering the irradiation effect based on the thermal activation theory is established.The evolutions of screw dislocations,edge dislocations,and stacking fault tetrahedrals(SFTs)(induced by irradiation)are included into the model.The interactions between dislocations and irradiation-induced SFTs are also considered.The constitutive model is numerically implemented on the ABAQUS platform through UMAT subroutine and applied to study the irradiation effect on the mechanical behavior of pure copper.The mechanical properties of single and polycrystalline copper are studied,and the simulation results show that the constitutive model can properly predict the mechanical behavior of irradiated pure copper.Especially for polycrystalline copper,the simulation results are in good agreement with the experimental data.