Displacement Vector Research Articles

Concrete Armour Unit Breakwater Physical Model Monitoring With 3D Modeling Tools

SEABIM® is a patented scan to BIM process that can generate a reliable and complete 3D model of a rubble-mound breakwater precast protective layer (Xbloc®, Accropode™, Core-loc™ etc.). Using computer vision techniques, the known 3D shape of the Concrete Armour Units (CAU) is detected in a high-resolution point cloud, which allows to obtain the position and orientation of each unit relative to the cloud. By superimposing 3D models produced from different scans at different dates, the movement of each block can then be quantified and represented as a vector. For natural riprap armours, a point cloud segmentation algorithm is applied. Each rock is identified, and then geometric characteristics of the armour can be computed (apparent rock diameters, placement density), and the movement between each scan can be monitored. This tool has been applied to numerous real-scale projects since 2019, both in the monitoring of newly built infrastructures and in the asset management phase (Aberdeen South Harbour, Calais port 2015, Nouvelle Route du Littoral etc...). Breakwater physical models in laboratories, in wave flumes or basins, also require monitoring of the movement of reduced-scale blocks between each wave series. With SEABIM® this movement can be computed with accuracy for each CAU, allowing for a consistent evaluation of the armour layer response to wave loads. The automatic block detection process can be launched on the point cloud from a survey of the initial physical model, obtained either by photogrammetry or LiDAR. Then a new survey is done after each wave series in order to create a sequence of 3D models. Those models can be compared by computing the displacement vectors of the centers of gravity of the units. The vectors can be visualized using color graded arrows within the 3D model. All the information can then be exported to a spreadsheet for further analysis. Using this technique, settlement, rocking and sliding effects can be identified. In this paper we will discuss in detail the adaptation of the survey techniques from real scale to laboratory projects (difference in scanning technique, scale effect on point cloud density etc.). Then we will explain the 3D modeling process and the applications of the 3D models for engineering analysis purposes.

Investigation of Tunnel Collapse Mechanisms through Numerical Analysis

In this study, tunnel excavations were introduced in two typical ground conditions found in unstable urban areas, and the corresponding tunnel collapse mechanisms were investigated using numerical analysis. The case studies of tunnel collapse in weathered rock and uncemented sandy soil layers revealed that a plastic area extended across the entire tunnel section from the tunnel shoulder to the ground surface. The development of the displacement vector indicated the possibility of a tunnel ceiling depression occurring around the plastic area. The numerical analysis results suggest that tunnel excavation, which accelerates bedrock weathering or induces ground relaxation following groundwater introduction, can lead to the collapse of uncemented sandy soils.

Analysis of lateral bearing mechanism for offshore monopiles in sand and development of an analytical method

Revealing lateral soil bearing mechanism around an offshore monopile is critical for understanding and solving soil-pile interactions. In this paper, the lateral bearing mechanism around an offshore monopile in sand is analyzed and a corresponding analytical method is established. Three soil bearing mechanisms, including wedge-type flow, triaxial extension and volume compression, near mudline under working load level are proposed by analyzing the distribution of the Lode's angle (θL) and soil displacement vectors. A corresponding analytical method to calculate the lateral response of the monopile is developed based on the strain wedge (SW) model. The proposed model is verified through three pile tests. Finally, the effect of the parameters of Duncan-Chang model on the monopiles in dry and saturated sand is discussed.

Motion correction and super-resolution for multi-slice cardiac magnetic resonance imaging via an end-to-end deep learning approach

Accurate reconstruction of a high-resolution 3D volume of the heart is critical for comprehensive cardiac assessments. However, cardiac magnetic resonance (CMR) data is usually acquired as a stack of 2D short-axis (SAX) slices, which suffers from the inter-slice misalignment due to cardiac motion and data sparsity from large gaps between SAX slices. Therefore, we aim to propose an end-to-end deep learning (DL) model to address these two challenges simultaneously, employing specific model components for each challenge. The objective is to reconstruct a high-resolution 3D volume of the heart (VHR) from acquired CMR SAX slices (VLR). We define the transformation from VLR to VHR as a sequential process of motion correction and super-resolution. Accordingly, our DL model incorporates two distinct components. The first component conducts motion correction by predicting displacement vectors to re-position each SAX slice accurately. The second component takes the motion-corrected SAX slices from the first component and performs the super-resolution to fill the data gaps. These two components operate in a sequential way, and the entire model is trained end-to-end. Our model significantly reduced inter-slice misalignment from originally 3.33±0.74 mm to 1.36±0.63 mm and generated accurate high resolution 3D volumes with Dice of 0.974±0.010 for left ventricle (LV) and 0.938±0.017 for myocardium in a simulation dataset. When compared to the LAX contours in a real-world dataset, our model achieved Dice of 0.945±0.023 for LV and 0.786±0.060 for myocardium. In both datasets, our model with specific components for motion correction and super-resolution significantly enhance the performance compared to the model without such design considerations. The codes for our model are available at https://github.com/zhennongchen/CMR_MC_SR_End2End.

Shear failure analysis of slip zone soil with different coarse particle shapes: Visualized shear test and PIV technology

Coarse particle shape in slip zone soil influences the mesoscopic structure of the soil, which in turn affects soil shear strength and failure behavior. In order to investigate the effect of particle shape on the shear characteristics of coarse–fine-grained mixed slip zone soil, three types of coarse particles (spheroidal, rounded, and angular) were selected for mixing and matching, and a total of 10 sets of medium-scale shear tests were designed for this paper. To quantify the shear deformation and failure process of slip zone soils, particle image velocimetry (PIV) technology and the hanging hammer method were used to obtain mesoscopic data of the soil (displacement vector data of soil particles and elevation data of the shear failure surface), which were used to calculate shear band thickness, shear dilatation, and roughness coefficient of the shear failure surface. The results indicate that coarse particle shape can considerably affect the macroscopic mechanical properties (internal friction angle and shear strength) and mesoscopic deformation characteristics (shear band thickness, shear dilatation, and shear surface morphology) of soils. Angular coarse particles have higher interlocking strength than spheroidal and rounded coarse particles, allowing angular coarse-grained slip zone soils to develop large shear band thickness and rough shear failure surfaces. In addition, mesoscopic damage analysis suggests that the damage rate of slip zone soils decreases with increasing coarse particle shape complexity. These findings enhance comprehension of the failure characteristics of soil-rock mixture slopes and serve as a good reference for the stability analysis of similar slopes.

First-principles calculations to investigate optical, phonon and electronic properties of quaternary sulfides SrRECuS3 (RE = La, Nd, Tm)

The structure and properties of three layered heterometallic quaternary sulfides SrLaCuS3, SrNdCuS3 and SrTmCuS3 were studied for the first time using first-principles calculations in the stoichiometric and nonstoichiometric approximations. The applied DFT-based computations were performed using a hybrid functional with the contribution of nonlocal exchange in the Hartree-Fock formalism. It was revealed that the nonstoichiometry of SrLaCuS3 and SrNdCuS3 must be considered for modeling phonon spectra, elastic properties and band gaps. The wavenumbers and types of the Raman and “silent” modes at the Г-point were determined. From the analysis of displacement vectors, the degree of participation of ions in each mode was determined. The elastic constants and elastic moduli of the reported sulfides were calculated.

Rayleigh wave propagation in centrosymmetric materials with micro-stiffness, flexoelectric and micro-inertia effects

A theoretical investigation of Rayleigh waves propagation in polarized media has been carried out using a reformulated flexoelectric theory for isotropic dielectrics with micro-inertia effect. Within this non-classical theory, the internal energy density is the functional of the strain tensor, dilatation gradient, deviatoric part of stretch gradient and rotation gradient tensors, polarization vector, and polarization gradient. The obtained system of governing equations additionally contains three material length-scale parameters to account the micro-stiffness effect, one material constant to capture the micro-inertia effect, two flexoelectric constants to describe the flexoelectric effect and three length scale parameters related to the polarization gradient. To solve the coupled governing equations, the method of Lamé-type potentials for mechanical displacement and electric polarization vectors is used. The influences of various factors such as micro-stiffness, flexoelectricity, electric quadrupoles and micro-inertia effects on the phase velocity of the Rayleigh waves in a homogeneous isotropic half-space are studied. It is found that above effects become significant with the increase of the wavenumber. This study can be important for the investigation of high frequency surface acoustic waves in dielectric materials.

Evaluation and mitigation of deformable image registration uncertainties for MRI-guided adaptive radiotherapy.

We evaluate the performance of a deformable image registration (DIR) software package in registering abdominal magnetic resonance images (MRIs) and then develop a mechanical modeling method to mitigate detected DIR uncertainties. Three evaluation metrics, namely mean displacement to agreement (MDA), DICE similarity coefficient (DSC), and standard deviation of Jacobian determinants (STD-JD), are used to assess the multi-modality (MM), contour-consistency (CC), and image-intensity (II)-based DIR algorithms in the MIM software package, as well as an in-house developed, contour matching-based finite element method (CM-FEM). Furthermore, we develop a hybrid FEM registration technique to modify the displacement vector field of each MIM registration. The MIM and FEM registrations were evaluated on MRIs obtained from 10 abdominal cancer patients. One-tailed Wilcoxon-Mann-Whitney (WMW) tests were conducted to compare the MIM registrations with their FEM modifications. For the registrations performed with the MIM-CC, MIM-MM, MIM-II, and CM-FEM algorithms, their average MDAs are 0.62±0.27, 2.39±1.30, 3.07±2.42, 1.04±0.72mm, and average DSCs are 0.94±0.03, 0.80±0.12, 0.77±0.15, 0.90±0.11, respectively. The p-values of the WMW tests between the MIM registrations and their FEM modifications are less than 0.0084 for STD-JDs and greater than 0.87 for MDA and DSC. Among the three MIM DIR algorithms, MIM-CC shows the smallest errors in terms of MDA and DSC but exhibits significant Jacobian uncertainties in the interior regions of abdominal organs. The hybrid FEM technique effectively mitigates the Jacobian uncertainties in these regions.

Multi-criteria parametric optimization of the displacement and weight of a shell of minimal surface on a circular contour consisting of two inclined ellipses under thermal and power loading with consideration of geometric nonlinearity

Thin shells are well suited to optimal design problems, and the finite element method and gradient descent method make it possible to solve inverse problems in structural and applied mechanics. The calculation process takes into account geometric nonlinearity. In modern numerical studies, a formulation taking into account geometric nonlinearity in the finite element method is used, namely, a stepwise loading procedure. It becomes necessary to take into account the relations between displacement vectors and their derivatives and strain increments. These relations help to determine the stiffness matrix of the finite element at each loading step, which leads to a qualitative study of real displacements in the structure. The geometrically nonlinear calculation of a shell of minimum surface on a circular contour consisting of two inclined ellipses in a multi-criteria parametric optimization is performed by the finite element method. The finite element method is a universal variational method that is focused on solving the most complex problems of elasticity theory and applied and structural mechanics using separate calculation complexes. Stiffness matrix - for the whole body is formed on the basis of the finite element stiffness matrix. The system of solving equations of the finite element method is formed using the Lagrange's variational principle, according to which the total potential energy P of a finite element model of a body in a state of stability and equilibrium has a minimum value. As part of the sensitivity analysis, gradients of the design variables of the structure, displacements in the form of partial derivatives of these characteristics along the design variables, and shell thickness are calculated. The sensitivity information serves as the basis for building an optimal design algorithm using the gradient descent method of the objective function. Standard multi-criteria parametric optimization allows for an average of 10% steel savings, with geometric nonlinearity increasing to 20%. At this study site, 24.4% of sheet steel was saved, which is a significant relative saving. This methodology of the authors shows high results for innovative design and production of steel thin shells in Ukraine and around the world, and also makes it possible to apply several types of optimization to one research object.

Facial Painting and 3D Stereophotogrammetric Analysis of Facial Dynamics: A Reliable Anatomical Educational Method.

(1) Background: Accurate knowledge of the dynamic anatomy of facial muscles is crucial for the use of functional and aesthetic botulinum toxin injections. We studied the reliability and relevance of facial painting as a pedagogic tool for the dynamic anatomy of facial muscles. (2) Methods: Different facial expressions were performed by a female model after a professional makeup artist applied makeup to the various facial muscles on her left hemiface. A 3D photograph was taken at the beginning and end of each movement using the VECTRA H2 Imaging System device. Cutaneous movements were visualized using displacement vectors. The correlation between the theoretical and dynamic positions of the makeup-muscle was assessed by two facial anatomy experts, thanks to a correlation scale. (3) Results: The overall average score for the 11 analyzed muscles or muscle groups was 3.36 out of 4, indicating a "strong" to "very strong" estimated correlation. There was a moderate agreement between Evaluator 1 and Evaluator 2 (ICC: 0.64; 95%CI: [0.244; 0.852]; p-value: 0.005). (4) Conclusions: The educational model with facial makeup provides an indirect but nonetheless precise and reliable representation of all facial muscles on the skin's surface. It is presented as a reliable and reproducible method, which exhibits great potential as a teaching tool.

A mathematical model of the gravitational potential of the planet taking into account tidal deformations

Objectives. This paper investigates the gravitational potential of a viscoelastic planet moving in the gravitational field of a massive attracting center (star), a satellite and one or more other planets moving in Keplerian elliptical orbits relative to the attracting center. Celestial bodies other than a viscoelastic planet are modeled by material points. Within the framework of the linear model of the theory of viscoelasticity, the problem of finding the vector of elastic displacement has been resolved. Traditionally, a solid body model is used to determine the Earth’s gravitational field, while tidal deformations are taken into account in the form of small corrections to the coefficients of the geopotential model. In this work, the viscoelastic ball model is used to take into account tidal effects. The relevance of the research topic is associated with high-precision forecasting of the movement of artificial satellites of the Earth, high-precision measurement of the Earth’s gravitational field.Methods. In this study the asymptotic and analytical methods developed by V.G. Vilke are used for mechanical systems containing viscoelastic elements of high rigidity, as well as methods of classical mechanics, mathematical analysis. The graphs were plotted using the Octave mathematical package.Results. After resolving the quasi-static problem of elasticity theory by calculating triple integrals over a spherical area, a formula for the gravitational potential of a deformable planet was obtained. In addition, the gravitational potential of the Earth was also calculated taking into account solid-state tidal effects from the Moon, Sun, and Venus at an external point. Graphs were constructed to show the dependence of the Earth’s gravitational potential on time.Conclusions. The theoretical and numerical results established herein show that the main contribution to the gravitational potential of the Earth is made by the Moon and the Sun. The influence of other planets in the solar system is small. The value of the gravitational potential at the outer point of the Earth, taking into account tidal effects, depends both on the position of the point in the moving coordinate system and on the relative position of celestial bodies.

A novel vector field analysis for quantitative structure changes after macular epiretinal membrane surgery.

The aim of this study was to introduce novel vector field analysis for the quantitative measurement of retinal displacement after epiretinal membrane (ERM) removal. We developed a novel framework to measure retinal displacement from retinal fundus images as follows: (1) rigid registration of preoperative retinal fundus images in reference to postoperative retinal fundus images, (2) extraction of retinal vessel segmentation masks from these retinal fundus images, (3) non-rigid registration of preoperative vessel masks in reference to postoperative vessel masks, and (4) calculation of the transformation matrix required for non-rigid registration for each pixel. These pixel-wise vector field results were summarized according to predefined 24 sectors after standardization. We applied this framework to 20 patients who underwent ERM removal to obtain their retinal displacement vector fields between retinal fundus images taken preoperatively and at postoperative 1, 4, 10, and 22months. The mean direction of displacement vectors was in the nasal direction. The mean standardized magnitudes of retinal displacement between preoperative and postoperative 1month, postoperative 1 and 4, 4 and 10, and 10 and 22months were 38.6, 14.9, 7.6, and 5.4, respectively. In conclusion, the proposed method provides a computerized, reproducible, and scalable way to analyze structural changes in the retina with a powerful visualization tool. Retinal structural changes were mostly concentrated in the early postoperative period and tended to move nasally.

Insights on structure-obstacle-soil interactions during the installation of monopiles for OWTs

ABSTRACT In this paper, a new three-dimensional FE model considering the structure-object-soil interaction is proposed to simulate the contact process during monopile installation. Using this model, the displacement vector diagrams and plastic strain contours of the obstacle and surrounding soil are obtained and analysed in detail to explain the soil failure mechanism and its effect on the installation resistance and moment profiles. Furthermore, the influences of the obstacle shape, size and contact position on the installation resistance of the monopile foundation are investigated. Results show that for obstacles with a dimensionless size (defined as the obstacle size to the monopile diameter) of 0.33, the resultant moment values exceed 1500kN·m, and the vertical resistance in the presence of obstacles is 0.88-1.53 times larger than that in the absence of obstacles. Further evaluation should be conducted to evaluate the effects of the resultant moments and vertical resistance on the actual installation process.

A model that predicts a real-time tumour surface using intra-treatment skin surface and end-of-expiration and end-of-inhalation planning CT images.

To develop a mapping model between skin surface motion and internal tumour motion and deformation using end-of-exhalation (EOE) and end-of-inhalation (EOI) 3D CT images for tracking lung tumours during respiration. Before treatment, skin and tumour surfaces were segmented and reconstructed from the EOE and the EOI 3D CT images. A non-rigid registration algorithm was used to register the EOE skin and tumour surfaces to the EOI, resulting in a displacement vector field that was then used to construct a mapping model. During treatment, the EOE skin surface was registered to the real-time, yielding a real-time skin surface displacement vector field. Using the mapping model generated, the input of a real-time skin surface can be used to calculate the real-time tumour surface. The proposed method was validated with and without simulated noise on 4D CT images from 15 patients at Léon Bérard Cancer Center and the 4D-lung dataset. The average centre position error, dice similarity coefficient (DSC), 95%-Hausdorff distance and mean distance to agreement of the tumour surfaces were 1.29 mm, 0.924, 2.76 mm, and 1.13 mm without simulated noise, respectively. With simulated noise, these values were 1.33 mm, 0.920, 2.79 mm, and 1.15 mm, respectively. A patient-specific model was proposed and validated that was constructed using only EOE and EOI 3D CT images and real-time skin surface images to predict internal tumour motion and deformation during respiratory motion. The proposed method achieves comparable accuracy to state-of-the-art methods with fewer pre-treatment planning CT images, which holds potential for application in precise image-guided radiation therapy.

Retrieval of sea ice drift in the Fram Strait based on data from Chinese satellite HaiYang (HY-1D)

Abstract. Melting of sea ice in the Arctic has accelerated due to global warming. The Fram Strait (FS) serves as a crucial pathway for sea ice export from the Arctic to the North Atlantic Ocean. Monitoring sea ice drift (SID) in the FS provides insight into how Arctic sea ice responds to the climate change. The SID has been retrieved from Sentinel-1 synthetic aperture radar (SAR), Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS), and Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and further exploration is needed for the retrieval of SID using optical imagery. In this paper, we retrieve SID in the FS using the Chinese HaiYang1-D (HY-1D) satellite equipped with the Coastal Zone Imager (CZI). A multi-template matching technique is employed to calculate the cross-correlation, and subpixel estimation is used to locate displacement vectors from the cross-correlation matrix. The dataset covering March to May 2021 was divided into hourly and daily intervals for analysis, and validation was performed using Copernicus Marine Environment Monitoring Service (CMEMS) SAR-based product and International Arctic Buoy Programme (IABP) buoy. A comparison with the CMEMS SID product revealed a high correlation with the daily interval dataset; however, due to the spatial and temporal variability in the sea ice motion, differences are observed with the hourly interval dataset. Additionally, validation with the IABP buoys yielded a velocity bias of −0.005 m s−1 and RMSE of 0.031 m s−1 for the daily interval dataset, along with a flow direction bias of 0.002 rad and RMSE of 0.009 rad, respectively. For the hourly interval dataset, the velocity bias was negligible (0 m s−1), with a RMSE of 0.036 m s−1, while the flow direction bias was 0.003 rad, with a RMSE of 0.010 rad. In addition, during the validation with buoys, we found that the accuracy of retrieving the SID flow direction is distinctly interrelated with the sea ice displacement.

Derivation of 3D Coseismic Displacement Field from Integrated Azimuth and LOS Displacements for the 2018 Hualien Earthquake

A 3D surface deformation field for an earthquake can aid in understanding fault behaviors and earthquake mechanisms. However, SAR-based 3D surface deformation estimates are often limited by insufficient observations and hampered by various error sources. In this study, we demonstrate the derivation of a 3D coseismic displacement field from different InSAR processing algorithms. The azimuth displacements from Multiple Aperture Interferometry (MAI) and Pixel Offset Tracking (POT) were integrated to ensure reliable displacements at low coherent areas. The 3D displacement field was inverted pixel-by-pixel by Line-of-Sight (LOS) displacement and integrated azimuth displacement. The results showed that MAI and POT could compensate for the weaknesses of each algorithm. Also, pixels with less than three sets of observations showed higher noise levels. Such noisy pixels were removed by a denoising criterion proposed herein. For the vertical direction, the proportion of pixels inverted with two sets of azimuth and one set of LOS displacements was 26.1%. After denoising, the proportion dropped to 2.4% due to the insufficiency of LOS displacements. This shows that the viewing angle influences the overall performance of 3D surface displacement inversion. Implementing various displacement vectors should reduce such limitations.

Planar defects in α-Ga2O3 thin films produced by HVPE

The defect structure of α-phase gallium oxide thin films was investigated using transmission electron microscopy (TEM). Epitaxial Ga2O3 films were grown via halide vapor-phase epitaxy on c-plane sapphire substrates. TEM analysis revealed a high density of extended planar defects within the films, primarily located along prismatic planes of {112¯0} type. Displacement vectors were determined using the invisibility criterion for stacking faults. The study encompassed both planar and cross-sectional views of the films. It is hypothesized that these defects form due to the motion of edge partial dislocations with the 13⟨11¯00⟩ Burgers vector. Various mechanisms of their formation have been explored.

Investigating micro-cutting mechanism of glow discharge polymer based on material properties and removal behaviors analysis

Glow discharge polymer (GDP) is the unique artificial material for the target balls in Inertial Confinement Fusion tests, while its practical micromachining such as the fabrication of microstructure on material surfaces keeps challenging due to its unclear micro-cutting mechanism. Hence, this paper seeks to investigate the micro-cutting mechanism of GDP from the perspective of the material flow and cutting energy. To achieve it, the energy conservation equation of three cutting modes triggered by different ratios of uncut chip thickness to the tool cutting edge radius (RTS) was established based on cutting deformation behaviors. Meanwhile, the diamond cutting tests and the FEM simulation at different RTS were developed. The experimental observation of cutting forces and specific cutting forces verified the evolution of three cutting modes, including shearing, shearing-ploughing and ploughing in the micro-cutting of GDP with the decrease of RTS. Next, from the change of the node displacement vector observed from the simulated results, it can be seen that the real-time material flow behavior during micro-cutting of GDP varies obviously with the evolution of cutting modes. Besides, the fracture toughness Gc, and the energy dissipation of different cutting modes were analyzed. The proportion of the energy spent on material fracture (Gc=9.95 N/mm) is the largest one in shearing and shearing-ploughing modes, while in ploughing mode, the material plastic deformation consumed the most energy. The above results reveal the specific material properties and removal behaviors of GDP and contribute to optimizing the machining strategies for the practical micromachining of microstructures on material surfaces.

Fractional-order velocity-dilatation-rotation viscoelastic wave equation and numerical solution based on a constant-Q model

Viscoelastic wave equations based on the constant- Q (CQ) model can accurately describe the amplitude dissipation and phase distortion of waves in anelastic media. However, only three velocity or displacement components can be obtained directly by solving such equations. Starting from the time-domain second-order displacement viscoelastic wave equation, we derived the decoupled P- and S-wave displacement vector viscoelastic wave equation by using the polarization difference of P- and S-wave propagation in isotropic media. The equation can be transformed into the velocity-dilatation-rotation viscoelastic wave equation containing the first-order temporal derivative and fractional Laplacian operators, which can be solved directly by using the staggered-grid finite-difference and pseudospectral methods. We use the low-rank decomposition method to approximate the derived mixed space-wavenumber domain fractional Laplacian operators for modeling wave propagation in heterogeneous attenuating media. We also demonstrated the precision of our equation by comparing the numerical solutions with the analytical solutions. Furthermore, compared with the conventional velocity-stress viscoelastic wave equation, experimental results demonstrate that our equation can separate the pure P and S waves from the mixed wavefield during wavefield continuation. In addition, it can be separated into an equation containing predominantly an amplitude attenuation or a phase distortion term.

An alternative potential method for mixed steady‐state elastic oscillation problems

We consider an alternative approach to investigate three‐dimensional exterior mixed boundary value problems (BVPs) for the steady‐state oscillation equations of the elasticity theory for isotropic bodies. The unbounded domain occupied by an elastic body, , has a compact boundary surface , which is divided into two disjoint parts, the Dirichlet part and the Neumann part , where the displacement vector (the Dirichlet‐type condition) and the stress vector (the Neumann‐type condition) are prescribed, respectively. Our new approach is based on the classical potential method and has several essential advantages compared with the existing approaches. We look for a solution to the mixed BVP in the form of a linear combination of the single‐layer and double‐layer potentials with densities supported on the Dirichlet and Neumann parts of the boundary, respectively. This approach reduces the mixed BVP under consideration to a system of boundary integral equations, which contain neither extensions of the Dirichlet or Neumann data nor the Steklov–Poincaré‐type operator involving the inverse of a special boundary integral operator, which is not available explicitly for arbitrary boundary surface. Moreover, the right‐hand sides of the resulting boundary integral equations system are vector functions coinciding with the given Dirichlet and Neumann data of the problem in question. We show that the corresponding matrix integral operator is bounded and coercive in the appropriate ‐based Bessel potential spaces. Consequently, the operator is invertible, which implies unconditional unique solvability of the mixed BVP in the class of vector functions belonging to the Sobolev space and satisfying the Sommerfeld–Kupradze radiation conditions at infinity. We also show that the obtained matrix boundary integral operator is invertible in the ‐based Besov spaces and prove that under appropriate boundary data a solution to the mixed BVP possesses ‐Hölder continuity property in the closed domain with , where is an arbitrarily small number.