JP2014203242A - Method for creating simulation model and simulation method - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a method for creating a simulation model capable of adjusting interactions of constituting objects, and a simulation method.SOLUTION: A method for creating a simulation model is a method for creating a composite model in which a computer can analyze a composite with at least one type of a second object of a different material characteristic existing in a first object. The method has: a step for determining a first object area obtained by reproducing the area of the first object and a second object area obtained by reproducing the area of the second object; and a step for setting a fixed area in which the first object area is fixed to the second object area about a boundary area between the first object area and the second object area, and a non-fixed area in which the first object area is not fixed to the second object area.

Description

  The present invention relates to a method for creating a simulation model of a composite in which at least one second object having different material properties is present in the first object, and a simulation method, and more particularly to an interaction between objects constituting the composite The present invention relates to a simulation model creation method and a simulation method capable of adjusting the parameters.

Recently, various dynamic simulations of composite materials using the finite element method have been proposed.
Among them, for example, a lot of modeling and analysis methods have been proposed for compound simulation that evaluates mechanical properties as a compound by creating a model in which fillers are arranged in a polymer and performing deformation analysis ( (See Patent Documents 1 and 2).

Patent Document 1 discloses a technique for completely fixing a filler and a polymer as a technique related to modeling.
In Patent Document 1, a first region that reproduces a region of a matrix material of a composite material in which particles are dispersed in a matrix material and a second region that reproduces a region of dispersed particles are defined, and the first region is defined as A mesh-divided model is formed, representative points belonging to the second region are determined, and force is applied to the nodes on the boundary according to the relative displacement of the nodes on the boundary of the model facing the second region with respect to the representative point A composite material model is created by adding force adjusting means to the model.

Patent Document 2 discloses a modeling method in which an interfacial phase is provided and the interaction between the two is adjusted by the material properties.
The rubber material model of Patent Document 2 includes a matrix model that models a rubber matrix, a filler model that models a filler, and an interface model that forms an interface between the matrix model and the filler model. This interface model continuously surrounds the filler model and has a thickness and is defined with viscoelastic properties different from the matrix model.

JP 2009-193290 A Japanese Patent No. 3668238

  In both the composite material model of Patent Document 1 and the rubber material model of Patent Document 2, the interface between the matrix and the filler is fixed. As described above, the conventional model is not necessarily sufficient because it may not match the actual phenomenon.

  An object of the present invention is to provide a simulation model creation method and a simulation method capable of solving the problems based on the prior art and adjusting the interaction of objects constituting a complex.

  To achieve the above object, according to a first aspect of the present invention, there is provided a composite model capable of being analyzed by a computer for a composite in which at least one second object having different material properties exists in the first object. A method for creating a simulation model for creating a first object region that reproduces a region of a first object and a second object region that reproduces a region of a second object; The first object region and the second object region are fixed to each other at the boundary region between the first object region and the second object region, and the first object region and the second object region are fixed to each other. The present invention provides a method for creating a simulation model characterized by having a step of setting a non-fixed region that is not.

For example, contact is defined in the non-stick region. Further, for example, a peeling condition is set in at least a part of the fixing region. A force depending on the distance from the boundary region may be defined in the non-fixed region.
Further, the elastic modulus of the second object region may be set higher than that of the first object region. Furthermore, viscoelasticity may be defined in at least the first object region.

  According to a second aspect of the present invention, there is provided a simulation method for simulating a computer-analysable model using a computer, wherein a deformation simulation is performed on a complex model created by the simulation model creation method of the present invention. To provide a simulation method characterized by calculating the mechanical response of the composite.

  According to the present invention, in the composite model of a composite in which at least one type of second object having different material characteristics exists in the first object, a fixed region and a non-fixed region are provided at the material boundary, thereby providing a composite The interaction between the objects constituting the body can be adjusted.

It is a schematic diagram which shows the simulation apparatus for implementing the creation method and simulation method of the simulation model of embodiment of this invention. (A) is a schematic diagram showing a first example of the complex model, (b) is a schematic diagram showing a second example of the complex model, and (c) is a schematic diagram of the complex model. 3 is a schematic diagram illustrating a fourth example of a complex model, and (e) is a schematic diagram illustrating a fifth example of the complex model. FIG. It is a flowchart which shows the preparation method of the simulation model of embodiment of this invention in order of a process. (A) is a schematic diagram which shows the complex model of Example 1, (b) is a schematic diagram which shows the complex model of Example 2, (c) is a complex model of Comparative Example 1. It is a schematic diagram which shows.

Hereinafter, a simulation model creation method and a simulation method according to the present invention will be described in detail based on preferred embodiments shown in the accompanying drawings.
FIG. 1 is a schematic diagram showing a simulation model creation method and a simulation apparatus for carrying out the simulation method according to the embodiment of the present invention.

The simulation apparatus 10 (hereinafter, referred to as a processing apparatus 10) shown in FIG. 1 has the following apparatus configuration including, for example, hardware such as a computer and its peripheral devices. The simulation model creation method and simulation method of the present invention are executed using the processing apparatus 10.
The processing apparatus 10 creates a composite model that can be analyzed by a computer for a composite in which at least one second object having different material properties from the first object is present in the first object. . Moreover, the processing apparatus 10 performs a deformation | transformation simulation on predetermined | prescribed simulation conditions about a complex model, and calculates the dynamic response of a complex.
The first object constituting the composite is, for example, a base material, and the second object is, for example, at least one kind of object having a material property different from that of the base material existing in the base material. is there.

The processing device 10 includes a processing unit 12, an input unit 14, and a display unit 16. The processing unit 12 includes a condition setting unit 20, a model generation unit 22, a calculation unit 24, a dynamic characteristic calculation unit 26, a memory 28, a display control unit 30, and a control unit 32. In addition, although not shown, it has a ROM and the like.
The processing unit 12 is controlled by the control unit 32. In the processing unit 12, the condition setting unit 20, the model generation unit 22, the calculation unit 24, and the dynamic characteristic calculation unit 26 are connected to the memory 28, and the condition setting unit 20, the model generation unit 22, the calculation unit 24, and the dynamic characteristic are connected. Data of the calculation unit 26 is stored in the memory 28.

  The input unit 14 is various input devices for inputting various information such as a mouse and a keyboard in accordance with an operator instruction. The display unit 16 displays, for example, the instruction input contents of the operator, the composite model, the progress and results of the simulation, and various known displays are used. The display unit 16 also includes a device such as a printer for displaying various types of information on an output medium.

  The processing apparatus 10 executes the program (computer software) stored in the ROM or the like by the control unit 32, so that the condition setting unit 20, the model generation unit 22, the calculation unit 24, and the display control unit 26 are functional. To form. As described above, the processing apparatus 10 may be configured by a computer in which each part functions by executing a program, or may be a dedicated apparatus in which each part is configured by a dedicated circuit.

The condition setting unit 20 receives information for creating a simulation model and information for deformation simulation performed on the simulation model via the input unit 14 and the like, and various types of information are stored in the memory 28. The Based on the information input to the condition setting unit 20, a simulation model is created by the model generation unit 22, and a deformation simulation is performed by the calculation unit 24 using, for example, a finite element method.
Specifically, the operator inputs an instruction to the condition setting unit 20 using the input unit 14 while looking at the input screen displayed on the display unit 16.

  Information for creating a simulation model is used in the finite element method other than information on a fixed region and a non-fixed region in a boundary region between a first object (base material) and a second object (object) described later. Known information for creating various known calculation models is included. For example, information on a region to be modeled, a material constant of a first object (base material), a size of the first object (base material), and a second object (in the first object (base material)) ( Size, shape, number and position information of the object), material constant of the second object (object), and a fixing region of a boundary region between the first object (base material) and the second object (object) described later And non-fixed area information. The first object (base material) is not particularly limited, but is, for example, rubber or polymer, and the object is not particularly limited, but is, for example, a filler. The type of the second object (object) is not limited to one type, and may be a plurality of types.

In addition, when the second object (object) is a filler, the size and shape information of the filler particle includes, for example, the position of the center point of the circle and the radius of the circle when the shape of the filler particle is a circular shape. Contains information. In the case of an elliptical shape, information on the major axis and minor axis, the position of the center point of the ellipse, and the orientation of the major axis or minor axis is included.
The material constant is a value such as an elastic modulus (Young's modulus), Poisson's ratio, viscoelasticity, and the like. As the elastic modulus, for example, the slope of the elastic region in the stress-strain or elongation curve (corresponding to the longitudinal elastic modulus of the elastic body) is used. In the case of a non-linear material, the elastic modulus is increased when the strain or elongation is small. The viscoelasticity is defined using, for example, a generalized Maxwell model.

  In the fixed region and the non-fixed region at the boundary region between the first object (base material) and the second object (object), the fixed region is, for example, that relative displacement does not occur or is in a node sharing state. Defined. Further, the non-fixed region is defined by relative displacement occurring or a state where no node is shared at the free boundary. The fixed area is, for example, information on a range sharing a node, and the information on a non-fixed area is information on a range not sharing a node at a free boundary.

  Further, for example, contact may be defined in a region other than the fixed region in the boundary region, that is, a non-fixed region. This contact is defined as a non-sticking area that does not bite into the second object area that reproduces the area of the second object, and if the second object area is an object, it does not bite into the object area that reproduces the object. Is done. For example, there is a boundary condition in which the node of the non-fixed area does not fall within the second object area (object area).

In addition to this, for example, a peeling condition can be set in at least a part of the fixing region. The peeling condition can be defined by, for example, the magnitude of the force acting on the node of the fixing region. It is assumed that a predetermined force value is set for the elastic modulus of the second object region (object region), and peeling occurs when the force acting on the node exceeds this force value.
Furthermore, a force depending on the distance from the boundary region may be defined in a region other than the fixed region in the boundary region, that is, a non-fixed region. The force depending on the distance includes both a force approaching the second object region (object region) and a force moving away from the second object region (object region). In the non-fixed area, the force acting as the distance from the surface of the second object area (object area) increases, and the force acting as the distance increases, conversely, the second object area (object area) ) The closer the distance from the surface, the smaller the applied force, and the longer the distance from the surface, the greater the applied force. Regarding peeling, even after satisfying the peeling conditions, either a force approaching the second object region (object region) or a force moving away from the second object region (object region) is applied. Also good.

  Information for deformation simulation includes known simulation conditions of the finite element method. For example, it is a simulation condition such as an external force for performing a simulation and a method of giving a forced displacement.

The model generation unit 22 calls information for creating the above-described simulation model from the memory 28, sets a range to be modeled using this, and in the first object region that reproduces the first object, A second object region that reproduces the second object is automatically set, and the first object region and the second object region are mesh-divided by automatic mesh generation means (not shown), and the material constant and the fixed region of the boundary region And a non-fixed region information is added to generate an executable complex model.
For example, if the first object is a base material and the second object is an object, the model generation unit 22 automatically sets the object region that reproduces the object in the base material region that reproduces the base material. Then, the base material region and the object region are divided into meshes by an automatic mesh generation means, and a composite model that can be executed is generated by adding material constants and information on the fixed region and the non-bonded region of the boundary region.
The complex model and its model creation method by the model generation unit 22 will be described later.

The generation of the complex model means that a file in which information on the position of nodes, the number of nodes, the number of elements, and the like and the material constants of each element are created on a computer is created. As a result, the complex model is constituted by numerical data that can be analyzed by the processing apparatus 10.
The complex model can be numerically analyzed by a computer (processing device 10). The possibility of numerical analysis means that deformation calculation can be performed on each element or the entire system by a numerical analysis method such as a finite element method, a finite volume method, a difference method, or a boundary element method. Specifically, for each element, a node coordinate value, an element shape, a material characteristic, and the like in the coordinate system are defined. For each element, for example, a triangular or quadrilateral element as a two-dimensional plane, and a three-dimensional element, for example, a tetrahedral element or a hexahedral element is preferably used.

  The calculation unit 24 performs a simulation calculation using a known finite element method by applying, for example, a forced displacement or an external force to a node of the complex model using the simulation conditions called from the memory 28 to the complex model. It is a part to do. Information on the deformed complex model obtained by the simulation operation is stored in the memory 28.

The mechanical property calculation unit 26 calls the operation result of the deformed composite model stored in the memory 28, and calculates values of parameters of the mechanical behavior of the entire composite model, for example, stress or strain distribution, loss energy, etc. Calculate the calculated mechanical properties. Further, when the mechanical characteristic is expressed as a linear characteristic, for example, a value such as Young's modulus or shear rigidity is calculated, and when it is expressed as a non-linear characteristic, for example, the value of each parameter of the superelastic potential is calculated.
The composite model before the simulation calculation, the composite model after the simulation calculation, the distribution of stress or strain acting on each finite element, the calculated value of the mechanical characteristics, and the like are output to the display unit 16 via the display control unit 30. .

  The display control unit 30 causes the display unit 16 to display various information stored in the memory 28, the composite model, simulation results, and the like.

Next, the simulation model of this embodiment will be described.
2A is a schematic diagram illustrating a first example of a complex model, FIG. 2B is a schematic diagram illustrating a second example of the complex model, and FIG. 2C is a complex model. (D) is a schematic diagram showing a fourth example of the complex model, and (e) is a schematic diagram showing a fifth example of the complex model. is there.
2A to 2E show composite models 40 to 40d that can be analyzed by a computer, but illustration of meshes and the like is omitted.
In the composite models 40 to 40d, the simulation model creation method, and the simulation method described below, the first material is described using the base material as an example, and the second object is described using the object as an example.

  As shown in FIG. 2A, the composite model 40 includes an object region 44 having a material characteristic different from that of the base material region 42 in the base material region 42. The material property is, for example, an elastic modulus. With respect to the boundary region 46 between the base material region 42 and the object region 44, a fixing region 48 where the base material region 42 and the object region 44 are fixed, and a base material region 42 and the object region 44 are fixed. There is no non-stick region 50. The boundary region 46 includes a fixed region 48 and a non-fixed region 50.

The number of fixing regions 48 is not limited to one, and may be a plurality of locations. For example, it may be two places like a complex model 40a shown in FIG.
The composite models 40 and 40a shown in FIGS. 2A and 2B are obtained by modeling a polymer in which Si particles are dispersed, for example. The Si particles are locally bonded to the polymer by, for example, a silane coupling agent, not the entire particle.

Further, the number of object regions 44 is not limited to one, and may be plural. For example, there may be two object regions 44 as in the complex model 40b shown in FIG. In this case, one object region 44 has one fixing region 48 in the boundary region 46, and the remaining object region 44 has two fixing regions 48 in the boundary region 46.
Furthermore, the type of the object region 44 existing in the base material region 42 is not limited to one, and may be a plurality of types. For example, two types of object regions 44 and object regions 52 may exist as in the complex model 40c illustrated in FIG. In this case, for the object region 44, there is one fixing region 48 in the boundary region 46, and in the object region 52, the boundary region 54 is all the fixing region 48. The object region 44 and the object region 52 have different material constants such as elastic modulus, for example.

Further, the object region 44 may be covered with another object region 56 as in the complex model 40d shown in FIG. In this case, one fixed region 48 exists in the boundary region 47 between the object region 44 and the object region 56. The boundary region 58 between the object region 56 and the base material region 42 is a fixed region 48.
In this way, by creating composite models of various forms, the interaction between the base material region 42 and the object regions 44, 52, 56 and the interaction between the object region 44 and the object region 56 can be adjusted. . Thereby, it is possible to deal with an event that could not be dealt with conventionally. In the boundary region, the number and ratio of the fixed region 48 and the non-fixed region 50 can be appropriately set according to the analysis target. The fixed area 48 and the non-fixed area 50 may be extremely small areas or dots.

  In the composite model, the elastic modulus of the object region 44 may be set higher than that of the base material region 42, and conversely, the elastic modulus may be lowered. For example, the elastic modulus of the object region 44 is at least 100 times the elastic modulus of the base material region 42. In addition, viscoelasticity may be defined in at least the base material region 42 out of the base material region 42 and the object region 44.

Next, a method for creating a simulation model according to the embodiment of the present invention will be described.
FIG. 3 is a flowchart showing a simulation model creation method according to the embodiment of the present invention in the order of steps.
A method for creating a simulation model will be described using an example of modeling a compound in which a polymer is used as a base material and one kind of filler particle is present as an object.

First, prior to creating the complex model, various information is set in the condition setting unit 20 based on the information input by the operator through the input unit 14 regarding the conditions for creating the complex model and the simulation conditions. And stored in the memory 28.
Specifically, the conditions for creating the composite model include materials such as the range of the region to be modeled as described above, material constants such as polymer elastic modulus, polymer viscoelasticity, and elastic modulus of filler particles. Information on constants, shape of filler particles, size, and position information. In addition to this, the number of fixed regions and non-fixed regions, peeling conditions in the fixed regions, contact conditions in the non-fixed regions, and force conditions depending on the distance from the boundary region in the non-fixed regions, etc. This is information on the fixing area.

Next, a condition for creating a complex model is called from the memory 28 to the model generation unit 22, and the model generation unit 22 sets a range of an area to be modeled (hereinafter referred to as a modeling area) (step). S10).
Next, in the modeling region, a filler region corresponding to the filler particles is set (step S12), and a polymer region corresponding to the polymer is set (step S14).
In step S12 and step S14, specifically, a filler region is secured in a predetermined modeling region based on the shape, size information, and position information of the filler particles, and the surrounding region is defined as a polymer region. Is set. For example, as shown in FIG. 2A, the object region 44 is determined based on the shape, size information, and position information of the filler particles. In the base material region 42, a region surrounding the object region 44 is set as a polymer region. A boundary region is set at the boundary between the filler region and the polymer region.

A known method can be used as a method for setting the boundary region. For example, the extraction of the node at the boundary between the filler region and the polymer region is performed when the filler region is represented by a circle, and whether or not the node to be examined is located on an arc using the center point and radius of the circle. Is determined, nodes located on the arc are extracted.
In the composite model, as is well known, the mechanical deformation behavior of the simulation model is determined by the relationship between length, force, and rigidity, and therefore it is not always necessary to adjust the size of the filler particles to the actual size.

Next, the model generation unit 22 performs mesh division on the set filler region (step S16). Further, the model generation unit 22 performs mesh division on the set polymer region (step S18).
The model generation unit 22 assigns the material constants of the polymer and filler particles called from the memory 28 to each element divided into meshes.
In step S16 and step S18, the mesh division performed on the base material region 42 and the object region 44 may be performed automatically, or may be performed by a predetermined method. Good.

Next, in the model generation unit 22, a fixing region is set for the boundary region between the filler region and the polymer region based on the information on the fixing region called from the memory 28. Then, fixing is set in this fixing region (step S20). The setting of the fixation is, for example, to share a node between the filler region and the polymer region.
Next, in the model generation unit 22, based on the information on the non-fixed area called from the memory 28, an interaction is set for the non-fixed area other than the fixed area set in step S20 (step S22). Examples of the interaction to be set include the contact described above, and the force depending on the distance from the boundary region to the non-fixed region described above. Note that the above-described peeling conditions can be set for the fixing region. In this way, a complex model, that is, a compound model is created.

Hereinafter, the simulation method of this embodiment will be described.
The complex model created as described above is subjected to simulation calculation in the calculation unit 24 based on the simulation conditions set by the condition setting unit 20.
The simulation operation is performed using a known finite element method. Since the simulation conditions used for the simulation calculation are already stored in the memory 28, the simulation conditions are called from the memory 28 and the simulation calculation is executed. The simulation conditions include, for example, information such as a forced displacement given to the complex model, an external force value, and a place to give.
The simulation operation is repeated until the solution converges. When the solution converges, the simulation calculation ends, and the deformed complex model is stored in the memory 28.

Then, the mechanical characteristic calculation unit 26 calls the result of the complex model deformed by the simulation calculation, and calculates the value of the mechanical characteristic of the compound reproduced by the complex model from the displacement of each node and the acting force. The calculated value is stored in the memory 28. The mechanical characteristics include, for example, various parameters when expressed by the apparent elastic modulus, shear rigidity, loss energy, and superelastic potential of the entire compound.
The calculated value is called from the memory 28 by the display control unit 30 and output to the display unit 16 for display on the screen. In addition, since the display unit 16 includes a printer, it can also output to a paper medium. The simulation method implemented by the processing apparatus 10 is performed as described above.

In the conventional model, the boundary region is a fixed region. However, in the simulation model creation method of this embodiment, at least one second object (object) having different material characteristics in the first object (base material). For a complex model of a complex in which there is a sticking region and a non-sticking region are provided in the boundary region. Thereby, the complex model which can adjust the interaction between the objects which comprise a complex can be obtained. In addition, a composite model in which the boundary region is a fixed region and the boundary region is a mixture of a fixed region and a non-fixed region can be used, and a composite model in accordance with an actual phenomenon can be obtained.
Further, in the simulation method, an analysis result in accordance with an actual phenomenon can be obtained by using the composite as described above.

  Note that the simulation model creation method of the present embodiment is not limited to the compound of polymer and filler particles exemplified as a composite. For example, as a composite, at least one of vulcanized rubber particles, resin balloons, glass beads, talc particles, polyethylene particles, polystyrene particles, acrylic particles, and carboxymethyl cellulose particles in the base material of the rubber composition. It may contain seeds.

  The present invention is basically configured as described above. The simulation model creation method and simulation method of the present invention have been described in detail above. However, the present invention is not limited to the above-described embodiment, and various improvements or modifications may be made without departing from the spirit of the present invention. Of course it is good.

Hereinafter, the method for creating a simulation model of the present invention will be described more specifically. In addition, this invention is not limited to the example shown below.
In this example, the composites in which the first object (base material) is a polymer and the second object (object) is filler particles are the examples 1, 2 and 4 shown in FIGS. A composite model of Comparative Example 1 was created. Then, each composite model was stretched using a finite element method, and loss energy and elastic modulus were obtained. The results are shown in Table 1 below. In Table 1, the loss energy and elastic modulus of Examples 1 and 2 are shown with the loss energy and elastic modulus of Comparative Example 1 being 100%, respectively.
The stretching condition of each complex model by the finite element method was uniaxial stretching under the periodic boundary condition. The uniaxial stretching strain was 10%.

In the composite model, the size of the base material is 50 elements × 50 elements, and the filler particles are represented by a circle having a radius of 10 elements. As material properties, the base material had an elastic modulus of 2 MPa and a Poisson's ratio of 0.49. The filler particles had an elastic modulus of 1000 MPa and a Poisson's ratio of 0.3. In addition, the fixed area was shared by nodes.
Note that the composite models of Examples 1 and 2 and Comparative Example 1 were created using the simulation model creation method described above, and thus detailed description of the creation method is omitted.

As shown in FIG. 4A, the composite model 60a of Example 1 has one filler particle 64 in the polymer material 62, and seven fixing regions in the boundary region 66 between the polymer material 62 and the filler particle 64. 68 is provided. In the boundary area 66, the areas other than the fixed area 68 are non-fixed areas.
In addition, the composite model 60b of Example 2 shown in FIG. 4B has one filler particle 64 in the polymer material 62, and three fixed regions in the boundary region 66 between the polymer material 62 and the filler particle 64. 68 is provided. In the boundary area 66, the areas other than the fixed area 68 are non-fixed areas.
Also, the composite model 60c of Comparative Example 1 shown in FIG. 4C has one filler particle 64 in the polymer material 62, and the boundary region 66 between the polymer material 62 and the filler particle 64 is entirely fixed region 68. It is a thing.

As shown in Table 1 above, when Comparative Example 1 was taken as 100%, Example 1 and Example 2 having a non-adhered region had low loss energy and elastic modulus. It can be seen that by changing the number of fixing regions as in Example 1 and Example 2, the loss energy and elastic modulus change, and the interaction between the filler particles and the base material can be adjusted.
In the case where the base material is a polymer and the filler particles are Si particles, the Si particles are not fixed to the entire region of the polymer, and the composite models of Examples 1 and 2 are partially fixed. In Examples 1 and 2, the loss energy and the elastic modulus are reduced as compared with Comparative Example 1, and this shows the same tendency as the result obtained when Si particles are present in the polymer. Comparative Example 1 corresponds to a composite model in which the base material is a polymer and the filler particles are carbon black particles. In Comparative Example 1, the case where Si particles are present in the polymer cannot be reproduced. Is clear.

10 Simulation device (processing device)
DESCRIPTION OF SYMBOLS 12 Processing part 14 Input part 16 Display part 20 Condition setting part 22 Model production | generation part 24 Calculation part 26 Mechanical property calculation part 28 Memory 30 Display control part 32 Control part 40, 40a-40d Composite body model 42 Base material area | regions 44, 52, 56 Object region 46, 47, 58, 66 Boundary region 48, 68 Adhering region 50 Non-adhering region 60a Composite model of Example 1 60b Composite model of Example 2 60c Composite model of Comparative example 1 62 Polymer region 64 Filler region

Claims (7)

A simulation model creation method for creating a composite model that can be analyzed by a computer for a composite in which at least one second object having different material properties exists in the first object,
Defining a first object region that reproduces the region of the first object and a second object region that reproduces the region of the second object;
A fixed region in which the first object region and the second object region are fixed with respect to a boundary region between the first object region and the second object region; the first object region; A method for creating a simulation model, comprising a step of setting a non-fixed area where the two object areas are not fixed.   The method for creating a simulation model according to claim 1, wherein contact is defined in the non-fixed region.   The method for creating a simulation model according to claim 1, wherein a peeling condition is set in at least a part of the fixed region.   The simulation model creation method according to claim 1, wherein a force depending on a distance from the boundary region is defined in the non-fixed region.   5. The simulation model creation method according to claim 1, wherein the second object region is set to have a higher elastic modulus than the first object region.   The simulation model creation method according to claim 1, wherein viscoelasticity is defined at least in the first object region. A simulation method using a computer to simulate a computer-analysable model,
A simulation method characterized by calculating a mechanical response of a complex by performing a deformation simulation on the complex model created by the simulation model creating method according to claim 1.