3D Nonlinear Parallel FEM for Analyzing Dynamic Response of the Large- Scale Saturated Soil Layers-Civil Structures Interaction Problem (Part I: Formulation and its Numerical Solution)

Jafril Tanjung

doi:10.5614/jts.2010.17.2.1

Authors

Jafril Tanjung Civil Engineering Department, Engineering Faculty, Andalas University, Padang

DOI:

https://doi.org/10.5614/jts.2010.17.2.1

Keywords:

3D nonlinear FEM, Parallel computation, Domain decomposition method, Saturated soil layers-civil Structure interaction problem.

Abstract

Abstract. In this analytical study a 3D nonlinear parallel FEM algorithm was developed to overcome the physical limitation of a single processor computer when analysing the dynamic response of interaction between a large-scale saturated soil layers and civil structures. The parallel algorithm was derived based on the Domain Decomposition Method (DDM) such that a boundary value problem can be converted into an interface problem. The DDM was used to partite a whole of analytical domain into several non-overlapping subdomains. One feature of the current algorithm is most parts of the computation processes are performed in subdomain basis without interprocessor communication. The interprocessor communication is only take place when solving the interface problem. An iterative
solver of the Conjugate Gradient (CG) algorithm incorporating with the interprocessor communication based on hypercube networking was used to solve it interface problem.

Abstrak. Dalam studi analitik ini, algoritma paralel elemen hingga non-linier dikembangkan untuk mengatasi kelemahan dari kemampuan komputer berprosesor tunggal dalam menganalisis tanggap dinamik pada masalah interkasi antara lapisan tanah jemuh dan struktur sipil dalam skala besar. Algoritma paralel dijabarkan berdasarkan motode dekomposisi domain untuk mengubah masalah nilai batas ke masalah interfes. Metode dekomposisi domain digunakan untuk membagi keseluruhan domain yang dinalisis menjadi beberapa subdomain yang tidak saling overlap. Keungulan dari algiritma yang dikembangkan ini adalah sebagian besar proses perhitungan dilakukan dalam basis subdomain tanpa adanya komunikasi antar prosesor. Komunikasi antar prosesor hanya dilakukan ketika menyelesaikan masalah interfes. Metode iteratif Conjugate Gradient, yang digabungkan dengan komunikasi antar prosesor berbasiskan jejaring hypercube, digunakan untuk menyelesaikan masalah interfes tersebut.

References

Baker, L., Smith, B.J., 1987, Parallel Programming, New York: McGraw-Hill.

Crouch, R.S. and Wolf, J.P., 1994, Unified 3D Critical State Bounding Surface Plasticity Model for Soils incorporating Continues Plasticity Loading Under Cyclic Path, Constitutive Relations, Int. Journal. Num. Anal. Met., Geomech. Vol. 18, pp, 735-758.

Dafalias, Y.F., Bounding Surface Plasticity I, 1986, Mathematical Foundation and Hypoplasticity, Journal of Eng. Mech. Geomech, Vol. 112, No. 12, pp, 966-987.

Farhat, C. and Roux, F.X., 1991, A Method of Finite Element Tearing and Interconnecting and Its Parallel Solution Algorithm, Int. Journal for Num. Methods Engg., Vol. 32, pp, 1205-1227.

Ghaboussi, A.M. et al., 1973, Seismic Analysis of Earth Dam Reservoir System, J. Soil Mech. And Found. Div., ASCE., Vol. 99, No. SM10, pp. 849-862.

Gill, P.E. and Murray, W., 1974, Numerical Methods for Constrained Optimization, New York: Academic Press.

Hughes, T.J.R., 1996, The Finite Element Method Linear Static and Dynamic Finite Element Analysis, New Jersey: Prentice Hall.

Kawamura, M. and Tanjung, J., 2000, Parallel FEM Analysis for Dynamic Responses of Saturated Soil Layer, Journal of Structural Engineering. JSCE, Vol. 46A, pp. 703-712.