Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at diffe...Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.展开更多
The failure behavior of the precast prestressed concrete T girder was investigated by destructive test and finite-element analysis,and the mid-span deflection,girder stiffness and the variation of the cross section st...The failure behavior of the precast prestressed concrete T girder was investigated by destructive test and finite-element analysis,and the mid-span deflection,girder stiffness and the variation of the cross section strain in the loading process were obtained,and the mechanical properties,mechanical behavior,elastic and plastic behavior and ultimate bearing capacity of T girder with large span were revealed.Furthermore,the relationship between the beam stiffness degradation,the neutral axis in cross-section,steel yielding and concrete cracking are investigated and analyzed.A method was proposed to predict the residual bearing capacity of a bridge based on the variation of the position of the cross section strain distribution and the section neutral axis,which provided a theoretical basis for predicting the stiffness detection and carrying capacity assessment of prestressed concrete beam.展开更多
With a geometrical model of porous material, a 3D finite-element analysis on the rolling process of spring steel60Si2Mn in the semi-solid state is carried out using software MARC. In terms of flat and groove rolling c...With a geometrical model of porous material, a 3D finite-element analysis on the rolling process of spring steel60Si2Mn in the semi-solid state is carried out using software MARC. In terms of flat and groove rolling conditions,stress field and strain field are studied. The simulation results show that the rigid-viscoplastic model can accuratelydescribe the semi-solid metal rolling process. Semi-solid slurry has the characteristics of low flow stress and goodfluidity. During groove rolling, distribution of stress and strain on the cross-section of deformation zone is moreuniform than that during flat rolling. The results of simulation are in good agreement with the experiment data, andshow that semi-solid material fits for groove rolling.展开更多
The finite-element modeling and simulations of the intra-body communication(IBC)were investigated to provide a theoretical basis for biomedical monitoring.A finite-element model for the whole human body was developed ...The finite-element modeling and simulations of the intra-body communication(IBC)were investigated to provide a theoretical basis for biomedical monitoring.A finite-element model for the whole human body was developed to simulate the IBC.The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths,and their attenuations were calculated.Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC,while the potential generally decreased along the axis of the body parts.Both signal attenuations of the two types IBC increased with increasing signal transmission distance,and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC.The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.展开更多
In this paper,we introduce new stable mixed finite elements of any order on polytopal mesh for solving second-order elliptic problem.We establish optimal order error estimates for velocity and super convergence for pr...In this paper,we introduce new stable mixed finite elements of any order on polytopal mesh for solving second-order elliptic problem.We establish optimal order error estimates for velocity and super convergence for pressure.Numerical experiments are conducted for our mixed elements of different orders on 2D and 3D spaces that confirm the theory.展开更多
In a full-arch implant rehabilitation ad modum Branemark, the distribution of stress and strain in mandibular bone is influenced by the type, number and position of implants used. In particular, the biomechanical beha...In a full-arch implant rehabilitation ad modum Branemark, the distribution of stress and strain in mandibular bone is influenced by the type, number and position of implants used. In particular, the biomechanical behaviour of the bone structure after complete osseointegration depends on the load transferred to the bone by each fixture. In this study, a finite-element analysis of two models was performed. Models of an all-on-four configuration and a six-implant configuration were compared in a worst-case scenario. A new V parameter is presented to aid the quantitative and comparative analysis of the all-on-four and six-implant configurations. The influence of orthotropy was also investigated, and a geometric change in the all-on-four configuration is presented.展开更多
A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series ofcorresponding experiments were made. The result...A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series ofcorresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of tem- perature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESE casting and structure of casting is fine and compact.展开更多
The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a sp...The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.展开更多
A two dimensional plane-stress finite-element type of analysis is pre-sented to predict the behaviour of geogrids embedded in sand under pullout loadingconditions.In the analysis the interactions between soil and geog...A two dimensional plane-stress finite-element type of analysis is pre-sented to predict the behaviour of geogrids embedded in sand under pullout loadingconditions.In the analysis the interactions between soil and geogrid are simulatedby non-linear springs.The stiffnesses of the springs can be determined from simpletests in a specially designed pullout box.The proposed finite element (FE) analysis isapplied to interpret test results from a large scale pullout test rig.The predicted be-haviour of the geogrid under pullout load agrees well with the observed data includingthe load-displacement properties,the displacement distribution along the longitudinaldirection and the mobilisation of the frictional and bearing resistance.展开更多
Background:The reverse obliquity inter-trochanteric fracture is a distinct fracture pattern that is mechanically different from most inter-trochanteric fractures and the optional treatment of it is still controversial...Background:The reverse obliquity inter-trochanteric fracture is a distinct fracture pattern that is mechanically different from most inter-trochanteric fractures and the optional treatment of it is still controversial.The purpose of this study was to compare differences in the efficacy of a novel nail(medial support nail[MSN-II])and proximal femoral nail anti-rotation(PFNA-II)in the treatment of reverse obliquity inter-trochanteric fractures(Arbeitsgemeinschaft fur Osteosynthesfrogen/Orthopedic Trauma Association[AO/OTA]31-A3.1)using finite-element analysis.Methods:Modeling software was used to establish a three-dimensional model of MSN-II and PFNA-II and an A3.1 inter-trochanteric fracture model.Abaqus software was used to implement different force loads to compare finite-element biomechanical parameters such as the maximum stress in implant and the displacement of fracture site.Results:The femoral stress,implant stress and fracture site displacement of MSN-II was less than that of PFNA-II.The results indicated that the maximal femoral stress was 581 MPa for PFNA-II and 443 MPa for the MSN-II.The maximum stress values in the PFNA-II and MSN-II models were 291 and 241 MPa,respectively.The maximal displacements of the fracture site were 1.47 and 1.16 mm in the PFNA-II and MSN-II models,respectively.Conclusions:Compared with PFNA-II for inter-trochanteric fracture(AO/OTA 31-A3.1),MSN-II which was designed with a triangular stability structure can provide better biomechanical stability.The MSN-II may be a feasible option for the treatment of reverse obliquity inter-trochanteric fracture.展开更多
Under the same conditions of external force, simulations on differences of deformational energy in structural zones, which have different deformational behavior, and that on the distribution of the differences have be...Under the same conditions of external force, simulations on differences of deformational energy in structural zones, which have different deformational behavior, and that on the distribution of the differences have been carried out by means of finite-element method. Shear deformational energy U_w is higher than volume deformational energy U_T by about one order of magnitude, and deformationat energy U_B( = U_w + U_T) and U_w show a trend to be larger→largest→small in quantity in structural zones of different deformational properties, which correspond to compressive, shear, tensile zones, respectively, but U_T shows a trend to be large→small gradually. This has a considerable significance in the study of tectonic heat and mineral liquid migration in association with the research on tectonic additional hydrostatic pressure.展开更多
The accumulated large amount of satellite magnetic data strengthens our capability of resolving the electrical conductivity of Earth’s mantle.To invert these satellite magnetic data,accurate and efficient forward mod...The accumulated large amount of satellite magnetic data strengthens our capability of resolving the electrical conductivity of Earth’s mantle.To invert these satellite magnetic data,accurate and efficient forward modeling solvers are needed.In this study,a new finite-element based forward modeling solver is developed to accurately and efficiently compute the induced electromagnetic field for a realistic 3D Earth.Firstly,the nodal-based finite element method with linear shape function on tetrahedral grid is used to assemble the final system of linear equations for the magnetic vector potential and electric scalar potential.The FGMRES solver with algebraic multigrid(AMG)preconditioner is used to quickly solve the final system of linear equations.The weighted moving least-square method is employed to accurately recover the electromagnetic field from the numerical solutions of magnetic vector and electric scalar potentials.Furthermore,a local mesh refinement technique is employed to improve the accuracy of the estimated electromagnetic field.At the end,two synthetic models are used to verify the accuracy and efficiency of our newly developed forward modeling solver.A realistic 3D Earth model is used to simulate the induced magnetic field at 450 and 200 km altitudes which are the planned flying altitudes of Macao’s geomagnetic satellites.The simulation indicates that(1)the amplitude of the mantle-induced magnetic field can reach 10–30 nT at 450 km altitude,which is 10–30%of the primary magnetic field.The induced magnetic field at 200 km altitude has larger amplitudes.These mantleinduced magnetic fields can be measured by Macao geomagnetic satellites;(2)the amplitude of the ocean-induced magnetic field can reach 5–30 nT at satellite altitudes,which needs to be carefully considered in the interpretation of satellite magnetic data.We are confident that our newly developed forward modeling solver will become a key tool for interpreting satellite magnetic data.展开更多
Accurate simulations of metamaterial devices are very important in the analysis of their electromagnetic properties.However,it is very difficult to make full-wave simulations of three-dimensional(3D)metamaterial devic...Accurate simulations of metamaterial devices are very important in the analysis of their electromagnetic properties.However,it is very difficult to make full-wave simulations of three-dimensional(3D)metamaterial devices due to the huge memory requirements and long computing time.In this paper,we present an efficient finiteelement method(FEM)to analyze 3D axisymmetric electromagnetic devices designed by the transformation-optics approach,such as invisibility cloaks and concentrators.In the proposed method,we use the edge-based vector basis functions to expand the transverse field components,and the node-based scalar basis functions to expand the angular component.The FEM mesh is truncated with a cylindrical perfectly matched layer.We have applied the method to investigate the scattering from spherical and ellipsoidal invisibility cloaks and circularly cylindrical concentrators,in which the permittivity and permeability are both inhomogeneous and anisotropic.Numerical results are presented to show the validity and efficiency of the method.展开更多
This paper presents the design,fabrication,and characterization of cantilever-type resonators with a novel stacked structure.Aluminum nitride is adopted as the material for both the structural layer and the piezoelect...This paper presents the design,fabrication,and characterization of cantilever-type resonators with a novel stacked structure.Aluminum nitride is adopted as the material for both the structural layer and the piezoelectric layer;this simplifies the fabrication process and improves the quality factor of the resonator.Both in-plane and out-of-planeflexural modes were investigated.The effect of the structural dimensions and electrode patterns on the resonator’s performance were also studied.Finite-element simulations and experiments examining anchor loss and thermoelastic damping,which are the main loss mechanisms affecting the quality factor of these resonators,were carried out.The optimal structural dimensions and electrode patterns of the cantilever-type resonators are presented.A quality factor of 7922 with a motional impedance of 88.52 kΩand a quality factor of 8851 with a motional impedance of 67.03 kΩwere achieved for the in-plane and out-of-planeflexural-mode resonators,respectively.The proposed resonator design will contribute to the development of high-performance devices such as accelerometers,gyroscopes,and pressure sensors.展开更多
Soils are not necessarily uniform and may present linearly varied or layered characteristics,for example the backfilled soils behind rigid retaining walls.In the presence of large lateral thrust imposed by arch bridge...Soils are not necessarily uniform and may present linearly varied or layered characteristics,for example the backfilled soils behind rigid retaining walls.In the presence of large lateral thrust imposed by arch bridge,passive soil failure is possible.A reliable prediction of passive earth pressure for the design of such wall is challenging in complicated soil strata,when adopting the conventional limit analysis method.In order to overcome the challenge for generating a kinematically admissible velocity field and a statically allowable stress field,finite element method is incorporated into limit analysis,forming finiteelement upper-bound(FEUB)and finite-element lower-bound(FELB)methods.Pseudo-static,original and modified pseudo-dynamic approaches are adopted to represent seismic acceleration inputs.After generating feasible velocity and stress fields within discretized elements based on specific criteria,FEUB and FELB formulations of seismic passive earth pressure(coefficient K_(P))can be derived from work rate balance equation and stress equilibrium.Resorting to an interior point algorithm,optimal upper and lower bound solutions are obtained.The proposed FEUB and FELB procedures are well validated by limit equilibrium as well as lower-bound and kinematic analyses.Parametric studies are carried out to investigate the effects of influential factors on seismic K_(P).Notably,true solution of K_(P) is well estimated based on less than 5%difference between FEUB and FELB solutions under such complex scenarios.展开更多
Regarding quality inspection of technologically important nanocomposite hard coatings based on Ti,B,Si,C,and N and bioceramics such as hydroxyapatite that are used in small-scale high-precision devices and bio-implant...Regarding quality inspection of technologically important nanocomposite hard coatings based on Ti,B,Si,C,and N and bioceramics such as hydroxyapatite that are used in small-scale high-precision devices and bio-implants,it is essential to study the failure mechanisms associated with nanoindentation,such as fracture,delamination,and chipping.The stress imposed by the indenter can affect the fracture morphology and the interfacial fracture energy,depending on indenter shape,substrate type,crystallographic properties,pre-existing flaws,internal microcracks,and pre-strain.Reported here are finite-element-based fracture studies that provide insights into the different cracking mechanisms related to the aforementioned failure process,showing that the fracture morphology is affected by the interaction of different cracking events.The interfacial fracture energy,toughness,and residual stress are calculated using existing models with minor adjustments,and it is found that increasing the indenter sharpness improves the shear stress distribution,making the coating more prone to separation.Depending on the prevailing type of stress,the stress distribution beneath the depression results in either crack formation or a dislocation pile-up leading to strain hardening.Different forms of resistances resulting from the indentation process are found to affect the tip–sample conduction,and because of its stronger induced plasticity than that of a Berkovich indenter tip,a sharper cube-corner tip produces more resistance.展开更多
Better torque performance and higher reliability have long been the focus of research for slotted limited-angle torque motors(LATMs),which are primarily used to position first-stage valves in electrohydraulic servosys...Better torque performance and higher reliability have long been the focus of research for slotted limited-angle torque motors(LATMs),which are primarily used to position first-stage valves in electrohydraulic servosystems.This paper presents a high reliability axial-flux slotted LATM with quasi-Halbach array for torque performance improvement including constant torque range(CTR)and output torque.Firstly,the structure with two sets of windings and the operation principle of the proposed slotted LATM is analyzed.Secondly,a brief design procedure is presented,the structure selections of open slot and double-stator single-rotor(DSSR)interior rotor with surface mounted quasi-Halbach permanent magnet(PM)array are illustrated,and the geometric parameters are optimized to obtain the optimal design of the proposed slotted LATM.Then,3-D finite-element method(FEM)is employed to compare the proposed slotted LATM with the conventional surface mounted PM slotted LATM in terms of cogging torque,no-load back EMF,and output torque,and the results show that the proposed LATM with quasi-Halbach array has a 10%improvement in output torque and a 25%improvement in CTR.Meanwhile,the flux linkages and torque performance of the two sets of windings under various conditions verify good magnetic isolation.Finally,prototypes of two different rotor types are manufactured and a series of experiments are performed to validate the analysis.展开更多
As members of doubly salient magnetless linear machines,linear variable flux reluctance(LVFR)and wound field flux reversal(LWFFR)machines inherit the merits of conventional magnetless linear machines such as low cost,...As members of doubly salient magnetless linear machines,linear variable flux reluctance(LVFR)and wound field flux reversal(LWFFR)machines inherit the merits of conventional magnetless linear machines such as low cost,high flux adjustment capability and high reliability.Furthermore,like linear switched reluctance machine,they have a very simple and compact long secondary,which are very attractive for long stroke applications.However,low force capability is their major defect.To solve this issue,new LVFR and LWFFR machine topologies were proposed in recent work,while lacking studies on their force improvement mechanism and further force evaluation.In this paper,LVFR and LWFFR machines with improved force performance are comparatively studied with the emphasis on their force capabilities.The operation principle of the two machines is analyzed based on magnetic field harmonics produced by flux modulation.Contributions of air-gap flux density harmonic components to no-load back electromagnetic forces of the two machines are analyzed and the average force equation is derived.Moreover,force capabilities of the both machines are investigated by means of the time-stepping finite-element analysis to verify the theoretical analysis.展开更多
The hidden water-bearing structures near the roadway tunnelling face are very likely to cause water seepage accidents in coal mines.Currently,transient electromagnetic(EM)technology has be-come an important method to ...The hidden water-bearing structures near the roadway tunnelling face are very likely to cause water seepage accidents in coal mines.Currently,transient electromagnetic(EM)technology has be-come an important method to detect water damage in advance of roadway excavation.In this paper,the time-domain finite element algorithm based on unstructured tetrahedron grids is used to accurate-ly simulate the geological body in front of the roadway excavation face and analyze its response.The authors detect the distance between the roadway excavation face and the low-resistivity water-bearing body,the resistivity difference between the low-resistivity body and surrounding rock,and the influence of the size of the low-resistivity body on the transient EM response.Furthermore,the common types of low-resistivity bodies in the roadway drivage process are used for modeling to analyze the attenuation of the detected EM response when there are low-resistivity bodies in front of the roadway.The research in this paper can help effectively detecting the water-bearing low-resistivity body in front of the roadway drivage and lay a foundation for reducing the risk of water seepage accidents.展开更多
Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularl...Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularly deep learning(DL),applied and relevant to computational mechanics(solid,fluids,finite-element technology)are reviewed in detail.Both hybrid and pure machine learning(ML)methods are discussed.Hybrid methods combine traditional PDE discretizations with ML methods either(1)to help model complex nonlinear constitutive relations,(2)to nonlinearly reduce the model order for efficient simulation(turbulence),or(3)to accelerate the simulation by predicting certain components in the traditional integration methods.Here,methods(1)and(2)relied on Long-Short-Term Memory(LSTM)architecture,with method(3)relying on convolutional neural networks.Pure ML methods to solve(nonlinear)PDEs are represented by Physics-Informed Neural network(PINN)methods,which could be combined with attention mechanism to address discontinuous solutions.Both LSTM and attention architectures,together with modern and generalized classic optimizers to include stochasticity for DL networks,are extensively reviewed.Kernel machines,including Gaussian processes,are provided to sufficient depth for more advanced works such as shallow networks with infinite width.Not only addressing experts,readers are assumed familiar with computational mechanics,but not with DL,whose concepts and applications are built up from the basics,aiming at bringing first-time learners quickly to the forefront of research.History and limitations of AI are recounted and discussed,with particular attention at pointing out misstatements or misconceptions of the classics,even in well-known references.Positioning and pointing control of a large-deformable beam is given as an example.展开更多
文摘Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.
基金the support from the Program for Natural Science Foundation of Zhejiang Province(LY16E080006)National Natural Science Foundation of China(51378240)+1 种基金2015 Jiangsu provincial building energy saving and construction industry science and technology project2016 Jiangsu provincial construction industry modernization base project.
文摘The failure behavior of the precast prestressed concrete T girder was investigated by destructive test and finite-element analysis,and the mid-span deflection,girder stiffness and the variation of the cross section strain in the loading process were obtained,and the mechanical properties,mechanical behavior,elastic and plastic behavior and ultimate bearing capacity of T girder with large span were revealed.Furthermore,the relationship between the beam stiffness degradation,the neutral axis in cross-section,steel yielding and concrete cracking are investigated and analyzed.A method was proposed to predict the residual bearing capacity of a bridge based on the variation of the position of the cross section strain distribution and the section neutral axis,which provided a theoretical basis for predicting the stiffness detection and carrying capacity assessment of prestressed concrete beam.
基金This project is supported by the National Natural Science Foundation of China under grant No. 50174003 and No. 59995440.
文摘With a geometrical model of porous material, a 3D finite-element analysis on the rolling process of spring steel60Si2Mn in the semi-solid state is carried out using software MARC. In terms of flat and groove rolling conditions,stress field and strain field are studied. The simulation results show that the rigid-viscoplastic model can accuratelydescribe the semi-solid metal rolling process. Semi-solid slurry has the characteristics of low flow stress and goodfluidity. During groove rolling, distribution of stress and strain on the cross-section of deformation zone is moreuniform than that during flat rolling. The results of simulation are in good agreement with the experiment data, andshow that semi-solid material fits for groove rolling.
基金Supported by the National Natural Science Foundation of China(60801050)the Excellent Talent Fund of Beijing(2011)Excellent Young Scholars Research Fund of Beijing Institute ofTechnology(2012)
文摘The finite-element modeling and simulations of the intra-body communication(IBC)were investigated to provide a theoretical basis for biomedical monitoring.A finite-element model for the whole human body was developed to simulate the IBC.The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths,and their attenuations were calculated.Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC,while the potential generally decreased along the axis of the body parts.Both signal attenuations of the two types IBC increased with increasing signal transmission distance,and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC.The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.
基金supported in part by the National Science Foundation Grant DMS-1620016supported in parts by HKSAR grant Q81Q and JRI of The Hong Kong Polytechnic University.
文摘In this paper,we introduce new stable mixed finite elements of any order on polytopal mesh for solving second-order elliptic problem.We establish optimal order error estimates for velocity and super convergence for pressure.Numerical experiments are conducted for our mixed elements of different orders on 2D and 3D spaces that confirm the theory.
文摘In a full-arch implant rehabilitation ad modum Branemark, the distribution of stress and strain in mandibular bone is influenced by the type, number and position of implants used. In particular, the biomechanical behaviour of the bone structure after complete osseointegration depends on the load transferred to the bone by each fixture. In this study, a finite-element analysis of two models was performed. Models of an all-on-four configuration and a six-implant configuration were compared in a worst-case scenario. A new V parameter is presented to aid the quantitative and comparative analysis of the all-on-four and six-implant configurations. The influence of orthotropy was also investigated, and a geometric change in the all-on-four configuration is presented.
文摘A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series ofcorresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of tem- perature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESE casting and structure of casting is fine and compact.
文摘The numerical analytic research approach of stress-strain state of anisotropic composite finite element area with different boundary conditions on the surface, is represented below. The problem is solved by using a spatial model of the elasticity theory. Differential equation system in partial derivatives reduces to one-dimensional problem using spline collocation method in two coordinate directions. Boundary problem for the system of ordinary higher-order differential equation is solved by using the stable numerical technique of discrete orthogonalization.
文摘A two dimensional plane-stress finite-element type of analysis is pre-sented to predict the behaviour of geogrids embedded in sand under pullout loadingconditions.In the analysis the interactions between soil and geogrid are simulatedby non-linear springs.The stiffnesses of the springs can be determined from simpletests in a specially designed pullout box.The proposed finite element (FE) analysis isapplied to interpret test results from a large scale pullout test rig.The predicted be-haviour of the geogrid under pullout load agrees well with the observed data includingthe load-displacement properties,the displacement distribution along the longitudinaldirection and the mobilisation of the frictional and bearing resistance.
基金This work was supported by a grant from the Capital Health Research and Development of Special Grants(No.2016-1-5012)。
文摘Background:The reverse obliquity inter-trochanteric fracture is a distinct fracture pattern that is mechanically different from most inter-trochanteric fractures and the optional treatment of it is still controversial.The purpose of this study was to compare differences in the efficacy of a novel nail(medial support nail[MSN-II])and proximal femoral nail anti-rotation(PFNA-II)in the treatment of reverse obliquity inter-trochanteric fractures(Arbeitsgemeinschaft fur Osteosynthesfrogen/Orthopedic Trauma Association[AO/OTA]31-A3.1)using finite-element analysis.Methods:Modeling software was used to establish a three-dimensional model of MSN-II and PFNA-II and an A3.1 inter-trochanteric fracture model.Abaqus software was used to implement different force loads to compare finite-element biomechanical parameters such as the maximum stress in implant and the displacement of fracture site.Results:The femoral stress,implant stress and fracture site displacement of MSN-II was less than that of PFNA-II.The results indicated that the maximal femoral stress was 581 MPa for PFNA-II and 443 MPa for the MSN-II.The maximum stress values in the PFNA-II and MSN-II models were 291 and 241 MPa,respectively.The maximal displacements of the fracture site were 1.47 and 1.16 mm in the PFNA-II and MSN-II models,respectively.Conclusions:Compared with PFNA-II for inter-trochanteric fracture(AO/OTA 31-A3.1),MSN-II which was designed with a triangular stability structure can provide better biomechanical stability.The MSN-II may be a feasible option for the treatment of reverse obliquity inter-trochanteric fracture.
文摘Under the same conditions of external force, simulations on differences of deformational energy in structural zones, which have different deformational behavior, and that on the distribution of the differences have been carried out by means of finite-element method. Shear deformational energy U_w is higher than volume deformational energy U_T by about one order of magnitude, and deformationat energy U_B( = U_w + U_T) and U_w show a trend to be larger→largest→small in quantity in structural zones of different deformational properties, which correspond to compressive, shear, tensile zones, respectively, but U_T shows a trend to be large→small gradually. This has a considerable significance in the study of tectonic heat and mineral liquid migration in association with the research on tectonic additional hydrostatic pressure.
基金supported by the National Natural Science Foundation of China(Grant Nos.72088101,41922027,41830107,41811530010)Innovation-Driven Project of Central South University(Grant No.2020CX0012)+1 种基金the National Natural Science Foundation of Hunan Province of China(Grant No.2019JJ20032)Macao Foundation and the pre-research project on Civil Aerospace Technologies funded by China’s National Space Administration(Grant Nos.D020308,D020303).
文摘The accumulated large amount of satellite magnetic data strengthens our capability of resolving the electrical conductivity of Earth’s mantle.To invert these satellite magnetic data,accurate and efficient forward modeling solvers are needed.In this study,a new finite-element based forward modeling solver is developed to accurately and efficiently compute the induced electromagnetic field for a realistic 3D Earth.Firstly,the nodal-based finite element method with linear shape function on tetrahedral grid is used to assemble the final system of linear equations for the magnetic vector potential and electric scalar potential.The FGMRES solver with algebraic multigrid(AMG)preconditioner is used to quickly solve the final system of linear equations.The weighted moving least-square method is employed to accurately recover the electromagnetic field from the numerical solutions of magnetic vector and electric scalar potentials.Furthermore,a local mesh refinement technique is employed to improve the accuracy of the estimated electromagnetic field.At the end,two synthetic models are used to verify the accuracy and efficiency of our newly developed forward modeling solver.A realistic 3D Earth model is used to simulate the induced magnetic field at 450 and 200 km altitudes which are the planned flying altitudes of Macao’s geomagnetic satellites.The simulation indicates that(1)the amplitude of the mantle-induced magnetic field can reach 10–30 nT at 450 km altitude,which is 10–30%of the primary magnetic field.The induced magnetic field at 200 km altitude has larger amplitudes.These mantleinduced magnetic fields can be measured by Macao geomagnetic satellites;(2)the amplitude of the ocean-induced magnetic field can reach 5–30 nT at satellite altitudes,which needs to be carefully considered in the interpretation of satellite magnetic data.We are confident that our newly developed forward modeling solver will become a key tool for interpreting satellite magnetic data.
基金supported in part by a Major Project of the National Science Foundation of China under Grant Nos.60990320 and 60990324in part by the National Science Foundation of China under Grant Nos.60871016,60802001,60921063 and 60901011+1 种基金in part by the Natural Science Foundation of Jiangsu Province under Grant No.BK2008031in part by the 111 Project under Grant No.111-2-05.
文摘Accurate simulations of metamaterial devices are very important in the analysis of their electromagnetic properties.However,it is very difficult to make full-wave simulations of three-dimensional(3D)metamaterial devices due to the huge memory requirements and long computing time.In this paper,we present an efficient finiteelement method(FEM)to analyze 3D axisymmetric electromagnetic devices designed by the transformation-optics approach,such as invisibility cloaks and concentrators.In the proposed method,we use the edge-based vector basis functions to expand the transverse field components,and the node-based scalar basis functions to expand the angular component.The FEM mesh is truncated with a cylindrical perfectly matched layer.We have applied the method to investigate the scattering from spherical and ellipsoidal invisibility cloaks and circularly cylindrical concentrators,in which the permittivity and permeability are both inhomogeneous and anisotropic.Numerical results are presented to show the validity and efficiency of the method.
基金supported in part by the National Key Research and Development Program of China(Grant No.2020YFB2008800)in part by the Nanchang Institute for Microtechnology of Tianjin University.
文摘This paper presents the design,fabrication,and characterization of cantilever-type resonators with a novel stacked structure.Aluminum nitride is adopted as the material for both the structural layer and the piezoelectric layer;this simplifies the fabrication process and improves the quality factor of the resonator.Both in-plane and out-of-planeflexural modes were investigated.The effect of the structural dimensions and electrode patterns on the resonator’s performance were also studied.Finite-element simulations and experiments examining anchor loss and thermoelastic damping,which are the main loss mechanisms affecting the quality factor of these resonators,were carried out.The optimal structural dimensions and electrode patterns of the cantilever-type resonators are presented.A quality factor of 7922 with a motional impedance of 88.52 kΩand a quality factor of 8851 with a motional impedance of 67.03 kΩwere achieved for the in-plane and out-of-planeflexural-mode resonators,respectively.The proposed resonator design will contribute to the development of high-performance devices such as accelerometers,gyroscopes,and pressure sensors.
基金The research was financially supported by National Natural Science Foundation of China(Grant Nos.52108302 and 52009046)Fundamental Research Funds for the Central Universities of Hua-qiao University(Grant No.ZQN-914).
文摘Soils are not necessarily uniform and may present linearly varied or layered characteristics,for example the backfilled soils behind rigid retaining walls.In the presence of large lateral thrust imposed by arch bridge,passive soil failure is possible.A reliable prediction of passive earth pressure for the design of such wall is challenging in complicated soil strata,when adopting the conventional limit analysis method.In order to overcome the challenge for generating a kinematically admissible velocity field and a statically allowable stress field,finite element method is incorporated into limit analysis,forming finiteelement upper-bound(FEUB)and finite-element lower-bound(FELB)methods.Pseudo-static,original and modified pseudo-dynamic approaches are adopted to represent seismic acceleration inputs.After generating feasible velocity and stress fields within discretized elements based on specific criteria,FEUB and FELB formulations of seismic passive earth pressure(coefficient K_(P))can be derived from work rate balance equation and stress equilibrium.Resorting to an interior point algorithm,optimal upper and lower bound solutions are obtained.The proposed FEUB and FELB procedures are well validated by limit equilibrium as well as lower-bound and kinematic analyses.Parametric studies are carried out to investigate the effects of influential factors on seismic K_(P).Notably,true solution of K_(P) is well estimated based on less than 5%difference between FEUB and FELB solutions under such complex scenarios.
文摘Regarding quality inspection of technologically important nanocomposite hard coatings based on Ti,B,Si,C,and N and bioceramics such as hydroxyapatite that are used in small-scale high-precision devices and bio-implants,it is essential to study the failure mechanisms associated with nanoindentation,such as fracture,delamination,and chipping.The stress imposed by the indenter can affect the fracture morphology and the interfacial fracture energy,depending on indenter shape,substrate type,crystallographic properties,pre-existing flaws,internal microcracks,and pre-strain.Reported here are finite-element-based fracture studies that provide insights into the different cracking mechanisms related to the aforementioned failure process,showing that the fracture morphology is affected by the interaction of different cracking events.The interfacial fracture energy,toughness,and residual stress are calculated using existing models with minor adjustments,and it is found that increasing the indenter sharpness improves the shear stress distribution,making the coating more prone to separation.Depending on the prevailing type of stress,the stress distribution beneath the depression results in either crack formation or a dislocation pile-up leading to strain hardening.Different forms of resistances resulting from the indentation process are found to affect the tip–sample conduction,and because of its stronger induced plasticity than that of a Berkovich indenter tip,a sharper cube-corner tip produces more resistance.
基金supported in part by the National Nature Science Foundation of China(NSFC)under Project 52122705。
文摘Better torque performance and higher reliability have long been the focus of research for slotted limited-angle torque motors(LATMs),which are primarily used to position first-stage valves in electrohydraulic servosystems.This paper presents a high reliability axial-flux slotted LATM with quasi-Halbach array for torque performance improvement including constant torque range(CTR)and output torque.Firstly,the structure with two sets of windings and the operation principle of the proposed slotted LATM is analyzed.Secondly,a brief design procedure is presented,the structure selections of open slot and double-stator single-rotor(DSSR)interior rotor with surface mounted quasi-Halbach permanent magnet(PM)array are illustrated,and the geometric parameters are optimized to obtain the optimal design of the proposed slotted LATM.Then,3-D finite-element method(FEM)is employed to compare the proposed slotted LATM with the conventional surface mounted PM slotted LATM in terms of cogging torque,no-load back EMF,and output torque,and the results show that the proposed LATM with quasi-Halbach array has a 10%improvement in output torque and a 25%improvement in CTR.Meanwhile,the flux linkages and torque performance of the two sets of windings under various conditions verify good magnetic isolation.Finally,prototypes of two different rotor types are manufactured and a series of experiments are performed to validate the analysis.
基金supported in part by the National Natural Science Foundation of China under Grant 51977099 and Grant 52177044in part by the Hong Kong Scholars Program under Grant XJ2019031+2 种基金in part by the China Postdoctoral Science Foundation under Grant 2019T120395in part by the Natural Science Foundation of Jiangsu Higher Education Institutions under Grant 21KJA470004in part by the Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘As members of doubly salient magnetless linear machines,linear variable flux reluctance(LVFR)and wound field flux reversal(LWFFR)machines inherit the merits of conventional magnetless linear machines such as low cost,high flux adjustment capability and high reliability.Furthermore,like linear switched reluctance machine,they have a very simple and compact long secondary,which are very attractive for long stroke applications.However,low force capability is their major defect.To solve this issue,new LVFR and LWFFR machine topologies were proposed in recent work,while lacking studies on their force improvement mechanism and further force evaluation.In this paper,LVFR and LWFFR machines with improved force performance are comparatively studied with the emphasis on their force capabilities.The operation principle of the two machines is analyzed based on magnetic field harmonics produced by flux modulation.Contributions of air-gap flux density harmonic components to no-load back electromagnetic forces of the two machines are analyzed and the average force equation is derived.Moreover,force capabilities of the both machines are investigated by means of the time-stepping finite-element analysis to verify the theoretical analysis.
文摘The hidden water-bearing structures near the roadway tunnelling face are very likely to cause water seepage accidents in coal mines.Currently,transient electromagnetic(EM)technology has be-come an important method to detect water damage in advance of roadway excavation.In this paper,the time-domain finite element algorithm based on unstructured tetrahedron grids is used to accurate-ly simulate the geological body in front of the roadway excavation face and analyze its response.The authors detect the distance between the roadway excavation face and the low-resistivity water-bearing body,the resistivity difference between the low-resistivity body and surrounding rock,and the influence of the size of the low-resistivity body on the transient EM response.Furthermore,the common types of low-resistivity bodies in the roadway drivage process are used for modeling to analyze the attenuation of the detected EM response when there are low-resistivity bodies in front of the roadway.The research in this paper can help effectively detecting the water-bearing low-resistivity body in front of the roadway drivage and lay a foundation for reducing the risk of water seepage accidents.
文摘Three recent breakthroughs due to AI in arts and science serve as motivation:An award winning digital image,protein folding,fast matrix multiplication.Many recent developments in artificial neural networks,particularly deep learning(DL),applied and relevant to computational mechanics(solid,fluids,finite-element technology)are reviewed in detail.Both hybrid and pure machine learning(ML)methods are discussed.Hybrid methods combine traditional PDE discretizations with ML methods either(1)to help model complex nonlinear constitutive relations,(2)to nonlinearly reduce the model order for efficient simulation(turbulence),or(3)to accelerate the simulation by predicting certain components in the traditional integration methods.Here,methods(1)and(2)relied on Long-Short-Term Memory(LSTM)architecture,with method(3)relying on convolutional neural networks.Pure ML methods to solve(nonlinear)PDEs are represented by Physics-Informed Neural network(PINN)methods,which could be combined with attention mechanism to address discontinuous solutions.Both LSTM and attention architectures,together with modern and generalized classic optimizers to include stochasticity for DL networks,are extensively reviewed.Kernel machines,including Gaussian processes,are provided to sufficient depth for more advanced works such as shallow networks with infinite width.Not only addressing experts,readers are assumed familiar with computational mechanics,but not with DL,whose concepts and applications are built up from the basics,aiming at bringing first-time learners quickly to the forefront of research.History and limitations of AI are recounted and discussed,with particular attention at pointing out misstatements or misconceptions of the classics,even in well-known references.Positioning and pointing control of a large-deformable beam is given as an example.
