An official website of the United States government
We present a new method to obtain dynamic body force at virtual interfaces to reconstruct shear wave motions induced by a source outside a truncated computational domain. Specifically, a partial differential equation (PDE)-constrained optimization method is used to minimize the misfit between measured motions at a limited number of sensors on the ground surface and their counterparts reconstructed from optimized forces. Numerical results show that the optimized forces accurately reconstruct the targeted ground motions in the surface and the interior of the domain. The proposed optimization framework yields a particular force vector among other valid solutions allowed by the domain reduction method (DRM). Per this optimized or inverted force vector, the reconstructed wave field is identical to its reference counterpart in the domain of interest but may differ in the exterior domain from the reference one. However, we remark that the inverted solution is valid and introduce a simple post-process that can modify the solution to achieve an alternative force vector corresponding to the reference wave field. We also study the desired sensor spacing to accurately reconstruct the wave responses for a given dominant frequency of interest. We remark that the presented method is omnidirectionally applicable in terms of the incident angle of an incoming wave and is effective for any given material heterogeneity and geometry of layering of a reduced domain. The presented inversion method requires information on the wave speeds and dimensions of only a reduced domain. Namely, it does not need any informa- tion on the geophysical profile of an enlarged domain or a seismic source profile outside a reduced domain. Thus, the computational cost of the method is compact even though it leads to the high-fidelity reconstruction of wave re- sponse in the reduced domain, allowing for studying and predicting ground and structural responses using real seismic measurements.
more »
« less
Convolutional neural network for identifying effective seismic force at a DRM layer for rapid reconstruction of SH ground motions
Abstract We introduce a novel data‐informed convolutional neural network (CNN) approach that utilizes sparse ground motion measurements to accurately identify effective seismic forces in a truncated two‐dimensional (2D) domain. Namely, this paper presents the first prototype of a CNN capable of inferring domain reduction method (DRM) forces, equivalent to incident waves, across all nodes in the DRM layer. It achieves this from sparse measurement data in a multidimensional setting, even in the presence of incoherent incident waves. The method is applied to shear (SH) waves propagating into a domain truncated by a wave‐absorbing boundary condition (WABC). By effectively training the CNN using input‐layer features (surface sensor measurements) and output‐layer features (effective forces at a DRM layer), we achieve significant reductions in processing time compared to PDE‐constrained optimization methods. The numerical experiments demonstrate the method's effectiveness and robustness in identifying effective forces, equivalent to incoherent incident waves, at a DRM layer.
more »
« less
- Award ID(s):
- 2053694
- PAR ID:
- 10477518
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Earthquake Engineering & Structural Dynamics
- Volume:
- 53
- Issue:
- 2
- ISSN:
- 0098-8847
- Format(s):
- Medium: X Size: p. 894-923
- Size(s):
- p. 894-923
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Consider the elastic scattering of a plane or point incident wave by an unbounded and rigid rough surface. The angular spectrum representation (ASR) for the time-harmonic Navier equation is derived in three dimensions. The ASR is utilized as a radiation condition to the elastic rough surface scattering problem. The uniqueness is proved through a Rellich-type identity for surfaces given by uniformly Lipschitz functions. In the case of flat surfaces with local perturbations, an equivalent variational formulation is deduced in a truncated bounded domain and the existence of solutions are shown for general incoming waves. The main ingredient of the proof is the radiating behavior of the Green tensor to the first boundary value problem of the Navier equation in a half-space.more » « less
-
We present a new convolutional neural network (CNN)-based element-wise classification method to detect a random number of voids in a 2D plain strain solid subjected to elastodynamics. We consider that an elastic wave source excites the solid including a random number of voids, and wave responses are measured by sensors placed around the solid. We present a CNN for resolving the inverse problem, which is formulated as an element-wise classification problem. The CNN is trained to classify each element into a regular or void element from measured wave signals. To this end, we generate training data consisting of input-layer features (i.e., measured wave signals at sensors) and output-layer features (i.e., element types of all elements). When the training data are generated, we utilize the level-set method to avoid an expensive re-meshing process, which is otherwise needed for each different configuration of voids. We also analyze how effectively the CNN performs on blind test data from a non-level-set wave solver that explicitly models the boundary of voids using an unstructured, fine mesh. Numerical results show that the suggested approach can detect the locations, shapes, and sizes of multiple elliptical and circular voids in the 2D solid domain in the test data set as well as a blind test data set.more » « less
-
Sparsification of neural networks is one of the effective complexity reduction methods to improve efficiency and generalizability. Binarized activation offers an additional computational saving for inference. Due to vanishing gradient issue in training networks with binarized activation, coarse gradient (a.k.a. straight through estimator) is adopted in practice. In this paper, we study the problem of coarse gradient descent (CGD) learning of a one hidden layer convolutional neural network (CNN) with binarized activation function and sparse weights. It is known that when the input data is Gaussian distributed, no-overlap one hidden layer CNN with ReLU activation and general weight can be learned by GD in polynomial time at high probability in regression problems with ground truth. We propose a relaxed variable splitting method integrating thresholding and coarse gradient descent. The sparsity in network weight is realized through thresholding during the CGD training process. We prove that under thresholding of L1, L0, and transformed-L1 penalties, no-overlap binary activation CNN can be learned with high probability, and the iterative weights converge to a global limit which is a transformation of the true weight under a novel sparsifying operation. We found explicit error estimates of sparse weights from the true weights.more » « less
-
In this paper, we develop a sparse grid stochastic collocation method to improve the computational efficiency in handling the steady Stokes-Darcy model with random hydraulic conductivity. To represent the random hydraulic conductivity, the truncated Karhunen-Loève expansion is used. For the discrete form in probability space, we adopt the stochastic collocation method and then use the Smolyak sparse grid method to improve the efficiency. For the uncoupled deterministic subproblems at collocation nodes, we apply the general coupled finite element method. Numerical experiment results are presented to illustrate the features of this method, such as the sample size, convergence, and randomness transmission through the interface.more » « less
