A Newton-CG Augmented Lagrangian-Based Finite-Element Method for Three-Dimensional Restoration of Geological Models
Publication: International Journal of Geomechanics
Volume 17, Issue 9
Abstract
In this paper, a Newton conjugate gradient (Newton-CG) augmented Lagrangian-based finite-element method is proposed to solve elastic contact problems without friction. The contact problem is discretized by the finite-element method. Based on a minimum theorem, the contact problem is reformulated as an optimization problem with inequality constraints. Then, the Newton-CG augmented Lagrangian method is applied to solve the constrained optimization problem. Numerical experiments were performed to validate the efficiency of the proposed method.
Get full access to this article
View all available purchase options and get full access to this article.
Acknowledgments
The authors sincerely appreciate the anonymous referees and the editor for their very helpful comments and suggestions, which greatly improved the quality of this paper. This work is supported by the Natural Science Foundation of Zhejiang Province of China under Grant LQ13A010006, the National Natural Science Foundation of China under Grant 11201382, the Youth Fund of Humanities and Social Sciences of the Ministry of Education under Grant 12YJC910008, the project of the science and technology department of Sichuan Province under Grant 2012ZR0154, and the National Major Scientific and Technological Special Project during the Thirteenth Five-year Plan Period 2016ZX05003-001 and 2017ZX05008-001.
References
Annavarapu, C., Hautefeuille, M., and Dolbow, J. E. (2013). “A Nitsche stabilized finite element method for frictional sliding on embedded interfaces. Part II: Intersecting interfaces.” Comput. Methods Appl. Mech. Eng., 267, 318–341.
Chabrand, P., Dubois, F., and Raous, M. (1998). “Various numerical methods for solving unilateral contact problems with friction.” Math. Comput. Modell., 28(4–8), 97–108.
Cheng, Y. M., Chen, W. S., and Zhang, Y. H. (2006). “New approach to determine three-dimensional contacts in blocks system: Penetration edges method.” Int. J. Geomech., 303–310.
Chouly, F., Hild, P., and Renard, Y. (2015). “Symmetric and non-symmetric variants of Nitsche’s method for contact problems in elasticity: Theory and numerical experiments.” Math. Comput., 84(293), 1089–1112.
Christensen, P. W., Klarbring, A., Pang, J. S., and Strömberg, N. (1998). “Formulation and comparison of algorithms for frictional contact problems.” Int. J. Numer. Methods Eng., 42(1), 145–173.
Dahlstrom, C. D. A. (1969). “Balanced cross sections.” Can. J. Earth Sci., 6(4), 743–757.
De Santi, M. R., Campos, J. L., and Martha, L. F. (2002). “A finite element approach for geological section reconstruction.” Proc., 21st GOCAD Meeting, Univ. of Nancy, Nancy, France.
Dezfooli, M. S., Abolmaali, A., and Razavi, M. (2015). “Coupled nonlinear finite-element analysis of soil-steel pipe structure interaction.” Int. J. Geomech., 04014032.
Fazel, M., Pong, T. K., Sun, D., and Tseng, P. (2013). “Hankel matrix rank minimization with applications to system identification and realization.” SIAM J. Matrix Anal. Appl., 34(3), 946–977.
Fichera, G. (1972). “Boundary value problems of elasticity.” Encyclopedia of physics, Vol.VIa/2, C. Truesdell, ed., Springer, Berlin.
Freund, R. W., and Nachtigal, N. M. (1994). “A new Krylov subspace method for symmetric indefinite linear systems.” Rep. ORNL/TM-12754, Oak Ridge National Laboratory, Oak Ridge, TN.
Hestenes, M. R. (1969). “Multiplier and gradient methods.” J. Optim. Theory Appl., 4(5), 303–320.
Hild, P., and Renard, Y. (2010). “A stabilized Lagrange multiplier method for the finite element approximation of contact problems in elastostatics.” Numer. Math., 115(1), 101–129.
Hiriart-Urruty, J.-B., Strodiot, J.-J., and Nguyen, V. H. (1984). “Generalized Hessian matrix and second-order optimality conditions for problems with C1,1 data.” Appl. Math. Opt., 11(1), 43–56.
Kikuchi, N., and Oden, J. T. (1988). Contact problems in elasticity: A study of variational inequalities and finite element methods, SIAM, Philadelphia.
Kinderleher, D., and Stampacchia, G. (1980). An introduction to variational inequalities and their application, Academic, New York.
Liu, Y.-J., Sun, D., and Toh, K.-C. (2012). “An implementable proximal point algorithmic framework for nuclear norm minimization.” Math. Program., 133(1), 399–436.
Maerten, F., and Maerten, L. (2015). “On a method for reducing interpretation uncertainty of poorly imaged seismic horizons and faults using geomechanically based restoration technique.” Interpretation, 3(4), SAA105–SAA116.
Maerten, L., and Maerten, F. (2006). “Chronologic modeling of faulted and fractured reservoirs using geomechanically based restoration: Technique and industry applications.” AAPG Bull., 90(8), 1201–1226.
MATLAB R2011a [Computer software]. MathWorks, Natick, MA.
Muron, P., and Mallet, J. L. (2003). “3D balanced unfolding: The tetrahedral approach.” 23rd GOCAD Meeting, Univ. of Nancy, Nancy, France.
Nitsche, J. (1971). “Über ein variationsprinzip zur lösung von dirichlet-problemen bei verwendung von teilräumen, die keinen randbedingungen unterworfen sind.” Abh. Math. Sem. Univ. Hamburg, 36(1), 9–15.
Powell, M. J. D. (1969). “A method for nonlinear constraints in minimization problems.” Optimization, R. Fletcher, ed., Academic, London.
Rockafellar, R. T. (1976). “Augmented Lagrangians and applications of the proximal point algorithm in convex programming.” Math. Oper. Res., 1(2), 97–116.
Schöberl, J. (1998). “Solving the Signorini problem on the basis of domain decomposition techniques.” Comput., 60(4), 323–344.
Schöberl, J. (2001). “Efficient contact solvers based on domain decomposition techniques.” Comput. Math. Appl., 42(8–9), 1217–1228.
Sun, D., Sun, J., and Zhang, L. (2008). “The rate of convergence of the augmented Lagrangian method for nonlinear semidefinite programming.” Math. Program., 114(2), 349–391.
Suppe, J. (1983). “Geometry and kinematics of fault-bend folding.” Am. J. Sci., 283(7), 684–721.
Tang, P., Wang, C., and Dai, X. (2016). “A majorized Newton-CG augmented Lagrangian-based finite element method for 3D restoration of geological models.” Comput. Geosci., 89, 200–206.
TetGen [Computer software]. Hang Si Research Group, Numerical Mathematics and Scientific Computing, Weierstrass Institute for Applied Analysis and Stochastics (WIAS), Berlin.
Wang, C., Sun, D., and Toh, K.-C. (2010). “Solving log-determinant optimization problems by a Newton-CG primal proximal point algorithm.” SIAM J. Optim., 20(6), 2994–3013.
Wriggers, P. (2006). Computational contact mechanics, 2nd Ed., Springer, Berlin.
Yang, J., Sun, D., and Toh, K.-C. (2013). “A proximal point algorithm for log-determinant optimization with group lasso regularization.” SIAM J. Optim., 23(2), 857–893.
Zhao, X.-Y., Sun, D., and Toh, K.-C. (2010). “A Newton-CG augmented Lagrangian method for semidefinite programming.” SIAM J. Optim., 20(4), 1737–1765.
Information & Authors
Information
Published In
International Journal of Geomechanics
Volume 17 • Issue 9 • September 2017
Copyright
© 2017 American Society of Civil Engineers.
History
Received: Apr 13, 2015
Accepted: Feb 3, 2017
Published online: Apr 28, 2017
Published in print: Sep 1, 2017
Discussion open until: Sep 28, 2017
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.