JP7199196B2 - Analysis device, analysis method, program and storage medium - Google Patents

Description

The present invention relates to an analysis device, analysis method, program and storage medium.

Conventionally, an analyzing apparatus and an analyzing method for analyzing the state quantity of an object to be analyzed are known (see, for example, Patent Document 1).

The above Patent Document 1 discloses a resin flow analysis method for obtaining flow behavior such as resin pressure, temperature, speed and shear rate in the resin flow process during plastic injection molding. In Patent Document 1, first, shape data of an object to be analyzed, molding condition data, pressure temperature boundary condition data, and the like are read. Next, the number of analysis steps is set for dividing the period from when the resin starts to flow from the inlet until it is finally filled into a plurality of analysis steps. Next, initial calculations such as initial viscosity are performed based on the set number of analysis steps. Next, the pressure distribution in the region filled with resin is calculated. Next, the velocity distribution in the entire filled area filled with all resin is calculated.

After that, in Patent Document 1, the temperature of the filled resin is calculated. In temperature calculation, first, various initial analysis parameters are set. A heat conduction, heat capacity, and heat flux matrix for each element is then formed. A matrix of elements for the entire fill area is then created. Next, the temperature distribution over the entire filling area is calculated. Here, an energy equation is used in calculating the temperature distribution. The energy equation includes a partial differential term of temperature T with respect to time t, a partial differential term of temperature T with respect to position (x, y), and the like.

Japanese Patent Application Laid-Open No. 2004-155005

Here, in the conventional flow analysis as described in Patent Literature 1, it is assumed that the flow progresses for a minute time (hereinafter referred to as Δt). In this case, the elapsed time is Δt (same for all regions) in all regions where the flow has progressed. For example, in a modeled object that is modeled by laminating materials such as resin, such as a 3D printer, the elapsed time when one layer is modeled is Δt. However, in practice, the elapsed time Δt of the portion to be shaped first differs from the elapsed time Δt of the portion to be shaped last. That is, in one modeled layer, the elapsed time Δt differs depending on the location.

In order to reflect the elapsed time for each position in the modeled layer in the flow analysis, it is necessary to divide the analysis steps so as to correspond to each position. For example, in order to form one layer, it is necessary to divide the path of a tool such as a nozzle that discharges material by Δt, and calculate the analysis steps for the number of divisions. This means that a relatively large number of analysis steps are required to calculate the temperature distribution for one built layer. Therefore, there is a problem that the analysis time becomes long.

SUMMARY OF THE INVENTION The present invention has been made to solve the problems described above, and one object of the present invention is to provide an analysis apparatus, an analysis method, a program, and a storage device capable of suppressing an analysis time from becoming long. It is to provide the medium.

In order to achieve the above object, an analysis device according to a first aspect of the present invention is an analysis device for analyzing state quantities of an analysis target, wherein a geometric model of the analysis target is mesh-divided into Calculation model creation unit that performs at least one of creating a calculation model of and reading the calculation model of the meshed analysis object from the outside, and the nodes that make up the mesh of the calculation model of the analysis object Alternatively, a time determination unit that performs at least one of determining the time of an element and reading the time from the outside, dividing the analysis of the state quantity of the analysis object into a plurality of steps, and dividing from the outside an analysis division unit that performs at least one of reading a plurality of steps of the analyzed analysis; and in the diffusion equation, the time of the time term at the node or element that constitutes the mesh of the computational model of the analysis object when the diffusion equation is discretized depends on the position of the mesh in the analysis object It is configured.

In the analysis apparatus according to the first aspect, as described above, in the diffusion equation, when the diffusion equation is discretized, the time of the time term in the node or element that constitutes the mesh of the calculation model of the analysis target is It is configured to depend on the position of the mesh on the object. As a result, since the time of the time term in the nodes or elements constituting the mesh of the computational model of the object to be analyzed depends on the position of the mesh, the elapsed time for each position is reflected in advance in the diffusion equation. As a result, for example, even if the generation time of the object to be analyzed differs within the same analysis step, it is not necessary to divide the analysis step so as to correspond to each location. As a result, it is possible to prevent the analysis time from becoming long.

In the analysis device according to the first aspect, the analysis object is preferably configured to be shaped by a programmed device, and in the discretized diffusion equation, the time of the time term at the node or element is , nodes or elements are configured to correspond to the time at which they are built. By configuring in this way, it is possible to suppress the analysis time from becoming long in the analysis of the state quantity of the object to be analyzed that is shaped by the programmed device.

In this case, preferably, in the discretized diffusion equation, the time of the time term at the node or element is configured to increase with the time that has elapsed since the portion corresponding to the node or element was shaped. . With this configuration, the time of the time term in the node or element corresponds to the time that has elapsed since the object was actually formed, so the state quantity of the object to be analyzed can be properly analyzed.

In the analysis apparatus in which the analysis object is formed by a programmed device, the state quantity of the analysis object preferably includes the temperature of the analysis object when the analysis object is formed. According to this configuration, it is possible to suppress the analysis time from becoming long in the analysis of the temperature of the analysis object when the analysis object is formed.

In this case, preferably, the thermal deformation analysis unit is further provided for analyzing the thermal deformation of the analysis object when the analysis object is shaped using the temperature of the analysis object analyzed by the state quantity analysis unit as a load condition. . With this configuration, it is possible to suppress the lengthening of the analysis time for the temperature of the analysis target when the analysis target is formed. time can be shortened.

In the analysis device according to the first aspect, preferably, in the discretized diffusion equation, the time of the time term at the node or element is set for each node or element. With this configuration, the time of the time term is set in relatively detail, so the state quantity of the object to be analyzed can be analyzed with high accuracy.

In the analysis device according to the first aspect, preferably, the initial value, the diffusion coefficient, and the boundary conditions of the state quantity of the object to be analyzed are determined, and the initial values, the diffusion coefficient, and the boundary conditions are read from the outside. and a condition creating unit that performs at least one of According to this configuration, the condition creation unit can easily (automatically without manual input) determine and read at least one of the initial value of the state quantity of the object to be analyzed, the diffusion coefficient, and the boundary condition. It can be performed.

In this case, the condition creation unit is preferably configured to determine the diffusion coefficient according to the site of the object to be analyzed. With this configuration, the state quantity of the object to be analyzed can be analyzed more accurately by, for example, making the diffusion coefficient different between the molding direction in which the object to be analyzed is molded and the direction perpendicular to the molding direction. be able to. Also, by setting the diffusion coefficient according to the time difference between adjacent elements or nodes, the state quantity of the object to be analyzed can be analyzed with higher accuracy.

In the analysis device according to the first aspect, the analysis object preferably includes a laminated article formed by laminating materials constituting the analysis object. With this configuration, unlike the conventional configuration in which a plurality of analysis steps need to be calculated to analyze the state quantity for one layer, the state quantity for one layer can be obtained by calculating one analysis step. can be analyzed. In addition, state quantities for a plurality of layers can be analyzed by one analysis step. There is no problem even if the analysis of one layer is divided into a plurality of analysis steps.

An analysis method according to a second aspect of the present invention is a computer-executed analysis method for analyzing state quantities of an analysis object, the analysis object being calculated by dividing a geometric model of the analysis object into meshes. The process of creating a model and/or loading a computational model of a meshed analysis object from the outside, and the time of the nodes or elements that make up the mesh of the computational model of the analysis object. a process of at least one of determining and reading the time from the outside, dividing the analysis of the state quantity of the analysis object into a plurality of steps, and the plurality of steps of the externally divided analysis and a process of analyzing the state quantity of the analysis object based on the diffusion equation for each of the plurality of divided steps, wherein the diffusion equation is The time of the time term in the nodes or elements constituting the mesh of the computational model of the object to be analyzed when discretized is configured to depend on the position of the mesh in the object to be analyzed.

In the analysis method according to the second aspect, as described above, in the diffusion equation, when the diffusion equation is discretized, the time of the time term in the nodes or elements that constitute the mesh of the calculation model of the analysis target is It is configured to depend on the position of the mesh on the object. As a result, since the time of the time term in the nodes or elements constituting the mesh of the computational model of the object to be analyzed depends on the position of the mesh, the elapsed time for each position is reflected in advance in the diffusion equation. This eliminates the need to divide the analysis step to correspond to each location. As a result, it is possible to provide an analysis method capable of suppressing an increase in analysis time.

A program according to a third aspect of the present invention is a program for causing a computer to analyze an object to be analyzed using an analysis method for analyzing state quantities of the object to be analyzed, the program comprising: a process of at least one of creating a computational model of the analysis object by doing, reading the computational model of the meshed analysis object from the outside, and generating the mesh of the computational model of the analysis object a process of at least one of determining the time of a constituent node or element and reading the time from the outside; dividing the analysis of the state quantity of the object to be analyzed into a plurality of steps; a process of reading a plurality of divided analysis steps; and a process of analyzing the state quantity of the object to be analyzed based on the diffusion equation for each of the plurality of divided steps. , in the diffusion equation, when the diffusion equation is discretized, the time of the time term at the node or element that constitutes the mesh of the computational model of the analysis object is configured to depend on the position of the mesh in the analysis object .

In the program according to the third aspect, as described above, in the diffusion equation, when the diffusion equation is discretized, the time of the time term at the node or element that constitutes the mesh of the calculation model of the analysis object is is configured to depend on the position of the mesh in As a result, since the time of the time term in the nodes or elements constituting the mesh of the computational model of the object to be analyzed depends on the position of the mesh, the elapsed time for each position is reflected in advance in the diffusion equation. This eliminates the need to divide the analysis step to correspond to each location. As a result, it is possible to provide a program capable of suppressing an increase in analysis time.

A storage medium according to a fourth aspect of the present invention stores a program for causing a computer to analyze the state quantity of an analysis target using an analysis method for analyzing the state quantity of the analysis target, and is computer-readable. A storage medium, at least of creating a calculation model of the analysis object by dividing the geometric model of the analysis object into meshes, and reading the mesh-divided calculation model of the analysis object from the outside. a process of at least one of: determining the time of nodes or elements that make up the mesh of the computational model of the object to be analyzed; reading the time from the outside; and the state of the object to be analyzed. A process of performing at least one of dividing the analysis of the quantity into a plurality of steps and reading the plurality of steps of the divided analysis from the outside; and a process of analyzing the state quantity of the analysis object, and in the diffusion equation, the time of the time term in the node or element that constitutes the mesh of the calculation model of the analysis object when the diffusion equation is discretized is the analysis A program is stored that is configured to be dependent on the position of the mesh on the object.

In the storage medium according to the fourth aspect, as described above, in the diffusion equation, when the diffusion equation is discretized, the time of the time term in the nodes or elements that constitute the mesh of the calculation model of the analysis target is It is configured to depend on the position of the mesh on the object. As a result, since the time of the time term in the nodes or elements constituting the mesh of the computational model of the object to be analyzed depends on the position of the mesh, the elapsed time for each position is reflected in advance in the diffusion equation. This eliminates the need to divide the analysis step to correspond to each location. As a result, it is possible to provide a storage medium capable of suppressing an increase in analysis time.

According to the present invention, as described above, it is possible to suppress the analysis time from becoming long.

It is a figure which shows the structure of the analysis apparatus by one Embodiment of this invention. FIG. 4 is a diagram showing an analysis object that has been divided into meshes; It is a figure for demonstrating the diffusion equation by the one-dimensional difference method. It is a flow diagram for explaining the analysis method according to one embodiment of the present invention. FIG. 10 is a diagram showing a toolpath; It is a figure which shows an example of an analysis object. It is a figure which shows the temperature (for 1 layer) of the analyzed analysis object. FIG. 10 is a diagram showing thermal deformation (for one layer) of the analyzed object to be analyzed; It is a figure which shows the displacement value of the analyzed analysis object. It is a figure which shows the difference between the analyzed analysis object and CAD data. FIG. 1 is a diagram (1) for explaining a simple calculation model; FIG. 2 is a diagram (2) for explaining a simple calculation model; It is a figure which shows the temperature of the analyzed analysis object (one dimension).

Embodiments embodying the present invention will be described below with reference to the drawings.

[This embodiment]
(Configuration of analysis device)
The configuration of an analysis device 100 according to this embodiment will be described with reference to FIG. Note that the analysis device 100 is configured to analyze the state quantity of the analysis object 210 formed by the programmed device 200 . For example, analysis device 100 is configured by a computer. Also programmed device 200 is, for example, a 3D printer. Further, the analysis object 210 is a laminated product formed by laminating materials (resin, metal sheath) constituting the analysis object 210 . For example, the analysis device 100 is a laminated article made of resin that is modeled by a 3D printer. Further, in the present embodiment, the state quantity of the analysis object 210 is the temperature of the analysis object 210 when the analysis object 210 is formed.

As shown in FIG. 1 , the analysis device 100 includes a calculation model creation section 10 . The calculation model creation unit 10 performs one of creating a calculation model of the analysis object 210 by dividing the analysis object 210 into meshes and reading the mesh-divided calculation model of the analysis object 210 from the outside. configured to do at least one. For example, pre-created three-dimensional CAD data of the analysis object 210 is read from the storage device 220 . Then, based on the read three-dimensional CAD data of the analysis object 210, the analysis object 210 is divided into meshes. Alternatively, data of the mesh-divided analysis object 210 is read from the storage device 220 . For example, a cylindrical analysis object 210 is divided into minute elements E, as shown in FIG. Also, the vertex of the minute element E is called a node N.

Moreover, as shown in FIG. The condition generating unit 20 determines the initial value of the state quantity, the diffusion coefficient, and the boundary conditions of the analysis object 210, and reads the initial values, the diffusion coefficient, and the boundary conditions from the outside (storage device 220). and at least one of The initial value of the state quantity of the analysis object 210 is the initial value of the temperature of the analysis object 210 . Also, the diffusion coefficient is a diffusion coefficient in a diffusion equation to be described later. Boundary conditions are temperature boundary conditions at the time of analysis (during simulation).

Moreover, the condition creating unit 20 is configured to determine the diffusion coefficient according to the part of the analysis object 210 . For example, the condition creating unit 20 determines different diffusion coefficients in the modeling direction in which the analysis object 210 is modeled and in the direction orthogonal to the modeling direction.

Here, the coefficient D corresponds to the diffusion coefficient in the diffusion equation represented by Equation (1) described later. Also, when the diffusion equation is applied to heat transfer, the diffusion equation applied to heat transfer is generally called the thermal diffusion equation. Here, the coefficient D is written as the coefficient α. The coefficient α is defined by the formula α=k/(c×ρ). α represents thermal diffusivity, k represents thermal conductivity, c represents specific heat, and ρ represents density. The condition generating unit 20 reads a diffusion coefficient (coefficient α) or reads physical properties such as thermal conductivity, specific heat, and density from the outside (storage device 220), and determines (calculates) the coefficient α from the read physical properties.

The analysis device 100 also includes a time determination unit 30 . The time determination unit 30 performs at least one of determining the time of the node N or the element E that constitutes the mesh of the computational model of the analysis object 210, and reading the time from the outside (storage device 220). is configured as

The analysis device 100 also includes an analysis division unit 40 . The analysis dividing unit 40 performs at least one of dividing the analysis of the state quantity of the analysis object 210 into a plurality of steps and reading the divided analysis steps from the outside (storage device 220). conduct. For example, in FIG. 3, which will be described later, the vertical axis represents the n-th step. Also, one step of analysis may correspond to one layer of the analysis object 210, which is a laminated product, or may correspond to a plurality of layers. By having one step of analysis correspond to multiple layers of the analysis object 210, the analysis of the analysis object 210 can be performed more quickly.

なお、上記の式（１）において、「ｕ」は、状態量（温度）を表している。「ｔ」は、時間を表している。「Ｄ」は拡散係数を表している。「Δ」は、ラプラシアンを表している。「Ｓ」は、発生量（境界から受ける熱など）を表している。なお、上記の式（１）では移流項を省略している一方、移流項を含めても問題はない。また、たとえば、有限要素法を用いて、上記の拡散方程式に基づいて、解析対象物２１０の状態量（温度）がシミュレーションされる。 The analysis device 100 also includes a state quantity analysis unit 50 . The state quantity analysis unit 50 is configured to analyze the state quantity (temperature) of the analysis object 210 based on the diffusion equation shown in the following formula (1) for each of the plurality of divided steps. .

In the above formula (1), "u" represents a state quantity (temperature). "t" represents time. "D" represents the diffusion coefficient. “Δ” represents the Laplacian. "S" represents the generated amount (such as heat received from the boundary). Although the advection term is omitted in the above equation (1), there is no problem even if the advection term is included. Also, for example, the state quantity (temperature) of the object to be analyzed 210 is simulated using the finite element method based on the above diffusion equation.

Here, in this embodiment, in the diffusion equation, the time of the time term at the node N or the element E that constitutes the mesh of the calculation model of the analysis object 210 when the diffusion equation is discretized is the mesh of the analysis object 210 is configured to depend on the position of In the conventional general diffusion equation, the state quantity u and the generated quantity S depend on the position of the mesh, while the time term (t in the above equation (1)) does not depend on the position of the mesh. That is, the state quantity u is expressed as u(t, x, y, z), and the generated quantity S is expressed as S(t, x, y, z). Time t, on the other hand, is an independent variable. In contrast, in this embodiment, the time term t is expressed as t(x, y, z).

Specifically, in this embodiment, the time (t) in the time term corresponds to the time at which the portion corresponding to node N or element E is built. Also, the time of the time term at node N or element E increases with the amount of time that has elapsed since the portion corresponding to node N or element E was formed. Also, the time of the time term at node N or element E is set for each node N or element E. FIG.

なお、上記の式（２）において、「θ」は、状態量（本実施形態では温度）を表している。また、「ｔ」および「ｘ」は、それぞれ、時間およびｘ軸（１次元）の座標を表している。また、「α」は、温度拡散率を表している。また、βは、０以上１以下の係数である。 Referring to FIG. 3, a description will be given based on the following equation (2) obtained by discretizing the diffusion equation of equation (1) above by the one-dimensional finite difference method.

In the above equation (2), "θ" represents a state quantity (temperature in this embodiment). Also, "t" and "x" represent time and x-axis (one-dimensional) coordinates, respectively. Also, "α" represents the thermal diffusivity. β is a coefficient of 0 or more and 1 or less.

For example, as shown in FIG. 3, a one-dimensional analysis object 210 model is used for explanation. The horizontal axis x in FIG. 3 represents the position of the computational grid. Also, i represents the number of the node N. Also, the vertical axis in FIG. 3 represents time. Also, n represents the n-th step when the analysis of the state quantity of the analysis object 210 is divided into a plurality of steps. FIG. 3 shows that the starting time of modeling in step _n is tn, and the finishing time of modeling in step n is tn ₊₁ . Let t _i be the time at which node i (computational grid i) is formed. In this case, in step n, the time of the part to be modeled first is t _i =t _n , so the time δt that has elapsed since this part was modeled is δt=t _n+1 −t _n . Also, since the last part to be modeled is t _i =t _n+1 , the time δt that has elapsed since this part was modeled is δt=t _n+1 −t _n+1 =0. In this way, the time (δt) of the time term increases in accordance with the time that has elapsed since the part (site) corresponding to the node N was formed. Note that although the node N has been described in FIG.

The analysis device 100 also includes a thermal deformation analysis unit 60 as shown in FIG. The thermal deformation analysis unit 60 uses the temperature of the analysis object 210 analyzed by the state quantity analysis unit 50 as a load condition to analyze the thermal deformation of the analysis object 210 when the analysis object 210 is formed. It is configured. Specifically, thermal deformation is analyzed based on the constitutive equation of an elastic body (stress=elastic modulus×strain). Discretizing and analyzing the constitutive equation of an elastic body is generally called structural analysis. Further, in the analysis of thermal deformation of the object 210 to be analyzed, time in the time term is treated as an independent variable (variable that does not depend on the position), unlike the temperature analysis described above.

Note that the calculation model creation unit 10, the condition creation unit 20, the time determination unit 30, the analysis division unit 40, the state quantity analysis unit 50, and the thermal deformation analysis unit 60 are stored in the storage medium 70 inside the analysis device 100. It is composed of a program 80 (software).

Next, an analysis method for analyzing the state quantity (temperature) of the analysis object 210 will be described with reference to FIG.

First, in step (process) S1, a computational model of the meshed analysis object 210 is created (and/or read from the outside).

Also, in step S2, the initial values of the state quantities of the object to be analyzed 210, the diffusion coefficients, and the boundary conditions are determined (and/or read from the outside).

Next, in step S3, modeling data is read from the outside. Note that the modeling data includes a material discharge portion and a trajectory P (tool path, see FIG. 5) of a laser beam when modeling the layers of the analysis object 210, and the like. For example, the modeling data is NC data specifying the three-dimensional coordinate (x, y, z) values along which the dispensing unit moves, the moving speed of the dispensing unit, etc. in an NC machine tool that moves the dispensing unit. .

Next, in step S4, the steps of the split analysis (calculation steps) are determined (and/or read from outside). Let N be the number of steps.

Next, in step S5, calculations in each step are defined. Specifically, at each step, the times of the nodes N or elements E that make up the mesh of the computational model of the analysis object 210 are determined (and/or read from the outside).

Next, in step S6, a counter i is set to 0.

Next, in step S7, the state quantity (temperature) of the analysis object 210 is analyzed based on the diffusion equation for each of the plurality of divided steps. Here, in the diffusion equation, the time (t) of the time term at node N or element E is configured to depend on the position of the mesh on the analysis object 210, so that it is divided into 1 A three-step analysis is performed.

For example, in a modeled product 230 including a model material 231 and a support material 232 as shown in FIG. 6, the state quantity (temperature) of the analysis object 210 is analyzed after one layer is generated. FIG. 7 shows that the part built last becomes relatively hot. This analysis for one layer is performed in one step.

Next, in step S8, thermal deformation of the analysis object 210 when the analysis object 210 is formed is analyzed using the analyzed temperature of the analysis object 210 as a load condition. For example, as shown in FIG. 8, the thermal load (thermal deformation) corresponding to the temperature change is calculated in one step of the modeled product 230 .

Next, in step S9, 1 is added to the counter i.

Next, in step S10, it is determined whether or not the counter i is smaller than the number N of steps. In the case of no in step S10, the process returns to step S7 and the analysis of the next step is performed. In step S10, if yes, the analysis is terminated.

Upon completion of the analysis, three-dimensional displacement values for the model 230 are determined, as shown in FIG. Also, as shown in FIG. 10, the difference between the modeled product 230 and the CAD data is obtained.

Next, with reference to FIGS. 11 and 12, the state quantity analysis method of this embodiment will be described using a simple calculation model.

A simple computational model is a one-dimensional analysis object 240, as shown in FIG. In a simple calculation model, element 1 is provided between node 1 and node 2, and element 2 is provided between node 2 and node 3. Also, an element 3 is provided between the node 3 and the node 4 . Further, as shown in FIG. 12, at time=0, the node 1 is formed in the analysis object 240 . Also, at time=1, element 1 between node 1 and node 2 is being built. Also, at time=2, elements 1 and 2, between node 1 and node 3, have been built. Also, at time=3, element 1, element 2 and element 3 between node 1 and node 4 are built.

Here, β=0 in the above equation (2). As a result, Equation (2) becomes an equation that is analyzed by an implicit method. Note that the implicit method is a method of solving simultaneous equations including unknowns, instead of using only the knowns of the current step to determine the unknowns of a new step. Also, θ _i,n =0, δx=1, β _i =0.1, and S _i =0.1. Here, in this embodiment, the time (δt) of the time term at the node N or the element E when the diffusion equation is discretized is configured to depend on the position of the mesh on the analysis object 240 . That is, δt ₁ =3, δt ₂ =2, δt ₃ =1, δt ₄ =0. The subscript i is 1, 2, 3 or 4. As a result, the following formulas (3-1) to (3-4) are obtained.

Solving the simultaneous equations of the above equations (3-1) to (3-4) yields θ _1,n+1 =0.223, θ _2,n+1 =0.189, and θ _3,n+1 =0.099. That is, there is a difference in the state quantity θ (temperature) at each position. That is, at each position, the state quantity (temperature) is affected by time.

On the other hand, when the time (δt) of the time term at the node N or the element E when discretizing the diffusion equation does not depend on the position of the mesh on the analysis target 240 as in the conventional case, δt ₁ =3, δt ₂ =3, δt ₃ =3, δt ₄ =3. In this case, the following formulas (4-1) to (4-4) are obtained.

Solving the simultaneous equations of the above equations (4-1) to (4-4) yields θ _1,n+1 =0.241, θ _2,n+1 =0.286, θ _3,n+1 =0.286, θ _{4 , n+1} =0.241. That is, the state quantity θ (temperature) becomes line symmetrical. That is, at each position, the state quantity (temperature) does not reflect the influence of time. Then, in order to reflect the effect of time on the state quantity (temperature) at each position, it is necessary to solve the diffusion equation (in multiple steps) at each of times 0, 1, 2 and 3. Therefore, the conventional method requires a long analysis time. On the other hand, in the present embodiment, it is possible to calculate the state quantity at a plurality of positions by one calculation (one step), so it is possible to shorten the analysis time.

(Effect of this embodiment)
Next, the effects of this embodiment will be described.

In this embodiment, as described above, in the diffusion equation, the time of the time term at the node N or the element E that constitutes the mesh of the calculation model of the analysis target 210 when the diffusion equation is discretized is the analysis target 210 is configured to depend on the position of the mesh in As a result, since the time of the time term at the node N or the element E that constitutes the mesh of the calculation model of the analysis object 210 depends on the position of the mesh, the elapsed time for each position is reflected in the diffusion equation in advance. . Thereby, for example, even if the generation time of the analysis object 210 differs within the same analysis step, it is not necessary to divide the analysis step so as to correspond to each location. As a result, it is possible to prevent the analysis time from becoming long.

Further, in this embodiment, as described above, the analysis object 210 is configured to be shaped by the programmed device 200, and in the discretized diffusion equation, the time term at the node N or the element E is configured to correspond to the time at which the portion corresponding to node N or element E is built. As a result, it is possible to suppress the analysis time from becoming long in the analysis of the state quantity of the analysis object 210 formed by the programmed device 200 .

Further, in this embodiment, as described above, in the discretized diffusion equation, the time of the time term at the node N or the element E is the time that has elapsed since the portion corresponding to the node N or the element E was formed. configured to grow accordingly. As a result, the time of the time term at the node N or the element E corresponds to the time that has elapsed since the model was actually formed, so the state quantity of the analysis object 210 can be properly analyzed.

Further, in the present embodiment, as described above, the state quantity of the analysis object 210 includes the temperature of the analysis object 210 when the analysis object 210 is formed. As a result, it is possible to suppress an increase in analysis time in the analysis of the temperature of the analysis object 210 when the analysis object 210 is formed.

Further, in the present embodiment, as described above, the temperature of the analysis object 210 analyzed by the state quantity analysis unit 50 is used as the load condition, and the thermal deformation of the analysis object 210 when the analysis object 210 is formed is A thermal deformation analysis unit 60 is provided for analysis. As a result, the analysis time of the temperature of the analysis object 210 when the analysis object 210 is molded is suppressed, so the total time required for the analysis of the temperature of the analysis object 210 and the analysis of thermal deformation is reduced. can save time.

Further, in this embodiment, the time of the time term at the node N or the element E is set for each node N or the element E in the discretized diffusion equation, as described above. As a result, the time of the time term is set relatively in detail, so the state quantity of the analysis object 210 can be analyzed with high accuracy.

Further, in the present embodiment, as described above, the initial value, the diffusion coefficient, and the boundary condition of the state quantity of the object to be analyzed 210 are determined, and the initial value, the diffusion coefficient, and the boundary condition are read from the outside. A condition generating unit 20 is provided for performing at least one of the above. As a result, the condition generating unit 20 easily (automatically without manual input) determines and reads at least one of the initial value of the state quantity of the object to be analyzed 210, the diffusion coefficient, and the boundary condition. be able to.

Further, in the present embodiment, as described above, the condition creation unit 20 is configured to determine the diffusion coefficient according to the site of the analysis target object 210 . As a result, for example, the state quantity of the analysis target object 210 can be analyzed with higher accuracy by making the diffusion coefficient different between the modeling direction in which the analysis target object 210 is shaped and the direction perpendicular to the modeling direction. can. Also, by setting the diffusion coefficient according to the time difference between adjacent elements or nodes, the state quantity of the analysis object 210 can be analyzed with higher accuracy.

Further, in the present embodiment, as described above, the analysis object 210 includes a laminated product formed by laminating the materials forming the analysis object 210 . As a result, unlike the conventional configuration that requires multiple analysis steps to analyze the state quantity of one layer, it is possible to analyze the state quantity of one layer by calculating one analysis step. can. In addition, state quantities for a plurality of layers can be analyzed by one analysis step. There is no problem even if the analysis of one layer is divided into a plurality of analysis steps.

[Example]
An example in which the analysis method of this embodiment is applied to the finite element method (the number of elements is 11) will be described with reference to FIG. The analysis target object 250 of the embodiment is formed in order from the upper right corner to the lower left corner. In the example shown in FIG. 13, the initial temperature was 200°C, the ambient temperature was 100°C, and the heat transfer coefficient was 100 W/m ² /°C. Also, the specific gravity, specific heat and thermal conductivity of the material were set to 1.05×10 ³ kg/m ³ , 1.3×10 ³ J/kg/°C and 0.2 W/m/°C, respectively. As a result, it was reproduced that the temperature of the early shaped elements was low and the temperature of the late shaped elements was high. By using the analysis method of this embodiment, the analysis of 11 elements can be performed in one step.

[Modification]
It should be noted that the embodiments and examples disclosed this time should be considered illustrative in all respects and not restrictive. The scope of the present invention is indicated by the scope of the claims rather than the description of the above-described embodiments and examples, and includes all changes (modifications) within the meaning and scope equivalent to the scope of the claims.

For example, in the above-described embodiments, an example in which the object to be analyzed is a modeled product (laminated product) that is modeled by a programmed device has been described, but the present invention is not limited to this. The present invention can also be applied to analysis of things other than shaped products, such as flow.

Further, in the above embodiment, an example in which the temperature of the object to be analyzed is analyzed has been shown, but the present invention is not limited to this. The present invention can also be applied to the analysis of magnetic fields other than temperature.

Also, in the above embodiment, the time of the time term at the node or element is set for each node or element, but the present invention is not limited to this. For example, the time of the time term may be set for each of multiple nodes or multiple elements.

Further, in the above embodiment, the calculation model creation unit, the condition creation unit, the time determination unit, the analysis division unit, the state quantity analysis unit, and the thermal deformation analysis unit are executed by a program (software) stored in the storage medium of the analysis device. Although a configured example is shown, the present invention is not limited to this. For example, the above program may be stored in a storage medium provided separately from the analysis device.

10 calculation model creation unit 20 condition creation unit 30 time determination unit 40 analysis division unit 50 state quantity analysis unit 60 thermal deformation analysis unit 70 storage medium 80 program 100 analysis device 200 device 210, 230, 240, 250 analysis object E element N node

Claims (12)

An analysis device for analyzing the state quantity of an analysis object,
at least one of creating a calculation model of the analysis object by dividing the shape model of the analysis object into meshes, and reading the calculation model of the mesh-divided analysis object from the outside. A calculation model creation unit that performs
a time determination unit that performs at least one of determining the time of nodes or elements that constitute the mesh of the computational model of the analysis object and reading the time from the outside;
an analysis dividing unit that performs at least one of dividing the analysis of the state quantity of the analysis object into a plurality of steps and reading the plurality of steps of the divided analysis from the outside;
a state quantity analysis unit that analyzes the state quantity of the analysis object based on a diffusion equation for each of the plurality of divided steps;
In the diffusion equation, when the diffusion equation is discretized, the time of the time term in the node or the element that constitutes the mesh of the calculation model of the analysis object depends on the position of the mesh in the analysis object. An analysis device configured to The analysis object is configured to be shaped by a programmed device,
2. The method of claim 1, wherein in the discretized diffusion equation, the time of the time term at the node or element is configured to correspond to the time at which the portion corresponding to the node or element is shaped. analysis equipment. In the discretized diffusion equation, the time of the time term at the node or the element is configured to increase according to the time that has elapsed since the part corresponding to the node or the element was formed. The analysis device according to claim 2. The analysis apparatus according to claim 2 or 3, wherein the state quantity of the analysis object includes the temperature of the analysis object when the analysis object is formed. Further comprising a thermal deformation analysis unit that analyzes thermal deformation of the analysis object when the analysis object is shaped using the temperature of the analysis object analyzed by the state quantity analysis unit as a load condition. Item 5. The analysis device according to item 4. 5. The analysis apparatus according to claim 1, wherein in said discretized diffusion equation, the time of the time term at said node or said element is set for each said node or said element. At least one of determining an initial value of the state quantity of the object to be analyzed, the diffusion coefficient, and the boundary condition, and reading the initial value, the diffusion coefficient, and the boundary condition from the outside. The analysis device according to any one of claims 1 to 6, further comprising a condition generation unit. 8. The analyzing apparatus according to claim 7, wherein said condition creating unit is configured to determine said diffusion coefficient according to a site of said object to be analyzed. The analyzing apparatus according to any one of claims 1 to 8, wherein the analysis object includes a laminated product formed by laminating materials constituting the analysis object. A computer-executed analysis method for analyzing state quantities of an object to be analyzed, comprising:
at least one of creating a calculation model of the analysis object by dividing the shape model of the analysis object into meshes, and reading the calculation model of the mesh-divided analysis object from the outside. process to do and
a process of at least one of determining times of nodes or elements that make up a mesh of the computational model of the object to be analyzed and reading the times externally;
a process of performing at least one of dividing the analysis of the state quantity of the analysis object into a plurality of steps and reading the divided plurality of steps of the analysis from the outside;
A process of analyzing the state quantity of the analysis object based on a diffusion equation for each of the plurality of divided steps,
In the diffusion equation, when the diffusion equation is discretized, the time of the time term in the node or the element that constitutes the mesh of the calculation model of the analysis object depends on the position of the mesh in the analysis object. A parsing method that is configured to A program for causing a computer to analyze the object to be analyzed using an analysis method for analyzing the state quantity of the object to be analyzed,
at least one of creating a calculation model of the analysis object by dividing the shape model of the analysis object into meshes, and reading the calculation model of the mesh-divided analysis object from the outside. process to do and
a process of at least one of determining times of nodes or elements that make up a mesh of the computational model of the object to be analyzed and reading the times externally;
a process of performing at least one of dividing the analysis of the state quantity of the analysis object into a plurality of steps and reading the divided plurality of steps of the analysis from the outside;
A process of analyzing the state quantity of the analysis object based on a diffusion equation for each of the plurality of divided steps,
In the diffusion equation, when the diffusion equation is discretized, the time of the time term in the node or the element that constitutes the mesh of the calculation model of the analysis object depends on the position of the mesh in the analysis object. A program that is configured to A computer-readable storage medium storing a program for causing a computer to analyze the state quantity of the analysis target using an analysis method for analyzing the state quantity of the analysis target,
at least one of creating a calculation model of the analysis object by dividing the shape model of the analysis object into meshes, and reading the calculation model of the mesh-divided analysis object from the outside. process to do and
a process of at least one of determining times of nodes or elements that make up a mesh of the computational model of the object to be analyzed and reading the times externally;
a process of performing at least one of dividing the analysis of the state quantity of the analysis object into a plurality of steps and reading the divided plurality of steps of the analysis from the outside;
A process of analyzing the state quantity of the analysis object based on a diffusion equation for each of the plurality of divided steps,
In the diffusion equation, when the diffusion equation is discretized, the time of the time term in the node or the element that constitutes the mesh of the calculation model of the analysis object depends on the position of the mesh in the analysis object. A storage medium storing a program configured to

Images