skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on May 5, 2026

Title: A Direct Solution to the Self-Sensing Inverse Problem Via the Primal-Dual Interior Point Method
Abstract Materials that are self-sensing via the piezoresistive effect have been widely explored for embedded sensing in aerospace structural composites, pressure sensing in touch pads, diagnostics in biomedical implant technology, and many other applications. In this approach, changes in electrical transport of the material are used as an indicator of stress/strain, pressure, or damage. However, engineers and other users of self-sensing materials are typically not directly interested in the electrical state of the material. Rather, they want to know the underlying mechanical state of the material that gives rise to an observed electrical change. Recovering material condition from electrical observations is referred to as the self-sensing inverse problem (SSIP). Prior work by the author has shown that the SSIP can be solved using an electrical impedance tomography (EIT)-generated conductivity map as an input, but this is undesirable because it requires solving a second inverse problem (i.e., the SSIP) on top of the EIT inverse problem. To that end, a direct formulation for the SSIP is herein presented. In this approach, voltage-current data is directly inverted to find the displacement field without using EIT as an intermediate step. Additionally, the direct SSIP is solved within the primal-dual interior point (PDIPM) framework such that an ℓ1-norm can be used on the regularization term, which promotes sparsity in the solution space. This approach is applied to a representative piezoresistive nanocomposite using the Laplace matrix as regularization to promote a spatially smooth solution. From the displacement field, strains and stresses are calculated and compared to a commercial finite element solution with good accuracy.  more » « less
Award ID(s):
2239039
PAR ID:
10616946
Author(s) / Creator(s):
Publisher / Repository:
American Society of Mechanical Engineers
Date Published:
ISBN:
978-0-7918-8875-9
Format(s):
Medium: X
Location:
Houston, Texas, USA
Sponsoring Org:
National Science Foundation
More Like this
  1. We present iSoft, a single volume soft sensor capable of sensing real-time continuous contact and unidirectional stretching. We propose a low-cost and an easy way to fabricate such piezoresistive elastomer-based soft sensors for instant interactions. We employ an electrical impedance tomography (EIT) technique to estimate changes of resistance distribution on the sensor caused by fingertip contact. To compensate for the rebound elasticity of the elastomer and achieve real-time continuous contact sensing, we apply a dynamic baseline update for EIT. The baseline updates are triggered by fingertip contact and movement detections. Further, we support unidirectional stretching sensing using a model-based approach which works separately with continuous contact sensing. We also provide a software toolkit for users to design and deploy personalized interfaces with customized sensors. Through a series of experiments and evaluations, we validate the performance of contact and stretching sensing. Through example applications, we show the variety of examples enabled by iSoft. 
    more » « less
  2. A comparative study is presented to solve the inverse problem in elasticity for the shear modulus (stiffness) distribution utilizing two constitutive equations: (1) linear elasticity assuming small strain theory, and (2) finite elasticity with a hyperelastic neo-Hookean material model. Assuming that a material undergoes large deformations and material nonlinearity is assumed negligible, the inverse solution using (2) is anticipated to yield better results than (1). Given the fact that solving a linear elastic model is significantly faster than a nonlinear model and more robust numerically, we posed the following question: How accurately could we map the shear modulus distribution with a linear elastic model using small strain theory for a specimen undergoing large deformations? To this end, experimental displacement data of a silicone composite sample containing two stiff inclusions of different sizes under uniaxial displacement controlled extension were acquired using a digital image correlation system. The silicone based composite was modeled both as a linear elastic solid under infinitesimal strains and as a neo-Hookean hyperelastic solid that takes into account geometrically nonlinear finite deformations. We observed that the mapped shear modulus contrast, determined by solving an inverse problem, between inclusion and background was higher for the linear elastic model as compared to that of the hyperelastic one. A similar trend was observed for simulated experiments, where synthetically computed displacement data were produced and the inverse problem solved using both, the linear elastic model and the neo-Hookean material model. In addition, it was observed that the inverse problem solution was inclusion size-sensitive. Consequently, an 1-D model was introduced to broaden our understanding of this issue. This 1-D analysis revealed that by using a linear elastic approach, the overestimation of the shear modulus contrast between inclusion and background increases with the increase of external loads and target shear modulus contrast. Finally, this investigation provides valuable information on the validity of the assumption for utilizing linear elasticity in solving inverse problems for the spatial distribution of shear modulus associated with soft solids undergoing large deformations. Thus, this work could be of importance to characterize mechanical property variations of polymer based materials such as rubbers or in elasticity imaging of tissues for pathology. 
    more » « less
  3. Nano-indentation is a promising method to identify the constitutive parameters of soft materials, including soft tissues. Especially when materials are very small and heterogeneous, nano-indentation allows mechanical interrogation where traditional methods may fail. However, because nano-indentation does not yield a homogeneous deformation field, interpreting the resulting load–displacement curves is non-trivial and most investigators resort to simplified approaches based on the Hertzian solution. Unfortunately, for small samples and large indentation depths, these solutions are inaccurate. We set out to use machine learning to provide an alternative strategy. We first used the finite element method to create a large synthetic data set. We then used these data to train neural networks to inversely identify material parameters from load–displacement curves. To this end, we took two different approaches. First, we learned the indentation forward problem, which we then applied within an iterative framework to identify material parameters. Second, we learned the inverse problem of directly identifying material parameters. We show that both approaches are effective at identifying the parameters of the neo-Hookean and Gent models. Specifically, when applied to synthetic data, our approaches are accurate even for small sample sizes and at deep indentation. Additionally, our approaches are fast, especially compared to the inverse finite element approach. Finally, our approaches worked on unseen experimental data from thin mouse brain samples. Here, our approaches proved robust to experimental noise across over 1000 samples. By providing open access to our data and code, we hope to support others that conduct nano-indentation on soft materials. 
    more » « less
  4. A path integral Monte Carlo method (PIMC) based on a Feynman-Kac formula for the Laplace equation with mixed boundary conditions is proposed to solve the forward problem of the electrical impedance tomography (EIT). The forward problem is an important part of iterative algorithms of the inverse EIT problem, and the proposed PIMC provides a local solution to find the potentials and currents on individual electrodes. Improved techniques are proposed to compute with better accuracy both the local time of reflecting Brownian motions (RBMs) and the Feynman-Kac formula for mixed boundary problems of the Laplace equation. Accurate voltage-to-current maps on the electrodes of a model 3-D EIT problem with eight electrodes are obtained by solving a mixed boundary problem with the proposed PIMC method. 
    more » « less
  5. Regularized learning problems in Banach spaces, which often minimize the sum of a data fidelity term in one Banach norm and a regularization term in another Banach norm, is challenging to solve. We construct a direct sum space based on the Banach spaces for the fidelity term and the regularization term and recast the objective function as the norm of a quotient space of the direct sum space. We then express the original regularized problem as an optimization problem in the dual space of the direct sum space. It is to find the maximum of a linear function on a convex polytope, which may be solved by linear programming. A solution of the original problem is then obtained by using related extremal properties of norming functionals from a solution of the dual problem. Numerical experiments demonstrate that the proposed duality approach is effective for solving the regularization learning problems. 
    more » « less