Force-displacement Research Articles

Displacement and Internal Force Response of Mechanically Connected Precast Piles Subjected to Horizontal Load Based on the m-Method

Mechanically connected precast piles are a type of precast piles that utilise snap-type mechanical connectors to restrain the pile ends of two identical or different precast piles at the top and bottom so as to quickly realise the purpose of the connection. However, the gap problem in the connectors of mechanically connected piles can lead to uneven and uniform deformation of the piles under horizontal loading, resulting in additional displacements and rotation angles of the piles at the connection. Solving the problem of calculating the internal force response of discontinuous deformed piles is a prerequisite for promoting and applying mechanically connected precast piles. Firstly, the theoretical derivation of mechanically connected piles with fixed constraints at the pile bottom is carried out. Secondly, the pile response equations of mechanically connected piles are established, and the theoretical solutions of pile displacement and internal force response of mechanically connected piles under horizontal loading are derived. Thirdly, the pile-soil model of the test pile is established using ABAQUS software (ABAQUS 2016) in combination with the design data of the test pile. The numerical simulation displacements and angles of rotation are compared with the test results. Finally, the theoretical and numerical simulation displacements and internal forces of the ordinary pile and the mechanically connected pile are compared. The relative errors of the displacements and angles of rotation of the established pile-soil model are less than 10%, indicating that the established model has good accuracy. The relative errors of the theoretical and numerical simulation displacements and internal forces of the mechanically connected pile are less than 10%, proving the correctness of the theoretical calculation by the m-method. This study can provide effective theoretical support and methodological guidance for the displacement and internal force response of discontinuous piles.

Comprehensive design strategy of the NSD–TMD controlled system with closed-form frequency-amplitude solutions

A negative stiffness device (NSD) utilizing pre-compressed helical springs has been integrated with diverse vibration control techniques. Recently, the incorporation of NSD has been proposed to enhance the efficacy of a tuned mass damper (TMD) controlled system. This study introduces two design strategies for the NSD–TMD system, which functions either as a TMD with an adjustable operating bandwidth or as a pure non-linear energy sink. First, the non-linear force of NSD was expressed in a Taylor series approximation form incorporating a negative linear stiffness and a positive cubic stiffness term. Second, the equations of motion were formulated, and closed-form solutions were derived using the complexification–averaging method. The accuracy of the Taylor series approximation was substantiated within a specified range through a comparison of force–displacement curves. Third, the accuracy of the dynamic closed-form solutions was affirmed by comparing the amplitude–frequency curves with the numerical simulations. For a case study, an actual timber pedestrian bridge with decaying natural frequency was selected, and the two design strategies were examined individually. The vibration control performance of the TMD and NSD–TMD was also comparatively analyzed. The assessment focused on the peak force transmissibility value and the effective operating frequency bandwidth. Results highlighted a notable enhancement in operational performance, particularly for primary structures exhibiting diminishing modal frequencies. In a TMD-controlled system, the efficacy of vibration control deteriorates and may even amplify the dynamic response when a modal frequency shift occurs in the primary structure. Conversely, both design strategies of NSD–TMD can achieve lower transmissibility and broaden the operating bandwidth.

Numerical analysis by comparing the Drucker-Prager and Mohr-Coulomb behaviors to optimize the best sites for constructing a new tunnel adjacent to the old one

The tunnels are one of the most significant infrastructures for facilitating mobility in mountainous terrain. It has several benefits, including shortening distances, time, and ease of transportation. To make it easier, developed countries create several tunnels close to each other. Engineers face the necessity of building these new tunnels near existing ones; this infrastructure must be safe. In this paper, for this numerical study, we used the OptumG2 software to analyze the position of a new tunnel created in six locations: bottom, top, left (two places at 12 m and 32 m), and right (two places at 12 m and 32 m) by the finite element method. This study was analyzed using two criteria: the equivalent of Drucker Prager and Mohr Coulomb. We are interested in the total field of the displacement, the vertical displacement of the model, the displacement, bending moment, shear, and normal force of the concrete lining. The findings of this numerical analysis using the Drucker Prager and Mohr Coulomb behaviors to determine the ideal locations for building a new tunnel next to an existing one show that the new bottom tunnel might have a significant influence on the stability of the existing tunnel. It was also observed that the maximum displacement present when we add a new tunnel at the bottom meets the Drucker-Prager criteria.

Cyclic performance of precast steel beam-to-CECFST column exterior sub-assemblages with dry connections: Experimental studies and mechanism analyses

Combining advantages of both reinforced concrete (RC) and concrete-filled steel tube (CFST), concrete-encased concrete-filled steel tube (CECFST) members demonstrate superior structural performance over pure RC or CFST counterparts. Nevertheless, current research in this domain has primarily focused on member and connection behaviour of CECFST columns, with no attention given to frame-level performance. This significantly impedes applications of this form of high-performance composite members into actual engineering practice. To bridge this research gap, a novel precast steel beam-to-CECFST column composite system is proposed in this study. This is followed by a comprehensive experimental programme, where seven precast steel beam-to-CECFST column exterior sub-assemblages were tested under cyclic loads. In the test specimens, three types of connections, including beam-to-beam (B-B), column-to-column (C-C) and beam-to-column (B-C), were incorporated to establish a precast steel beam-to-CECFST column sub-assemblage. Horizontal displacement versus shear force relationship and failure mode of the test specimens are reported, followed by a detailed analysis of test results and findings. Furthermore, a comprehensive mechanism analysis is presented, with emphasis on energy dissipation behaviour and deformation mechanism of the proposed precast composite system. The test results show that the proposed precast composite system has comparable load-carrying capacity and energy dissipation capacity in comparison with conventional cast-in-place counterparts. This indicates that despite using precast beam-to-beam and column-to-column connections, structural performance of the composite frame system would not be significantly affected, showing that this form of precast steel beam-to-CECFST column composite system has considerable potential to be adopted in actual engineering practice.

Study on Vibration Reduction Effect of the Building Structure Equipped with Intermediate Column–Lever Viscous Damper

Generally speaking, the traditional lever amplification damping system is installed between adjacent columns in a building, which occupies a significant amount of space in the building. In contrast to amplification devices in different forms, the damper displacement of the intermediate column damper system is smaller, and the vibration reduction efficiency is lower. In light of these drawbacks, this study proposes a new amplification device for energy dissipation and vibration reduction, which is based on an intermediate column–lever mechanism with a viscous damper (CLVD). Initially, a specific simplified mechanical model of CLVD is derived. Subsequently, an equivalent Kelvin mechanical model of CLVD is derived to intuitively reflect CLVD’s damping and stiffness effect. The damping ratio added by CLVDs to the structure is calculated according to that model; the additional damping ratio and additional stiffness are utilized to calculate the displacement ratio Rd and shear force ratio Rv of the structure with CLVDs to the structure without CLVDs. Rd and Rv are introduced to evaluate the vibration reduction effect of the structure with CLVDs, and the effects of various parameters (such as intermediate column position, beam’s bending line stiffness, lever amplification factor, damping coefficient, and earthquake intensity) on Rd and Rv are analyzed. The results indicate that when the ratio of the distance from the intermediate column to the edge column to the span of the beam is 0.5, CLVD owns the optimal vibration reduction effect. Increasing the beam’s bending line stiffness is beneficial for CLVD to control structural displacement and shear force; when the leverage amplification factor is too large, the CLVD provides the structure with stiffness as the main factor, followed by damping. Additionally, when the ratio of the displacement amplification factor to the geometric amplification factor satisfies fd/γ = 1/21−0.5α, the CLVD has the optimal displacement control effect on the structure. After that, measures are provided to optimize the CLVD in different situations in order to effectively control the inter-story displacement and the story shear force of the structure. Consequently, a nine-story frame is taken as an example to elaborate the application of CLVDs in the design for energy dissipation and vibration reduction. The results reveal that the CLVD scheme adopting the proposed optimization method can effectively enhance the displacement amplification ability of CLVDs, resulting in an additional damping ratio of up to 12%. At the same time, the inter-story displacement was reduced by almost 40% under fortification earthquakes. Through the research in this study, designers can obtain a new choice in structural vibration reduction design.

The method for calculating phase angle between exciter force of vibration exciter and roller displacement

Introduction. In the process of soil compaction it is important to have information about the current density of the layer, as it enables to quickly adjust the load on the compacted material. The field methods of compaction quality assessment do not cope with this task, as they make point estimation within the pavement area. Therefore, continuous compaction monitoring systems installed on vibratory road rollers are becoming increasingly common. The systems developed by BOMAG and AMMANN require, among other things, the phase angle between the exciter force and the roller movement to calculate the compaction quality index. The phase angle is determined by the unbalance position sensor, which is very labor-intensive. In addition, continuous compaction monitoring systems include an accelerometer. The purpose of this paper is to develop an indirect method for calculating the phase angle from accelerometer readings.The method of research. In order to achieve the purpose of the work, a roller-soil single-mass model in a typical mode for vibratory rollers (periodic loss of contact) has been studied. As a result of modeling it has been found that the reaction of the compacted material has the main influence on the vertical component of a roller acceleration and practically does not affect the horizontal component. This is confirmed by the experimental data.Results. The phase angle can be determined by mutual correlation of the horizontal and vertical acceleration signals of the roller obtained with the accelerometer.Conclusion. The study proposes a new method of calculating the phase angle between the exciter force and the roller displacement, which eliminates the direct measurement of this angle. The calculation of the angle is based on the readings of a two-axis accelerometer installed on the road roller. The proposed method enables to simplify the system of continuous compaction control and reduce the labor intensity of phase angle measurement.

Low-velocity impact analysis and multi-objective optimization of hybrid carbon/basalt fibre reinforced composite laminate

Three-dimensional in-plane progressive damage model was developed for carbon/basalt fibre interlayer hybrid composites in this study, and cohesive element model was employed to simulate interlaminar damage. Numerical simulation of low-velocity progressive impact damage of hybrid composite laminate was carried out, and its accuracy was verified by physical impact test, with the maximum error of the maximum contact force, maximum displacement and energy absorption not exceeding 6.1% and the maximum error in damage area size within 4.8%. The effects of lay-up structure and angle on the impact properties of hybrid laminate were investigated through simulation. Finally, real number coding strategy was proposed for the parametric representation of the lay-up material and angle of each layer. Multi-objective optimization of hybrid laminate based on Kriging surrogate model and hybrid algorithm composed of elitist non-dominated sorting genetic (NSGA2) algorithm and multi-objective particle swarm (MPSO) algorithm was conducted. The optimal trade-off solution was selected from the Pareto solutions using technique for ordering preferences by similarity to ideal solution (TOPSIS) coupled with principal component analysis (PCA) methods. The coefficient of variance between peak force and maximum displacement of the optimized laminate was 49.3%, 44.9% and 31.4% lower than those of the pure CFRP laminate, pure BFRP laminate, and hybrid laminate before optimization, respectively, indicating that the balance degree of impact displacement and force has been substantially improved. The progressive impact damage simulation and parametric optimization design methodology proposed in this study provide useful guidance for the design of hybrid composite laminate.

Experimental and numerical study of underwater laser transmission welding

This work demonstrates the feasibility of the underwater laser transmission welding (LTW) of polycarbonate sheets. The experiments are performed to investigate the effect of water depth, laser power, beam diameter, and scan speed on the weld strength, weld dimensions, and bond morphology. The results of mechanical properties and bond morphology of the weld joint in underwater are compared with air environment. The lap shear tests are performed to obtain the shear force–displacement of the welded joints. A three-dimensional (3-D) numerical model is developed to estimate the interface temperature in the underwater condition, which is validated by comparing the numerical weld width and weld depth with their corresponding experimental value. The bond morphology of the weldments is analyzed by the optical microscope and scanning electron microscopy (SEM). The results show that a maximum mean breaking force of 567.5 N is obtained at the 20 mm water depth which is 6 % higher than air due to lower thermal degradation and wider bubble-less region inside the welded region.

Study on prestress distribution and structural performance of heptagonal six-five-strut alternated cable dome with inner hole

The cable dome structure is a highly efficient prestressed space structure with reasonable force, attractive shape, and large spanning capacity. It is increasingly widely used in stadiums, exhibition halls, terminals, and other large-scale public buildings and the roof structure of various transport hub projects. To solve the problem of weak overall structural stability of the traditional cable dome structure, a new type of prestressed space structure is proposed - heptagonal six-five-strut alternated cable dome with an inner hole. The composition of the structural system was given and a theoretical study of the distribution of structural prestress was carried out; based on the symmetry and periodicity conditions, the initial prestress distribution of the structure was accurately calculated using the node balance method, and the formulae for the geometric parameters of the structural cable-strut were derived; the numerical model was established by using the finite element software ANSYS 2022 R1, and the force characteristics of the structure under full-span uniform load, half span uniform live load, wind load, and temperature were studied. The effects of structural prestress, overall shape, and support stiffness, totaling three categories and six parameters, on the structural mechanical performance of the internal force of the member, node displacement, and support reaction force were analysed in detail and the sensitivity of the parameters is assessed quantitatively, accordingly, the values of the main design parameters such as thick-span ratio are given as suggestions. The analysis results show that the new structural system has a high load-carrying capacity and deformation resistance, and its proposal provides a new scheme and new ideas for the selection and design of the cable dome.

Chaotic Vortex-Induced Vibrations of Rigid Cylinders with Nonlinear Snapping Support

Vortex-Induced Vibration (VIV) is a complex fluid–structure interaction in offshore structures. Traditionally, this phenomenon is considered periodic; however, many of its signals show chaotic behavior. The basic model already employed by other researchers is a rigid circular cylinder with linear springs and dampers. In this work, nonlinear snapping support is used to model nonlinearity in the system. To numerically simulate the flow, Reynolds-Averaged Navier–Stokes (RANS) equations for two-dimensional incompressible unsteady flows are applied. The degree of nonlinearity of the system can be changed by manipulating [Formula: see text], which is one of the geometric properties of the spring and takes values between 0 and 1. The 0–1 test, Poincaré section, and Fast Fourier Transform are used to analyze the cylinder and lift force behavior. Also, the Hilbert transform is applied to the signals, and the phase shift between displacement and lift force is obtained. The results show that the system behavior consists of branches: branch I and branch II. The large amplitudes occur in branch II. It is found that chaos emerges at the beginning of branch II, regardless of the value of [Formula: see text]. By raising the [Formula: see text] value, the span of branch II becomes more expansive, and its first point is placed at lower reduced velocities. Also, the wake dynamics becomes more regular at the end of branch I and more irregular at the beginning of branch II with the increase in [Formula: see text]. When the cylinder displacement signal is chaotic, the lift force behavior is also chaotic, but not vice versa.

Suture materials in maxillofacial surgery – A comprehensive review

By holding tissues together to speed up the healing process, surgical sutures are used to secure and restore active or trauma-induced wounds. In dentistry, a variety of sutures are employed. Both absorbable and nonabsorbable suture materials are part of the primary category. A well-fitting suture prevents the displacement forces caused by external agents, functional movements, and muscle insertions. The potential of novel suture materials to resemble wounds and cure wounds is driving up their manufacture in recent times. In order to help the clinician obscure the best possible soft tissue care, this review provides an update on the features and guiding principles of suture selection.

Experimental, analytical and numerical study of reinforced concrete beams impacted by deformable bodies at low velocities

In this work, the structural behavior of Reinforced Concrete (RC) beams struck by rigid and deformable impactors at low velocities is studied. An experimental campaign consisting of three-point bending, crushing, and drop-weight tests is carried out to characterize and compare the quasi-static and dynamic behaviors of the colliding bodies. Three repetitions per configuration are tested to obtain reliable data and to check the consistency and repeatability of the proposed experiments. The results are accurately reproduced using Non-Linear Finite Element (NLFE) models to obtain insights regarding the influence of strain rates in the experiments. Furthermore, a simple analytical method that considers the deformability of both bodies through their corresponding force–displacement relationships is presented. The proposed model is observed to predict with a reasonable accuracy the contact force profiles and the deformations of both bodies in almost all of the studied scenarios, while being orders of magnitude faster than the NLFE simulations.

Development of an inertia-driven resonant piezoelectric stack pump based on the flexible support structure

This paper proposes an inertia-driven resonant piezoelectric stack pump based on a flexible support structure to solve the problem that the piezoelectric stack cannot effectively drive the diaphragm pump to transport liquid due to too small output displacement and too high resonant frequency when one end is fixed. Under the inertial force generated by the vibration of the piezoelectric stack’s mass center during its deformation, the whole piezoelectric stack will vibrate with the flexible support structure; and a large displacement and inertial force can be achieved to drive the pump at the resonant frequency. Piezoelectric pumps are designed with a diaphragm pump and a piezoelectric stack based on the flexible support structure. The piezoelectric vibrator includes a piezoelectric stack, a preloading component and a flexible support plate. A fixed support plate and three flexible support plates with different stiffnesses were fabricated and assembled with the same piezoelectric stack and diaphragm pump respectively to construct four piezoelectric pump prototypes with different resonant frequencies. The temperature rise characteristics of the piezoelectric stack were experimentally studied to determine the safe range of the driving voltage and frequency. Then the output performances of the piezoelectric pumps were tested. Under a sinusoidal driving voltage of 100 Vpp, the piezoelectric pump based on the fixed support structure cannot pump water, while the piezoelectric pumps based on the flexible support structure achieved the maximum flow rates of 89.0 ml min−1, 123.4 ml min−1 and 197.4 ml min−1 at the resonant frequencies of 262 Hz, 297 Hz and 354 Hz, and the maximum backpressures of 4.4 kPa, 7.5 kPa and 11.0 kPa at 266 Hz, 309 Hz and 365 Hz.

Fracture behaviors of the 1500 MPa-grade anti-oxidation hot-stamped steel: Measurement, numerical simulation and experimental validation

Ultra-high strength hot-stamped steel has been increasingly and widely used for the improvement of the collision safety and lightweight level. The 1500 MPa-grade anti-oxidation hot-stamped steel (AHS-1500), as a newly developed hot-stamped steel, is promising in the automotive industry due to its advantages of anti-oxidation, high strength, high hardenability, and no coating. It is necessary to figure out its fracture failure behaviors and further predict the fracture risk for promoting the large-scale application. In this paper, the fracture behaviors of the AHS-1500 steel plate are investigated by experiments and simulation. The fracture tests are conducted under static and dynamic loads by considering various stress states, namely uniaxial tension, notched static tension, shear, punching, center-hole static tension, and VDA bending. The fracture failure behaviors are obtained. Moreover, an Swift+Hockett-Sherby combined hardening model is adopted to extrapolate the post-necking flow stress to obtain the stress-strain relation at large strains. A fracture failure model is further established based on the GISSMO fracture failure criterion to accurately describe the failure behaviors under various stress states. The nonlinear fracture failure strain and the necking instability strain are determined with the consideration of the influence of element size and strain rate. Finally, the three-point bending experiment of a cap beam part is conducted and simulated to verify the effectiveness of the proposed model. The simulated results match the experimental ones very well, with the absolute relative errors of only 1.4 % and 4.6 % for the peak force and peak displacement, respectively. This research is favourable for contributing to high-precision numerical analysis for the design of automotive structural components and providing technical support for large-scale promotion of the AHS-1500 steel.

A rate-dependent damage mechanics model for predicting plasticity and ductile fracture behavior of sheet metals at high strain rates

Uniaxial tensile tests were performed on H340 steel sheets at different strain rates (10-4 to 103s-1) and temperatures (–30 °C and 280 °C). Oscillation-free forces were measured during high-speed tests at strain rates of up to 1000 s-1 using specimens with different stress states. Material hardening curve, strain rate sensitivity, temperature effects, and Taylor-Quinney coefficient for adiabatic temperature calculations were determined in experiments. Digital image correlation (DIC) technique was employed to measure displacement, deformation, and local strain fields. Meanwhile, the temperature fields of specimen gauge section were measured with a high–speed thermal camera in the uniaxial tensile tests at various strain rates. Damage and fracture-related parameters were calibrated and validated using porosity measurements on SEM micrographs and in combination with finite element (FE) simulations. A rate– and temperature–dependent plasticity and damage mechanics model (e2MBW) was proposed and calibrated to predict the plasticity and fracture behavior of H340 under different loading speeds. The study demonstrated good agreement in the overall experimental and simulated force–displacement responses and local strain evolution across all fracture specimens at loading speeds from 0.005 mm/s to 10000 mm/s.

Effect of heat treatment on AlSi10Mg composite 3D printed energy absorption structures

To enhance the energy absorption characteristics of the energy-absorbing structure in metal alloy 3D printing, various heat treatment conditions were studied on the 3D printed hourglass multi-cellular energy-absorbing structure (HTMEAS). The results reveal that as the annealing temperature rises, the structural morphology of the α-Al+Si phase in AlSi10Mg alloy undergoes alterations, particularly when the annealing temperature remains below 300°C. After heat treatment, the HTMEAS can contract layer by layer along the gradient and undergo controllable and orderly deformation. The maximum energy absorption, displacement, and peak force reached 758 J, 30 mm, and 30 kN, respectively.These values represent a significant enhancement of 146.8 %, 157.8 %, and 50 % compared to those without heat treatment, showcasing commendable mechanical properties. This advancement establishes a solid theoretical framework and assurance for the future development and refinement of 3D printing energy-absorbing structures.

Unloading properties based elastoplastic continuous contact force model: Restitution and stiffness characterization

A nonlinear contact force model is developed for spherical elastoplastic impacts. The nonlinear contact force model is constructed based on an existing piecewise model and the other based on the force–displacement diagram. The energetic impact law of the coefficient of restitution is used to solve for the unloading restitution stiffness. An integration of the dynamic equation of motion and proper use of initial conditions is performed to solve for the maximum deformation. Furthermore, the contact force law from the force–displacement diagram is used at maximum deformation to obtain a physical expression of the final indentation in terms of the maximum deformation and compression and restitution stiffness. In addition, the principle of the unloading linear stiffness is used and two equivalent closed form solutions are derived for the unloading linear stiffness. Moreover, a closed form solution of the maximum force is derived based on the unloading stiffness closed form solutions. The final and unloading indentation properties of the new model are investigated by performing numerical simulations for all coefficient of restitution values ranging from elastic to plastic impacts.

Elastic deformation analysis of rotational wood-dowel welding joint system based on the variational method

Abstract This study aimed to investigate the elastic deformation of rotational wood-dowel welding joints by using the variational method, as there was no efficient tool to access the deformation behavior of welding joints perpendicular to axis orientation so far. Based on the variational method, the displacement control equation of axially-loaded dowel-welded wood joints was derived. By incorporating the boundary conditions and using an iterative approach, the parameters required in the control equation were determined. The analytical expression ϕ(r) for the deformation function of the wood substrate under elastic conditions was proposed, thereby determining the displacement and axial force of the single wood-dowel. The elasticity modulus was the primary influence parameter on the function ϕ(r). However, the welding depth had a more critical impact, and reasonable combination of parameters could lead to a better interaction between the wood dowel and substrate. The comparison between calculated results and test trials collected from the previous studies indicated that the elastic solution method could accurately estimate the pull-out ultimate bearing capacity and deformation characteristics of the welding joints. This work can provide unique perspective to understand the technology of rotational wood-dowel welding joints, and a good reference suggestion for the manufacture of multi-dowel welding joints.

Numerical and artificial neural network analysis of tunnel liner design parameters – a comparative study

ABSTRACT This paper proposes the development of an analytical-based Artificial Neural Network (ANN) model for predicting the concrete segmental lining parameters and comparing it with the finite element method, Rocscience tool. To achieve this aim, the dataset was prepared for twelve different loading conditions using the analytical equations by Muir Wood and Curtis. Using this dataset, an ANN model was developed for predicting the segmental lining parameters, viz. displacement, bending moment, axial and shear force. Furthermore, the Rocscience tool was used to analyse the aforementioned design parameters, and the results of these two methods were compared and summarised. It is found that the developed ANN model shows reasonable results of liner design parameters with the Rocscience program. The efficiency of both techniques was evaluated using the statistical measures viz. Mean Absolute Error, Root Mean Square Error, and Coefficient of Efficiency. It is observed that the ANN model is best fitted for the prediction of liner shear and axial force, while it could not find the optimised solution for liner displacement and bending moment. The finding of this study showed that the developed ANN model can be applied to predict the tunnel liner design parameters instead of the alternative complicated and time-consuming techniques.

Phase-field fracture modelling of piezoelectric quasicrystals

In this paper, an electro-elastic phase field model is developed for brittle fracture in piezoelectric quasicrystals. The elastic field of quasicrystals is composite of a phonon field and a phason field. The phonon field, phason field and electric field are coupled with each other, and the phonon and electric loadings trigger the variation of both the phonon and phason elastic energies, which then together contribute to crack evolution. The cracks initiation and propagation in piezoelectric quasicrystals under quasi-static mechanical loading and electric loading can be well predicted. Several numerical examples are performed to evaluate the effects of electric field and anisotropic fracture toughness. Simulation results show that the electric field influences the initial force, peak force, fracture displacement and crack path. For the single-edge notched tensile test, the effects of electric field force on the initial force, peak force and fracture displacement are obvious. For the single-edge notched shear test, the electric field force has no effect on the initial force, but it affects the peak force, fracture displacement and crack path. For the porous plate tensile test, a high-intensity electric field can deflect crack substantially resulting in passing through different holes. Further simulation indicates that the anisotropic fracture toughness has a significant effect on peak force, fracture displacement, and crack path. The present model can be used to simulate various fracture problems of piezoelectric quasicrystals.