JP2015099573A - Wire harness analysis method, wire harness analyzing apparatus, and program - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a wire harness analysis method configured to take into consideration an internal state of a member constituting a wire harness or the work environment in assembling the wire harness on a vehicle.SOLUTION: A wire harness analysis method includes: referring to a physical property value for each of separate elements of a wire harness (step S71); correcting the physical property value on the basis of both or one of an internal state of a member constituting the wire harness to be analyzed and the work environment in assembling the wire harness on a vehicle (step S72); and calculating a shape of the wire harness, on the basis of the corrected physical property value (step S73).

Description

  The present invention relates to a wire harness analysis device, a wire harness analysis method, and a program for analyzing a wire harness.

  Before manufacturing a real wire harness, a virtual wire harness is modeled on a computer, and the virtual wire harness is studied at the design stage. Examples of this type of simulation are disclosed in Patent Documents 1 to 4.

JP 2003-141949 A JP 2004-362542 A JP 2005-85191 A JP 2011-22594 A

  In the simulation methods described in Patent Documents 1 to 4, a plurality of electric wires constituting the wire harness, a protective member (tape, corrugated tube, protector, etc.) for protecting the electric wires, and various types of bundling the plurality of electric wires. Fixtures (clamps, clips, grommets, fastening bands, tapes, etc.) are regarded as one continuous body, and this wire harness is discretized into several elements. Various physical property values measured using a real wire harness are assigned to each element. For example, in Patent Document 2, length l, sectional area A, sectional secondary moment I, sectional secondary pole moment J, longitudinal elastic modulus E, and transverse elastic modulus G are assigned as physical property values.

  However, as described above, in the approximate model that captures the electric wire, the protective member, and the fixture together, the influence of the electric wire, the protective member, and the fixture constituting the wire harness, for example, how to bundle the electric wires, The influence of the tightening degree, the rigidity of the protective member, etc. is not reflected in the analysis result. For example, the wire harness has a portion wound with tape, and the performance of the wire harness changes depending on the tension applied by the tape to the portion wound with tape. As described above, the internal state of the members constituting the wire harness (how to bundle the electric wires arranged in each section, the degree to which the fixture tightens the electric wires, the rigidity of the protective member, the state of the members constituting the wire harness, The analysis result that does not take into account the situation)) cannot accurately perform a study assuming a real wire harness. On the other hand, it is conceivable that each electric wire, protective member, and fixture are regarded as objects of analysis, and each electric wire, protective member, and fixture are discretized into several elements. However, this method requires a large amount of computation. End up.

  In the above approximate model, the work environment when the wire harness is routed in the vehicle is not taken into consideration. An example of a working environment where the wire harness is routed in the vehicle is the outside air temperature. When the outside air temperature is high and when the outside air temperature is low, the rigidity of the wire harness is different, and when the wire harness is assembled to the vehicle, the load required for the worker to assemble the wire harness is different. In the analysis result in which the work environment in which the wire harness is assembled in this manner (the wire harness or the environment around the vehicle when the wire harness is assembled in the vehicle is referred to as the work environment) is not taken into consideration, the actual wire harness is used. The assumed study cannot be performed accurately.

  The present invention has been made in view of the above-described circumstances. The purpose of the present invention is to adopt an approximate model that captures a wire harness as one continuous body, while the internal state of members constituting the wire harness or the wire harness is An object of the present invention is to provide a wire harness analysis method, a wire harness analysis device, and a program that can take into consideration the work environment when assembled to a vehicle.

In order to achieve the above-described object, the wire harness analysis method according to the present invention is characterized by the following (1) to (10).
(1) A wire harness analysis method for calculating a shape of a wire harness that is sequentially deformed using a computer,
A step of referring to a physical property value for each small region of the wire harness that has been discretized;
A correction step of correcting the physical property value based on one or both of an internal situation of a member constituting the wire harness to be analyzed and a work environment when the wire harness to be analyzed is assembled to a vehicle When,
Based on the physical property values corrected in the correction step, a calculation step for calculating the shape of the wire harness;
Be provided.
(2) A wire harness analysis method configured as described in (1) above,
In the correction step, for each element of the wire harness discretized by an element, the physical property value is corrected based on one or both of the internal situation and the work environment.
(3) A wire harness analysis method configured as described in (2) above,
The internal state includes a ratio of the electric wire occupying on a cross section of the wire harness to the cross sectional area of the wire harness.
about.
(4) A wire harness analysis method configured as described in (2) or (3) above,
The internal situation includes a variation in the angle in the axial direction of the electric wire with respect to the axial direction of the wire harness at an arbitrary position in the longitudinal direction of the wire harness.
about.
(5) A wire harness analysis method having any one of the constitutions (2) to (4),
The internal situation includes a tension acting on the band-shaped member when a portion of the wire harness is wound by the band-shaped member, or a tightening force exerted on the position by the band-shaped member.
about.
(6) The wire harness analysis method having any one of the configurations (2) to (5),
The internal situation includes the rigidity of a protective member that covers a portion of the wire harness.
about.
(7) The wire harness analysis method having any one of the configurations (2) to (6),
The internal situation includes the amount of wire movement allowed within the wire bundle,
about.
(8) The wire harness analysis method having any one of the configurations (2) to (7),
The internal situation includes the ratio of the short axis to the long axis of the cross section of the wire bundle at an arbitrary position in the longitudinal direction of the wire harness.
about.
(9) The wire harness analysis method having any one of the above configurations (2) to (8),
The internal situation includes a numerical value for evaluating a location where a thick wire in the wire bundle is located,
about.
(10) The wire harness analysis method having any one of the above configurations (2) to (9),
The external environment includes outside air temperature or humidity under a situation where the wire harness is assembled to a vehicle.
about.

In order to achieve the above-described object, the wire harness analyzer according to the present invention is characterized by the following (11).
(11) A wire harness analyzer that calculates the shape of a wire harness that is sequentially deformed,
A physical property value for each small region of the wire harness that has been discretized is recorded, and an internal situation of a member constituting the wire harness to be analyzed and the wire harness to be analyzed are assembled to a vehicle. A recording section in which one or both of the work environments are recorded;
For each small region of the wire harness that has been discretized, either the internal state of the members constituting the wire harness to be analyzed and the work environment when the wire harness to be analyzed is assembled to a vehicle Or a calculation unit that corrects a physical property value based on both, and calculates the shape of the wire harness based on the corrected physical property value;
Be provided.

In order to achieve the above-described object, a program according to the present invention is characterized by the following (12).
(12) A program for causing a computer to execute each step of the wire harness analysis method according to any one of (1) to (10).

According to the wire harness analysis method having the configuration of (1) to (10), the wire harness analysis apparatus having the configuration of (11), and the program having the configuration of (12), the wire harness is captured as one continuous body. While adopting the approximate model, it is possible to consider the internal conditions of the members constituting the wire harness or the work environment when the wire harness is assembled to the vehicle. For this reason, since an analysis closer to the actual wire harness can be performed, it is possible to examine the work by the worker according to the work environment in consideration of the analysis result.
Or, the manufacturing result of the wire harness is managed by utilizing the analysis result, for example, the tension of the tape when the tape is wound around a predetermined portion of the wire harness is controlled to a predetermined value, so that the operator feels uncomfortable. In addition, it is possible to manufacture a wire harness having stable characteristics with no variation in the accuracy of work performed by an operator.

  According to the wire harness analysis method, the wire harness analysis apparatus, and the program of the present invention, the internal situation of the members constituting the wire harness or the wire harness is assembled to the vehicle while adopting the approximate model that captures the wire harness as one continuous body. The working environment can be taken into account.

  The present invention has been briefly described above. Further, the details of the present invention will be further clarified by reading through a mode for carrying out the invention described below (hereinafter referred to as “embodiment”) with reference to the accompanying drawings. .

FIG. 1 (A) is a diagram showing the appearance of the wire harness, FIG. 1 (B) is a diagram showing a state where the wire harness of FIG. 1 (A) is discretized, and FIG. It is the figure which represented the wire harness of FIG. 1 (A) with the beam element and the node. FIG. 2 is a diagram for explaining the degree of freedom in the wire harness represented by beam elements and nodes. FIG. 3 is a cross-sectional view of the wire harness cut along a plane perpendicular to the axial direction. 4A is a cross-sectional view when the wire harness is cut along a plane passing through the axial direction, and FIG. 4B is a cross-sectional view as viewed in the direction of the line IVB-IVB in FIG. It is. FIG. 5 is a cross-sectional view of the wire harness cut along a plane perpendicular to the axial direction. FIG. 6 is a block diagram illustrating a hardware configuration of the wire harness analyzer according to the embodiment of the present invention. FIG. 7 is a flowchart illustrating a procedure of an analysis method performed by the wire harness analyzer according to the embodiment of the present invention.

  Specific embodiments relating to the present invention will be described below with reference to the drawings.

[Wire harness analysis algorithm]
Various methods have been proposed for an analysis method for calculating the shape of the wire harness. For example, Japanese Patent Application Laid-Open No. 2004-362542, filed by the applicant of the present application and already published as an application, describes in detail. In the embodiment of the present invention, a mode in which a wire harness is modeled using the technique described in Japanese Patent Application Laid-Open No. 2004-362542 will be described. Hereinafter, a method for modeling a wire harness according to an embodiment of the present invention will be described. FIG. 1 (A) is a diagram showing the appearance of the wire harness, FIG. 1 (B) is a diagram showing a state where the wire harness of FIG. 1 (A) is discretized, and FIG. It is the figure which represented the wire harness of FIG. 1 (A) with the beam element and the node.

  First, in the embodiment of the present invention, the shape of the wire harness is calculated using a finite element method. The numerical analysis method of the present invention is not limited to the one using the finite element method.

  In the present embodiment, first, the wire harness is discretized. As shown in FIG. 1A, the wire harness 1 in which a plurality of electric wires 11 are bundled by a protective member such as a tape 12 can be regarded as a continuous body. Next, as shown in FIG. 1B, such a wire harness 1 is divided (discretized) into several beam elements C1, C2, C3,. The wire harness can be regarded as a connection of a finite number of beam elements.

  Therefore, as shown in FIG. 1C, the wire harness is represented as a plurality of beam elements C1, C2, C3,... Coupled by a plurality of nodes N1, N2, N3,. The physical properties required for the beam elements are as follows.

Length l (see Fig. 1 (B))
Cross section A (See Fig. 1 (B))
Sectional secondary moment I
Cross section secondary pole moment J
Longitudinal elastic modulus E
Transverse elastic modulus G
Although not directly expressed in the physical property values, density ρ, Poisson's ratio ν, and the like are also used to obtain them.

  As shown in FIG. 2, each beam element C (C1, C2, C3,...) Has two nodes α (α1, α2, α3,...). In the three-dimensional space, the node α has three translation components (x component, y component, z component) of the translation vector and three rotation components (u component, v component, w component) of the rotation axis vector. , Has a total of 6 degrees of freedom.

In the figure,
Ui (x): x-component of translation vector of i-th element Ui (y): y-component of translation vector of i-th element Ui (z): z-component of translation vector of i-th element θi (u): i-th element U component θi (v) of the rotation vector of i: v component of the rotation vector of the i-th element θi (w): w component αi of the rotation vector of the i-th element indicates the node of the element.

  In this manner, a finite element model of a wire harness that is divided into an arbitrary number of elements is created.

Further, in structural mechanics with large deformation such as a wire harness, the equilibrium equation of the finite element method generally takes the form of the following equation.
([K] + [KG]) {x} = {F} (1)
Here, [K]: Overall stiffness matrix, [KG]: Overall geometric stiffness matrix, {x}:
Displacement vector, {F}: Load vector (also called force vector)

  Numerical analysis of the mathematical formula (1) (for details of the numerical analysis, refer to Japanese Patent Application Laid-Open No. 2004-362542), an unknown displacement vector {x} is obtained to obtain the shape of the target wire harness. Can be calculated. In Japanese Patent Application Laid-Open No. 2004-362542, the wire harness is assumed to be an elastic body in the numerical analysis, but the present invention can also be applied to the numerical analysis assumed to be an elastic-plastic body.

[Calculation method of physical property values]
In the conventional analysis method, the measured physical property value is directly assigned to the element. In the embodiment of the present invention, the measured physical property value is not assigned to an element as it is, but the internal state of a member constituting a wire harness to be analyzed or the wire harness to be analyzed is assembled to a vehicle. The physical property value corrected in consideration of the working environment is assigned to the element.

[Correction taking into account the internal conditions of the members constituting the wire harness]
In a real wire harness, even if two sections are composed of the same wire, protective member, and fixture, how to bundle the wires arranged in each section, the degree to which the fixture tightens the wire, Depending on the situation of the members constituting the wire harness such as rigidity (this situation is referred to as “internal situation of the wire harness” in this specification), the physical property values in the two sections will be different. The present inventor does not consider the internal state of the wire harness in the conventional analysis method in which the wire harness is constructed by the approximate model in which the electric wire, the protective member, and the fixture are gathered together. The knowledge that the performance of the calculated wire harness and the performance of the actual wire harness to be modeled would be shifted.

Among the internal conditions of the wire harness, the following factors can be cited as factors that greatly affect the physical property values.
(A) Ratio of the cross-sectional area of the wire occupying the cross-section with the wire harness (wire occupancy rate)
(B) Variation in angle in the axial direction of the electric wire with respect to the axial direction of the wire harness (C) When a portion having the wire harness is wound by the strip member, the tension acting on the strip member, or the strip member is applied to the portion Tightening force exerted (D) Rigidity of protective member covering a portion of the wire harness (E) Distance that the wire can move in the cross section of the wire harness (F) Ratio of short axis to long axis in the cross section of the wire bundle (G ) Distance between the center of the wire bundle and the center of the thick wire in the wire bundle in the cross section of the wire bundle

  Hereinafter, a method for calculating a physical property value in consideration of the factors (A) to (G) will be described in detail.

[A. Method of reflecting wire occupancy ratio to physical property values]
First, the wire occupation ratio will be described. FIG. 3 is a cross-sectional view of the wire harness cut along a plane perpendicular to the axial direction. The cross-sectional view of the wire harness shown in FIG. 3 shows a state in which three electric wires W1, W2, and W3 having different wire diameters are bundled and the outer circumferences of these electric wires W1, W2, and W3 are covered with the protective member P. . In such a cross section of the wire harness, the wire occupation ratio O [%] is calculated as follows.
O = S1 / S0 × 100
S1: Total cross-sectional area of all electric wires W1, W2, W3 S0: Cross-sectional area inside protective member

  In a real wire harness, the wire occupation ratio O often takes a value between 55% and 75%. The electric wire occupancy O can be calculated by measuring the wire harness to be modeled using a caliper, a measure, a weight scale or the like.

  In the wire harness, when the wire occupation rate O of the electric wires accommodated in a certain protective member P is different, the wire occupation rate O is larger (in other words, the wire is “dense”) smaller (in other words, Bending rigidity EI and other various rigidity (shaft rigidity EA, torsional rigidity GJ, shear rigidity GA) are larger than those of “sparser” electric wires. Based on this point, the wire occupancy O (O1, O2, O3,...) Is assigned to each element C (C1, C2, C3,...), And the wire occupancy factor Ok is set for each numerical range of the wire occupancy O. Create a defined correspondence table. Table 1 shows an example of the correspondence table.

  The wire occupancy coefficient Ok is associated with each of the four numerical ranges of the wire occupancy O (O <55%, 55% ≦ O <65%, 65% ≦ O <75%, 75% ≦ O), − Numerical values of 0.15, 0, 0.6, and 1.2 are set. The electric wire occupation ratio coefficient Ok is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the physical property value is assigned to the element C by the electrical wire occupancy coefficient Ok so that the physical property value increases if the electrical wire occupancy O is large, and the physical property value decreases if the electrical wire occupancy O is small. By correcting the physical property value, the wire occupation ratio O is reflected in the physical property value. Here, the electric wire occupancy rate O is divided into four numerical ranges, but the numerical value range can be further divided, or the electric wire occupancy factor Ok can be expressed as a function of the electric wire occupancy O. A specific calculation formula for correcting the physical property value by the electric wire occupancy coefficient Ok will be described in detail later.

[B. Method of reflecting angle variation in physical property values]
First, the variation in angle will be described. 4A is a cross-sectional view when the wire harness is cut along a plane passing through the axial direction, and FIG. 4B is a cross-sectional view as viewed in the direction of the line IVB-IVB in FIG. It is. The cross-sectional view of the wire harness shown in FIG. 4 shows a state in which the electric wires W1, W2, and W3 are bundled and the outer circumferences of these electric wires W1, W2, and W3 are covered with the protective member P. At an arbitrary position in the longitudinal direction of the wire harness shown in FIG. 4, the extending direction of the electric wire W1 and the axial direction of the protective member P (the normal direction in the cross section of the protective member P shown by the one-dot chain line in FIG. 4) And the electric wire W1 are arranged so that the angle ω1 (where ω1 is a numerical value satisfying 0 ≦ ω1 <90 °). Further, the electric wire W2 is arranged so that the extending direction of the electric wire W2 and the axial direction of the protective member P form an angle ω2 (where ω2 is a numerical value satisfying 0 ≦ ω2 <90 °). Further, the electric wire W3 is arranged so as to make an angle ω3 (where ω3 is a numerical value satisfying 0 ≦ ω3 <90 °) with respect to the extending direction of the electric wire W3 and the axial direction of the protective member P. As FIG. 4 shows, the electric wire which comprises a wire harness is not necessarily arrange | positioned in parallel with the axial direction of a wire harness, but is arrange | positioned so that it may cross | intersect another electric wire. In such a cross section at an arbitrary position in the longitudinal direction of the wire harness, the angular variation R is calculated as the root mean square of the angle ωn as follows.
R (z) = {(ω1 2 + ω2 2 + ω3 2, ... + ωn 2, ... + ωN 2) / N} 1/2
ωn: Angle n: 1, 2, 3,... N natural number that the wire n makes with respect to the axial direction of the wire harness z: Variable that defines the position of the wire harness in the longitudinal direction

  In the actual wire harness, the angle variation R can be calculated by measuring the wire harness to be modeled using a caliper, a measure, a weight scale, or the like.

  In the wire harness, when the angle variation R of the electric wires accommodated in a certain protective member P is different, the wire having a smaller angle variation R (in other words, more wires extending in the same direction) is larger. Bending rigidity EI and other various rigidity (axial rigidity EA, torsional rigidity GJ, shear rigidity GA) are larger than those (in other words, those having fewer wires extending in the same direction). Based on this point, an angle variation R (R1, R2, R3,...) Is assigned to each element C (C1, C2, C3,...), And an angle variation coefficient Rk is set for each numerical range of the angle variation R. Create a defined correspondence table. Table 2 shows an example of the correspondence table.

  Angle variation coefficient Rk corresponds to each of four numerical ranges of angle variation R (0 ° ≦ R <2 °, 2 ° ≦ R <4 °, 4 ° ≦ R <6 °, 6 ° ≦ R). In addition, numerical values of 0, -0.05, -0.1, and -0.15 are determined. The angle variation coefficient Rk is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the property value is assigned to the element C by the angle variation coefficient Rk so that the physical property value increases when the angle variation R is small, and the physical property value decreases when the angle variation R is large. By correcting the physical property value, the angle variation R is reflected in the physical property value. Here, the angle variation R is divided into four numerical ranges. However, the numerical range may be further divided, or the angle variation coefficient Rk may be expressed as a function of the angle variation R. A specific calculation formula for correcting the physical property value by the angle variation coefficient Rk will be described in detail later.

[C. Method of reflecting the tension acting on the belt-like member on the physical property value]
First, the strip member will be described. As for the electric wire which comprises a wire harness, the outer periphery of an electric wire is tied with the fastening band, and the state which was bundled is maintained. Or the branch line branched from the trunk line of the electric wire is wound on a sheet or tape, so that the state where the branch lines are bundled is maintained. Alternatively, the bundled state is maintained by accommodating the electric wires in the hollow cylindrical tube. The band-shaped member of the present invention corresponds to a fastening band, a sheet, a tape, or a tube for maintaining a state in which the electric wires are bundled.

  Next, the tension acting on the band member and the tightening force exerted by the band member will be described. FIG. 5 is a cross-sectional view of the wire harness cut along a plane perpendicular to the axial direction. The cross-sectional view of the wire harness shown in FIG. 5 shows a state in which three electric wires W1, W2, and W3 having different wire diameters are bundled and the outer circumferences of these electric wires W1, W2, and W3 are bound by the fastening band B. . During the tightening operation, a tension T acts on the tightening band B in a circular tangential direction surrounding the electric wires W1, W2, and W3, and the tightening band B contacts the electric wires W1, W2, and W3 in the process of pulling the tightening band B. . At this time, according to the tension T acting on the tightening band B, a circular vertical normal force surrounding the electric wires W1, W2, W3 acts on the tightening band B from the electric wires W1, W2, W3. A force that resists the vertical drag acts on the electric wires W1, W2, and W3 as the tightening force U.

  When attaching the fastening band to the electric wire, a bundling tool such as a band gun is used. The tightening strength by this tightening band is specified in the product specification. For this reason, the tension T can be specified by referring to the tightening strength specified in the product specification.

  In the wire harness, when the tension T and the tightening force U acting on the tightening band B are different, the one having the larger tension T and the tightening force U (in other words, the one in which the electric wire is “dense”) is smaller (in other words, Bending rigidity EI and other various rigidity (shaft rigidity EA, torsional rigidity GJ, shear rigidity GA) are larger than those of “sparser” electric wires. In consideration of this point, each element C (C1, C2, C3,...) Is assigned a tension T (T1, T2, T3,...) And a tightening force (U1, U2, U3,...). A correspondence table that defines the tension coefficient Tk and the tightening coefficient Uk for each numerical value range of U is created. Table 3 shows an example of the correspondence table. Table 3 shows a correspondence table in which the tension coefficient Tk is determined for each numerical value range of the tension T.

  The tension coefficient Tk is associated with each of the four numerical ranges of the tension T (T <1N, 1N ≦ T <3N, 3N ≦ T <5N, 5N ≦ T), −0.05, 0, 0.05, A numerical value of 0.1 is set. The tension coefficient Tk is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the physical property value assigned to the element C is corrected by the tension coefficient Tk so that the physical property value increases when the tension T is large, and the physical property value decreases when the tension T is small. Thus, the tension T is reflected in the physical property value. Here, the tension T is divided into four numerical ranges, but the numerical ranges can be further divided, or the tension coefficient Tk can be expressed as a function of the tension T. A specific calculation formula for correcting the physical property value by the tension coefficient Tk will be described in detail later.

[D. Method of reflecting the rigidity of the protective member to the physical property value]
First, there are a tape, a sheet, a corrugated tube, and the like as a protective member for protecting an electric wire. In the present embodiment, a case in which the rigidity of these protective members increases in the order of tape, sheet, and corrugated tube will be described as an example, and a method of reflecting the rigidity of the protective member on the physical property values will be described.

  In the wire harness, when the type of the protective member for protecting the electric wire is different, the bending rigidity EI and other various rigidity (axial rigidity EA, torsional rigidity GJ, shear rigidity GA) increase as the portion is covered by the highly rigid protective member. . Based on this point, each element C (C1, C2, C3,...) Is assigned an identifier V (V1, V2, V3,...) That identifies the type of protection member, and a protection member coefficient Vk is determined for each identifier V. Create a corresponding table. Table 4 shows an example of the correspondence table.

  The protection member coefficient Vk is set to a numerical value of 0, 0.1, or 0.3 in association with each identifier V (tape, sheet, corrugated tube) that identifies the type of protection member. The protection member coefficient Vk is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the rigidity of the protective member is reflected in the physical property value by correcting the physical property value assigned to the element C by the protective member coefficient Vk determined for each tape, sheet, and corrugated tube. A specific calculation formula for correcting the physical property value by the protective member coefficient Vk will be described in detail later.

[E. Method of reflecting the amount of wire movement]
First, the movement amount of the electric wire will be described. In the actual wire harness, the wire bundle often takes a value between 55% and 75%. This can be considered that the electric wire which comprises an electric wire bundle can move the remaining space which is not occupied by the electric wire. The movement amount D of the electric wire is calculated as an average of distances when the electric wires move in the space in a cross section with the wire harness. The movement amount D [mm] of the electric wire is calculated as follows.
D = {(D1 ² + D2 ² + D3 ² ,... + Dn ² ,... + DN ² ) / N} ^1/2
Dn: Distance that electric wire Wn can move n: n: 1, 2, 3, ... N natural number

  In an actual wire harness, the amount of movement D of the electric wire can be calculated by measuring the wire harness to be modeled using a caliper, a measure, a weight scale, or the like.

  In the wire harness, when the movement amount D of the electric wire accommodated in a certain protective member P is different, the bending rigidity EI and other various rigidity (axial rigidity EA, torsion rigidity GJ, Shear rigidity GA) increases. Based on this point, the wire movement amount D (D1, D2, D3,...) Is assigned to each element C (C1, C2, C3,...), And the wire movement amount for each numerical range of the wire movement amount D. A correspondence table defining the coefficient Dk is created. Table 5 shows an example of the correspondence table.

  The movement amount coefficient Dk of the electric wire corresponds to each of four numerical ranges of the movement amount D of the electric wire (0 ≦ D <1 mm, 1 mm ≦ D <3 mm, 3 mm ≦ D <5 mm, 5 mm ≦ D). Numerical values of 0.15, -0.3, and -0.45 are determined. The movement amount coefficient Dk of the electric wire is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the element C is given by the moving amount coefficient Dk of the electric wire so that the physical property value becomes larger if the moving amount D of the electric wire is small and the physical property value becomes smaller if the moving amount D of the other electric wire is large. By correcting the assigned physical property value, the movement amount D of the electric wire is reflected in the physical property value. In this example, the movement amount D of the electric wire is divided into four numerical ranges. However, the numerical value range can be further divided, or the movement amount coefficient Dk of the electric wire can be expressed as a function of the movement amount D of the electric wire. is there. A specific calculation formula for correcting the physical property value by the movement amount coefficient Dk of the electric wire will be described in detail later.

[F. Method of reflecting the ratio of the short axis to the long axis in the cross section of the wire bundle]
In the actual wire harness, the cross section of the wire bundle is not necessarily a perfect circle, and there is a portion that is an ellipse. In the actual wire harness, the ratio of the short axis to the long axis (hereinafter simply referred to as the short and long axis ratio of the electric wire bundle) F in the electric wire bundle having an elliptical cross section can be calculated by the following equation.
F = Ls / Ll
Ls: Length of short axis Ll: Length of long axis

  The short / major axis ratio F of the wire bundle can be calculated by measuring the wire harness to be modeled using a caliper, a measure, a weight scale or the like.

  In a wire harness, when the cross-sectional area of a wire bundle accommodated in a certain protective member P is constant and the minor / major axis ratio F is different, the wire bundle having a larger minor / major axis ratio F (in other words, an electric wire bundle) ) Has a larger bending rigidity EI and a smaller torsional rigidity GJ than a smaller one (in other words, a cross section of the wire bundle is away from the perfect circle). Based on this point, the short and long axis ratio F (F1, F2, F3,...) Of the wire bundle is assigned to each element C (C1, C2, C3,...), And the numerical range of the short and long axis ratio F of the wire bundle is assigned. A correspondence table (in the case of bending rigidity) in which the short / long axis ratio coefficient Fk of the wire bundle is determined for each is created. Table 6 shows an example of the correspondence table.

  The short and long axis ratio coefficient Fk of the wire bundle is determined by four numerical ranges of the short and long axis ratio F of the wire bundle (F <0.33, 0.33 ≦ F <0.5, 0.5 ≦ F <1, F = 1) Numerical values of 0.07, 0.15, 0.35, and 1 are determined in association with each. This short and long axis ratio coefficient Fk of the wire bundle is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the shortness of the electric wire bundle is such that the physical property value increases when the short / long axis ratio F of the electric wire bundle is large, and the physical property value decreases when the short / major axis ratio F of the electric wire bundle is small. By correcting the physical property value assigned to the element C by the long axis ratio coefficient Fk, the short major axis ratio F of the wire bundle is reflected in the physical property value. Here, the minor axis ratio F of the electric wire bundle is divided into four numerical ranges, but the numerical range is further divided finely, or the minor axis ratio of the electric wire bundle as a function of the minor axis ratio F of the electric wire bundle. It is also possible to express the coefficient Fk. A specific calculation formula for correcting the physical property value by the minor axis ratio coefficient Fk of the wire bundle will be described in detail later.

[G. Method of reflecting the distance between the center of the wire bundle and the center of the thick wire in the wire bundle in the cross section of the wire bundle]
In the cross section of the wire bundle, the rigidity of the wire bundle varies depending on the position of the thick wire in the wire bundle. The position of the thick electric wire in the cross section of the electric wire bundle can be evaluated by the following formula when the radius of the electric wire bundle is R and the radius of the thick electric wire in the electric wire bundle is r.
A = 1-r / R
R: radius of wire bundle r: radius of thick wire

  For example, A is 0 when the center of the wire bundle matches the center of the thick wire. On the other hand, when the thick electric wire is located on the outermost periphery of the electric wire bundle, A = 1−r / R. For example, 0.5A indicates that a thick electric wire in the wire bundle is arranged at an exactly middle point from the center of the wire bundle to the outermost periphery. Hereinafter, a numerical value mA (where m is 0 <m ≦ 1) that evaluates a location where a thick wire in the wire bundle is located is referred to as a thick material position evaluation value L.

  In the actual wire harness, the thick position evaluation value L can be calculated by measuring the wire harness to be modeled using a caliper, a measure, a weight scale, or the like.

  In the wire harness, when the thick position evaluation value L of the electric wire accommodated in a certain protection member P is different, the thick position evaluation value L is larger (in other words, the thick electric wire is positioned on the outer periphery of the electric wire bundle) Bending rigidity EI and torsional rigidity GJ are larger than those of a small one (in other words, a thick electric wire is located at the center of the electric wire bundle). In consideration of this point, a thick object position evaluation value L (L1, L2, L3,...) Is assigned to each element C (C1, C2, C3,...), And a thick object is assigned for each numerical range of the thick object position evaluation value L. A correspondence table that defines the position evaluation value coefficient Lk is created. Table 7 shows an example of the correspondence table.

  The thick object position evaluation value coefficient Lk has four numerical ranges of the thick object position evaluation value L (0, 0 <L ≦ 0.3A, 0.3A <L ≦ 0.6A, 0.6A <L ≦ A), respectively. The numerical values 1, 1.1, 1.2, and 1.3 are determined in association with. The thick object position evaluation value coefficient Lk is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the thick body position evaluation value coefficient Lk is such that the physical property value decreases when the thick object position evaluation value L is small, and the physical property value increases when the thick object position evaluation value coefficient Lk is large. By correcting the physical property value assigned to the element C by the above, the thick object position evaluation value L is reflected in the physical property value. Here, the thick object position evaluation value L is divided into four numerical ranges. However, the numerical value range is further divided, or the thick object position evaluation value coefficient Lk is expressed as a function of the thick object position evaluation value L. Is also possible. A specific calculation formula for correcting the physical property value by the thick position evaluation value coefficient Lk will be described in detail later.

[Correction taking into account the work environment when the wire harness is assembled to the vehicle]
An example of a working environment where the wire harness is routed in the vehicle is the outside air temperature. In the case where the outside air temperature is high, it is considered that the rigidity of the electric wire, the fixture, and the protective member constituting the wire harness is reduced due to softening of the resin or the like. Thus, depending on the working environment in which the wire harness is assembled to the vehicle (this environment is referred to as “the working environment of the wire harness” in this specification), the physical property values of the wire harness are different. The present inventor has found that the conventional analysis method does not consider the work environment of the wire harness, and therefore cannot accurately perform a study assuming a real wire harness. Less than,
(H) A method for calculating a physical property value in consideration of factors of the wire harness and the vehicle outside temperature and humidity when the wire harness is assembled to the vehicle will be described in detail.

[H. Method to reflect outside temperature and humidity]
When the outside air temperature is high or the humidity is high, the rigidity of the wires, fixtures, and protective members constituting the wire harness decreases. In other words, the bending rigidity EI and other various rigidity increase as the outside air temperature increases or the humidity increases. (Axial rigidity EA, torsional rigidity GJ, shear rigidity GA) is reduced. Based on this point, the outside air temperature K (K1, K2, K3,...) And the humidity H (H1, H2, H3,...) Are assigned to each element C (C1, C2, C3,...) A correspondence table that defines the work environment coefficient Ek for each numerical value range of the humidity H is created. Table 8 shows an example of the correspondence table.

  The working environment coefficient Ek includes three numerical ranges of outside air temperature K (0 ° C. ≦ K <12 ° C., 12 ° C. ≦ K <24 ° C., 24 ° C. ≦ K <36) and three numerical ranges of humidity H (40% ≦ H <55%, 55% ≦ H <70%, 70% ≦ H), the numerical value of the work environment coefficient Ek is determined in association with each. The work environment coefficient Ek is a numerical value for correcting the physical property value assigned to the element C. In the embodiment of the present invention, the physical property value decreases when the outside air temperature K is large, the physical property value increases when the outside air temperature K is small, and the physical property value decreases when the humidity H is large. On the other hand, the outside air temperature K and the humidity H are reflected in the physical property value by correcting the physical property value assigned to the element C by the work environment coefficient Ek so that the physical property value becomes large if the humidity H is small. Here, the outside air temperature K and the humidity H are divided into three numerical ranges, but the numerical range can be further divided, or the work environment coefficient Ek can be expressed as a function with respect to the outside air temperature K and the humidity H. is there. A specific calculation formula for correcting the physical property value by the work environment coefficient Ek will be described in detail later.

[Calculation formula for correcting physical property values]
Subsequently, the physical property values are determined by various coefficients described in the above items [Correction considering internal conditions of members constituting the wire harness] and Item [Correction considering the work environment when the wire harness is assembled to the vehicle]. An example of the calculation formula to be corrected will be described. Below, the bending rigidity EI is mentioned as a physical-property value of an electric wire bundle, and the case where this bending rigidity EI is correct | amended is demonstrated. The corrected bending rigidity EI ₁ is calculated by the following calculation formula.
EI ₁ = (1 + Ok + Rk + Tk + Dk) × Fk ₀₁ × Ek ₀₁ × EI ₀₁
EI ₁ : Bending stiffness after correction Ok: Electric wire occupancy factor Rk: Angle variation coefficient Tk: Tension coefficient Dk: Electric wire movement coefficient Fk ₀₁ : Short and long axis ratio coefficient Ek ₀₁ : Electric wire bundle working environment Coefficient EI ₀₁ : Bending rigidity of wire bundle (standard state)
The reference state is a state in which various coefficients shown in Tables 1 to 8 are the following values.
Electric wire occupancy factor Ok = 0
Angle variation coefficient Rk = 0
Tension coefficient Tk = 0
Protection member coefficient Vk = 0
Wire movement coefficient Dk = 0
Short / long axis ratio coefficient Fk = 1
Work environment coefficient Ek = 1
Thick position evaluation value coefficient Lk = 1
The bending rigidity EI ₀₁ of the wire bundle is a physical property value when the internal state of the wire harness and the work environment of the wire harness are in the standard state.

For example, for an element C, the wire occupancy O is 70%, the angle variation R is 2 °, the tension T is 3N, the wire movement distance D is 3 mm, the short / long axis ratio F of the wire bundle is 1, the outside air temperature K Is 20 ° C. and humidity H is 60%, the various coefficients specified with reference to the various correspondence tables are:
Electric wire occupancy factor Ok = 0.6
Angle variation coefficient Rk = −0.05
Tension coefficient Tk = 0.05
Wire movement coefficient Dk = −0.3
Short and long axis ratio coefficient Fk ₀₁ = 1 of the wire bundle
Work environment factor Ek ₀₁ = 1 for the wire bundle
become. Therefore, the corrected bending stiffness EI ₁ is
EI ₁ = (1 + 0.6−0.05 + 0.05−0.3) × 1 × 1 × EI ₀₁
= (1 + 0.3) × 1 × 1 × EI ₀₁
= 1.3EI ₀₁
Is calculated. In this calculation formula, the numerical value “0.3” multiplied by the bending rigidity EI ₀₁ of the wire bundle corresponds to the correction amount. Further, the factors that correct the bending rigidity EI ₀₁ of the wire bundle upward are the wire occupancy O and the tension T, and the degree of correction of the wire occupancy O upward is large. On the other hand, the factors for correcting the bending stiffness EI ₀₁ of the wire bundle downward are the angle variation coefficient Rk and the wire movement amount coefficient Dk.

Furthermore, the equation for correcting the bending rigidity of the thick wire and the protective member in consideration of the thick wire and the protective member is as follows.
EI ₁ = (1 + Ok + Rk + Tk + Dk) × Fk ₀₁ × Ek ₀₁ × EI ₀₁
+ Ek ₀₂ × Lk × EI ₀₂
+ (1 + Vk) × Ek ₀₃ × Fk ₀₃ × EI ₀₃
Ek _02: working environment factor thick wire Lk: Futomono location estimate coefficients _{EI 02:} thick wire bending stiffness Vk: protective member Factor _{Ek 03:} working environment coefficient of the protective member _{Fk 03:} short length axis ratio factor of the protection member EI ₀₃ : Bending rigidity of protective member

For example, in an element C, the wire occupancy O is 70%, the angle variation R is 2 °, the tension T is 3N, the wire movement amount D is 3 mm, the short and long axis ratio F of the wire bundle is 0.7, thick The object position evaluation value L is 0.2 A, the minor axis ratio F of the protective member is 0.7, the identifier V of the protective member is “sheet”, the outside air temperature K is 20 ° C., and the humidity H is 60%. Then, various coefficients specified by referring to various correspondence tables are
Electric wire occupancy factor Ok = 0.6
Angle variation coefficient Rk = −0.05
Tension coefficient Tk = 0.05
Wire movement coefficient Dk = −0.3
Short length axis ratio factor of the wire bundle _Fk 01 = 0.35
Work environment factor Ek ₀₁ = 1 for the wire bundle
Work environment factor Ek ₀₂ = 1 for thick wires
Thick position evaluation value coefficient Lk = 1.1
Protection member coefficient Vk = 0.1
Work environment coefficient Ek ₀₃ = 1 of protective member
Protective member short / long axis ratio coefficient Fk ₀₃ = 0.35
become. Therefore, the corrected bending stiffness EI ₁ is
EI ₁ = (1 + 0.6−0.05 + 0.05−0.3) × 0.35 × 1 × EI ₀₁
+ 1 × 1.1 × EI ₀₂
+ (1 + 0.1) × 1 × 0.35 × EI ₀₃
= (1 + 0.3) × 0.35 × 1 × EI ₀₁ + 1.1 × EI ₀₂ + 0.385 × EI ₀₃
= 0.455EI ₀₁ + 1.1 × EI ₀₂ + 0.385 × EI ₀₃
Is calculated. In this calculation formula, the bending rigidity EI ₀₂ of the thick electric wire is obtained by the numerical value “1.1” multiplied by the bending rigidity EI _{02 of the} thick electric wire and the numerical value “0.385” multiplied by the bending rigidity EI ₀₃ of the protective member. And the bending rigidity EI _{03 of the} protective member is corrected. The numerical value obtained by adding the first term relating to the bending stiffness EI ₀₁ of the wire bundle, the second term relating to the bending stiffness EI ₀₂ of the thick wire, and the third term relating to the bending stiffness EI ₀₃ of the protective member is the corrected bending stiffness EI. Calculated as ₁ .

  As described above, the physical property value assigned to the element C is corrected. Here, the form for correcting the bending rigidity EI has been described, but the same applies to correction of other various rigidity (axial rigidity EA, torsional rigidity GJ, shear rigidity GA).

Note that the calculation formula described in this item [Calculation formula for correcting physical property values] is an example, and various calculation formulas for correcting physical property values are conceivable. Of course, by appropriately combining the various factors described in the item [Correction taking into account the internal conditions of members constituting the wire harness] and Item [Correction taking into account the work environment when the wire harness is assembled to the vehicle] A calculation formula for correcting the value can be obtained. Various calculation formulas for the purpose of reflecting one or both of the internal state of the wire harness and the work environment of the wire harness in the physical property values can be applied to the present invention. This calculation formula can be generalized by the following function.
M1 = M0 (Zi, Zo)
M0: physical property value before correction M1: physical property value after correction Zi: variable obtained by quantifying the internal state of the wire harness Zo: variable obtained by quantifying the work environment of the wire harness

[Hardware Configuration of Wire Harness Analyzer of Embodiment of the Present Invention]
FIG. 6 is a block diagram illustrating a hardware configuration of the wire harness analyzer according to the embodiment of the present invention. The wire harness analysis device according to the embodiment of the present invention includes an input unit 611, a database unit 612, a program storage unit 613, a data storage unit 614, a display unit 615, and a processing unit 616. When the wire harness analyzer of the present invention is configured by a general-purpose PC, for example, the input unit 611 is realized by various input interfaces such as a keyboard, a mouse, and a numeric keypad, and the database unit 612 and the program storage unit 613 are a hard disk drive (HDD). The data storage unit 614 is realized by a RAM (Random Access Memory), the display unit 615 is realized by various output devices such as a CRT display and a liquid crystal display, and the processing unit 616 is executed by a CPU (Central Processing Unit). Realized. The database unit 612 stores information (including physical property values) for modeling the wire harness, and the digitized internal state and work environment of the wire harness. The program storage unit 613 stores a program for causing the processing unit 616 to execute the calculation algorithm described in the above [Wire harness analysis algorithm] and [Physical property value calculation method]. The data storage unit 614 records data input and output from the processing unit 616 executing the calculation algorithm described in [Wire harness analysis algorithm] and [Physical property value calculation method].

[Analysis method by the wire harness analyzer of the embodiment of the present invention]
Then, the analysis method by the wire harness analyzer of embodiment of this invention is demonstrated with reference to FIG. FIG. 7 is a flowchart illustrating a procedure of an analysis method performed by the wire harness analyzer according to the embodiment of the present invention.

  First, the processing unit 616 accesses the database unit 612 to refer to the finite element model of the wire harness divided into an arbitrary number of elements (step S71). At this time, the physical property value before correction assigned to each element is read. Subsequently, the processing unit 616 performs a calculation for correcting the physical property value as described in the above [Method for calculating physical property value] (step S72), and writes the corrected physical property value in the data storage unit 614. Then, the processing unit 616 refers to the corrected physical property value stored in the data storage unit 614, and calculates the shape of the wire harness as described in the above [Wire harness analysis algorithm] (step S73).

[Summary]
As described above, according to the wire harness analysis method, the wire harness analysis apparatus, and the program according to the embodiment of the present invention, while adopting the approximate model that captures the wire harness as one continuous body, It is possible to consider the external environment when the wire harness is assembled to the vehicle. For this reason, since an analysis closer to the actual wire harness can be performed, it is possible to examine the work by the worker according to the work environment in consideration of the analysis result.
In addition, the manufacturing requirements of the wire harness are managed using the analysis result, for example, the tension of the tape when the tape is wound around a predetermined portion of the wire harness is controlled to a predetermined value, so that the operator feels uncomfortable. In addition, it is possible to manufacture a wire harness having stable characteristics with no variation in the accuracy of work performed by an operator.

Here, the features of the embodiments of the wire harness analysis method, the wire harness analysis device, and the program according to the present invention described above are briefly summarized and listed in the following [1] to [12], respectively.
[1] A wire harness analysis method for calculating a shape of a wire harness that is sequentially deformed using a computer,
A reference step (S71) for referring to a physical property value for each small region of the discrete wire harness;
A correction step of correcting the physical property value based on one or both of an internal situation of a member constituting the wire harness to be analyzed and a work environment when the wire harness to be analyzed is assembled to a vehicle (S72),
A calculation step (S73) for calculating the shape of the wire harness based on the physical property value corrected in the correction step;
A wire harness analysis method comprising:
[2] In the correction step, for each element of the wire harness discretized by an element, the physical property value is corrected based on one or both of the internal situation and the work environment. 1] The wire harness analysis method according to [1].
[3] The internal state includes a ratio of an electric wire occupying on a cross section of the wire harness to a cross sectional area of the wire harness.
The wire harness analysis method according to [2], wherein
[4] The internal situation includes variations in the angle of the wire in the axial direction with respect to the axial direction of the wire harness at an arbitrary position in the longitudinal direction of the wire harness.
The wire harness analysis method according to [2] or [3], wherein
[5] The internal state includes a tension acting on the band-shaped member when a portion of the wire harness is wound by the band-shaped member, or a tightening force exerted on the portion by the band-shaped member.
The wire harness analysis method according to any one of [2] to [4], wherein:
[6] The internal situation includes the rigidity of a protective member that covers a portion of the wire harness.
The wire harness analysis method according to any one of [2] to [5], wherein
[7] The internal situation includes the amount of movement of the electric wire allowed within the electric wire bundle.
The wire harness analysis method according to any one of [2] to [6], wherein
[8] The internal situation includes the ratio of the short axis to the long axis of the cross section of the wire bundle at an arbitrary position in the longitudinal direction of the wire harness.
The wire harness analysis method according to any one of [2] to [7], wherein
[9] The internal situation includes a numerical value for evaluating a location where a thick wire in the wire bundle is located.
The wire harness analysis method according to any one of [2] to [8], wherein:
[10] The external environment includes outside air temperature or humidity under a situation where the wire harness is assembled to a vehicle.
The wire harness analysis method according to any one of [2] to [9], wherein:
[11] A wire harness analyzer for calculating the shape of a wire harness that is sequentially deformed,
A physical property value for each small region of the wire harness that has been discretized is recorded, and an internal situation of a member constituting the wire harness to be analyzed and the wire harness to be analyzed are assembled to a vehicle. A recording unit (database unit 612) in which one or both of the work environments are recorded;
For each small region of the wire harness that has been discretized, either the internal state of the members constituting the wire harness to be analyzed and the work environment when the wire harness to be analyzed is assembled to a vehicle Or a calculation unit (processing unit 616) that corrects the physical property value based on both, and calculates the shape of the wire harness based on the corrected physical property value;
A wire harness analyzing apparatus comprising:
[12] A program for causing a computer to execute each step of the wire harness analysis method according to any one of [1] to [10].

611 Input unit 612 Database unit 613 Program storage unit 614 Data storage unit 615 Display unit 616 Processing unit

Claims (12)

A wire harness analysis method for calculating a shape of a wire harness that is sequentially deformed using a computer,
A step of referring to a physical property value for each small region of the wire harness that has been discretized;
A correction step of correcting the physical property value based on one or both of an internal situation of a member constituting the wire harness to be analyzed and a work environment when the wire harness to be analyzed is assembled to a vehicle When,
Based on the physical property values corrected in the correction step, a calculation step for calculating the shape of the wire harness;
A wire harness analysis method comprising: The physical property value is corrected based on one or both of the internal situation and the work environment for each element of the wire harness discretized by the element in the correction step. The wire harness analysis method as described. The internal state includes a ratio of the electric wire occupying on a cross section of the wire harness to the cross sectional area of the wire harness.
The wire harness analysis method according to claim 2. The internal situation includes a variation in the angle in the axial direction of the electric wire with respect to the axial direction of the wire harness at an arbitrary position in the longitudinal direction of the wire harness.
The wire harness analysis method according to claim 2 or 3, wherein The internal situation includes a tension acting on the band-shaped member when a portion of the wire harness is wound by the band-shaped member, or a tightening force exerted on the position by the band-shaped member.
The wire harness analysis method according to any one of claims 2 to 4, wherein the wire harness analysis method is provided. The internal situation includes the rigidity of a protective member that covers a portion of the wire harness.
The wire harness analysis method according to any one of claims 2 to 5, characterized in that: The internal situation includes the amount of wire movement allowed within the wire bundle,
The wire harness analysis method according to any one of claims 2 to 6, characterized in that: The internal situation includes the ratio of the short axis to the long axis of the cross section of the wire bundle at an arbitrary position in the longitudinal direction of the wire harness.
The wire harness analysis method according to any one of claims 2 to 7, wherein the wire harness analysis method is provided. The internal situation includes a numerical value for evaluating a location where a thick wire in the wire bundle is located,
The wire harness analysis method according to any one of claims 2 to 8, wherein the wire harness analysis method is provided. The working environment includes outside air temperature or humidity under the condition that the wire harness is assembled to a vehicle.
The wire harness analysis method according to any one of claims 2 to 9, wherein the wire harness analysis method is provided. A wire harness analyzer for calculating the shape of a wire harness that is sequentially deformed,
A physical property value for each small region of the wire harness that has been discretized is recorded, and an internal situation of a member constituting the wire harness to be analyzed and the wire harness to be analyzed are assembled to a vehicle. A recording section in which one or both of the work environments are recorded;
For each small region of the wire harness that has been discretized, either the internal state of the members constituting the wire harness to be analyzed and the work environment when the wire harness to be analyzed is assembled to a vehicle Or a calculation unit that corrects a physical property value based on both, and calculates the shape of the wire harness based on the corrected physical property value;
A wire harness analyzing apparatus comprising: The program for making a computer perform each step of the wire harness analysis method of any one of Claim 1 to 10.

Images