• Journal of Korea Game Society
• /
• v.18 no.5
• /
• pp.123-132
• /
• 2018

Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.8 no.9
• /
• pp.3231-3249
• /
• 2014

Deniable encryption, introduced in 1997 by Canetti, Dwork, Naor, and Ostrovsky, guarantees that the sender or the receiver of a secret message is able to "fake" the message encrypted in a specific ciphertext in the presence of a coercing adversary, without the adversary detecting that he was not given the real message. Sender - side deniable encryption scheme is considered to be one of the classification of deniable encryption technique which defined as resilient against coercing the sender. M. H. Ibrahim presented a sender - side deniable encryption scheme which based on public key and uncertainty of Jacobi Symbol [6]. This scheme has several problems; (1) it can't be able to derive the fake message $M_f$ that belongs to a valid message set, (2) it is not secure against Quadratic Residue Problem (QRP), and (3) the decryption process is very slow because it is based dramatically on square root computation until reach the message as a Quadratic Non Residue (QNR). The first problem is solved by J. Howlader and S. Basu's scheme [7]; they presented a sender side encryption scheme that allows the sender to present a fake message $M_f$ from a valid message set, but it still suffers from the last two mentioned problems. In this paper we present a new sender-side deniable public-key encryption scheme with fast decryption by which the sender is able to lie about the encrypted message to a coercer and hence escape coercion. While the receiver is able to decrypt for the true message, the sender has the ability to open a fake message of his choice to the coercer which, when verified, gives the same ciphertext as the true message. Compared with both Ibrahim's scheme and J. Howlader and S. Basu's scheme, our scheme enjoys nice two features which solved the mentioned problems: (1) It is semantically secure against Quadratic Residue Problem; (2) It is as fast, in the decryption process, as other schemes. Finally, applying the proposed deniable encryption, we originally give a coercion resistant internet voting model without physical assumptions.

• Structural Engineering and Mechanics
• /
• v.19 no.4
• /
• pp.425-440
• /
• 2005

In this paper, the behavior of a crack between two half-planes of functionally graded materials subjected to arbitrary tractions is resolved using a somewhat different approach, named the Schmidt method. To make the analysis tractable, it is assumed that the Poisson's ratios of the mediums are constants and the shear modulus vary exponentially with coordinate parallel to the crack. By use of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations in which the unknown variables are the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surfaces are expanded in a series of Jacobi polynomials. This process is quite different from those adopted in previous works. Numerical examples are provided to show the effect of the crack length and the parameters describing the functionally graded materials upon the stress intensity factor of the crack. It can be shown that the results of the present paper are the same as ones of the same problem that was solved by the singular integral equation method. As a special case, when the material properties are not continuous through the crack line, an approximate solution of the interface crack problem is also given under the assumption that the effect of the crack surface interference very near the crack tips is negligible. It is found that the stress singularities of the present interface crack solution are the same as ones of the ordinary crack in homogenous materials.

• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.

• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.73-83
• /
• 2022

This paper is concerned with an analysis of a surface edge crack emanated from a sharp contact edge. For a geometrical model, a square wedge is in contact with a half plane whose materials are identical, and a surface perpendicular crack initiated from the contact edge exists in the half plane. To analyze this crack problem, it is necessary to evaluate the stress field on the crack line which are induced by the contact tractions and pseudo-dislocations that simulate the crack, using the Bueckner principle. In this Part I, the stress filed in the half plane due to the contact is re-summarized using an asymptotic analysis method, which has been published before by the author. Further focus is given to the stress field in the half plane due to a pseudo-edge dislocation, which will provide a stress solution due to a crack (i.e. a continuous distribution of edge dislocations) later, using the Burgers vector. Essential result of the present work is the corrective functions which modify the stress field of an infinite domain to apply for the present one which has free surfaces, and thus the infiniteness is no longer preserved. Numerical methods and coordinate normalization are used, which was developed for an edge crack problem, using the Gauss-Jacobi integration formula. The convergence of the corrective functions are investigated here. Features of the corrective functions and their application to a crack problem will be given in Part II.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.280-290
• /
• 2002

Three-dimensional (3-D) electromagnetic (EM) modeling algorithm has been developed using finite element method (FEM) to acquire more efficient interpretation techniques of EM data. When FEM based on nodal elements is applied to EM problem, spurious solutions, so called 'vector parasite', are occurred due to the discontinuity of normal electric fields and may lead the completely erroneous results. Among the methods curing the spurious problem, this study adopts vector element of which basis function has the amplitude and direction. To reduce computational cost and required core memory, complex bi-conjugate gradient (CBCG) method is applied to solving complex symmetric matrix of FEM and point Jacobi method is used to accelerate convergence rate. To verify the developed 3-D EM modeling algorithm, its electric and magnetic field for a layered-earth model are compared with those of layered-earth solution. As we expected, the vector based FEM developed in this study does not cause ny vector parasite problem, while conventional nodal based FEM causes lots of errors due to the discontinuity of field variables. For testing the applicability to high frequencies 100 MHz is used as an operating frequency for the layer structure. Modeled fields calculated from developed code are also well matched with the layered-earth ones for a model with dielectric anomaly as well as conductive anomaly. In a vertical electric dipole source case, however, the discontinuity of field variables causes the conventional nodal based FEM to include a lot of errors due to the vector parasite. Even for the case, the vector based FEM gave almost the same results as the layered-earth solution. The magnetic fields induced by a dielectric anomaly at high frequencies show unique behaviors different from those by a conductive anomaly. Since our 3-D EM modeling code can reflect the effect from a dielectric anomaly as well as a conductive anomaly, it may be a groundwork not only to apply high frequency EM method to the field survey but also to analyze the fold data obtained by high frequency EM method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.9-16
• /
• 2014

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The "Hamilton-Jacobi(H-J)" equation and computationally robust numerical technique of "up-wind scheme" lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H-J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes are not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Journal of Radiation Protection and Research
• /
• v.16 no.2
• /
• pp.41-54
• /
• 1991

This paper aimed to analyse dose sensitivity to the controllable parameters of indoor radon $(^{222}Rn)$ and its decay products (Rn-D) by applying the input~output linear system theory. Physical behaviors of $^{222}Rn\;&\;Rn-D$ were analyzed in terms of $(^{222}Rn)$ gas -generation, -migation and -infiltration to indoor environments, and the performance output-function, i. e. mean dose equivalent to Tracho-Bronchial (TB) lung region, was assessed to the following extented ranges of the controllable paramenters; a) the ventilation rate $constant({\lambda}_v)\;:\;0{\sim}50[h^{-l}].\;b)$ the attachment rate $constant({\lambda}_a)\;:\;0{\sim}500[h^{-l}].\;c)$ the unattached-deposition rate constant (${\lambda}^u_d)\;:\;0-50[h-l]$. A linear input-output model was reconstructed from the original models in literatures, as follows, which was modified into the matrices consisting of 111 nodal equations; a) indoor $^{222}Rn\;&\;Rn-D$ Behaviour; Jacobi-Porstendoerfer-Bruno model.