• 설비공학논문집
• /
• 제6권3호
• /
• pp.247-258
• /
• 1994

In this study, experimental investigations to find cold energy storage performance have been made for two different temperatures at condenser. Temperatures at inlet and outlet of condenser were measured to calculate global heat transfer coefficient of direct contact method in our cold energy storage system. Also storage performance by direct contact method was compared with that of Ice-On-Coil type ice storage which was calculated by analytic solution. Results show that, in the case of $-8.0^{\circ}C$ at condenser inlet, heat transfer coefficient of direct contact method is 3.25 times higher than that of conventional method and COP of system is improved by using R141b as refrigerant which produces gas hydrate and has higher phase change temperature than $0.0^{\circ}C$.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1994년도 추계학술대회 논문집
• /
• pp.837-842
• /
• 1994

The boundary element method was used for studying singularities of an interface crack with contact zones. Because the contact zone size was extremely small in a tention field, a large number of Gaussian points were is used for numerical integration of the Kernels. In order to estimate the contact zone size, iteration method was used. The interface crack tips with contact zones showed no oscillatory behavior and J-integral values across the interface were conserved.

• Coupled systems mechanics
• /
• 제7권5호
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.551-555
• /
• 2010

Electrowetting is a versatile tool to handle tiny droplets and forms a backbone of digital microfluidics. Numerical analysis is necessary to fully understand the dynamics of electrowetting, especially in designing electrowetting-based devices, such as liquid lenses and reflective displays. We developed a numerical method to analyze the general contact-line problems, incorporating dynamic contact angle models. The method is based on the conservative level set method to capture the interface of two fluids without loss of mass. We applied the method to the analysis of spreading process of a sessile droplet for step input voltages and oscillation of the droplet for alternating input voltages in electrowetting. The result was compared with experimental data. It is shown that contact line friction significantly affects the contact line motion and the oscillation amplitude. The pinning process of contact line was well represented by including the hysteresis effect in the contact angle models. In level set method, in the mean time, material properties are made to change smoothly across an interface of two materials with different properties by introducing an interpolation or smoothing scheme. So far, the weighted arithmetic mean (WAM) method has been exclusively adopted in level set method, without complete assessment for its validity. We viscosity, thermal conductivity, electrical conductivity, and permittivity, can be an alternative. I.e., the WHM gives more accurate results than the WAM method in certain circumstances. The interpolation scheme should be selected considering various characteristics including type of property, ratio of property of two fluids, geometry of interface, and so on.

• Computers and Concrete
• /
• 제27권3호
• /
• pp.199-210
• /
• 2021

The aim of this paper was to examine the continuous and discontinuous contact problems between the functionally graded (FG) layer pressed with a uniformly distributed load and homogeneous half plane using an analytical method and FEM. The FG layer is made of non-homogeneous material with an isotropic stress-strain law with exponentially varying properties. It is assumed that the contact at the FG layer-half plane interface is frictionless, and only the normal tractions can be transmitted along the contacted regions. The body force of the FG layer is considered in the study. The FG layer was positioned on the homogeneous half plane without any bonds. Thus, if the external load was smaller than a certain critical value, the contact between the FG layer and half plane would be continuous. However, when the external load exceeded the critical value, there was a separation between the FG layer and half plane on the finite region, as discontinuous contact. Therefore, there have been some steps taken in this study. Firstly, an analytical solution for continuous and discontinuous contact cases of the problem has been realized using the theory of elasticity and Fourier integral transform techniques. Then, the problem modeled and two-dimensional analysis was carried out by using ANSYS package program based on FEM. Numerical results for initial separation distance and contact stress distributions between the FG layer and homogeneous half plane for continuous contact case; the start and end points of separation and contact stress distributions between the FG layer and homogeneous half plane for discontinuous contact case were provided for various dimensionless quantities including material inhomogeneity, distributed load width, the shear module ratio and load factor for both methods. The results obtained using FEM were compared with the results found using analytical formulation. It was found that the results obtained from analytical formulation were in perfect agreement with the FEM study.

• 한국표면공학회지
• /
• 제35권2호
• /
• pp.129-137
• /
• 2002

The purpose of solderability assessment is to predict the effectiveness of soldering process. It is important for companies pursuing zero defects manufacturing because poor solderability is the major cause of two third of soldering failures. The most versatile solderability method is wetting balance method. However, there exist so many indices for wettability in the wetting balance test e.g. time to reach 2/3 values of maximum wetting force, tine to reach zero wetting force, maximum withdrawal force. In this study, three solderability assessment methods, which were the maximum withdrawal force, the wetting balance and the dynamic contact angle (DCA), were evaluated by comparing each other. The wetting balance technique measures the solderability by recording the forces exerted from the specimen after being dipped into the molten solder. Then the force at equilibrium state can be used to calculate a contact angle, which is known as static contact angles. The DCA measures contact angles occurred during advancing and withdrawing of the specimen and the contact angles are known as dynamic contact angles. The maximum withdrawal force uses the maximum force during withdrawal movement and then a contact angle can be calculated. In this study, the maximum withdrawal force method was found to be an objective index for measuring the solderability and the experiment results indicated good agreement between the maximum withdrawal force and the wetting balance method.

• 소성∙가공
• /
• 제33권4호
• /
• pp.248-254
• /
• 2024

The triply periodic minimal surface (TPMS) structure is characterized by a high surface-to-volume (S/V) ratio and the separated internal structure for flow. Combining the different TPMS structures can provide unique flow and strength characteristics. This paper investigates the effects of dimension, density and arrangement of the unit cell of the TPMS on contact and flow areas of combined TPMS structures. Several representative TPMS structures, including primitive, gyroid and diamond structures, are adopted to design gradient and heterogeneous types TPMS structures. The estimation method of contact and flow areas using an image processing technique is proposed. Python software is used to predict contact and flow area. The influence of the combination method of TPMS on contact and flow areas in the contact surface of combined TPMS structures with different shapes is investigated. Based on the results of the investigation, an appropriate combination method of TPMS structures is discussed.

• Structural Engineering and Mechanics
• /
• 제72권5호
• /
• pp.557-568
• /
• 2019

Accurate modeling of contact interface in bolted joints is crucial in predicting the dynamic behavior for bolted assemblies under external load. This paper presents a contact pressure distribution based non-uniform virtual material method to describe the joint interface of assembly structure, which is connected by sparsely distributed multi-bolts. Firstly, the contact pressure distribution of bolted joints is obtained by the nonlinear static analysis in the finite element software ANSYS. The contact surface around bolt hole is divided into several sub-layers, and contact pressure in each sub-layer is thought to be evenly. Then, considering multi-asperity contact at the micro perspective, the relationship between contact pressure and interfacial virtual material parameters for each sub-layer is established by using the fractal contact theory. Finally, an experimental platform for the dynamic characteristics testing of a beam lap structure with double-bolted joint is constructed to validate the efficiency of proposed method. It is found that the theoretical results are in good agreement with experimental results by impact response in both time- and frequency-domain, and the relative errors of the first four natural frequencies are less than 1%. Furthermore, the presented model is used to examine the effect of rough contact surface on dynamic characteristics of bolted joint.

• 한국생산제조학회지
• /
• 제23권5호
• /
• pp.478-484
• /
• 2014

Ultrasonic atomic force microscopy(Ultrasonic-AFM) can be used to obtain images of the elastic properties of a subsurface and to evaluate the elastic properties by measuring the contact resonance frequency. When a tip is in contact with the sample, it is necessary to understand the cantilever behavior and the tip-sample interaction for the quantitative and reliable analysis. Therefore, precise analysis models that can accurately simulate the tip-sample contact are required; these can serve as good references for predicting the contact resonance frequency. In this study, modal analyses of the first four modes were performed to calculate the contact resonance frequency by using a spring model, and the deformed shapes of the cantilever were visualized at each mode. We presented the contact characteristics of the cantilever with a variety of contact conditions by applying the contact area, contact material thickness, and material properties as the parameters for the FEM analysis.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 2000년도 추계학술대회 논문집
• /
• pp.517-522
• /
• 2000

When a body including a crack inside is subjected to the compressive forces, the crack is closed and sliding occurs on the crack surfaces. In this work, a subsurface crack subjected to a static or moving compressive load is analyzed with the finite element method considering friction on the crack surface. The friction on the crack surface is assumed to follow the Coulomb friction law. A numerical method based on the finite element method and iterative method is applied in this work. And the result is compared with the frictional contact of crack by ANSYS using contact 12 element. The numerical results of two methods are compared with the wellknown analytical solutions, and the accuracy of iterative method is checked..