Importance of Three-Dimensional Piezoelectric Coupling Modeling in Quantitative Analysis of Piezoelectric Actuators

This paper demonstrates the importance of three-dimensional(3-D)piezoelectric coupling in the electromechanical behavior of piezoelectric devices using three-dimensional finite element analyses based on weak and strong coupling models for a thin cantilevered piezoelectric bimorph actuator.It is found that there is a significant difference between the strong and weak coupling solutions given by coupling direct and inverse piezoelectric effects(i.e.,piezoelectric coupling effect).In addition,there is significant longitudinal bending caused by the constraint of the inverse piezoelectric effect in the width direction at the fixed end(i.e.,3-D effect).Hence,modeling of these effects or 3-D piezoelectric coupling modeling is an electromechanical basis for the piezoelectric devices,which contributes to the accurate prediction of their behavior.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.

Evaluation value of three-dimensional finite element model analysis for bone mineral density and bone metabolism activity in patients with osteoporosis

Objective: To study the evaluation value of three-dimensional finite element model analysis for bone mineral density (BMD) and bone metabolism activity in patients with osteoporosis. Methods: A total of 218 patients who were diagnosed with osteoporosis in the hospital between February 2014 and January 2017 were collected as observation group, and 100 healthy volunteers who received physical examination in the hospital during the same period were selected as normal control group. The femoral head of the two groups was analyzed by three-dimensional finite element model, and the femoral head BMD levels and serum bone metabolism index contents were measured. Pearson test was used to evaluate the evaluation value of femoral head three-dimensional finite element model for osteoporosis. Results: The cancellous bone and cortical bone Von Mises stress value of observation group were lower than those of normal control group, and femoral neck BMD value of observation group was lower than that of normal control group;serum bone metabolism index BGP content was lower than that of normal control group while NBAP, TRACP-5b and CTX-1 contents were higher than those of normal control group. Pearson test showed that the cancellous bone and cortical bone Von Mises stress value of patients with osteoporosis were directly correlated with BMD value and bone metabolism index contents. Conclusion: The three-dimensional finite element model analysis resultsof patients with osteoporosis can objectively reflect the femoral headBMD value and bone metabolism activity, and is a reliable way to evaluate the risk of long-term fractures.

RIGID FINITE ELEMENT AND LIMIT ANALYSIS

According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.

A pseudo-coupled numerical approach for stability analysis of frozen soil slopes based on finite element limit analysis method

To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on the computed thermal-hydro field, the stability of frozen soil slope is evaluated. Although the shear strength for frozen soil is very complicated and is usually represented by a nonlinear MC failure criterion, a simple linear MC yield criterion is utilized. In this method, the internal friction angle is expressed as a function of volumetric ice content and the cohesion is fitted as a simple bilinear expression of Tand volumetric water content. To assess slope stability, the limit analysis is employed in conjunction with the recently developed a-section search algorithm. A frozen soil slope example is used to examine the proposed pseudo-coupled numerical approach, and numerical studies validate its effectiveness. Based on numerical results, it is seen that slope stability may be remarkably influenced by warming air （or grotmd surface） temperature. With increasing ground surface temperature, slope stability indicated by FOS may reduce to 1.0, implying that wanning air temperature could be a trigger of frozen soil slope failure.

Upper-bound limit analysis based on the natural element method

The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic the- orem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nons- moothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.

Development of Finite Element Limiting Analysis Method and Its Applications in Geotechnical Engineering

The Finite Element Limiting Analysis Method(LELAM) has the advantage of combining a numerical analysis method with traditional limiting equilibrium methods.It is particularly applicable to the analysis and design of geotechnical engineering.In the early 20th century,FELAM has been developed vigorously in domestic geotechnical engineering over international common finite element procedures.It has made great achievements in basic theory research and computational precision,thus broadening the application fields in practical projects.In order to gradually make innovations in geotechnical design methods,some of our research results are presented,mainly including geotechnical safety factor definitions,the principles for use of the method concerned,the overall failure criterion,the deduction and selection of the yield criterion,and the measurement to improve the computational precision,etc..The application field has been broadened from two-dimensional to three-dimensional,from soil slope to jointed rock slope and foundation,from stable seepage to non-stable seepage,from slope and foundation to tunnel.This method has also been used in search of many hidden sliding surfaces of complex landslides,conducting the structural support design considering the interaction between the soil and the structure,and computing simulation foundation bearing plates load tests,etc..

PLASTIC LIMIT ANALYSIS OF INCOMPATIBLE FINITE ELEMENT METHOD

This paper describes an incompatible finite element model satisfying the consistency condition of energy to solve the numerical precision problem of finite element solution in perfectly plastic analysis. In this paper the reason and criterion of the application of the model to plastic limit analysis are discussed, and an algorithm of computing plastic limit load is given.

Stochastic analysis of excavation-induced wall deflection and box culvert settlement considering spatial variability of soil stiffness

In this study, a three-dimensional (3D) finite element modelling (FEM) analysis is carried out to investigate the effects of soil spatial variability on the response of retaining walls and an adjacent box culvert due to a braced excavation. The spatial variability of soil stiffness is modelled using a variogram and calibrated by high-quality experimental data. Multiple random field samples (RFSs) of soil stiffness are generated using geostatistical analysis and mapped onto a finite element mesh for stochastic analysis of excavation-induced structural responses by Monte Carlo simulation. It is found that the spatial variability of soil stiffness can be described by an exponential variogram, and the associated vertical correlation length is varied from 1.3 m to 1.6 m. It also reveals that the spatial variability of soil stiffness has a significant effect on the variations of retaining wall deflections and box culvert settlements. The ignorance of spatial variability in 3D FEM can result in an underestimation of lateral wall deflections and culvert settlements. Thus, the stochastic structural responses obtained from the 3D analysis could serve as an effective aid for probabilistic design and analysis of excavations.

Mechanical analysis of the femoral neck dynamic intersection system with different nail angles and clinical applications

BACKGROUND The femoral neck dynamic intersection system(FNS)is mechanically more stable than other internal fixation techniques.Current studies have confirmed that the structural design of FNS has good biomechanical properties in European and American populations.However,whether the suitability of the FNS's 130°main nail angle design for Asian populations has been thoroughly investigated remains unclear.AIM To compare the biomechanical stability differences among different main nail angles of the FNS in the treatment of femoral neck fractures in Asian populations.METHODS Computed tomography data of the femur of healthy adult male volunteers were imported into Mimics software to create a three-dimensional model of the femur.The model was adapted to the curve using Geomagic software and imported into Solidworks software to construct the Pauwels I femoral neck fracture model and design the FNS internal fixation model using different main nail angles.Afterward,the models were assembled with the FNS fracture model and meshed using the preprocessing Hypermesh software.Subsequently,they were imported into Abaqus software to analyze and evaluate the biomechanical effects of different angles of the FNS main nail on the treatment of femoral neck fractures.RESULTS The peak displacement of the proximal femur under different angles of FNS fixation under stress was 7.446 millimeters in the 120°group and 7.416 millimeters in the 125°group;in the 130°,135°,and 140°FNS fixation groups,the peak displacement was 7.324 millimeters,8.138 millimeters,and 8.246 millimeters,respectively.In the 120°and 125°FNS fixation groups,the maximum stresses were concentrated at the main nail and the anti-rotation screw,which intersected the fracture line of the femur neck,resulting in peak stresses of 200.7 MPa and 138.8 MPa,respectively.Peak stresses of 208.8 MPa,219.8 MPa,and 239.3 MPa were observed on the angular locking plate distal to the locking screw in the 130°,135°,and 140°fixation groups.CONCLUSION FNS has significant stress distribution properties,a minimal proximal femoral displacement,and an optimal stability for treating femoral neck fractures in Asian populations when performed with a 130°main nail angle.

Finite element analysis of slope stability by expanding the mobilized principal stress Mohr's circles-Development, encoding and validation

In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies.However,the application of finite element method(FEM)to slope stability as a strength reduction method(SRM)or as finite element limit analysis(FELA)is not always a success for the drawbacks that characterize both methods.To increase the performance of finite element analysis in this problem,a new approach is proposed in this paper.It consists in gradually expanding the mobilized stress Mohr’s circles until the soil failure occurs according to a prescribed non-convergence criterion.The present approach called stress deviator increasing method(SDIM)is considered rigorous for three main reasons.Firstly,it preserves the definition of the factor of safety(FOS)as the ratio of soil shear strength to the mobilized shear stress.Secondly,it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane,on which the shear strength takes place.Thirdly,by introducing the concept of equivalent stress loading,the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor.The new numerical procedure was encoded in a Fortran computer code called S^(4)DINA and verified by several examples.Comparisons with other numerical methods such as the SRM,gravity increasing method(GIM)or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

Limit load and failure mechanisms of a vertical Hoek-Brown rock slope

The problem considered in this short note is the limit load determination of a vertical rock slope.The classical limit theorem is employed with the use of adaptive finite elements and nonlinear programming to determine upper and lower bound limit loads of a Hoek-Brown vertical rock slope.The objective function of the mathematical programming problem is such as to optimize a boundary load,which is known as the limit load,resembling the ultimate bearing capacity of a strip footing.While focusing on the vertical slope,parametric studies are carried out for several dimensionless ratios such as the dimensionless footing distance ratio,the dimensionless height ratio,and the dimensionless rock strength ratio.A comprehensive set of design charts is presented,and failure envelopes shown with the results explained in terms of three identified failure mechanisms,i.e.the face,the toe,and the Prandtl-type failures.These novel results can be used with great confidence in design practice,in particularly noting that the current industry-based design procedures for the presented problem are rarely found.

A MATHEMATICAL PROGRAMMING ALGORITHM FOR LIMIT ANALYSIS

This paper deals with the limit analyses of perfect rigid-plastic continua.Based on the kinematic theorem of the limit analysis theory,a mathematical programming finite element formula for determining the upper bound load multiplier has been established,and an iteration algorithm proposed accordingly.In this algorithm the plastic and rigid zones are distinguished for every iteration step,and the goal function is modified gradually.The difficulties caused by the nonsmoothness of the goal function are over- come.Some examples solved by this algorithm are presented.

Lower Bound Limit Analysis of Anisotropic Soils

Previous approaches can only tackle anisotropic problems with cohesion varying with direction.A novel linearization of the Mohr-Coulomb yield criterion associated with plane strain problem has been achieved by simulating the Mohr’s circle with orientation lines inσ-τspace,which allows for lower bound solution of soils with cohesion and friction coefficient varying with direction.The finite element lower limit analysis formulation using the modified anisotropic yield criterion is then developed.Several examples are given to illustrate the capability and effectiveness of the proposed numerical procedure for computing rigorous lower bounds for anisotropic soils.

Upper Bound Limit Analysis of Anisotropic Soils

In this paper,a novel discretization method inσ-τspace is developed to calculate the upper bound limit loads and failure modes of anisotropic Mohr-Coulomb materials.To achieve this objective,the Mohr-Coulomb yield criterion is linearized inσ-τspace,which allows for upper bound solution of soils whose cohesion and friction coefficient varying with direction.The finite element upper limit analysis formulation using the modified anisotropic yield criterion is then developed.Several examples are given to illustrate the capability and effectiveness of the proposed numerical procedure for computing rigorous upper bounds for anisotropic soils.

Stress Distribution in Maxillary Central Incisors without Ferrules: A Finite Element Analysis of Post and Core Systems

Purpose: The purpose of this study was to identify optimal post and core materials for central incisors without ferrules using three-dimensional finite element analysis and three-point bending tests. Methods: Stress analyses were performed with six models: cast metal post and core (MP), composite resin core alone, straight fiber-reinforced post-composite resin core (FSR), tapered fiber-reinforced post-composite resin core, straight titanium post-composite resin core (TSR), and tapered titanium post-composite resin core (TTR). A 100-N load was applied to the lingual surface at a 45° angle to the long axis of the tooth. Maximum von Mises stress distributions were calculated with finite element analysis software. Five samples each of composite resin, straight fiber-reinforced post, straight titanium post, straight fiber-reinforced post and composite resin, and straight titanium post and composite resin were subjected to three-point bending tests, followed by analysis of variance and Tukey’s multiple comparison test. Results: Stress distribution was optimal on TTR. Maximum von Mises stress on the cervical side of the post was greatest in TSR (693 MPa) and TTR (556 MPa). Maximum stress on the apical side of the post was greatest in MP (110 MPa). Maximum stress in surrounding dentin was lowest in MP (203 MPa) and TTR (250 MPa). Gap distance was smallest in MP (0.09 mm) and largest in FSR (0.26 mm). Mean maximum three-point bending force was lowest in composite resin (26.9 N/mm) and highest in titanium post and composite resin (97.1 N/mm). Titanium post bending strength was consistently greater than that of the fiber-reinforced post (p < 0.01). Conclusion: These results revealed optimal stress distribution and high bending strength with the tapered titanium post and resin combination, suggesting that this combination can most effectively prevent root or post fracture in an anterior tooth without a ferrule.

Arc-length technique for nonlinear finite element analysis

Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ordinary solution techniques lead to instability near the limit points and also have problems in case of snap-through and snap-back. Thus they fail to predict the complete load-displacement response. The arc-length method serves the purpose well in principle, received wide acceptance in finite element analysis, and has been used extensively. However modifications to the basic idea are vital to meet the particular needs of the analysis. This paper reviews some of the recent developments of the method in the last two decades, with particular emphasis on nonlinear finite element analysis of reinforced concrete structures.

Current Limiting Characteristic Analysis of Magnetic-controlled Switcher Type Fault Current Limiter

A novel magnetic-controlled switcher type fault current limiter (FCL) based on the topology of the saturated iron core high temperature superconducting FCL is proposed. The magnetic field distribution of the FCL iron core is analyzed by FEA software ANSYS. The current limiting characteristic is investigated by both 3-D field-circuit coupled simulation and Matlab. The experiments on the 220 V/50 A test model show that the FCL can limit the fault current swiftly and effectively,and the FCL has the advantages of simple and reliable structure, flexible control strategy. The simulation and experimental results prove that the theoretical expectation and current limiting performance is satisfactory for practical use.

An efficient probabilistic design approach for tunnel face stability by inverse reliability analysis

In order to maintain the safety of underground constructions that significantly involve geo-material uncertainties,this paper delivers a new computation framework for conducting reliability-based design(RBD)of shallow tunnel face stability,utilizing a simplified inverse first-order reliability method(FORM).The limit state functions defining tunnel face stability are established for both collapse and blow-out modes of the tunnel face failure,respectively,and the deterministic results of the tunnel face support pressure are obtained through three-dimensional finite element limit analysis(FELA).Because the inverse reliability method can directly capture the design support pressure according to prescribed target reliability index,the computational cost for probabilistic design of tunnel face stability is greatly reduced.By comparison with Monte Carlo simulation results,the accuracy and feasibility of the proposed method are verified.Further,this study presents a series of reliability-based design charts for vividly understanding the limit support pressure on tunnel face in both cohesionless(sandy)soil and cohesive soil stratums,and their optimal support pressure ranges are highlighted.The results show that in the case of sandy soil stratum,the blowout failure of tunnel face is extremely unlikely,whereas the collapse is the only possible failure mode.The parametric study of various geotechnical uncertainties also reveals that ignoring the potential correlation between soil shear strength parameters will lead to over-designed support pressure,and the coefficient of variation of internal friction angle has a greater influence on the tunnel face failure probability than that of the cohesion.