EP2808812A1 - Analysevorrichtung und simulationsverfahren - Google Patents

Analysevorrichtung und simulationsverfahren Download PDF

Info

Publication number
EP2808812A1
EP2808812A1 EP12866680.7A EP12866680A EP2808812A1 EP 2808812 A1 EP2808812 A1 EP 2808812A1 EP 12866680 A EP12866680 A EP 12866680A EP 2808812 A1 EP2808812 A1 EP 2808812A1
Authority
EP
European Patent Office
Prior art keywords
unit
particle
cell
unit cell
analyzed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP12866680.7A
Other languages
English (en)
French (fr)
Other versions
EP2808812A4 (de
Inventor
Yoshitaka Ohnishi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sumitomo Heavy Industries Ltd
Original Assignee
Sumitomo Heavy Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sumitomo Heavy Industries Ltd filed Critical Sumitomo Heavy Industries Ltd
Publication of EP2808812A1 publication Critical patent/EP2808812A1/de
Publication of EP2808812A4 publication Critical patent/EP2808812A4/de
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Definitions

  • the present invention relates to an analysis device which analyzes a particle system, and a simulation method.
  • a molecular dynamics method (hereinafter, referred to as an MD method, and for example, see the related art) is known.
  • an MD method a molecular dynamics method
  • the motion of a particle can be physically handled more exactly compared to a finite element method (for example, see the related art) or the like.
  • the MD method if the number of particles increases, since the amount of calculation increases rapidly, only a comparatively small number of particles are handled in practical use. Accordingly, in the related art, there are many cases where the MD method is primarily used for the purpose, in which the shape of the object to be analyzed is less involved, such as prediction of physical properties of materials.
  • a period boundary condition is applied. Under the period boundary condition, a substantially cubic unit cell is defined, and the same space as the unit cell is repeatedly assumed in three orthogonal directions of x, y, and z (this is referred to as a replica cell). In addition, it is assumed that the exact same motion of an atom as the motion of an atom in the unit cell occurs in the respective replica cells. For example, if an atom in the unit cell exits out of the boundary surface of the unit cell, the same kind of atom having the same speed as the exited atom enters the corresponding position of the unit cell from an adjacent replica cell through the opposite boundary surface of the unit cell simultaneously.
  • an object to be analyzed spreads to a real structure, such as a gear, a motor, or a beam.
  • an analysis device including a cell acquisition unit configured to acquire a unit cell under a period boundary condition reflecting symmetry of an object to be analyzed, a particle system acquisition unit configured to acquire a system including a plurality of particles in the unit cell acquired by the cell acquisition unit, and a numerical computation unit configured to numerically compute a governing equation, which governs the motion of each particle of the system acquired by the particle system acquisition unit, under the period boundary condition.
  • the numerical computation unit moves the boundary of the unit cell to simulate the movement of the object to be analyzed during computation.
  • a simulation method in which, when analyzing an object having rotational symmetry using a renormalized molecular dynamics method, a period boundary condition reflecting rotational symmetry of the object is set.
  • a real structure has a shape which is not isotropic. Accordingly, the above-described instinctive period boundary condition for use in the MD method is unsuitable as a boundary condition which is applied when simulating a real structure using the RMD method.
  • the unit cell corresponding to the moving object it is possible to use the unit cell corresponding to the moving object to be analyzed.
  • any arbitrary combination of the above-described structural components or rearrangement of the structural components and the expressions of the invention among a device, a method, a system, a computer program, a recording medium having a computer program recorded thereon, and the like are also effective as the embodiments of the invention.
  • a period boundary condition reflecting rotational symmetry of the object is set.
  • the user generates a unit particle system and a unit cell in the following exemplary procedure as an initial condition for numerical computation in the analysis device according to the embodiment.
  • the user selects a real object having rotational symmetry, such as a gear or a disk, as an object to be analyzed.
  • the user generates an object region having a similar shape to the shape of the object to be analyzed in a three-dimensional virtual space using the analysis device or other computation devices.
  • the user disposes particles in the MD method, that is, a plurality of particles corresponding to a real-world atom or molecule in the generated object region.
  • the number of particles to be disposed is about an Avogadro number.
  • the user applies a conversion rule in the RMD method to a system including a plurality of particles disposed.
  • the number of particles is reduced to about millions by millions of times of renormalization.
  • the user defines a unit cell under a period boundary condition reflecting rotational symmetry of the object to be analyzed for the renormalized system.
  • the user defines a system including a plurality of particles in the defined unit cell as a unit particle system.
  • the unit particle system is a part of the renormalized system.
  • a unit cell may be defined for the system including a plurality of particles disposed, and a conversion rule may be applied to a system including particles in the defined unit cell.
  • Fig. 1 is a schematic view showing an example of a unit cell 300.
  • an object to be analyzed has a uniform disk shape.
  • an object region hereinafter, referred to as a renormalized object region
  • Particles 310 of the renormalized system are disposed in the renormalized object region 302.
  • the user defines a first plane 306 which includes a rotation axis 304 (an axis extending in a direction orthogonal to the paper surface of Fig. 1 ) as a reference of rotational symmetry of the renormalized object region 302, and a second plane 308 which includes the rotation axis 304 and is orthogonal to the first plane 306. Since the renormalized object region 302 is divided into four regions by the first plane 306 and the second plane 308, the user selects one region of the divided regions as the unit cell 300.
  • the user defines a system including a plurality of particles 310 in the selected unit cell 300 as a unit particle system.
  • the unit particle system is a subset of the renormalized system. With this meaning, it can be said that the definition of the unit particle system is an operation to cut the unit particle system from the renormalized system.
  • the unit cell 300 is a fan-shaped region having a predetermined thickness, and a central angle ⁇ thereof is 90 degrees.
  • the surface of the unit cell 300 includes a first boundary surface 312 which is a part of the first plane 306, and a second boundary surface 314 which is a part of the second plane 308, and the unit cell 300 is a region which is sandwiched between these boundary surfaces.
  • the first boundary surface 312 and the second boundary surface 314 are planes which extend from the rotation axis 304 to the boundary of the renormalized object region 302 outward in a radial direction.
  • the second boundary surface 314 as a part of the second plane 308, alternatively, a surface which is obtained by rotating the first boundary surface 312 by 90 degrees (an angle of 1/4 of 360 degrees) in a counterclockwise direction around the rotation axis 304 may be defined as the second boundary surface 314.
  • a corresponding particle 310b enters the unit cell 300 from a corresponding position on the second boundary surface 314 at a corresponding speed.
  • a particle runs outside through the second boundary surface 314.
  • Fig. 2 is a block diagram showing the functions and configuration of an analysis device 100 according to the embodiment.
  • respective blocks shown in the drawing can be realized by hardware, for example, an element, such as a central processing unit (CPU) of a computer, or a mechanical device, or can be realized by software, for example, a computer program or the like, functional blocks which are realized by collaboration of hardware and software are drawn. Accordingly, it is understood by those skilled in the art of contact with the specification that these functional blocks can be realized by combination of hardware and software in various ways.
  • the analysis device 100 is connected to an input device 102 and a display 104.
  • the input device 102 may be a device which receives an input related to processing executed on the analysis device 100 from the user, for example, a keyboard or a mouse.
  • the input device 102 may be configured to receive an input from a network, such as Internet, or a recording medium, such as a compact disk (CD) or a digital versatile disk (DVD).
  • a network such as Internet
  • a recording medium such as a compact disk (CD) or a digital versatile disk (DVD).
  • the analysis device 100 includes a unit cell acquisition unit 114, a unit particle system acquisition unit 116, a numerical computation unit 118, a display control unit 122, a unit particle information retaining unit 112, a unit cell information retaining unit 136, and a replica particle information retaining unit 124.
  • the unit cell acquisition unit 114 acquires the unit cell 300 defined by the user through the input device 102.
  • the unit cell acquisition unit 114 registers information, which expresses the acquired unit cell 300 in a virtual space, in the unit cell information retaining unit 136.
  • Information which expresses the unit cell 300 in a virtual space is, for example, functions representing a plane, a curved surface, a line, or a curve, or a combination thereof, and in particular, includes information which expresses the first boundary surface 312 and the second boundary surface 314.
  • the unit particle system acquisition unit 116 acquires the unit particle system defined by the user through the input device 102.
  • the unit particle system includes Q (where Q is a natural number) unit particles.
  • the unit particle system acquisition unit 116 registers a particle ID for specifying a unit particle of the unit particle system, the position vector (hereinafter, simply referred to as position) of the unit particle, and the speed vector (hereinafter, simply referred to as speed) of the unit particle in the unit particle information retaining unit 112 in association with one another.
  • the numerical computation unit 118 numerically computes a governing equation, which governs the motion of each unit particle of the unit particle system acquired by the unit particle system acquisition unit 116, under the period boundary condition.
  • the numerical computation unit 118 performs repetitive computation according to an equation of motion of discretized particles.
  • a condition that the object to be analyzed moves in the radial direction is imposed.
  • the numerical computation unit 118 moves the first boundary surface 312 and the second boundary surface 314 in the same direction during the computation (see a broken-line arrow of Fig. 1 ).
  • the direction of movement corresponds to the direction of movement of the object to be analyzed.
  • the numerical computation unit 118 updates the position and speed of each unit particle of the unit particle system based on the computation result.
  • the numerical computation unit 118 includes a replica particle generation unit 126, a force computation unit 128, a particle state computation unit 130, a state update unit 132, an end condition determination unit 134, and a boundary movement unit 120.
  • the duplicated particle is called a replica particle.
  • the replica particle generation unit 126 registers a particle ID for specifying a replica particle, the position of the replica particle, and the speed of the replica particle in the replica particle information retaining unit 124 in association with one another.
  • colored circles represent a unit particle and three replica particles corresponding to the unit particle.
  • the force computation unit 128 refers to the unit particle information retaining unit 112 and the replica particle information retaining unit 124, and for each unit particle of the unit particle system, computes a force applied to the unit particle based on the distance between particles. At this time, the force computation unit 128 may compute the distance between each unit particle and all replica particles to compute the force applied to the unit particle. In particular, the force computation unit 128 may compute the distance between each unit particle and replica particles, which are a replica of the unit particle to compute the force applied to the unit particle by the replica particles.
  • the force computation unit 128 decides unit particles or replica particles (hereinafter, referred to as near particle) which are at a distance from the i-th unit particle smaller than a predetermined cutoff distance. For each near particle, the force computation unit 128 computes the force applied to the i-th unit particle by the near particle based on the potential energy function between the near particle and the i-th unit particle and the distance between the near particle and the i-th unit particle. In particular, the force computation unit 128 calculates the force from the value of the gradient of the potential energy function in the value of the distance between the near particle and the i-th unit particle. The force computation unit 128 adds the force applied to the i-th unit particle by the near particle for all near particles to calculate the force applied to the i-th unit particle.
  • near particle unit particles or replica particles
  • the particle state computation unit 130 refers to the unit particle information retaining unit 112, and applies the force computed by the force computation unit 128 to the equation of motion of discretized particles for each unit particle to compute at least one of the position and speed of the unit particle. In this embodiment, the particle state computation unit 130 computes both the position and speed of the unit particle.
  • the particle state computation unit 130 computes the speed of the unit particle from the equation of motion of discretized particles including the force computed by the force computation unit 128.
  • the particle state computation unit 130 substitutes the force computed by the force computation unit 128 in the equation of motion of particles discretized using a predetermined minute time interval ⁇ t based on a predetermined numerical analysis method, such as a leapfrog method or an Euler method, for the i-th unit particle, thereby computing the speed of the unit particle.
  • a predetermined numerical analysis method such as a leapfrog method or an Euler method
  • the particle state computation unit 130 calculates the position of the unit particle based on the computed speed of the unit particle.
  • the particle state computation unit 130 applies the computed speed of the unit particle to a relationship expression of position and speed of discretized particles using the time interval ⁇ t based on a predetermined numerical analysis method for the i-th unit particle, thereby computing the position of the unit particle. In the computation, the position of the unit particle computed in the previous repetitive computation cycle is used.
  • the state update unit 132 updates the position and speed of each unit particle retained in the unit particle information retaining unit 112 to the position and speed computed by the particle state computation unit 130.
  • the state update unit 132 refers information of the first boundary surface 312 and information of the second boundary surface 314 retained in the unit cell information retaining unit 136, and for each unit particle, performs determination about whether or not the unit particle exits from the first boundary surface 312 or the second boundary surface 314 to the outside of the unit cell 300 as a result of updating the position.
  • the state update unit 132 rotates a unit particle, which is determined to exit to the outside of the unit cell 300, around the rotation axis 304 toward the unit cell 300 along with the direction of the speed.
  • the state update unit 132 rotates the position vector and the speed vector of a unit particle exited from the first boundary surface 312 to the outside of the unit cell 300 in Fig. 1 by 90 degrees in a counterclockwise direction, and rotates the position vector and speed vector of a unit particle exited from the second boundary surface 314 to the outside of the unit cell 300 by 90 degrees in a clockwise direction.
  • the state update unit 132 may add the position and speed computed by the particle state computation unit 130 in the unit particle information retaining unit 112, instead of updating the unit particle information retaining unit 112.
  • the end condition determination unit 134 performs determination about whether or not to end repetitive computation in the numerical computation unit 118.
  • An end condition for ending repetitive computation is, for example, that repetitive computation is performed a predetermined number of times, an end instruction is received from the outside, or the system reaches a normal state.
  • the end condition determination unit 134 ends repetitive computation in the numerical computation unit 118.
  • the boundary movement unit 120 moves the first boundary surface 312 and the second boundary surface 314 in the same direction.
  • the amount of single movement in this case is set to a unit movement distance which is a distance obtained by multiplying a set movement speed by the time interval ⁇ t.
  • the boundary movement unit 120 registers a boundary surface obtained by moving the first boundary surface 312 in parallel at the unit movement distance in the direction of movement in the unit cell information retaining unit 136 as a new first boundary surface 312.
  • Processing when a unit particle exits to the outside of the unit cell 300 with the movement of the boundary surface is equivalent to processing in the state update unit 132. After the boundary movement processing in the boundary movement unit 120, the process returns to the replica particle generation unit 126.
  • the display control unit 122 displays a mode of temporal development of the unit particle system or a state at a certain time on the display 104 based on the position and speed of each unit particle of the unit particle system retained in the unit particle information retaining unit 112.
  • Fig. 3 is a data structural diagram showing an example of the unit particle information retaining unit 112.
  • the unit particle information retaining unit 112 retains a particle ID, the position of the particle, and the speed of the particle in association with one another.
  • the replica particle information retaining unit 124 has the same data structure as shown in Fig. 3 .
  • an example of the retaining unit is a hard disk or a memory. It is understood by those skilled in the art of contact with the specification that the respective units can be realized by a CPU (not shown), a module of an installed application program, a module of a system program, a memory which temporarily stores the contents of data read from a hard disk, or the like based on the description of the specification.
  • Fig. 4 is a flowchart showing an example of a sequence of processing in the analysis device 100.
  • the unit cell acquisition unit 114 acquires information of the unit cell defined by the user (S202).
  • the unit particle system acquisition unit 116 acquires information of the unit particle system defined by the user (S204).
  • the replica particle generation unit 126 generates replica particles corresponding to each unit particle (S206).
  • the force computation unit 128 computes the force applied to each unit particle (S208).
  • the particle state computation unit 130 computes the position and speed of each unit particle based on the equation of motion of particles (S210).
  • the state update unit 132 updates the position and speed of each unit particle (S212).
  • the state update unit 132 performs determination about whether or not there is a unit particle outside the unit cell (S214).
  • the state update unit 132 rotates the position and speed of the exited particle (S216), and the process returns to Step S214.
  • the end condition determination unit 134 performs determination about whether or not the end condition is satisfied (S218).
  • the boundary movement unit 120 moves the boundary surfaces of the unit cell (S220), and the process returns to Step S206.
  • the display control unit 122 displays the computation result on the display 104 (S222).
  • the period boundary condition reflecting rotational symmetry is set.
  • the first boundary surface 312 and the second boundary surface 314 of the unit cell 300 move in the same direction to correspond to the movement of the object to be analyzed under the period boundary condition during numerical computation. Accordingly, even when the object to be analyzed moves, the period boundary condition can be applied in the same manner as when the object to be analyzed is stationary. That is, a particle system corresponding to the entire object to be analyzed in a moving state does not need to be computed, and only the particles in the unit cell 300 are to be computed, whereby it is possible to simulate the behavior of the object to be analyzed in the moving state. Therefore, it is possible to reduce the number of particles which are handled in numerical computation, and to reduce a calculation load.
  • the inventors have simulated how a workpiece is deformed by spinning using the analysis device 100 as a preliminary experiment.
  • Figs. 5A and 5B are schematic views showing a unit particle system 400 corresponding to a workpiece and a work particle model 402 corresponding to a work.
  • Fig. 5A is a perspective view
  • Fig. 5B is a top view.
  • the workpiece has a disk shape, and the central angle K of the unit cell 404 is 36.5 degrees.
  • the workpiece moves toward a work at a predetermined speed.
  • the same speed as the direction from the rotation axis 410 toward the center of the work particle model 402 is set.
  • Figs. 6A to 6D are side views showing a mode of temporal development of the unit particle system 400.
  • Figs. 7A to 7D are side views showing a mode of temporal development of the unit particle system 400.
  • Figs. 6A, 6B, 6C, 6D , 7A, 7B, 7C, and 7D are time-sequentially arranged in this order. According to the comparison and examination of the inventors, the simulation result describes that plastic deformation in spinning is satisfactory. In this way, if the analysis device 100 of this embodiment is used, it becomes possible to more accurately understand a phenomenon in an object to be analyzed while moving.
  • a cylindrical period boundary condition is often used in a finite element method or the like.
  • the MD method in general, since an object region is a region on nano scale, there is less need for devising a boundary condition with an object to be analyzed on macro scale in mind.
  • the inventors have found out that, if the MD method is developed to the RMD method, since the shape of an object to be analyzed can be one of important factors in a simulation, it is advantageous to use a period boundary condition reflecting symmetry of the object to be analyzed in the RMD method.
  • a cylindrical period boundary condition is applied, whereby it is possible to reduce a calculation load in numerical computation and to improve precision of the simulation result.
  • the invention is not limited thereto, and there is no limit to the dimension of the virtual space, such as a two-dimensional virtual space.
  • the invention is not limited thereto, and for example, when N is a natural number equal to or greater than 2, one of N divided regions of the renormalized object region 302 may be defined as a unit cell. In this case, a period boundary condition is also realized by continuous rotational symmetry of the renormalized object region 302.
  • the replica particle generation unit may rotate each unit particle of the unit particle system around the rotation axis 304 by an angle of ⁇ ', 2 ⁇ ', ..., and (N-1) ⁇ ' and may duplicate the unit particle to generate replica particles.
  • a unit cell may be a region sandwiched between a first boundary surface which extends in a radial direction from a rotation axis and a second boundary surface which is obtained by rotating the first boundary surface at an angle of 1/M (where M is a natural number equal to or greater than 2) of 360 degrees around the rotation axis.
  • a unit cell may be defined as a region sandwiched between a first boundary surface which extends in a radial direction from a rotation axis and a second boundary surface which is obtained by rotating the first boundary surface at an angle of 1/P of 360 degrees around the rotation axis.
  • the invention is not limited thereto.
  • a method for numerical analysis like a Verlet method, a method in which, when computing the position of the particle, the position of the particle may be computed directly from the force applied to the particle, and the speed of the particle may not be computed explicitly is known, and the technical idea of this embodiment may be applied to this method.

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Biophysics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Biotechnology (AREA)
  • Evolutionary Biology (AREA)
  • Medical Informatics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
EP12866680.7A 2012-01-26 2012-08-09 Analysevorrichtung und simulationsverfahren Withdrawn EP2808812A4 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2012013721A JP5844165B2 (ja) 2012-01-26 2012-01-26 解析装置
PCT/JP2012/005069 WO2013111204A1 (ja) 2012-01-26 2012-08-09 解析装置およびシミュレーション方法

Publications (2)

Publication Number Publication Date
EP2808812A1 true EP2808812A1 (de) 2014-12-03
EP2808812A4 EP2808812A4 (de) 2016-05-11

Family

ID=48872983

Family Applications (1)

Application Number Title Priority Date Filing Date
EP12866680.7A Withdrawn EP2808812A4 (de) 2012-01-26 2012-08-09 Analysevorrichtung und simulationsverfahren

Country Status (4)

Country Link
US (1) US20140309975A1 (de)
EP (1) EP2808812A4 (de)
JP (1) JP5844165B2 (de)
WO (1) WO2013111204A1 (de)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6261130B2 (ja) * 2014-06-04 2018-01-17 国立研究開発法人海洋研究開発機構 粒子シミュレーション装置、粒子シミュレーション方法及び粒子シミュレーションプログラム
EP3218833A2 (de) * 2014-11-14 2017-09-20 D.E. Shaw Research, LLC Unterdrückung der interaktion zwischen gebundenen partikeln
JP6458501B2 (ja) * 2015-01-06 2019-01-30 富士通株式会社 シミュレーションプログラム、シミュレーション方法、およびシミュレーション装置
JP6900118B2 (ja) * 2017-10-24 2021-07-07 住友重機械工業株式会社 シミュレーション方法、シミュレーション装置、及びプログラム
US20220114309A1 (en) * 2020-10-08 2022-04-14 Dassault Systèmes Americas Corp. Modeling and Simulating Material Microstructures
CN113126773B (zh) * 2021-05-08 2023-05-30 北京理工大学 一种基于虚拟现实技术的交互式分子模拟系统
CN115329647B (zh) * 2022-10-14 2022-12-27 四川轻化工大学 一种基于sph的筒形件旋压成形过程分析方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002243722A (ja) * 2001-02-16 2002-08-28 Hideaki Umeyama 薬物の分子運動解析方法およびそれを用いる薬物分子のスクリーニング方法
JP4726189B2 (ja) * 2004-12-16 2011-07-20 キヤノン株式会社 粒子挙動計算結果表示方法
JP4666357B2 (ja) * 2005-04-04 2011-04-06 住友重機械工業株式会社 シミュレーション方法
JP5052985B2 (ja) * 2007-07-31 2012-10-17 住友重機械工業株式会社 分子シミュレーション方法、分子シミュレーション装置、分子シミュレーションプログラム、及び該プログラムを記録した記録媒体
JP5320806B2 (ja) * 2008-04-24 2013-10-23 横浜ゴム株式会社 回転体のシミュレーション方法
JP5402274B2 (ja) * 2009-06-11 2014-01-29 横浜ゴム株式会社 回転体のシミュレーション方法、装置及びプログラム
US8912973B2 (en) * 2011-05-04 2014-12-16 The Penn State Research Foundation Anisotropic metamaterial gain-enhancing lens for antenna applications

Also Published As

Publication number Publication date
EP2808812A4 (de) 2016-05-11
JP2013152658A (ja) 2013-08-08
JP5844165B2 (ja) 2016-01-13
WO2013111204A1 (ja) 2013-08-01
US20140309975A1 (en) 2014-10-16

Similar Documents

Publication Publication Date Title
US20140309975A1 (en) Analysis device and simulation method
Okereke et al. Finite element applications
Krüger et al. The lattice Boltzmann method
Matouš et al. A review of predictive nonlinear theories for multiscale modeling of heterogeneous materials
Fritzen et al. Nonlinear reduced order homogenization of materials including cohesive interfaces
Sotomayor et al. Role of cell regularity and relative density on elasto-plastic compression response of random honeycombs generated using Voronoi diagrams
Cole et al. Validation of a commercial fluid-structure interaction solver with applications to air cushion vehicle flexible seals
Biancolini et al. Radial basis functions mesh morphing: A comparison between the bi-harmonic spline and the wendland c2 radial function
Seth et al. Development of a dual-handed haptic assembly system: SHARP
Tschisgale et al. A constraint-based collision model for Cosserat rods
JP6129193B2 (ja) 解析装置
Fuge et al. The MechProcessor: Helping novices design printable mechanisms across different printers
Su et al. Accelerating the force-coupling method for hydrodynamic interactions in periodic domains
Barreto et al. A guide to modeling the geotechnical behavior of soils using the discrete element method
WO2014132300A1 (ja) 解析装置
JP5523364B2 (ja) 解析装置
Wang et al. Haptic interaction with fluid based on smooth particles and finite elements
Carrick et al. Packaging optimization using the dynamic vector fields method
Ganzenmüller et al. Consistent temperature coupling with thermal fluctuations of smooth particle hydrodynamics and molecular dynamics
Moore et al. Advancements in multiresolution analysis
Lyu et al. Integrating peridynamics with material point method for elastoplastic material modeling
JP5777506B2 (ja) 解析装置およびシミュレーション方法
Blom et al. Partitioned Fluid–Structure–Acoustics Interaction on Distributed Data: Numerical Results and Visualization
Patro et al. Speeding up particle trajectory simulations under moving force fields using graphic processing units
Krüger et al. Numerical methods for fluids

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20140617

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

DAX Request for extension of the european patent (deleted)
RA4 Supplementary search report drawn up and despatched (corrected)

Effective date: 20160408

RIC1 Information provided on ipc code assigned before grant

Ipc: G06F 19/00 20110101AFI20160404BHEP

Ipc: G06F 17/50 20060101ALI20160404BHEP

17Q First examination report despatched

Effective date: 20190529

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN WITHDRAWN

18W Application withdrawn

Effective date: 20190808