• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.164-175
• /
• 2011

Using natural electromagnetic (EM) fields at low frequencies, magnetotelluric (MT) surveys can investigate conductivity structures of the deep subsurface and thus are used to explore geothermal energy resources and investigate proper sites for not only geological $CO_2$ sequestration but also enhanced geothermal system (EGS). Moreover, marine MT data can be used for better interpretation of marine controlled-source EM data. In the interpretation of MT data, MT modeling schemes are important. This study improves a three dimensional (3D) MT modeling algorithm which uses edge finite elements. The algorithm computes magnetic fields by solving an integral form of Faraday's law of induction based on a finite difference (FD) strategy. However, the FD strategy limits the algorithm in computing vertical magnetic fields for a topographic model. The improved algorithm solves the differential form of Faraday's law of induction by making derivatives of electric fields, which are represented as a sum of basis functions multiplied by corresponding weightings. In numerical tests, vertical magnetic fields for topographic models using the improved algorithm overcome the limitation of the old algorithm. This study recomputes induction vectors and tippers for a 3D hill and valley model which were used for computation of the responses using the old algorithm.

• Journal of the Korean GEO-environmental Society
• /
• v.22 no.1
• /
• pp.13-20
• /
• 2021

A suite of 3D large deformation finite element (FE) analyses was performed to investigate the load transfer mechanism and penetration resistance of spudcan foundations in heterogeneous soil profile consisting of sand and clay. The Elasto-Plastic models following Mohr-Coulomb and Tresca failure criteria were adopted for sand and clay, respectively. The accuracy of the numerical model was validated against centrifuge test measurements. The dense sand behavior with dilation is modeled using the non-associated flow rule. An investigation study consisting of key parameters, which includes variation in soil stratigraphy (sand-clay, sand-clay-sand), strength parameters of sand and clay (��' and su) and normalized height ratio of the sand layer (Hs/D) was conducted to assess the penetration behavior of spudcan. Based on calculated outputs, it was demonstrated that these parameters have a significant influence on the penetration resistance of spudcan. The calculated penetration resistance profiles are compared with the published (sand overlying clay) analytical model. It is confirmed that for the case of two-layer soil, the available theoretical model provides an accurate estimate of peak penetration resistance (qpeak). In the case of three-layer soil, the presence of a third stiff layer affects the penetration resistance profile due to the squeezing of the soil.

• The korean journal of orthodontics
• /
• v.36 no.6
• /
• pp.412-421
• /
• 2006

Objective: The result of finite element analysis depends on material properties, structural expression, density of element, and boundar or loading conditions. To represent proper elastic behavior, a finite element model was made using Hounsfield unit (HU) values in CT images. Methods: A 13 year 6 month old male was used as the subject. A 3 dimensional visualizing program, Mimics, was used to build a 3D object from the DICOM file which was acquired from the CT images. Model 1 was established by giving 24 material properties according to HU. Model 2 was constructed by the conventional method which provides 2 material properties. Protraction force of 500g was applied at a 45 degree downward angle from Frankfort horizontal (FH) plane. Results: Model 1 showed a more flexible response on the first premolar region which had more forward and downward movement of the maxillary anterior segment. Maxilla was bent on the sagittal plane and frontal plane. Model 2 revealed less movement in all directions. It moved downward on the anterior part and upward on the posterior part, which is clockwise rotation of the maxilla. Conclusion: These results signify that different outcomes of finite element analysis can occur according to the given material properties and it is recommended to use HU values for more accurate results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.10
• /
• pp.1644-1652
• /
• 2003

The monitor product is packed by cushioning materials because the monitor can be broken during transportation. However, the addition of the cushioning material increased the volume of the product. Therefore, it is required that the usage of cushioning material be minimized. In practice, design engineers have followed the ad hoc design with experiences of predecessors. Automation of the design process is very important for the reduction of engineering cost, and can be achieved by an excellent design process and software development. According to Axiomatic design, a design flow is defined and a software system is developed for automated design. At first, a basic model is defined. A user can modify the model from menus and design is carried out according to the input from the user. Finite element models are automatically generated based on the design. A nonlinear finite element analysis program called LS/DYNA3D is linked for the impact analysis. The process of Design of Experiments using orthogonal array is installed to minimize the maximum acceleration in drop test. Therefore, a new design can be proposed by the system. The program is designed according to the Independence Axiom of Axiomatic design. FRs and DPs of the software system are defined and decomposed by zigzagging process. Independent modules can be generated by analysis of the full design matrix and each module is coded as class in Object Oriented Programming (OOP). Design results are discussed.

• Interaction and multiscale mechanics
• /
• v.2 no.1
• /
• pp.45-68
• /
• 2009

An aim of the study is to develop an efficient numerical simulation technique that can handle the two-scale analysis of fluid permeation filters fabricated by the partial sintering technique of small spherical ceramics. A solid-fluid mixture homogenization method is introduced to predict the mechanical characters such as rigidity and permeability of the porous ceramic filters from the micro-scale geometry and configuration of partially-sintered particles. An extended finite element (X-FE) discretization technique based on the enriched interpolations of respective characteristic functions at fluid-solid interfaces is proposed for the non-interface-fitted mesh solution of the micro-scale analysis that needs non-slip condition at the interface between solid and fluid phases of the unit cell. The homogenization and localization performances of the proposed method are shown in a typical two-dimensional benchmark problem whose model has a hole in center. Three-dimensional applications to the body-centered cubic (BCC) and face-centered cubic (FCC) unit cell models are also shown in the paper. The 3D application is prepared toward the computer-aided optimal design of ceramic filters. The accuracy and stability of the X-FEM based method are comparable to those of the standard interface-fitted FEM, and are superior to those of the voxel type FEM that is often used in such complex micro geometry cases.

• Earthquakes and Structures
• /
• v.19 no.5
• /
• pp.331-345
• /
• 2020

The multi-story timber structure with high platform base is one of the important architectural types in the traditional Chinese buildings. To study the dynamic characteristics and seismic responses on this kind of traditional structure, the 3-D finite element models of Xi'an drum tower which included the high platform base, upper timber structure and whole structure was established considering the structural form and material performance parameters of the structure in this study. By the modal analysis, the main frequencies and mode shapes of this kind of traditional building were obtained and investigated. The three kinds of earthquake excitations included El-Centro wave, Taft wave and Lanzhou wave were separately imposed on the upper timber structure model and the overall structure model, and the seismic responses on the tops of columns were analyzed. The results of time history analysis show that the seismic response of the upper timber structure is obviously amplified by high platform base. After considering the effect of high platform base, the mean value on the lateral displacement increments of the top column in the overall structure is more than 20.478% and the increase of dynamic coefficients was all above 0.818 under the above three different earthquake excitations. Obviously, it shows that the existence of high platform base has a negative influence on the seismic responses of upper timber structure. And the high platform base will directly affect the safety of the upper timber structure. Therefore, the influence of high platform base on the dynamic response of its upper timber structure cannot be neglected.

• Steel and Composite Structures
• /
• v.8 no.5
• /
• pp.403-421
• /
• 2008

The design provisions for semi-rigid steel frames have been incorporated in codes of practice for steel structures. In order to do the same, it is necessary to know the experimental moment-relative rotation (M-${\theta}_r$) behaviour of beam-to-column connections. In spite of numerous publications and collection of several connection databases, there is no unified approach for the semi-rigid design of steel frames. Amongst the many connection models available, the Frye-Morris polynomial model, with its limitations reported in the literature, is simple to adopt at least for the linear design space. However this model requires more number of connection tests and regression analyses to make it a realistic prediction model. In this paper, 3D nonlinear finite element (FE) analysis of beam-column connection specimens, carried out using ABAQUS software, for evaluating the M-${\theta}_r$ behaviour of semi-rigid top and seat-angle (TSA) bolted connections are described. The finite element model is validated against experimental behaviour of the same connection with regard to their moment-rotation behaviour, stress distribution and mode of failure of the connections. The calibrated FE model is used to evaluate the performance of the Frye-Morris polynomial model. The results of the numerical parametric studies carried out using the validated FE model have been used in proposing modifications to the Frye-Morris model for TSA connection in terms of the powers of the size parameters.

• Journal of Dental Rehabilitation and Applied Science
• /
• v.29 no.2
• /
• pp.111-118
• /
• 2013

Purpose of present study is to investigate the effects of thread pitch in coronal portion in scalloped implant with 2 different connections on loading stress using 3 dimensional finite element analysis. Scalloped implant with 4 different thread pitches (0.4mm, 0.5mm, 0.6, and 0.7mm) in the coronal part was modeled with 2 different implant-abutment connections. Platform matching connection had the same implant and abutment diameter so that they were in flush contact at the periphery while platform mismatching connection had smaller abutment diameter than implant so that their connection was made away from periphery of implant-bone interface. Occlusal loading of 100N force was applied vertically and 30 degree obliquely to all 8 models and the maximum von Mises bone stress was identified. Loading stress as highly concentrated in cortical bone. Platform mismatching scalloped implant with small thread pitch (0.4mm) model had consistently lowest maximum von Mises bone stress in vertical and oblique loads. Platform matching model had lowest maximum von Mises bone stress with 0.6mm thread pitch in vertical load and with 0.4mm thread pitch in oblique load. Platform mismatching connection had important roles in reducing maximum von Mises bone stress. Scalloped implant with smaller coronal thread pitch showed trend of reducing maximum von Mises bone stress under load.

• Smart Structures and Systems
• /
• v.4 no.2
• /
• pp.221-232
• /
• 2008

This paper represents the final results of a research program sponsored by the European Commission through project WIND-CHIME ($\underline{W}$ide Range Non-$\underline{IN}$trusive $\underline{D}$evices toward $\underline{C}$onservation of $\underline{HI}$storical Monuments in the $\underline{ME}$diterranean Area), in which the possibility of using advanced seismic protection technologies to preserve historical monuments in the Mediterranean area is investigated. In the current research, the dynamic characteristics of two outstanding Mamluk-Style minarets, which similar minarets were reported to experience extensive damage during Dahshur 1992 earthquake, are investigated. The first minaret is the Qusun minaret (1337 A.D, 736 Hijri Date (H.D)) located in El-Suyuti cemetery on the southern side of the Salah El-Din citadel. The minaret is currently separated from the surrounding building and is directly resting on the ground (no vaults underneath). The total height of the minaret is 40.28 meters with a base rectangular shaft of about 5.42 ${\times}$ 5.20 m. The second minaret is the southern minaret of Al-Sultaniya (1340 A.D, 739 H.D). It is located about 30.0 meters from Qusun minaret, and it is now standing alone but it seems that it used to be attached to a huge unidentified structure. The style of the minaret and its size attribute it to the first half of the fourteenth century. The minaret total height is 36.69 meters and has a 4.48 ${\times}$ 4.48 m rectangular base. Field investigations were conducted to obtain: (a) geometrical description of the minarets, (b) material properties of the minarets' stones, and (c) soil conditions at the minarets' location. Ambient vibration tests were performed to determine the modal parameters of the minarets such as natural frequencies and mode shapes. A $1/16^{th}$ scale model of Qusun minaret was constructed at Cairo University Concrete Research Laboratory and tested under free vibration with and without SMA wire dampers. The contribution of SMA wire dampers to the structural damping coefficient was evaluated under different vertical loads and vibration amplitudes. Experimental results were used along with the field investigation data to develop a realistic 3-D finite element model that can be used for seismic risk evaluation of the minarets. Examining the updated finite element models under different seismic excitations indicated the vulnerability of such structures to earthquakes with medium to high a/v ratio. The use of SMA wire dampers was found feasible for reducing the seismic risk for this type of structures.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.