WO2014129125A1 - 信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム - Google Patents

信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム Download PDF

Info

Publication number
WO2014129125A1
WO2014129125A1 PCT/JP2014/000537 JP2014000537W WO2014129125A1 WO 2014129125 A1 WO2014129125 A1 WO 2014129125A1 JP 2014000537 W JP2014000537 W JP 2014000537W WO 2014129125 A1 WO2014129125 A1 WO 2014129125A1
Authority
WO
WIPO (PCT)
Prior art keywords
target device
stress
input
failure rate
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2014/000537
Other languages
English (en)
French (fr)
Inventor
耕太 大石
邦保 松本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyota Motor Corp filed Critical Toyota Motor Corp
Priority to CN201480009173.8A priority Critical patent/CN104995624B/zh
Priority to DE112014000935.0T priority patent/DE112014000935T5/de
Priority to US14/768,937 priority patent/US11085843B2/en
Priority to KR1020157022249A priority patent/KR101658748B1/ko
Publication of WO2014129125A1 publication Critical patent/WO2014129125A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/22—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring the force applied to control members, e.g. control members of vehicles, triggers
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Definitions

  • the present invention relates to a reliability design support device, a reliability design support method, and a reliability design support program.
  • a general risk assessment is performed by the following process, for example.
  • step (2) “(3) Hazard identification, risk level estimation given by hazard” is executed. This includes the process of identifying hazards caused by failures and behaviors caused by product use and misuse (3-1: identification of hazards), and the process of estimating the degree of risk given by hazards (3-2: risk level Estimate).
  • Patent Document 1 discloses a reliability design support device that calculates the failure rate of home appliances and the like installed and used indoors.
  • the reliability design support apparatus clarifies a reference value for occurrence of a failure, and calculates a failure rate by comparing the probability density distribution of the characteristic values of parts with the reference value.
  • failure rates are calculated using a stress-strength model.
  • the input stress in the market is defined as a distribution by measurement or the like, and the failure rate is calculated by collating the distribution with the probability density distribution of the component characteristic values.
  • the reliability design support apparatus described in Patent Literature 1 calculates a failure rate by clarifying a reference value for occurrence of a failure.
  • a reference value for occurrence of a defect cannot be determined, and thus the reliability design support apparatus cannot be used.
  • the input stress distribution is used as described above.
  • the input stress distribution is determined by factors such as the user's physical characteristics and usage (road surface conditions when using the product) in addition to design factors such as the shape and configuration of the parts constituting the product.
  • the general stress-strength model does not consider these factors.
  • the failure rate can be calculated with high accuracy by the stress-strength model.
  • the conventional technology has a problem that the product failure rate cannot be calculated accurately when a factor that changes according to the user and the environment has a large correlation with the input stress.
  • the present invention has been made in view of the above problems, and is a reliability design support that accurately calculates a product failure rate even when a factor that changes according to the user or the environment has a large correlation with input stress.
  • the main purpose is to provide a device, a reliability design support method, and a program.
  • One aspect of a reliability design support apparatus is an approximation to a data string in which a value of a target factor that changes according to a user and an environment and an input stress measurement value measured when the target apparatus is used are associated with each other.
  • Relational expression calculating means for calculating a relational expression by performing calculation, input stress calculating means for calculating an input stress input to the target device according to the relational expression and the index data of the target factor, and Failure rate calculation means for calculating a failure rate of the target device according to the input stress and the characteristic value distribution of the target device.
  • the index data is a probability density distribution related to the physical characteristics of the user
  • the input stress calculation means relates to the input stress based on the relational expression and the probability density distribution. It is preferable that a stress distribution as a distribution is calculated, and the failure rate calculation means calculates a failure rate of the target device from an overlap between the stress distribution and the characteristic value distribution.
  • the input stress calculation unit calculates a threshold line by substituting the worst value of the target factor into the relational expression
  • the failure rate calculation unit includes the characteristic value distribution and the threshold line. It is desirable to calculate the failure rate of the target device from the distribution region formed by
  • the input stress distribution adjusting means for adjusting the input stress calculated by the input stress calculating means by performing finite element analysis using the part shape data used for the target device is further provided. It is desirable to have.
  • the target factor is a physical feature of the user, a usage environment of the target device, an actual measurement value when the target device is used, or a characteristic value of the target device, and the index data is ,
  • the target factor is a physical feature of the user, the probability density distribution on the physical feature of the user, and when the target factor is a measured value when the target device is used, A plurality of values are actually measured data, and are the worst values that can be calculated from the properties of the target device when the target factor is the characteristic value of the target device.
  • One aspect of the reliability design support method is an approximation to a data string in which a value of a target factor that changes according to a user or an environment and an input stress measurement value measured when the target device is used are associated with each other.
  • a relational expression calculating step for calculating a relational expression by performing a calculation; an input stress calculating step for calculating an input stress input to the target device according to the relational expression and the index data of the target factor; and
  • a failure rate calculation step of calculating a failure rate of the target device according to the input stress and the characteristic value distribution of the target device is executed.
  • One aspect of the program according to the present invention is an approximate calculation for a data string in which a value of a target factor that changes according to a user or an environment and an input stress measurement value measured when the target device is used are associated with a computer.
  • a relational expression calculating step for calculating a relational expression, an input stress calculating step for calculating an input stress input to the target device according to the relational expression and the index data of the target factor, and the input A failure rate calculating step of calculating a failure rate of the target device according to the stress and the characteristic value distribution of the target device is executed.
  • a reliability design support apparatus a reliability design support method, and a program for accurately calculating a product failure rate even when a target factor that changes according to a user or an environment has a large correlation with input stress. Can be provided.
  • FIG. 1 is a block diagram showing a configuration of a reliability design support apparatus 1 according to a first exemplary embodiment.
  • FIG. 3 is a conceptual diagram showing an operation of the reliability design support apparatus 1 according to the first exemplary embodiment.
  • FIG. 3 is a conceptual diagram showing an operation of the reliability design support apparatus 1 according to the first exemplary embodiment.
  • FIG. 3 is a conceptual diagram showing an operation of the reliability design support apparatus 1 according to the first exemplary embodiment.
  • 3 is a flowchart showing a risk assessment method using the reliability design support device 1 according to the first exemplary embodiment;
  • 1 is a block diagram showing a hardware configuration of a reliability design support apparatus 1 according to a first exemplary embodiment.
  • FIG. 3 is a block diagram showing a configuration of a reliability design support apparatus 1 according to a second exemplary embodiment.
  • FIG. 1 is a block diagram showing the configuration of the reliability design support apparatus according to the present embodiment.
  • the reliability design support apparatus 1 is an apparatus that performs each process of risk assessment (the above (1) to (4)) at the time of product development.
  • the configuration and operation relating to “(2-1) Estimating failure rate” in the above-described process will be mainly described.
  • a device that is a target for which the reliability design support device 1 performs design support may be an arbitrary product.
  • the target device is an electrically assisted bicycle that is in close contact with the human body, a small mobility, a wearable robot, or the like.
  • These target devices are mainly devices that come into contact with the human body, and factors such as weight (for example, weight) and environment (such as road surface conditions in which the target device is used) that vary from user to user are dominant with respect to input stress. It is a device (very influential). For example, when the user's weight increases or decreases, the stress applied to the product changes greatly.
  • the reliability design support device 1 includes a relational expression calculating means 10, an input stress calculating means 20, a characteristic value distribution calculating means 30, and a failure rate calculating means 40. Details of each means will be described below.
  • the relational expression calculation means 10 reads the stress measurement value file.
  • the stress measurement value file is a file in which a data string in which a value of a certain factor (also referred to as a target factor in the following description) and a stress measurement value corresponding to the value of the target factor are associated is described. It is.
  • the target factor is a factor that can change depending on the user and the environment
  • the stress measurement value is a value measured when the target device is used.
  • the stress measurement file describes a data string in which the weight of each user's weight is associated with the stress measurement value when each user uses the target device.
  • the target factors described in the stress measurement file include, in addition to weight, physical characteristics such as height and shoe size, and the ratio between height and handle height of small mobile mobility.
  • the relational expression calculating means 10 calculates a relational expression indicating the relation between the value of the target factor and the magnitude of the stress by analyzing the data string of the stress measurement value file using a general approximation method.
  • the relational expression calculating means 10 calculates the relational expression using, for example, a response surface method as an approximation method.
  • the relational expression calculating means 10 may use multiple regression analysis or the like as an approximation method.
  • the relational expression calculation unit 10 supplies the calculated relational expression to the input stress calculation unit 20.
  • the input stress calculation means 20 reads the index data file in which the index data of the target factor is described.
  • the index data file describes the probability density distribution of the target factor.
  • the index data file describes graph information (or similar information) in which the horizontal axis is the weight and the vertical axis is the distribution ratio.
  • it is desirable to use known data with high reliability such as human body statistical data for the index data file.
  • the input stress calculation means 20 calculates an input stress distribution based on the relational expression calculated by the relational expression calculation means 10 and the index data file.
  • the worst value that can be logically calculated in addition to the probability density distribution as described above, it is possible to describe actually measured statistical data and the worst value that can be logically calculated.
  • measured statistical data for example, the distribution of the positions of both feet of the user when using small mobile mobility can be used.
  • the worst values that can be calculated logically include, for example, the maximum torque and maximum speed of the motor determined by the design of the target site, and the chemical substance concentration (value standardized by IEC (International Electrotechnical Commission) standards). To do.
  • the input stress distribution calculated by the input stress calculation means 20 is information indicating the stress magnitude on the horizontal axis and the occurrence frequency on the vertical axis.
  • the input stress calculation means 20 calculates the input stress distribution by performing analysis using, for example, the Monte Carlo method on the relational expression calculated by the relational expression calculation means 10 and the index data file.
  • the input stress calculation unit 20 supplies the calculated input stress distribution to the failure rate calculation unit 40. Note that the input stress calculated by the input stress calculation means 20 is not necessarily in the form of distribution, but may be a linear expression as shown in FIG.
  • the characteristic value distribution calculating unit 30 reads the characteristic value file.
  • the characteristic value file is a file in which characteristic values such as the strength of parts used in the target device are described.
  • the characteristic value file a plurality of corresponding parts of the target device are sampled, and measurement values obtained by measuring the magnitude of stress of each part are described. The stress is measured by, for example, a general tensile test.
  • the characteristic value distribution calculating unit 30 calculates a characteristic value distribution from the characteristic value file with the horizontal axis as stress and the vertical axis as frequency of occurrence.
  • the characteristic value distribution calculating unit 30 supplies the calculated characteristic value distribution to the failure rate calculating unit 40.
  • each means reads out the file format information, but is not necessarily limited to this, and is stored in, for example, a database. Each data may be read out.
  • the characteristic value distribution calculating unit 30 calculates the characteristic value distribution of the target device from the characteristic value file.
  • the present invention is not limited to this.
  • an input file describing a characteristic value distribution may be directly input to the reliability design support apparatus 1.
  • the failure rate calculation means 40 calculates the failure rate of the target device according to the overlap between the input stress distribution and the characteristic value distribution. Details of the failure rate calculation method will be described later with reference to FIG.
  • the failure rate calculation means 40 outputs the calculated failure rate to an arbitrary output means (for example, a display device) or a storage device (for example, a file system).
  • FIG. 2 is a conceptual diagram illustrating a first operation example of the reliability design support apparatus 1.
  • the target factor is the weight of the user.
  • the weight of the user is a different factor for each user.
  • the relational expression calculation means 10 reads all data strings in the stress measurement value file.
  • the data string for example, a user with body weight Xkg and a stress measurement value when a user with body weight Xkg uses the target device are described in association with each other.
  • the relational expression calculating means 10 calculates a relational expression (FIG. 2A) by using all data strings in the stress measurement value file as analysis targets of an approximation method such as a response surface method.
  • the input stress calculation means 20 reads an index data file indicating statistical information of the weight distribution as shown in FIG.
  • the index data file may be a file described as appropriate according to the application of the target device. For example, when the target device is used only by Japanese people, the probability density distribution of the weights of Japanese people may be used instead of the probability density distribution of the weights of all human beings.
  • the input stress calculation means 20 performs a known analysis (for example, analysis by the Monte Carlo method) using a relational expression (FIG. 2 (A)) and statistical information on the weight distribution (FIG. 2 (B)), thereby providing an input stress distribution ( FIG. 2C is calculated.
  • the input stress distribution is a distribution indicating the relationship between the magnitude of stress and the occurrence frequency as described above.
  • the characteristic value distribution calculating means 30 analyzes the strength of the parts used in the target product and calculates the characteristic value distribution (for example, intensity distribution) of the parts (FIG. 2D).
  • the component characteristic distribution is a distribution representing the relationship between the magnitude of stress and the frequency of failure as shown in FIG.
  • the failure rate calculation means 40 collates the input stress distribution (FIG. 2 (C), FIG. 2 (D)) with the characteristic value distribution (FIG. 2 (D)), and calculates the failure rate. Specifically, the failure rate calculation unit 40 calculates a ratio between the overlapping portion of both distributions and the entire characteristic value distribution of parts as a failure rate. That is, the failure rate calculation unit 40 calculates the failure rate by performing calculation using the total number of components as the denominator and the number of components related to the overlapping portion as the numerator.
  • FIG. 3 is a conceptual diagram illustrating a second operation example of the reliability design support apparatus 1.
  • the target factor is the temperature of the user's living environment where the target device is used.
  • the temperature of the living environment is also a factor that changes according to the use state of the user.
  • the relational expression calculation means 10 reads all data strings in the stress measurement value file.
  • the data string for example, a temperature and a stress measurement value when the target device is used at the temperature are described in association with each other.
  • the relational expression calculation means 10 calculates a relational expression (FIG. 3A) by using all data strings in the stress measurement value file as analysis targets of an approximation method such as a response surface method.
  • the input stress calculation means 20 reads the index data file. For example, in the index data file, a plurality of temperature data measured at the usage location of the target device are described. The input stress calculation means 20 extracts the maximum value (the maximum temperature, indicating the worst value) from the temperature data (FIG. 3B). The input stress calculation means 20 calculates a threshold straight line (FIG. 3C) indicating the input stress maximum value by substituting the extracted maximum value into the relational expression (FIG. 3A).
  • the failure rate calculation means 40 calculates a failure rate by collating the threshold line (FIGS. 3C and 3D) indicating the maximum input stress value and the characteristic value distribution (FIG. 3D). The failure rate calculation means 40 calculates the failure rate based on a range where the stress is smaller than the threshold line as shown in the figure.
  • FIG. 4 is a conceptual diagram illustrating a third operation example of the reliability design support apparatus 1.
  • the target device is small mobile mobility.
  • the target factor is a velocity distribution when the target device is used. Since the operation example is different from the first operation example described with reference to FIG. 2 only in the index data file, the contents of the index data file will be described below.
  • the input stress calculation means 20 performs a known analysis (for example, analysis by the Monte Carlo method) using a relational expression (FIG. 4A) and a velocity distribution (FIG. 4B) to thereby input stress distribution (FIG. C)) is calculated.
  • a method for risk assessment using the reliability design support apparatus 1 will be described with reference to FIG.
  • a designer assumes a user of a target device (or a device that may come into contact with the target device) that is a target of risk assessment (S1). Then, the designer clarifies the usage method and foreseeable misuse by the assumed user (or the person who may contact the target device) (S1).
  • the designer identifies the hazard that is caused by the behavior caused by the use and misuse of the target device (S3-1). Then, the designer estimates the degree of risk based on the influence of the hazard on the user (S3-2). Subsequently, the designer defines an allowable risk occurrence frequency, that is, an allowable failure rate, based on the estimated risk level, the planned number of target devices to be shipped, and the like (S4). For example, the designer defines the allowable failure rate so that there is a probability that none of the target devices will have a significant effect on the user.
  • the designer sets the reliability design support apparatus 1 described above (S2-1-a). Then, the designer prepares various files (index data file, stress measurement value file, characteristic value file) to be input to the reliability design support apparatus 1.
  • the reliability design support device 1 calculates a failure rate by the above-described method, and compares the calculated failure rate with an allowable failure rate. The designer appropriately adjusts the design of the target device until the calculated failure rate becomes smaller than the allowable failure rate.
  • the reliability design support apparatus 1 for handling this allowable failure rate may include failure rate evaluation means (not shown) in the subsequent stage of the failure rate calculation means 40 in addition to the configuration of FIG.
  • the failure rate evaluation unit compares the failure rate calculated by the failure rate calculation unit 40 with the allowable failure rate input by the designer.
  • the reliability design support apparatus 1 includes an input device 110, a CPU (Central Processing Unit) 120, a ROM (Read Only Memory) 130, a RAM (Random Access Memory) 140, and an output device 150 as hardware.
  • the configuration illustrated in FIG. 6 is an example, and other configurations may be used.
  • a storage device such as a USB (Universal Serial Bus) memory may be configured to be detachable.
  • the input device 110 is an information input device such as a keyboard and a mouse.
  • the ROM 130 is a storage device used for storage such as BIOS (Basic Input / Output System) and IPL (Initial Program Loader).
  • BIOS Basic Input / Output System
  • IPL Initial Program Loader
  • the RAM 140 is a memory that holds a program and is used as a working area of the CPU 120 that executes the program.
  • the CPU 120 reads and executes a program and controls output to the output device 150.
  • the output device 150 is a device that outputs arithmetic processing results and the like under the control of the CPU 120.
  • the output device 150 is, for example, a liquid crystal display device or a printer device.
  • the reliability design support device 1 uses the input stress (for example, FIG. 2C) and FIG. 3 (C)).
  • the reliability design support device 1 calculates the failure rate of the target device using this input stress. That is, the reliability design support apparatus 1 calculates the input stress considering the change of the factor value. Therefore, the reliability design support apparatus 1 can calculate an appropriate failure rate even when a change factor that changes depending on the user or the environment is dominant with respect to the input stress.
  • the reliability design support device 1 targets a human body support device (for example, an electrically assisted bicycle, a small mobile mobility, a wearable robot, etc.) in which physical characteristics such as weight and height have a great influence on input stress. Even in this case, the failure rate can be accurately calculated. This makes it possible to design a target device that achieves high reliability and low cost.
  • the reliability design support device 1 can calculate the failure rate in consideration of the user's physical characteristics (weight, height, etc.) as the target factors as shown in FIG. In this case, the reliability design support apparatus 1 can calculate an input stress distribution based on a statistical probability density distribution having a large parameter, that is, a highly reliable distribution. Therefore, the reliability design support apparatus 1 can calculate a highly reliable input stress distribution. Therefore, the reliability design support device 1 can accurately calculate the failure rate.
  • the reliability design support device 1 can also calculate the failure rate based on the maximum value (worst value) of the target factor as shown in FIG. Thereby, the reliability design support apparatus 1 can calculate the failure rate in consideration of the case where the target apparatus is used in the worst environment (for example, when used at the highest temperature). By evaluating this failure rate, safer product design can be realized.
  • the maximum value of the target factor is used.
  • the present invention is not necessarily limited thereto. For example, a failure rate may be calculated in consideration of the minimum value according to the nature of the factor.
  • the reliability design support apparatus 1 according to the present embodiment is characterized in that a failure rate is calculated in consideration of the shape characteristics of components used in the target apparatus.
  • the reliability design support apparatus 1 according to the present embodiment will be described below with respect to differences from the first embodiment. In the following description, the same processing units as those in the first embodiment are given the same names and the same reference numerals, and detailed descriptions thereof are omitted.
  • FIG. 7 is a block diagram showing a configuration of the reliability design support apparatus 1 according to the present embodiment.
  • the reliability design support apparatus 1 according to the present embodiment is configured to further include an input stress distribution adjusting unit 50 in addition to the configuration of FIG.
  • the input stress calculation means 20 calculates the stress distribution by the above method. At this time, the input stress calculation means 20 calculates the stress on the horizontal axis as a value in Newton units, calculates the input stress distribution with the vertical axis as the occurrence frequency, and inputs the input stress distribution to the input stress adjustment means 50. Supply.
  • a part shape file is input to the input stress adjusting means 50.
  • the part shape file is a file in which data indicating the shape of a part constituting the target device is stored.
  • the part shape file stores data (part shape data) in a format that can be used for finite element analysis described later.
  • the input stress adjustment means 50 performs finite element analysis based on the part shape file and the input stress distribution, and calculates the adjusted input stress distribution.
  • This input stress distribution is a distribution in which the stress on the horizontal axis is calculated using megapascal (Mpa) units, and the vertical axis is the frequency of occurrence.
  • the input stress adjustment means 50 supplies the failure rate calculation means 40 with the input stress distribution (stress distribution with the stress on the horizontal axis as the unit of megapascal (Mpa)) adjusted using the component shape file.
  • the failure rate calculation means 40 may calculate the failure rate by the same method as in the first embodiment.
  • the input stress distribution is calculated in consideration of the part shape used for the target device.
  • the reliability design support apparatus 1 concerning this Embodiment can implement
  • Each means (relational expression calculating means 10, input stress calculating means 20, characteristic value distribution calculating means 30, failure rate calculating means 40, input stress of the reliability design support apparatus 1 according to the first embodiment and the second embodiment described above.
  • the adjusting means 50) is executed by the CPU 120 in the form of a program.
  • the program can be stored and provided to a computer using various types of non-transitory computer readable media.
  • Non-transitory computer readable media include various types of tangible storage media (tangible storage medium).
  • non-transitory computer-readable media examples include magnetic recording media (eg flexible disks, magnetic tapes, hard disk drives), magneto-optical recording media (eg magneto-optical discs), CD-ROMs (Read Only Memory), CD-Rs, CD-R / W, semiconductor memory (for example, mask ROM, PROM (Programmable ROM), EPROM (Erasable ROM), flash ROM, RAM (random access memory)) are included.
  • the program may also be supplied to the computer by various types of temporary computer-readable media. Examples of transitory computer readable media include electrical signals, optical signals, and electromagnetic waves.
  • the temporary computer-readable medium can supply the program to the computer via a wired communication path such as an electric wire and an optical fiber, or a wireless communication path.

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)

Abstract

　関係式算出手段（１０）は、ユーザ又は環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する。入力ストレス算出手段（２０）は、関係式と、対象因子の指標データと、に応じて対象装置に入力される入力ストレスを算出する。故障率算出手段（４０）は、入力ストレスと、対象装置の特性値分布と、に応じて対象装置の故障率を算出する。

Description

信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム
　本発明は信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラムに関する。
　近年、製造業者に対する安全性要求が高まってきている。これに伴い、製品開発時におけるリスクアセスメントの実施が必須となってきている。一般的なリスクアセスメントは、例えば以下のようなプロセスで実行される。
　はじめに、「（１）製品の使用、誤使用の明確化」を実行する。これは、製品の使用するユーザや意図される使用シーンを特定し、予見可能な誤使用を明確化する工程である。続いて、「（２）故障率の見積もり、及び発生頻度の見積もり」を実行する。これは、製品の使用及び誤使用が製品を構成する部品に与える入力ストレス（負荷）を見積り、故障率を見積もる工程（２－１：故障率の見積もり）と、誤使用頻度や部品故障確率を考慮してリスクの発生頻度を見積もる工程（２－２：発生頻度の見積もり）と、から構成される。
　また、上記（２）の工程と並行して、「（３）ハザードの特定、ハザードが与えるリスクの程度の見積もり」を実行する。これは、製品使用及び誤使用により引き起こされる故障や挙動が引き起こすハザードを特定する工程（３－１：ハザードの特定）と、ハザードが与えるリスクの程度を見積もる工程（３－２：リスクの程度の見積もり）と、から構成される。
　そして、上記（２）及び（３）の工程後に「（４）リスクの評価」を実行する。これは、上記工程により算出したリスクの発生頻度、及びリスクの程度、の２軸で評価を行い、リスクの評価を行う工程である。上記した（１）～（４）の工程の内、故障率を見積もる工程（２－１）において高精度な故障率の見積もりを行うことが非常に重要である。
　特許文献１には、屋内に設置して利用される家電製品等の故障率を算出する信頼性設計支援装置が開示されている。当該信頼性設計支援装置は、不良発生の基準値を明確化し、部品の特性値の確率密度分布と当該基準値を照合することにより故障率を算出する。
　また、屋外環境で使用される製品の市場では、ストレス－ストレングスモデルを用いた故障率の算出が行われている。ストレス－ストレングスモデルを用いた手法では、市場における入力ストレスを計測等により分布として定義し、当該分布と部品特性値の確率密度分布を照合して故障率を算出する。
特開２００２－１４９７１４号公報
田村　泰彦、飯塚　悦功著、不具合に関する設計知識の運用に関する研究 : ストレス-ストレングスモデルによる知識獲得、Journal of the Japanese Society for Quality Control 31(1), P168-180
　しかしながら、上述の特許文献１の技術を用いた場合には以下の問題がある。特許文献１に記載の信頼性設計支援装置は、不良発生の基準値を明確化することにより故障率を算出している。しかし、入力ストレスのばらつきが大きい場合には不良発生の基準値を定めることができないため、当該信頼性設計支援装置を使用することができない。
　また上述の一般的なストレス－ストレングスモデルでは、上述のように入力ストレス分布を用いる。入力ストレス分布は、製品を構成する部品形状や構成等の設計因子に加え、ユーザの身体的特徴等や使われ方（製品を使用する際の路面状況等）の因子によって決定される。一般的なストレス－ストレングスモデルではこれらの因子に関する考慮が行われていない。
　たとえば自動車のように大きなエネルギーを生み出す商品の場合、ユーザの体重等の因子や使われ方の因子は入力ストレスにほとんど影響しない。そのため、自動車等の製品設計時には、ストレス－ストレングスモデルによって高精度に故障率を算出することができる。
　一方、ユーザの体重や路面状況等の因子、すなわちユーザや環境に応じて変化する因子が入力ストレス分布に大きな相関を持つ場合、一般的なストレス－ストレングスモデルでは導出された故障率の精度も低くなってしまう。
　すなわち、従来の技術では、ユーザや環境に応じて変化する因子が入力ストレスに大きな相関を持つ場合に製品故障率を精度良く算出できないという問題があった。
　本発明は、上記の問題に鑑みてなされたものであり、ユーザ又は環境に応じて変化する因子が入力ストレスに大きな相関を持つ場合であっても製品故障率を精度良く算出する信頼性設計支援装置、信頼性設計支援方法、及びプログラムを提供することを主たる目的とする。
　本発明にかかる信頼性設計支援装置の一態様は、ユーザや環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する関係式算出手段と、前記関係式と、前記対象因子の指標データと、に応じて前記対象装置に入力される入力ストレスを算出する入力ストレス算出手段と、前記入力ストレスと、前記対象装置の特性値分布と、に応じて前記対象装置の故障率を算出する故障率算出手段と、を有する、ものである。
　また、この一態様において、前記指標データは、前記ユーザの身体的特徴に関する確率密度分布であり、前記入力ストレス算出手段は、前記関係式と、前記確率密度分布と、を基に前記入力ストレスに関する分布であるストレス分布を算出し、前記故障率算出手段は、前記ストレス分布と前記特性値分布と、の重複から前記対象装置の故障率を算出する、ことが望ましい。
　また、この一態様において、前記入力ストレス算出手段は、前記対象因子の最悪値を前記関係式に代入することにより閾値直線を算出し、前記故障率算出手段は、前記特性値分布と前記閾値直線が形成する分布領域から前記対象装置の故障率を算出する、ことが望ましい。
　さらに、この一態様において、前記対象装置に使用される部品形状データを使用して有限要素解析を行うことにより、前記入力ストレス算出手段が算出した前記入力ストレスを調整する入力ストレス分布調整手段を更に有する、ことが望ましい。
　なお、この一態様において、前記対象因子は、前記ユーザの身体的特徴、前記対象装置の使用環境、前記対象装置の使用時の実測値、または前記対象装置の特性値であり、前記指標データは、前記対象因子が前記ユーザの身体的特徴である場合に、前記ユーザの身体的特徴にかかる確率密度分布であり、前記対象因子が前記対象装置の使用時の実測値である場合に、当該実測値が複数列挙された実測データであり、前記対象因子が前記対象装置の特性値である場合に、前記対象装置の性質から算出できる最悪値である。
　本発明にかかる信頼性設計支援方法の一態様は、ユーザや環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する関係式算出ステップと、前記関係式と、前記対象因子の指標データと、に応じて前記対象装置に入力される入力ストレスを算出する入力ストレス算出ステップと、前記入力ストレスと、前記対象装置の特性値分布と、に応じて前記対象装置の故障率を算出する故障率算出ステップと、を実行するものである。
　本発明にかかるプログラムの一態様は、コンピュータに、ユーザや環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する関係式算出ステップと、前記関係式と、前記対象因子の指標データと、に応じて前記対象装置に入力される入力ストレスを算出する入力ストレス算出ステップと、前記入力ストレスと、前記対象装置の特性値分布と、に応じて前記対象装置の故障率を算出する故障率算出ステップと、を実行させる、ものである。
　本発明では、ユーザや環境に応じて変化する対象因子が入力ストレスに大きな相関を持つ場合であっても製品故障率を精度良く算出する信頼性設計支援装置、信頼性設計支援方法、及びプログラムを提供することができる。
実施の形態１にかかる信頼性設計支援装置１の構成を示すブロック図である。 実施の形態１にかかる信頼性設計支援装置１の動作を示す概念図である。 実施の形態１にかかる信頼性設計支援装置１の動作を示す概念図である。 実施の形態１にかかる信頼性設計支援装置１の動作を示す概念図である。 実施の形態１にかかる信頼性設計支援装置１を用いたリスクアセスメント手法を示すフローチャートである。 実施の形態１にかかる信頼性設計支援装置１のハードウェア構成を示すブロック図である。 実施の形態２にかかる信頼性設計支援装置１の構成を示すブロック図である。
＜実施の形態１＞
　以下、図面を参照して本発明の実施の形態について説明する。図１は、本実施の形態にかかる信頼性設計支援装置の構成を示すブロック図である。信頼性設計支援装置１は、上述の製品開発時におけるリスクアセスメントの各プロセス（上述の（１）～（４））を行う装置である。以下の説明では、上述のプロセス内の「（２－１）故障率の見積もり」に関する構成及び動作について主に説明する。
　なお、信頼性設計支援装置１が設計支援を行う対象とする装置（対象装置）は、任意の製品でよい。例えば対象装置は、人体と密着する電動アシスト自転車、小型移動モビリティ、装着型ロボット等である。これらの対象装置は、主に人体と接触する装置であり、ユーザ毎に変化する因子（例えば体重）及び環境（対象装置が使用される路面状況等）の因子が入力ストレスに対して支配的な（非常に影響の大きい）装置である。例えば、ユーザの体重が増減すると、製品に与えられるストレスが大きく変化する。
　信頼性設計支援装置１は、関係式算出手段１０、入力ストレス算出手段２０、特性値分布算出手段３０、及び故障率算出手段４０を有する。以下、各手段の詳細について説明する。
　関係式算出手段１０は、ストレス計測値ファイルを読み込む。ストレス計測値ファイルとは、ある因子（以下の説明では、対象因子とも記載する。）の値と、当該対象因子の値に対応するストレス計測値と、が関連付けられたデータ列が記述されたファイルである。ここで対象因子はユーザや環境に応じて変化し得る因子であり、ストレス計測値は対象装置を使用した場合に計測した値である。例えばストレス計測ファイルには、各ユーザの体重の重さと、各ユーザが対象装置を使用した場合のストレス計測値と、を関連付けたデータ列が記述されている。なお、ストレス計測ファイルに記述される対象因子として、体重に加えて、身長、靴の大きさ等の身体的な特徴や、身長と小型移動モビリティのハンドルの高さとの比等が挙げられる。
　関係式算出手段１０は、ストレス計測値ファイルのデータ列を一般的な近似手法を用いて解析することにより、対象因子の値とストレスの大きさとの関係を示す関係式を算出する。関係式算出手段１０は、たとえば近似手法として応答曲面法を用いて関係式を算出する。なお関係式算出手段１０は、近似手法として重回帰分析等を用いてもよい。関係式算出手段１０は、算出した関係式を入力ストレス算出手段２０に供給する。
　入力ストレス算出手段２０は、対象因子の指標データが記述された指標データファイルを読み出す。例えば指標データファイルは、対象因子の確率密度分布が記述されている。この場合に指標データファイルには、横軸を体重として縦軸を分布割合とするグラフ情報（またはこれに類する情報）が記述されている。なお、指標データファイルには、人体統計データのように信頼性の高い既知のデータを用いることが望ましい。入力ストレス算出手段２０は、関係式算出手段１０が算出した関係式と、指標データファイルと、を基に入力ストレス分布を算出する。
　指標データファイルには、上述のように確率密度分布の他に、実測した統計データや論理的に算出できる最悪値を記述しておくことも可能である。実測した統計データとして、例えば小型移動モビリティ使用時におけるユーザの両足の位置の分布を用いることができる。また論理的に算出できる最悪値として、例えば対象須地の設計によって定まるモータの最大トルク、最高速度、使用環境の化学物質濃度（ＩＥＣ（International Electrotechnical Commission）規格によって標準化されている値）等が該当する。
　入力ストレス算出手段２０が算出する入力ストレス分布とは、横軸にストレスの大きさを取り、縦軸に発生頻度を示す情報である。入力ストレス算出手段２０は、関係式算出手段１０が算出した関係式と指標データファイルに対し、たとえばモンテカルロ法を用いた解析を行うことにより入力ストレス分布を算出する。入力ストレス算出手段２０は、算出した入力ストレス分布を故障率算出手段４０に供給する。なお入力ストレス算出手段２０が算出する入力ストレスは、必ずしも分布の形式ではなく、後述する図３に示すように直線式等であってもよい。
　特性値分布算出手段３０は、特性値ファイルを読み出す。特性値ファイルは、対象装置に使用される部品の強度等の特性値が記述されたファイルである。特性値ファイルには、対象装置の該当部品を複数サンプリングし、各部品のストレスの大きさを計測した計測値が記述されている。ストレスの計測は、例えば一般的な引っ張り試験等により行う。特性値分布算出手段３０は、特性値ファイルから横軸をストレスとし、縦軸を発生頻度とする特性値分布を算出する。特性値分布算出手段３０は、算出した特性値分布を故障率算出手段４０に供給する。
　なお上述の説明では、各手段（関係式算出手段１０、入力ストレス算出手段２０、特性値分布算出手段３０）はファイル形式の情報を読み出したが、必ずしもこれに限られず、例えばデータベースに格納された各データを読み出すようにしてもよい。
　また上述の説明では、特性値分布算出手段３０が特性値ファイルから対象装置の特性値分布を算出するものとしたが、必ずしもこれに限られない。例えば、信頼性設計支援装置１に対して特性値分布を記述した入力ファイルが直接入力される構成であってもよい。
　故障率算出手段４０は、入力ストレス分布と特性値分布との重複に応じて、対象装置の故障率を算出する。故障率の算出手法の詳細は、図２を参照して後述する。故障率算出手段４０は、算出した故障率を任意の出力手段（例えばディスプレイ装置）や記憶装置（例えばファイルシステム）に出力する。
　続いて図２を参照して、信頼性設計支援装置１の第１の動作例を説明する。図２は、信頼性設計支援装置１の第１の動作例を示す概念図である。以下の説明では、対象因子はユーザの体重であるものとする。ユーザの体重は、ユーザ毎に異なる因子である。
　関係式算出手段１０は、ストレス計測値ファイルの全データ列を読み出す。当該データ列には、例えば体重Ｘｋｇのユーザと、体重Ｘｋｇのユーザが対象装置を使用した場合のストレス計測値と、が関連付けられて記載されている。関係式算出手段１０は、ストレス計測値ファイルの全データ列を応答曲面法等の近似手法の解析対象とすることにより関係式（図２（Ａ））を算出する。
　また入力ストレス算出手段２０は、図２（Ｂ）に示すような体重分布の統計情報を示す指標データファイルを読み出す。ここで指標データファイルは、対象装置の用途に応じて適宜記述されたファイルであれば良い。例えば、対象装置が日本人のみしか使用しない場合、全人類の体重の確率密度分布を用いるのではなく、日本人の体重の確率密度分布を用いてもよい。
　入力ストレス算出手段２０は、関係式（図２（Ａ））と体重分布の統計情報（図２（Ｂ））を用いた既知の解析（例えばモンテカルロ法による解析）を行うことにより入力ストレス分布（図２（Ｃ））を算出する。入力ストレス分布とは、前述のようにストレスの大きさと発生頻度の関係を示す分布である。
　特性値分布算出手段３０は、対象商品に用いられる部品の強度分析を行い、部品の特性値分布（例えば強度分布）を算出する（図２（Ｄ））。ここで部品の特性分布とは、図２（Ｄ）に示すようにストレスの大きさと故障の発生頻度との関係を表す分布である。
　故障率算出手段４０は、入力ストレス分布（図２（Ｃ）、図２（Ｄ））と、特性値分布（図２（Ｄ））を照合し、故障率を算出する。詳細には故障率算出部４０は、両分布の重複部分と、部品の特性値分布の全体と、の比を故障率として算出する。すなわち故障率算出部４０は、部品の全体個数を分母とし、重複部分にかかる部品個数を分子とした計算を行うことにより故障率を算出する。
　続いて図３を参照して、信頼性設計支援装置１の第２の動作例を説明する。図３は、信頼性設計支援装置１の第２の動作例を示す概念図である。以下の説明では、対象因子は対象装置が使用されるユーザの住環境の温度であるものとする。住環境の温度もユーザの使用状態に応じて変化する因子である。
　関係式算出手段１０は、ストレス計測値ファイルの全データ列を読み出す。当該データ列には、例えば温度と、当該温度において対象装置を使用した場合のストレス計測値と、が関連付けられて記載されている。関係式算出手段１０は、ストレス計測値ファイルの全データ列を応答曲面法等の近似手法の解析対象とすることにより関係式（図３（Ａ））を算出する。
　入力ストレス算出手段２０は、指標データファイルを読み出す。例えば指標データファイルには、対象装置の使用場所で測定した温度データが複数記述されている。入力ストレス算出手段２０は、当該温度データの中で最大値（最高温度であり、最悪値を示す。）を抽出する（図３（Ｂ））。入力ストレス算出手段２０は、抽出した最大値を関係式（図３（Ａ））に代入することにより入力ストレス最大値を示す閾値直線（図３（Ｃ））を算出する。
　故障率算出手段４０は、入力ストレス最大値を示す閾値直線（図３（Ｃ）、図３（Ｄ））と、特性値分布（図３（Ｄ））を照合して故障率を算出する。故障率算出手段４０は、図示するように閾値直線よりもストレスが小さい範囲を基に故障率を算出する。
　続いて図４を参照して、信頼性設計支援装置１の第３の動作例を説明する。図４は、信頼性設計支援装置１の第３の動作例を示す概念図である。以下の説明では、対象装置が小型移動モビリティであるものとする。また対象因子は、対象装置の使用時の速度分布であるものとする。当該動作例は、図２を参照して説明した第１の動作例と指標データファイルのみが異なるため、指標データファイルの内容について以下に説明する。
　指標データファイルには、図４（Ｂ）に示すように対象装置の内部制御により上限値（最大速度）が制限される速度分布が記載されている。入力ストレス算出手段２０は、関係式（図４（Ａ））と速度分布（図４（Ｂ））を用いた既知の解析（例えばモンテカルロ法による解析）を行うことにより入力ストレス分布（図４（Ｃ））を算出する。
　続いて、図５を参照して信頼性設計支援装置１を用いたリスクアセスメントの一手法について説明する。はじめに設計者は、リスクアセスメントの対象となる対象装置のユーザ（または対象装置と接触する可能性があるもの）を想定する（Ｓ１）。そして設計者は、想定したユーザ（または対象装置と接触する可能性があるもの）による使用方法及び予見可能な誤使用を明確化する（Ｓ１）。
　次に設計者は、対象装置の使用及び誤使用により引き起こされる挙動が原因として生じるハザードを特定する（Ｓ３－１）。そして設計者は、ハザードがユーザに与える影響を基に、リスクの程度を見積もる（Ｓ３－２）。続いて設計者は、見積もったリスクの程度や対象装置の出荷予定台数等を基に、許容可能なリスクの発生頻度、すなわち許容故障率を定義する（Ｓ４）。例えば設計者は、対象装置の全出荷予定台数に対して１台もユーザに大きな影響を与えない確率となるように許容故障率を定義する。
　Ｓ３－１～Ｓ４の工程と並行して、設計者は上述の信頼性設計支援装置１の設定を行う（Ｓ２－１－ａ）。そして、設計者は信頼性設計支援装置１に入力する各種ファイル（指標データファイル、ストレス計測値ファイル、特性値ファイル）の準備を行う。
　信頼性設計支援装置１は、上述の手法により故障率を算出し、算出した故障率と許容故障率を比較する。算出した故障率が許容故障率よりも小さくなるまで、設計者は対象装置の設計を適宜調整する。なお、この許容故障率を扱うための信頼性設計支援装置１は、図１の構成に加えて故障率算出手段４０の後段に故障率評価手段（図示せず）を有していればよい。故障率評価手段は、故障率算出手段４０が算出した故障率と、設計者が入力する許容故障率と、を比較する。
　次に図６を参照して信頼性設計支援装置１のハードウェア構成について説明する。信頼性設計支援装置１は、ハードウェアとして入力装置１１０、ＣＰＵ（Central Processing Unit）１２０、ＲＯＭ（Read Only Memory）１３０、ＲＡＭ（Random Access Memory）１４０、及び出力装置１５０を有する。なお、図６に示す構成は一例であり、その他の構成であってもよい。また、ＵＳＢ（Universal Serial Bus）メモリ等の記憶装置を着脱可能に構成してもよい。
　入力装置１１０は、キーボード、マウス等の情報入力装置である。ＲＯＭ１３０は、ＢＩＯＳ（Basic Input/Output System）やＩＰＬ（Initial Program Loader）などの記憶に用いられる記憶装置である。ＲＡＭ１４０は、プログラムを保持し、プログラムを実行するＣＰＵ１２０のワーキングエリア等として用いられるメモリである。ＣＰＵ１２０は、プログラムの読み込み及び実行を行い、出力装置１５０への出力制御を行う。出力装置１５０は、ＣＰＵ１２０の制御によって演算処理結果等を出力する装置である。出力装置１５０は、例えば液晶ディスプレイ装置やプリンタ装置等である。
　続いて本実施の形態にかかる信頼性設計支援装置１の効果について説明する。信頼性設計支援装置１は、上述のようにユーザ毎または対象装置の使用環境に応じて変化する対象因子の確率分布や最大値／最小値を基に入力ストレス（例えば図２（Ｃ）、図３（Ｃ））を算出している。そして信頼性設計支援装置１は、この入力ストレスを用いて対象装置の故障率を算出している。すなわち信頼性設計支援装置１は、因子の値の変化を考慮した入力ストレスの算出を行っている。そのため信頼性設計支援装置１は、ユーザや環境によって変化する変化因子が入力ストレスに対して支配的な場合であっても、適切な故障率の算出が可能となる。特に信頼性設計支援装置１は、体重や身長等の身体的な特徴が入力ストレスに対して大きな影響を与える人体サポート機器（例えば電動アシスト自転車、小型移動モビリティ、装着型ロボット等）を対象とする場合であっても、その故障率を正確に算出することができる。これにより、高信頼かつ低コストを実現した対象装置の設計が可能となる。
　信頼性設計支援装置１は、図２に示すように対象因子としてユーザの身体的特徴（体重、身長等）を考慮して故障率を算出することができる。この場合に信頼性設計支援装置１は、母数が大きな統計的な確率密度分布、すなわち信頼性の高い分布を基に入力ストレス分布を算出することができる。よって信頼性設計支援装置１は、信頼性の高い入力ストレス分布を算出することができる。そのため信頼性設計支援装置１は、故障率を精度よく算出することが可能となる。
　信頼性設計支援装置１は、図３に示すように対象因子の最大値（最悪値）を基に故障率を算出することもできる。これにより信頼性設計支援装置１は、対象装置が最も劣悪な環境で使用された場合（例えば最も高温下で使用された場合）を考慮した故障率を算出することができる。この故障率を評価することにより、より安全な商品設計を実現することができる。なお、上述の説明では対象因子の最大値を用いたが必ずしもこれに限られず、例えば因子の性質に応じて最小値を考慮した故障率を算出してもよい。
＜実施の形態２＞
　本実施の形態にかかる信頼性設計支援装置１は、対象装置に使用される部品の形状特性を考慮した故障率算出を行うことを特徴とする。本実施の形態にかかる信頼性設計支援装置１について、実施の形態１と異なる点を以下に説明する。なお、以下の説明において実施の形態１と同一の処理部については、同一名称及び同一符号を付し、その詳細な説明は省略する。
　図７は、本実施の形態にかかる信頼性設計支援装置１の構成を示すブロック図である。本実施の形態にかかる信頼性設計支援装置１は、図１の構成に加えて入力ストレス分布調整手段５０を更に備える構成である。
　入力ストレス算出手段２０は、上述の手法でストレス分布を算出する。この際に、入力ストレス算出手段２０は、横軸のストレスをニュートン単位とする値で算出し、縦軸を発生頻度とする入力ストレス分布を算出し、当該入力ストレス分布を入力ストレス調整手段５０に供給する。
　入力ストレス調整手段５０には、上述の入力ストレス分布に加え、部品形状ファイルが入力される。部品形状ファイルは、対象装置を構成する部品の形状を示すデータが格納されたファイルである。当該部品形状ファイルは、後述する有限要素解析に使用できるフォーマットでデータ（部品形状データ）が格納されている。
　入力ストレス調整手段５０は、部品形状ファイルと、入力ストレス分布と、を基に有限要素解析を行い、調整済みの入力ストレス分布を算出する。この入力ストレス分布は、横軸のストレスをメガパスカル（Ｍｐａ）単位とする値で算出し、縦軸を発生頻度とする分布である。入力ストレス調整手段５０は、部品形状ファイルを用いて調整した入力ストレス分布（横軸のストレスをメガパスカル（Ｍｐａ）単位とするストレス分布）を故障率算出手段４０に供給する。故障率算出手段４０は、実施の形態１と同様の手法で故障率を算出すればよい。
　続いて本実施の形態にかかる信頼性設計支援装置１の効果について説明する。上述のように本実施の形態では、対象装置に使用される部品形状をも考慮した入力ストレス分布を算出している。これにより、本実施の形態にかかる信頼性設計支援装置１は、より精度の高い故障率の算出を実現することができる。
　上述の実施の形態１及び実施の形態２にかかる信頼性設計支援装置１の各手段（関係式算出手段１０、入力ストレス算出手段２０、特性値分布算出手段３０、故障率算出手段４０、入力ストレス調整手段５０）は、ＣＰＵ１２０によってプログラムの形式で実行される。プログラムは、様々なタイプの非一時的なコンピュータ可読媒体（non-transitory computer readable medium）を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体（tangible storage medium）を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体（例えばフレキシブルディスク、磁気テープ、ハードディスクドライブ）、光磁気記録媒体（例えば光磁気ディスク）、ＣＤ－ＲＯＭ（Read Only Memory）、ＣＤ－Ｒ、ＣＤ－Ｒ／Ｗ、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ（Programmable ROM）、ＥＰＲＯＭ（Erasable PROM）、フラッシュＲＯＭ、ＲＡＭ（random access memory））を含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体（transitory computer readable medium）によってコンピュータに供給されてもよい。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。
　以上、本発明を上記実施形態に即して説明したが、本発明は上記実施形態の構成にのみ限定されるものではなく、本願特許請求の範囲の請求項の発明の範囲内で当業者であればなし得る各種変形、修正、組み合わせを含むことは勿論である。
　この出願は、２０１３年２月２１日に出願された日本出願特願２０１３－０３１７４７を基礎とする優先権を主張し、その開示の全てをここに取り込む。
１　信頼性設計支援装置
１０　関係式算出手段
２０　入力ストレス算出手段
３０　特性値分布算出手段
４０　故障率算出手段
５０　入力ストレス調整手段
１１０　入力装置
１２０　ＣＰＵ
１３０　ＲＯＭ
１４０　ＲＡＭ
１５０　出力装置 

Claims (7)

  1. 　ユーザや環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する関係式算出手段と、
    　前記関係式と、前記対象因子の指標データと、に応じて前記対象装置に入力される入力ストレスを算出する入力ストレス算出手段と、
    　前記入力ストレスと、前記対象装置の特性値分布と、に応じて前記対象装置の故障率を算出する故障率算出手段と、
    　を有する、信頼性設計支援装置。
  2. 　前記指標データは、前記ユーザの身体的特徴に関する確率密度分布であり、
    　前記入力ストレス算出手段は、前記関係式と、前記確率密度分布と、を基に前記入力ストレスとしてストレス分布を算出し、
    　前記故障率算出手段は、前記ストレス分布と前記特性値分布と、の重複から前記対象装置の故障率を算出する、
    　ことを特徴とする請求項１に記載の信頼性設計支援装置。
  3. 　前記入力ストレス算出手段は、前記対象因子の最悪値を前記関係式に代入することにより閾値直線を算出し、
    　前記故障率算出手段は、前記特性値分布と前記閾値直線が形成する分布領域から前記対象装置の故障率を算出する、
    　ことを特徴とする請求項１に記載の信頼性設計支援装置。
  4. 　前記対象装置に使用される部品形状データを使用して有限要素解析を行うことにより、前記入力ストレス算出手段が算出した前記入力ストレスを調整する入力ストレス分布調整手段を更に有する、
    　ことを特徴とする請求項１～請求項３のいずれか１項に記載の信頼性設計支援装置。
  5. 　前記対象因子は、前記ユーザの身体的特徴、前記対象装置の使用環境、前記対象装置の使用時の実測値、または前記対象装置の特性値であり、
    　前記指標データは、
    　前記対象因子が前記ユーザの身体的特徴である場合に、前記ユーザの身体的特徴にかかる確率密度分布であり、
    　前記対象因子が前記対象装置の使用時の実測値である場合に、当該実測値が複数列挙された実測データであり、
    　前記対象因子が前記対象装置の特性値である場合に、前記対象装置の性質から算出できる最悪値である、請求項１～請求項４のいずれか１項に記載の信頼性設計支援装置。
  6. 　ユーザや環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する関係式算出ステップと、
    　前記関係式と、前記対象因子の指標データと、に応じて前記対象装置に入力される入力ストレスを算出する入力ストレス算出ステップと、
    　前記入力ストレスと、前記対象装置の特性値分布と、に応じて前記対象装置の故障率を算出する故障率算出ステップと、
    　を実行する信頼性設計支援方法。
  7. 　コンピュータに、
    　ユーザや環境に応じて変化する対象因子の値と、対象装置の使用時に計測した入力ストレス計測値と、を関連づけたデータ列に対して近似計算を行って関係式を算出する関係式算出ステップと、
    　前記関係式と、前記対象因子の指標データと、に応じて前記対象装置に入力される入力ストレスを算出する入力ストレス算出ステップと、
    　前記入力ストレスと、前記対象装置の特性値分布と、に応じて前記対象装置の故障率を算出する故障率算出ステップと、を実行させる、
    　信頼性設計支援プログラム。
PCT/JP2014/000537 2013-02-21 2014-01-31 信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム Ceased WO2014129125A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201480009173.8A CN104995624B (zh) 2013-02-21 2014-01-31 可靠性设计辅助装置以及可靠性设计辅助方法
DE112014000935.0T DE112014000935T5 (de) 2013-02-21 2014-01-31 Zuverlässigkeitsentwurf-Unterstützungsvorrichtung, Zuverlässigkeitsentwurf-Unterstützungsverfahren und Zuverlässigkeitsentwurf-Unterstützungsprogramm
US14/768,937 US11085843B2 (en) 2013-02-21 2014-01-31 Reliability design assistance device, reliability design assistance method, and reliability design assistance program
KR1020157022249A KR101658748B1 (ko) 2013-02-21 2014-01-31 신뢰성 설계 지원 장치 및 신뢰성 설계 지원 방법

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013031747A JP5761225B2 (ja) 2013-02-21 2013-02-21 信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム
JP2013-031747 2013-02-21

Publications (1)

Publication Number Publication Date
WO2014129125A1 true WO2014129125A1 (ja) 2014-08-28

Family

ID=51390932

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/000537 Ceased WO2014129125A1 (ja) 2013-02-21 2014-01-31 信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム

Country Status (6)

Country Link
US (1) US11085843B2 (ja)
JP (1) JP5761225B2 (ja)
KR (1) KR101658748B1 (ja)
CN (1) CN104995624B (ja)
DE (1) DE112014000935T5 (ja)
WO (1) WO2014129125A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101705249B1 (ko) * 2016-02-19 2017-02-10 울산과학기술원 유한요소법을 이용한 구조물 신뢰성 평가를 위한 인터페이스, 시스템 및 방법
CN110276093B (zh) * 2019-04-29 2023-06-02 北京圣涛平试验工程技术研究院有限责任公司 电线电缆可靠性评估方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002149714A (ja) * 2000-08-31 2002-05-24 Toshiba Corp 信頼性設計支援装置および信頼性設計支援方法およびプログラムを記録した媒体およびプログラム製品
JP2007108843A (ja) * 2005-10-11 2007-04-26 Toshiba Corp 半導体装置設計支援方法、半導体装置設計支援システム、半導体装置設計支援プログラム

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06195370A (ja) * 1992-12-22 1994-07-15 Mitsubishi Kasei Corp 数値データの広がりの評価方法
JPH06348683A (ja) * 1993-04-16 1994-12-22 Sony Corp 集積回路のシミュレーション方法
JP4613555B2 (ja) 2004-09-08 2011-01-19 サンスター株式会社 有限要素解析による歯ブラシ柄部の設計方法、及び該設計に用いるコンピュータ解析システム
US8290753B2 (en) * 2006-01-24 2012-10-16 Vextec Corporation Materials-based failure analysis in design of electronic devices, and prediction of operating life
DE102010024683A1 (de) 2010-06-23 2011-12-29 Siemens Aktiengesellschaft Verfahren zum rechnergestützten Ermitteln eines von einem technischen System ausgehenden Risikos für einen Nutzer
CN102589891B (zh) * 2012-03-02 2013-12-18 北京理工大学 一种车用电机驱动系统耐久寿命的估计方法
CN102663201B (zh) * 2012-04-24 2013-12-25 北京航空航天大学 一种考虑电气互连结构可靠性的电子产品可靠性预计方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002149714A (ja) * 2000-08-31 2002-05-24 Toshiba Corp 信頼性設計支援装置および信頼性設計支援方法およびプログラムを記録した媒体およびプログラム製品
JP2007108843A (ja) * 2005-10-11 2007-04-26 Toshiba Corp 半導体装置設計支援方法、半導体装置設計支援システム、半導体装置設計支援プログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TETSURO SEKI, JITSUMU NI YAKUDATSU SERIES 3 SHINRAISEI DATA KAISEKI NYUMON, 5 October 2011 (2011-10-05), pages 152 - 154 *

Also Published As

Publication number Publication date
CN104995624B (zh) 2018-05-08
DE112014000935T5 (de) 2015-11-26
JP2014164308A (ja) 2014-09-08
US20160003695A1 (en) 2016-01-07
US11085843B2 (en) 2021-08-10
CN104995624A (zh) 2015-10-21
KR101658748B1 (ko) 2016-09-21
JP5761225B2 (ja) 2015-08-12
KR20150108884A (ko) 2015-09-30

Similar Documents

Publication Publication Date Title
CN103323703B (zh) 一种电缆接头故障诊断方法
WO2025036056A1 (zh) 一种基于熵值算法的动力电池电压故障在线诊断处理方法
JP6141235B2 (ja) 時系列データにおける異常を検出する方法
CN103487682A (zh) 一种电压暂降扰动下敏感客户电能体验质量的预警方法
CN103827653B (zh) 检测异常的发生的方法、设备及计算机程序
CN105512812A (zh) 基于动态仿真模型的核电厂设备故障预警分析方法及系统
RU2016142123A (ru) Способ оценки нормальности или ненормальности измеренного значения физического параметра двигателя летательного аппарата
CN111340367B (zh) 基于多维度关联性数据分析的少油设备状态评估方法及系统
Karbalaei et al. A comparison of voltage collapse proximity indicators
CN104075825A (zh) 一种电力电缆光纤测温基准数据测量方法
CN108681319B (zh) 一种基于传递函数的变压器绕组故障识别方法
CN113177343A (zh) 一种输电线路弧垂计算方法、设备及计算机介质
JP5761225B2 (ja) 信頼性設計支援装置、信頼性設計支援方法、及び信頼性設計支援プログラム
JP2015052487A (ja) 材料強度評価装置
CN119989759A (zh) 一种基于数字孪生技术的电缆接头的状态评估方法和系统
CN104713730A (zh) 一种根据振动信号确定飞机发动机退化率的方法
CN115219916A (zh) 电池温升情况检测方法、装置和电子设备
JP7107386B2 (ja) 情報処理装置
JP2016512949A (ja) コンピュータ支援により送配電網のインピーダンスを求める方法
KR20170140753A (ko) 소프트웨어 안전성 분석 방법 및 장치
CN120086610A (zh) 一种基于数据分析的保护器成品质量检测方法及系统
CN105891645A (zh) 确定防滑刹车控制装置振动故障分布的方法
JP2013200245A (ja) 前処理方法
CN110567576A (zh) 一种变电站厂界噪声超标原因的确定方法及装置
CN103644884A (zh) 基于温度测量的电力电缆接头施工规范性检测方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14754057

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 20157022249

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 1120140009350

Country of ref document: DE

Ref document number: 112014000935

Country of ref document: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14754057

Country of ref document: EP

Kind code of ref document: A1