JP4936135B2 - Deformation simulation method of membrane structure - Google Patents

Deformation simulation method of membrane structure Download PDF

Info

Publication number
JP4936135B2
JP4936135B2 JP2007255902A JP2007255902A JP4936135B2 JP 4936135 B2 JP4936135 B2 JP 4936135B2 JP 2007255902 A JP2007255902 A JP 2007255902A JP 2007255902 A JP2007255902 A JP 2007255902A JP 4936135 B2 JP4936135 B2 JP 4936135B2
Authority
JP
Japan
Prior art keywords
membrane structure
film
membrane
data
fluid
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.)
Expired - Fee Related
Application number
JP2007255902A
Other languages
Japanese (ja)
Other versions
JP2009087025A (en
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.)
Yokohama Rubber Co Ltd
Original Assignee
Yokohama Rubber Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yokohama Rubber Co Ltd filed Critical Yokohama Rubber Co Ltd
Priority to JP2007255902A priority Critical patent/JP4936135B2/en
Publication of JP2009087025A publication Critical patent/JP2009087025A/en
Application granted granted Critical
Publication of JP4936135B2 publication Critical patent/JP4936135B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Description

本発明は、膜構造物の変形シミュレーション方法に関する。   The present invention relates to a deformation simulation method for a membrane structure.

超音波を用いて海底の地形を計測するソナーを備えた船舶がある。このような船舶においては、船体が海中に臨む箇所にソナーを収容する収容凹部を設けるとともに、この収容凹部を覆うように超音波が通過可能な膜から構成された膜構造物を取り付けてソナーの保護を図っている。
そして、膜構造物と収容凹部との間に所定の圧力で海水を満たすことにより膜構造物の形状を保っている。
この場合、船体が大きく上下に揺動すると、膜構造物がいったん海面の上に現れたのち、海中に没することになり、膜構造物が空中から海面に没する際に海水から圧力を受けて変形する。
もし、膜構造物が大きく変形してソナーに接触するとソナーの損傷が懸念されるため、膜構造物は、それが変形してもソナーに接触しない程度の強度または遊間、形状保持のための内圧を確保することが必要となる。
そのため、膜構造物の設計においては、膜構造物が周囲の流体からの圧力変動を受けた場合にどのように変形するかを、シミュレーションによって解析しておくことが必要となる。
従来から構造物の構造解析に際しては有限要素法が用いられている(特許文献1参照)。
ところが、上述したように膜構造物が周囲の海水などの流体からの圧力を受けて大きく変形する現象を解析する場合には、有限要素法による膜構造物の構造解析と、流体の挙動を再現する流体解析とを連動して行うことが、膜構造物の変形挙動を正確に再現する上で必要となる。
特開2006−175994
There are ships equipped with sonar that measures the topography of the seabed using ultrasound. In such a ship, an accommodation recess for accommodating sonar is provided at a location where the hull faces the sea, and a membrane structure composed of a membrane through which ultrasonic waves can pass is attached so as to cover the accommodation recess. We are trying to protect it.
And the shape of a membrane structure is maintained by filling seawater with a predetermined pressure between a membrane structure and an accommodation recessed part.
In this case, if the hull swings up and down greatly, the membrane structure will appear once on the sea surface and then will be submerged in the sea. When the membrane structure is submerged from the air to the sea surface, it will receive pressure from seawater. And deform.
If the membrane structure is greatly deformed and contacts the sonar, the sonar may be damaged. Therefore, the membrane structure is not strong enough to contact the sonar even if it is deformed. It is necessary to ensure.
Therefore, in the design of the membrane structure, it is necessary to analyze by simulation how the membrane structure deforms when it receives pressure fluctuations from the surrounding fluid.
Conventionally, the finite element method has been used for structural analysis of structures (see Patent Document 1).
However, as described above, when analyzing a phenomenon in which the membrane structure undergoes large deformation under pressure from a fluid such as surrounding seawater, the structure analysis of the membrane structure by the finite element method and the behavior of the fluid are reproduced. In order to accurately reproduce the deformation behavior of the membrane structure, it is necessary to perform the fluid analysis in conjunction with the fluid analysis.
JP 2006-175994 A

しかしながら、上述した構造解析と流体解析とを連動させて行う場合には、膨大な計算を行うことが必要となるため、解析に要するコストが多大なものとなり、現実的には不可能であった。
仮に、膜構造物の周囲の流体の挙動を考慮せず、流体から加わる圧力のみを考慮して有限要素法による膜構造物の構造解析を行った場合には、現実には起こりえない過剰な変形結果となり、膜構造物の変形挙動を正確に得ることができないものであった。
本発明は、このような事情に鑑みてなされたものであり、その目的は、コストを抑制しつつ膜構造物の変形挙動を正確に再現する上で有利な膜構造物の変形シミュレーション方法を提供することにある。
However, when the structural analysis and the fluid analysis described above are performed in conjunction with each other, it is necessary to perform an enormous amount of calculation, which increases the cost required for the analysis, which is practically impossible. .
If the structure of the membrane structure is analyzed by the finite element method, considering only the pressure applied from the fluid without considering the behavior of the fluid around the membrane structure, an excessive amount that cannot occur in reality will occur. As a result of deformation, the deformation behavior of the membrane structure could not be obtained accurately.
The present invention has been made in view of such circumstances, and an object thereof is to provide a membrane structure deformation simulation method that is advantageous in accurately reproducing the deformation behavior of the membrane structure while suppressing cost. There is to do.

上記目的を達成するために本発明は、流体中において膜で構成された膜構造物に圧力が作用した際に発生する前記膜構造物の変形挙動を有限要素法によりコンピュータを用いて解析する膜構造物の変形シミュレーション方法であって、前記コンピュータは、操作者による操作入力を受け付ける入力手段と、有限要素法による解析処理を行う解析手段と、データの出力を行う出力手段とを備え、操作者が前記入力手段に対して操作入力を行うことにより前記膜構造物の形状を示す形状データを前記解析手段に与える形状データ入力ステップと、前記解析手段が前記形状データに基づいて前記膜構造物を要素分割し複数の要素と複数の節点を設定することで構造データを生成する構造データ生成ステップと、操作者が前記入力手段に対して操作入力を行うことにより前記膜の材料特性を示す材料特性データを前記解析手段に設定する材料特性設定ステップと、操作者が前記入力手段に対して操作入力を行うことにより前記膜構造物の拘束条件を示す拘束条件データを前記解析手段に設定する拘束条件設定ステップと、操作者が前記入力手段に対して操作入力を行うことにより前記膜構造物に作用する圧力を示す圧力データを前記解析手段に設定する圧力設定ステップと、前記解析手段が、前記構造データ生成ステップで生成された構造データと前記材料特性設定ステップで設定された材料特性データと前記拘束条件設定ステップで設定された拘束条件データと前記圧力設定ステップで設定された圧力データとに基づいて有限要素法による計算を行い前記膜構造物の変形挙動を示すデータを前記出力手段を介して出力する出力ステップとを含み、前記材料特性設定ステップによって設定される前記材料特性データは前記膜の密度を含み、前記材料特性設定ステップは、前記膜を構成する材料によって特定される密度に、前記膜に接する流体の単位面積当たりの質量を付加質量として加算した値を前記膜の密度として設定する付加質量加算ステップを含むことを特徴とする。   In order to achieve the above object, the present invention provides a membrane for analyzing the deformation behavior of the membrane structure generated when a pressure is applied to the membrane structure composed of a membrane in a fluid using a finite element method. A structure deformation simulation method, wherein the computer includes an input unit that receives an operation input by an operator, an analysis unit that performs an analysis process using a finite element method, and an output unit that outputs data. A shape data input step for giving shape data indicating the shape of the membrane structure to the analysis means by performing an operation input to the input means, and the analysis means is configured to input the membrane structure based on the shape data. A structure data generation step for generating structure data by dividing an element and setting a plurality of elements and a plurality of nodes, and an operator inputs an operation to the input means. A material property setting step for setting the material property data indicating the material property of the film in the analysis unit by performing an operation input to the input unit by the operator to set the constraint condition of the film structure. A constraint condition setting step for setting the constraint condition data to be indicated in the analysis means, and pressure data indicating the pressure acting on the membrane structure when the operator performs an operation input to the input means. The pressure setting step, the analysis means includes the structure data generated in the structure data generation step, the material property data set in the material property setting step, the constraint condition data set in the constraint condition setting step, and the Data showing the deformation behavior of the membrane structure by calculation using the finite element method based on the pressure data set in the pressure setting step Output through the output means, the material property data set by the material property setting step includes the density of the film, and the material property setting step is specified by the material constituting the film And an additional mass adding step of setting a value obtained by adding the mass per unit area of the fluid in contact with the membrane as an additional mass to the density as the density of the membrane.

本発明によれば、膜構造物の膜の密度に膜に接する流体の単位面積当たりの質量を付加質量として加算して膜構造物の変形挙動を解析するようにしたので、流体中における膜構造物の変形挙動を正確に再現する上で有利となり、流体解析を必要としないのでコストの抑制を図る上でも有利となる。   According to the present invention, since the mass per unit area of the fluid in contact with the membrane is added to the membrane density of the membrane structure as an additional mass to analyze the deformation behavior of the membrane structure, the membrane structure in the fluid This is advantageous in accurately reproducing the deformation behavior of an object, and is also advantageous in reducing costs because fluid analysis is not required.

次に、本発明の実施の形態による膜構造物の変形シミュレーション方法について図面を参照して説明する。
図1は本発明方法を実行するために使用されるコンピュータ10の構成を示すブロック図、図2はコンピュータ10の機能ブロック図である。
まず、本発明方法を実行するコンピュータコンピュータ10について説明する。
図1に示すように、コンピュータ10は、CPU12と、不図示のインターフェース回路およびバスラインを介して接続されたROM14、RAM16、ハードディスク装置18、ディスク装置20、キーボード22、マウス24、ディスプレイ26、プリンタ28などを有している。
ROM14は制御プログラムなどを格納し、RAM16はワーキングエリアを提供するものである。
ハードディスク装置18は本発明に係る膜構造物の変形シミュレーション方法を実現するための解析プログラムを格納している。
ディスク装置20はCDやDVDなどの記録媒体に対してデータの記録および/または再生を行うものである。
キーボード22およびマウス24は、操作者による操作入力を受け付けるものであり、入力手段10A(図2)を構成するものである。
ディスプレイ26はデータを表示出力するものであり、プリンタ28はデータを印刷出力するものであり、ディスプレイ26およびプリンタ28によってデータを出力する出力手段10C(図2)が構成されている。
Next, a deformation simulation method for a membrane structure according to an embodiment of the present invention will be described with reference to the drawings.
FIG. 1 is a block diagram showing a configuration of a computer 10 used for executing the method of the present invention, and FIG. 2 is a functional block diagram of the computer 10.
First, the computer 10 that executes the method of the present invention will be described.
As shown in FIG. 1, a computer 10 includes a CPU 12, a ROM 14, a RAM 16, a hard disk device 18, a disk device 20, a keyboard 22, a mouse 24, a display 26, and a printer connected via an interface circuit (not shown) and a bus line. 28 and so on.
The ROM 14 stores a control program and the like, and the RAM 16 provides a working area.
The hard disk device 18 stores an analysis program for realizing the film structure deformation simulation method according to the present invention.
The disk device 20 records and / or reproduces data on a recording medium such as a CD or a DVD.
The keyboard 22 and the mouse 24 receive an operation input by the operator, and constitute the input means 10A (FIG. 2).
The display 26 displays and outputs data, and the printer 28 prints and outputs data. The display 26 and the printer 28 constitute output means 10C (FIG. 2).

図2に示すように、コンピュータ10は、機能的には、入力手段10A、解析手段10B、出力手段10Cを含んで構成されている。
入力手段10Aは、操作者による操作入力を受け付けることにより、膜構造物を有限要素法によって解析するにあたって必要なデータを入力するものであり、それらデータについては後述する。
解析手段10Bは、入力手段10Aによって入力されたデータに基づいて有限要素法による解析処理を行うものであり、図1に示すハードディスク装置18に格納されている解析プログラムがRAM16にロードされ、CPU12が解析プログラムに基づいて動作することで実現される。
出力手段10Cは、解析手段10Bによる計算結果から構成されるデータを出力するものである。
As shown in FIG. 2, the computer 10 functionally includes an input means 10A, an analysis means 10B, and an output means 10C.
The input means 10A inputs data necessary for analyzing the membrane structure by the finite element method by receiving an operation input by the operator, and these data will be described later.
The analysis means 10B performs analysis processing by the finite element method based on the data input by the input means 10A. The analysis program stored in the hard disk device 18 shown in FIG. This is realized by operating based on an analysis program.
The output means 10C outputs data composed of the calculation results by the analysis means 10B.

次に、膜構造物の変形シミュレーション方法の手順について説明する。
図3は本実施の形態の膜構造物30の変形シミュレーション方法を示すフローチャートチャート、図4(A)は膜構造物30の外面が海面上に臨んだ状態を模式的に示す断面説明図、(B)は膜構造物30の外面が海面下に水没した状態を模式的に示す断面説明図である。なお、図4(A)、(B)は、膜構造物体30が設けられた船体2をその前後方向に対して直交する平面で破断した状態で示している。
本実施の形態では、図4(A)に示すように、膜構造物30が、船舶に設けられた不図示のソナーを保護するためのカバーを構成するものである場合について説明する。
すなわち、超音波を用いて海底の地形を計測するソナーを備えた船舶において、船舶の船体2が海中に臨む箇所にソナーを収容する収容凹部4が設けられている。
膜構造物30は、膜32と、膜32の周囲に形成された取り付け部(不図示)とを備えている。
膜32は、超音波が通過可能で可撓性を有する材料、例えば、ゴムと、ゴムを補強する補強材料とによって構成されており、このような補強材料として繊維材料が採用可能である。
膜構造物30は、収容凹部4を覆うように前記取り付け部が収容凹部4の開口を形成する船体2部分に例えばねじなどにより取着され、これにより、膜構造物30と収容凹部4との間に閉鎖された内側空間が形成される。
Next, the procedure of the deformation simulation method for the membrane structure will be described.
FIG. 3 is a flowchart chart showing a deformation simulation method of the membrane structure 30 of the present embodiment, FIG. 4A is a cross-sectional explanatory view schematically showing a state in which the outer surface of the membrane structure 30 faces the sea surface, B) is a cross-sectional explanatory view schematically showing a state in which the outer surface of the membrane structure 30 is submerged under the sea surface. 4A and 4B show a state in which the hull 2 provided with the membrane structure object 30 is broken along a plane orthogonal to the front-rear direction.
In the present embodiment, as shown in FIG. 4A, a case where the membrane structure 30 constitutes a cover for protecting a sonar (not shown) provided in the ship will be described.
That is, in a ship provided with a sonar that measures the topography of the seabed using ultrasonic waves, an accommodation recess 4 that accommodates the sonar is provided at a location where the hull 2 of the ship faces the sea.
The membrane structure 30 includes a membrane 32 and an attachment portion (not shown) formed around the membrane 32.
The membrane 32 is made of a flexible material through which ultrasonic waves can pass, for example, rubber, and a reinforcing material that reinforces the rubber, and a fiber material can be adopted as such a reinforcing material.
The membrane structure 30 is attached to the hull 2 part where the mounting portion forms the opening of the accommodation recess 4 so as to cover the accommodation recess 4 by, for example, screws, and thereby the membrane structure 30 and the accommodation recess 4. A closed inner space is formed between them.

そして、膜構造物30は、内面32Aと収容凹部4との間に所定の圧力で海水を満たすことにより膜32が円筒面上を延在するように形状を保持している。
具体的には、図4(A)に示すように、膜32は、船体2の前後方向と並行する軸心Oを中心とした半径Rの円筒面上を延在しており、膜32の2つの面のうち一方の面が前記内側空間(収容凹部4)に臨む内面32Aを構成し、他方の面が膜32を挟んで前記内側空間の反対側に位置する外側空間に臨む外面32Bを構成している。
また、本実施の形態では、膜32の輪郭(膜32を平面視したときの形状)は矩形であり、図5に示すように、膜32は軸心Oに沿った長さLと、軸心Oと直交する幅Wとを有し、長さLは幅Wよりも大きな寸法であり、例えば、L>10Wである。
The membrane structure 30 holds the shape so that the membrane 32 extends on the cylindrical surface by filling seawater with a predetermined pressure between the inner surface 32 </ b> A and the housing recess 4.
Specifically, as shown in FIG. 4A, the membrane 32 extends on a cylindrical surface having a radius R centering on an axis O parallel to the longitudinal direction of the hull 2. One of the two surfaces constitutes an inner surface 32A facing the inner space (accommodating recess 4), and the other surface has an outer surface 32B facing the outer space located on the opposite side of the inner space with the membrane 32 in between. It is composed.
Further, in the present embodiment, the outline of the film 32 (the shape when the film 32 is viewed in plan view) is a rectangle, and the film 32 has a length L along the axis O and an axis as shown in FIG. It has a width W orthogonal to the center O, and the length L is larger than the width W, for example, L> 10W.

次に膜構造物の変形シミュレーション方法について図2および図3のフローチャートを参照しつつ詳しく説明する。
まず、操作者が入力手段10Aに対して操作入力を行うことにより、膜構造物30の形状を示す形状データD1(図2)を解析手段10Bに与える(ステップS10:形状データ入力ステップ)。
次に、解析手段10Bは形状データD1に基づいて膜構造物30を要素分割し複数の要素と複数の節点を設定することで構造データを生成する(ステップS12:構造データ生成ステップ)。
ここで、操作者は、膜構造物30を解析するにあたって最適な構造モデルを決定し、その構造モデルに対応した要素分割がなされるように入力手段10Aに対して操作入力を行うことで解析手段10Bに設定を行う。最適な構造モデルをどのように決定するかは任意であるが、例えば、最終的に得られたシミュレーション結果(解析結果)と、膜構造物30の変形量の実測データとを比較して決定するなど、従来公知の有限要素法の手法に基づいて行えばよい。
なお、本実施の形態では、膜32として、ゴムと補強材料とを用いているので形状データとしてそれらゴムおよび補強材料の双方の形状データを設定する必要がある。具体的には、補強材料が繊維材料であれば、補強材料を構成する繊維材料の方向や層数などを反映した形状データを設定する必要がある。
また、構造データもゴムおよび補強材料の双方の構造データを含むものとなる。
Next, the deformation simulation method of the membrane structure will be described in detail with reference to the flowcharts of FIGS.
First, when an operator inputs an operation to the input unit 10A, shape data D1 (FIG. 2) indicating the shape of the membrane structure 30 is given to the analysis unit 10B (step S10: shape data input step).
Next, the analysis unit 10B generates structure data by dividing the film structure 30 into elements based on the shape data D1 and setting a plurality of elements and a plurality of nodes (step S12: structure data generation step).
Here, the operator determines an optimum structural model for analyzing the membrane structure 30, and performs an operation input to the input means 10A so that the element division corresponding to the structural model is performed. Set to 10B. How to determine the optimum structural model is arbitrary, but for example, it is determined by comparing the finally obtained simulation result (analysis result) with the actual measurement data of the deformation amount of the membrane structure 30. For example, a conventional finite element method may be used.
In the present embodiment, since rubber and a reinforcing material are used as the film 32, it is necessary to set the shape data of both the rubber and the reinforcing material as the shape data. Specifically, if the reinforcing material is a fiber material, it is necessary to set shape data reflecting the direction and the number of layers of the fiber material constituting the reinforcing material.
The structural data also includes structural data of both rubber and reinforcing material.

次に、操作者は、入力手段10Aに対して操作入力を行うことにより膜32の材料特性を示す材料特性データD2(図2)を解析手段10Bに設定する(ステップS14:材料特性設定ステップ)。
材料特性としては、例えば、ヤング率、ポアソン比、密度(質量密度)などが挙げられる。
本実施の形態では、膜32として、ゴムと補強材料とを用いているので材料特性としてそれらゴムおよび補強材料の双方のデータを設定する。
Next, the operator sets the material property data D2 (FIG. 2) indicating the material properties of the film 32 in the analysis unit 10B by performing an operation input on the input unit 10A (step S14: material property setting step). .
Examples of material properties include Young's modulus, Poisson's ratio, density (mass density), and the like.
In this embodiment, since rubber and a reinforcing material are used as the film 32, data of both the rubber and the reinforcing material are set as material characteristics.

本実施の形態では、材料特性設定ステップは、ステップS16に示す付加質量加算ステップを含んでいる。
付加質量加算ステップは、膜32を構成する材料によって特定される密度に、膜32に接する流体(本実施の形態では海水)の単位面積当たりの質量を付加質量として加算し、加算された値を膜32の密度として設定する(ステップS16)。
In the present embodiment, the material property setting step includes an additional mass addition step shown in step S16.
In the additional mass addition step, the mass per unit area of the fluid (seawater in this embodiment) in contact with the membrane 32 is added to the density specified by the material constituting the membrane 32 as an additional mass, and the added value is obtained. The density of the film 32 is set (step S16).

以下、付加質量の加算方法について詳細に説明する。
すなわち、加算質量設定ステップおける膜32に接する流体は、その質量を容易に計算することができるように、膜32を囲む簡略化された3次元形状にモデル化されている。
本実施の形態では、モデル化された流体の形状が、円柱体(図4(A))、あるいは、円柱体をその軸心を含む平面で2分割した円柱分割体(図4(B))である場合について説明する。
まず、図4(A)の場合について説明すると、前述したように、膜32は軸心Oを中心とする円筒面上を延在し、したがって、断面で見ると、膜32は軸心Oを中心とした半径Rを有する円弧状に延在している。
したがって、膜32によって構成される円弧に対応する中心角をθとし、円弧の両端を結ぶ直線(弦)の長さを2aとすると、aは次の式(1)で示される。
a=Rsin(θ/2)……(1)
モデル化された流体は、長さ2aを直径とし(aを半径とし)、軸心Oと平行に延在する軸心を有する円柱体C0を呈している。
この際、この円柱体C0を軸心O方向から見た状態で円柱体C0の内部に膜32が円弧状に延在している。
この場合、膜32に接する流体の単位面積当たりの質量(円柱体C0の単位長さ当たりの質量)である付加質量m0は(1)式に基づいて(2)式で示されることになる。なお、(2)式は、例えば、機械工学便覧などに記載されている柱状物体が振動する際の付加質量の算出式に基づいている。
m0=ρπa=ρπRsin(θ/2)……(2)
ただしρは流体の密度。
すなわち、図4(A)に図示した状態は、膜32が全て流体中に没した状態であり、したがって、膜32の内面32Aおよび外面32Bの双方の全面に流体が接しており、流体の圧力が内面32Aおよび外面32Bの双方の全面にわたって均一に加わった状態を表していることになる。
そして、後述する有限要素法による膜32の変形挙動の解析計算に際して、膜32の全面にわたって均一に円柱体C0の付加質量m0が付加されたものとして扱う。
すなわち、ステップS16(付加質量加算ステップ)では、膜32の実際の密度に付加質量m0を加えたものを膜32の密度として設定することになる。
Hereinafter, a method for adding the additional mass will be described in detail.
That is, the fluid in contact with the membrane 32 in the additional mass setting step is modeled in a simplified three-dimensional shape surrounding the membrane 32 so that the mass can be easily calculated.
In the present embodiment, the shape of the modeled fluid is a cylindrical body (FIG. 4A), or a cylindrical divided body obtained by dividing the cylindrical body into two planes including its axis (FIG. 4B). The case where it is is demonstrated.
First, the case of FIG. 4A will be described. As described above, the film 32 extends on the cylindrical surface centered on the axis O, and therefore, when viewed in cross section, the film 32 has the axis O. It extends in a circular arc shape having a radius R at the center.
Therefore, when the central angle corresponding to the arc formed by the film 32 is θ and the length of the straight line (string) connecting both ends of the arc is 2a, a is expressed by the following equation (1).
a = Rsin (θ / 2) (1)
The modeled fluid has a cylindrical body C0 having a length 2a as a diameter (a being a radius) and an axis extending parallel to the axis O.
At this time, the film 32 extends in an arc shape inside the cylindrical body C0 when the cylindrical body C0 is viewed from the direction of the axis O.
In this case, the additional mass m0, which is the mass per unit area of the fluid in contact with the membrane 32 (the mass per unit length of the cylindrical body C0), is expressed by the formula (2) based on the formula (1). The formula (2) is based on, for example, a formula for calculating an additional mass when a columnar object vibrates described in a mechanical engineering manual or the like.
m0 = ρπa 2 = ρπR 2 sin 2 (θ / 2) (2)
Where ρ is the density of the fluid.
That is, the state illustrated in FIG. 4A is a state in which the membrane 32 is entirely submerged in the fluid. Therefore, the fluid is in contact with both the inner surface 32A and the outer surface 32B of the membrane 32, and the pressure of the fluid Represents a state in which both are uniformly applied to both the inner surface 32A and the outer surface 32B.
Then, in the analysis calculation of the deformation behavior of the film 32 by the finite element method described later, it is assumed that the additional mass m0 of the cylindrical body C0 is uniformly added over the entire surface of the film 32.
That is, in step S16 (additional mass addition step), the actual density of the film 32 plus the additional mass m0 is set as the density of the film 32.

次に、図4(B)の場合について説明すると、モデル化された流体は、円柱体C0をその軸心を含む平面で2分割した円柱分割体C1を呈している。
この際、分割円柱体C1を円柱体C0の軸心O方向から見た状態で分割円柱体C1の内部に膜32が円弧状に延在するように形成されている。
この場合、膜32に接する流体の単位面積当たりの質量(円柱分割体C1の単位長さ当たりの質量)である付加質量m1は(3)式で示され、(2)式で示されたものの1/2となる。
m1=(1/2)ρπa=(1/2)ρπRsin(θ/2)……(3)
すなわち、図4(B)に図示した状態は、膜32が全て流体の表面から空中に臨んだ状態であり、したがって、膜32の外面32Bは空中に面し膜32の内面32Aの全面のみに流体が接しており、流体の圧力が内面32Aの全面にわたって均一に加わった状態を表していることになる。
そして、この場合は、後述する有限要素法による膜32の変形挙動の解析計算に際して、膜32の全面にわたって均一に円柱分割体C1の付加質量m1が付加されたものとして扱う。
すなわち、ステップS16(付加質量加算ステップ)では、膜32の実際の密度に付加質量m1を加えたものを膜32の密度として設定することになる。
Next, the case of FIG. 4B will be described. The modeled fluid presents a cylindrical divided body C1 obtained by dividing the cylindrical body C0 into two by a plane including its axis.
At this time, the film 32 is formed so as to extend in an arc shape inside the divided cylinder C1 when the divided cylinder C1 is viewed from the direction of the axis O of the cylinder C0.
In this case, the additional mass m1, which is the mass per unit area of the fluid in contact with the membrane 32 (the mass per unit length of the cylindrical divided body C1), is expressed by the equation (3) and expressed by the equation (2). 1/2.
m1 = (1/2) ρπa 2 = (1/2) ρπR 2 sin 2 (θ / 2) (3)
That is, the state shown in FIG. 4B is a state in which the film 32 faces the air from the surface of the fluid. Therefore, the outer surface 32B of the film 32 faces the air and covers only the entire inner surface 32A of the film 32. This means that the fluid is in contact and the pressure of the fluid is uniformly applied over the entire inner surface 32A.
In this case, in the analysis calculation of the deformation behavior of the film 32 by the finite element method described later, the additional mass m1 of the cylindrical divided body C1 is treated as being added uniformly over the entire surface of the film 32.
That is, in step S16 (additional mass addition step), the actual density of the film 32 plus the additional mass m1 is set as the density of the film 32.

次に、操作者は、入力手段10Aに対して操作入力を行うことにより膜構造物30の構造データの各節点を拘束する拘束条件を示す拘束条件データD3(図2)を解析手段10Bに設定する(ステップS18:拘束条件設定ステップ)。
拘束条件は、各節点に対して該節点の持つ各自由度(6自由度)の固定あるいは自由を設定するものであり、従来公知の有限要素法の手法に基づいて設定される。
Next, the operator sets constraint condition data D3 (FIG. 2) indicating the constraint condition for restraining each node of the structure data of the membrane structure 30 to the analysis unit 10B by performing an operation input to the input unit 10A. (Step S18: restriction condition setting step).
The constraint condition is to set a fixed or free degree of freedom (6 degrees of freedom) of each node for each node, and is set based on a conventionally known finite element method.

次に、操作者が入力手段10Aに対して操作入力を行うことにより膜構造物30に作用する圧力を示す圧力データD4(図2)を解析手段10Bに設定する(ステップS20:圧力設定ステップ)。
圧力設定ステップによる圧力データD4の設定は、圧力が加わる膜32の位置と、該位置に加わる圧力の時間的変化とが設定されることでなされ、言い換えると、動的な圧力データが設定される。
圧力データD4をどのように設定するかは、膜構造物30と流体とがどのように動いているかに応じて決定されるものであり、膜32の内面32Aに加わる圧力を内圧、外面32Bに加わる圧力を外圧とした場合、外圧のみを変化させ、内圧を一定としてもよいし、外圧および内圧の双方を変化させてもよい。
Next, pressure data D4 (FIG. 2) indicating the pressure acting on the membrane structure 30 is set in the analysis means 10B when the operator performs an operation input to the input means 10A (step S20: pressure setting step). .
The setting of the pressure data D4 in the pressure setting step is performed by setting the position of the film 32 to which pressure is applied and the temporal change of the pressure applied to the position, in other words, dynamic pressure data is set. .
How to set the pressure data D4 is determined according to how the membrane structure 30 and the fluid are moving. The pressure applied to the inner surface 32A of the membrane 32 is applied to the inner pressure and the outer surface 32B. When the applied pressure is an external pressure, only the external pressure may be changed, the internal pressure may be constant, or both the external pressure and the internal pressure may be changed.

次に、解析手段10Bは、前記構造データ生成ステップで生成された構造データと前記材料特性設定ステップで設定された材料特性データD2と前記拘束条件設定ステップで設定された拘束条件データD3と前記圧力設定ステップで設定された圧力データD4とに基づいて有限要素法による計算を行い、膜構造物30の変形挙動を示すデータD10(図2)を出力手段10Cを介して出力する(ステップS22:出力ステップ)。
本実施の形態では、出力手段10Cから出力されるデータD10は節点の変位量を示す変位量データD10として得られる。
この変位量データD10は、例えば、膜構造物30を3次元的に示した画像などによって表現されるが、変位量データD10をどのような形態で出力し、表現するかは任意である。
なお、本実施の形態では、前記付加質量加算ステップによって付加質量m0が加算され膜32の外面32B全面が流体中に没している状態での解析結果と、付加質量m1が加算され膜32の外面32B全面が空中に臨んでいる状態での解析結果との双方が得られることになる。
したがって、膜32の変形挙動は、その全部が空中に臨む状態と、その全部が流体中に没する状態との2つの状態でのみ解析されることになるが、膜32の外面32Bの一部が空中に臨み、残りが流体中に没している状態は、前記2つの状態の中間の状態として扱うことが可能である。
Next, the analysis means 10B includes the structure data generated in the structure data generation step, the material property data D2 set in the material property setting step, the constraint condition data D3 set in the constraint condition setting step, and the pressure. Calculation based on the finite element method is performed based on the pressure data D4 set in the setting step, and data D10 (FIG. 2) indicating the deformation behavior of the membrane structure 30 is output via the output means 10C (step S22: output). Step).
In the present embodiment, the data D10 output from the output means 10C is obtained as the displacement amount data D10 indicating the displacement amount of the node.
The displacement amount data D10 is expressed by, for example, an image showing the film structure 30 three-dimensionally, but it is arbitrary in what form the displacement amount data D10 is output and expressed.
In the present embodiment, the additional mass m0 is added by the additional mass adding step, and the analysis result in a state where the entire outer surface 32B of the membrane 32 is submerged in the fluid, and the additional mass m1 is added to the membrane 32. Both of the analysis results in a state where the entire outer surface 32B faces the air are obtained.
Therefore, the deformation behavior of the film 32 is analyzed only in two states, a state in which the entire surface is exposed to the air and a state in which the entire surface is immersed in the fluid, but a part of the outer surface 32B of the film 32 is analyzed. The state in which the air enters the air and the rest is submerged in the fluid can be treated as an intermediate state between the two states.

そして、得られた解析結果(膜構造物30の変形挙動を示すデータ)に基づいて、流体から圧力を受けて膜32が変形した場合における膜32とソナーとの距離の接近度合い、干渉の有無などを確認することにより、膜構造物32の設計評価を行い(ステップS24:評価ステップ)、一連の処理を終了する。   Based on the obtained analysis results (data indicating the deformation behavior of the membrane structure 30), the degree of proximity of the distance between the membrane 32 and the sonar when the membrane 32 is deformed by receiving pressure from the fluid, and the presence or absence of interference By confirming the above, the design evaluation of the film structure 32 is performed (step S24: evaluation step), and the series of processes is terminated.

以上説明したように、本実施の形態によれば、膜構造物30の膜32を構成する材料によって特定される密度に、膜32に接する流体の単位面積当たりの質量を付加質量として加算し、加算された値を膜32の密度として設定するようにしたので、有限要素法により膜構造物30の変形挙動を解析するに際して、膜32に接する流体の圧力の影響が反映されるため、流体中における膜構造物30の変形挙動を正確に再現する上で有利となり、流体解析を必要としないのでコストの抑制を図る上でも有利となる。
また、有限要素法により膜構造物30の変形挙動を解析するに際して、膜32に接する流体の単位面積当たりの質量を付加質量として加算するだけでよいため、膜構造物30の形状や大きさを変化させながら繰り返して解析を行ったとしても解析に要する手間や時間がかからないため、設計の効率化を図る上で極めて有利となる。
As described above, according to the present embodiment, the mass per unit area of the fluid in contact with the membrane 32 is added as an additional mass to the density specified by the material constituting the membrane 32 of the membrane structure 30. Since the added value is set as the density of the membrane 32, the influence of the pressure of the fluid in contact with the membrane 32 is reflected when analyzing the deformation behavior of the membrane structure 30 by the finite element method. This is advantageous in accurately reproducing the deformation behavior of the membrane structure 30 in the above, and is not advantageous in reducing costs because fluid analysis is not required.
Further, when analyzing the deformation behavior of the membrane structure 30 by the finite element method, it is only necessary to add the mass per unit area of the fluid in contact with the membrane 32 as an additional mass, so that the shape and size of the membrane structure 30 can be changed. Even if the analysis is repeatedly performed while changing, it does not take time and effort required for the analysis, which is extremely advantageous in improving the efficiency of the design.

なお、本実施の形態では、膜構造物30が船体に設けられたソナーを保護するカバーを構成するものである場合について説明したが、本発明方法によって解析される膜構造物の形態や構成は任意である。
また、本実施の形態では、流体が海水や水である場合について説明したが、流体は水に限定されるものではなく、さまざまな液体あるいは気体であってもよい。
したがって、膜構造物としては、袋状を呈し内部に流体を蓄える膜構造物、例えば、空気を蓄えて港の岸壁と船舶の船体との衝突を防止する防舷材、あるいは、真水や液体燃料などの流体を蓄え船舶によって海中を曳航される袋体、あるいは、水が満たされたクッション、あるいは、気体(ガス)によって膨張するエアバッグなどが挙げられる。
また、膜構造物は袋状に構成されている必要は無く、例えば、衝立のように板状を呈していてもよい。
In the present embodiment, the case where the membrane structure 30 constitutes a cover for protecting the sonar provided on the hull has been described. However, the form and configuration of the membrane structure analyzed by the method of the present invention are as follows. Is optional.
Moreover, although the case where the fluid is seawater or water has been described in the present embodiment, the fluid is not limited to water, and may be various liquids or gases.
Therefore, as the membrane structure, a membrane structure that has a bag shape and stores fluid therein, for example, a fender that stores air and prevents collision between the harbor quay and the ship's hull, or fresh water or liquid fuel For example, a bag body that stores a fluid such as towed in the sea by a ship, a cushion filled with water, an airbag that inflates with gas (gas), and the like.
Further, the membrane structure does not have to be configured in a bag shape, and may have a plate shape like a partition, for example.

また、本実施の形態では、膜構造物の膜を構成する材料としてゴムおよび補強材料を用いた場合について説明したが、材料は任意である。
例えば、船舶の船底部分を構成する鋼板などで形成された壁部を膜構造物として扱い、本発明方法によって船底部分の変形挙動を解析することも可能である。
In the present embodiment, the case where rubber and a reinforcing material are used as the material constituting the membrane of the membrane structure has been described, but the material is arbitrary.
For example, it is also possible to treat a wall portion formed of a steel plate or the like constituting a ship bottom portion of a ship as a membrane structure and analyze the deformation behavior of the ship bottom portion by the method of the present invention.

また、本実施の形態では、膜構造物30の膜32の外面32Bが流体中に没した状態と、空気中に臨んだ状態との2つの状態に対応してそれぞれ付加質量を求めて膜32の密度に加算するようにしたが、これら2つの状態の中間の状態、すなわち、膜32の外面32Bの一部が流体中に没し、残りが空気中に臨んだ状態、言い換えると、流体が膜32の外面32Bの全体ではなく一部に接した状態に対応して付加重量を求めて膜32の密度に加算して解析を行うようにしてもよい。
さらには、この中間の状態を時間経過に応じて、すなわち、膜32が空中から流体中に移動する過程、あるいは、膜32が流体中から空中へ移動する過程を複数段階に分けて付加重量を求め、求められた複数段階の付加重量を膜32の密度にそれぞれ加算して解析を行うようにしてもよい。
このようにすると、膜32に接する流体が膜32に与える影響の時間的な変化を解析結果に正確に反映する上で有利となり、膜構造物30の変形挙動の再現をより正確に行う上で有利となる。
Further, in the present embodiment, the additional mass is obtained for each of the two states of the state where the outer surface 32B of the membrane 32 of the membrane structure 30 is submerged in the fluid and the state where the outer surface 32B is exposed to the air. However, in the state where the outer surface 32B of the membrane 32 is partially submerged in the fluid and the rest is exposed to the air, in other words, the fluid The analysis may be performed by obtaining the additional weight corresponding to a state in which the membrane 32 is in contact with a part rather than the entire outer surface 32B and adding it to the density of the membrane 32.
Furthermore, this intermediate state is increased over time, that is, the process in which the membrane 32 moves from the air to the fluid or the process in which the membrane 32 moves from the fluid to the air is divided into a plurality of stages to increase the additional weight. The obtained additional weights in a plurality of stages may be added to the density of the membrane 32 for analysis.
This is advantageous in accurately reflecting the temporal change in the influence of the fluid in contact with the film 32 on the analysis result 32 in the analysis result, and in reproducing the deformation behavior of the film structure 30 more accurately. It will be advantageous.

なお、本実施の形態では、膜32が円筒面上を延在し、膜32に接する流体が膜32を囲む3次元形状にモデル化され、モデル化された流体の3次元形状が、前記円筒面の軸心Oと平行に延在する軸心Oを有する円柱体C0を呈し、かつ、円柱体C0の軸心方向から見た状態で円柱体C0の内部に膜32が延在するように形成されている場合について説明した。
また、モデル化された流体の3次元形状は前記円筒面の軸心と平行に延在する軸心を有する円柱体をその軸心を含む平面で2分割した円柱分割体C1を呈し、かつ、分割円柱体C1を前記軸心方向から見た状態で円柱分割体C1の内部に膜32が延在するように形成されている場合について説明した。
すなわち、円柱体C0および円柱分割体C1の軸心の両端は、前記軸心と直交する平面で構成されている。
In the present embodiment, the membrane 32 extends on the cylindrical surface, the fluid in contact with the membrane 32 is modeled into a three-dimensional shape surrounding the membrane 32, and the modeled three-dimensional shape of the fluid is the cylinder. A cylindrical body C0 having an axial center O extending in parallel with the axial center O of the surface is presented, and the film 32 extends inside the cylindrical body C0 as viewed from the axial center direction of the cylindrical body C0. The case where it is formed has been described.
The modeled three-dimensional shape of the fluid presents a cylinder divided body C1 obtained by dividing a cylinder having an axis extending parallel to the axis of the cylindrical surface into two by a plane including the axis, and The case where the film 32 is formed so as to extend inside the columnar divided body C1 in a state where the divided columnar body C1 is viewed from the axial direction has been described.
That is, both ends of the axial center of the cylindrical body C0 and the cylindrical divided body C1 are configured by planes orthogonal to the axial center.

図5(A)は膜32および円柱分割体C1の一例を示す平面図、(B)は(A)のB矢視図、(C)は(B)のC矢視図である。なお、図面の簡略化を図るために、膜32は平面として描かれている。
図5に示すように、モデル化された円柱分割体C1によって膜32に加わる長さL方向における荷重の分布は長さL方向において均一となる。
しかしながら、実際には、膜32に加わる荷重のうち、膜32の長さL方向の両端近傍箇所に加わる荷重は、両端近傍を除く箇所に加わる荷重に比較して軽減することが想定され、言い換えると、流体の長さL方向における荷重の分布は長さL方向の両端に近づくにつれて次第に低下する分布となる。
例えば、膜32の長さLと幅Wが、例えば、L>10Wである場合には、膜32の長さ方向の両端近傍を除く箇所が占める面積が大半となるため、荷重分布が長さL方向において一様であるとして扱ってもよい。
しかしながら、幅Wに対して長さLがあまり小さくない場合、例えば、L>2Wといった場合には、膜32の長さ方向の両端近傍を除く箇所が占める面積に対して膜32の長さ方向の両端近傍の面積が占める割合が大きくなることから、荷重分布が長さL方向の両端に近づくにつれて次第に低下することを考慮する必要がある。
このような場合には、円柱分割体C1の軸心の両端を平面として扱うのではなく、流体の形状を、荷重分布が長さL方向の両端に近づくにつれて次第に低下する分布となるようにモデル化することが流体によって膜32に加わる荷重の分布を正確に反映して解析結果(膜構造物30の変形挙動を示すデータ)を得る上で好ましい。例えば、図5に示すように、円柱分割体C1の軸心の両端を平面ではなく球面としてモデル化することが好ましい。
以上、流体の形状を円柱分割体C1としてモデル化した場合について説明したが、流体の形状を円柱体C0としてモデル化した場合についても同様である。
FIG. 5A is a plan view showing an example of the membrane 32 and the cylindrical divided body C1, FIG. 5B is a view taken along the arrow B in FIG. 5A, and FIG. 5C is a view taken along the arrow C in FIG. In order to simplify the drawing, the film 32 is drawn as a plane.
As shown in FIG. 5, the distribution of the load in the length L direction applied to the film 32 by the modeled cylindrical divided body C1 is uniform in the length L direction.
However, in actuality, among the loads applied to the film 32, it is assumed that the load applied to the vicinity of both ends in the length L direction of the film 32 is reduced compared to the load applied to the positions excluding the vicinity of both ends. Then, the distribution of the load in the length L direction of the fluid is a distribution that gradually decreases as it approaches both ends in the length L direction.
For example, when the length L and the width W of the film 32 are, for example, L> 10 W, the area occupied by the portions excluding the vicinity of both ends in the length direction of the film 32 becomes most, so the load distribution is long. It may be treated as being uniform in the L direction.
However, when the length L is not so small with respect to the width W, for example, when L> 2W, the length direction of the film 32 with respect to the area occupied by the portions excluding the vicinity of both ends in the length direction of the film 32 It is necessary to consider that the load distribution gradually decreases as it approaches both ends in the length L direction because the area occupied by the area in the vicinity of both ends increases.
In such a case, instead of treating both ends of the axis of the cylindrical divided body C1 as a plane, the model of the fluid is a distribution that gradually decreases as the load distribution approaches both ends in the length L direction. It is preferable to obtain the analysis result (data indicating the deformation behavior of the membrane structure 30) by accurately reflecting the distribution of the load applied to the membrane 32 by the fluid. For example, as shown in FIG. 5, it is preferable to model both ends of the axial center of the cylindrical divided body C1 as a spherical surface instead of a plane.
Although the case where the shape of the fluid is modeled as the cylindrical divided body C1 has been described above, the same applies to the case where the shape of the fluid is modeled as the cylindrical body C0.

また、流体の3次元形状は円柱体C0あるいは円柱分割体C1に限定されるものではなく、膜32の形状に対応したさまざまな形状を用いることができる。
図6(A)は膜32および半球体C10の他の例を示す平面図、(B)は(A)のB矢視図、(C)は(B)のC矢視図である。
例えば、図6に示すように、膜32の長さLと高さWが等しい場合には、流体の荷重分布が長さL方向および高さW方向の両端に近づくにつれて低下することを考慮して、流体を球体、あるいは、半球体C10としてモデル化することが流体によって膜32に加わる荷重の分布を正確に反映して解析結果(膜構造物30の変形挙動を示すデータ)を得る上で好ましい。
また、膜32の形状は任意であり、膜32の輪郭(膜32を平面視したときの形状)は矩形に限定されるものではなく、三角形、多角形、円形、楕円形などさまざまな形状が採用可能であることは無論である。また、膜32は、円筒面上を延在する形状に限定されるものではなく、例えば、平面状を延在しても、あるいは、球面上を延在してもよい。
Further, the three-dimensional shape of the fluid is not limited to the cylindrical body C0 or the cylindrical divided body C1, and various shapes corresponding to the shape of the film 32 can be used.
6A is a plan view showing another example of the membrane 32 and the hemisphere C10, FIG. 6B is a view taken in the direction of arrow B in FIG. 6A, and FIG. 6C is a view taken in the direction of arrow C in FIG.
For example, as shown in FIG. 6, when the length L and the height W of the membrane 32 are equal, it is considered that the fluid load distribution decreases as it approaches both ends in the length L direction and the height W direction. Therefore, modeling the fluid as a sphere or hemisphere C10 accurately reflects the distribution of the load applied to the membrane 32 by the fluid to obtain an analysis result (data indicating the deformation behavior of the membrane structure 30). preferable.
The shape of the film 32 is arbitrary, and the outline of the film 32 (the shape when the film 32 is viewed in plan) is not limited to a rectangle, and various shapes such as a triangle, a polygon, a circle, and an ellipse can be used. Of course, it can be adopted. Further, the film 32 is not limited to the shape extending on the cylindrical surface, and may be, for example, a flat shape or a spherical surface.

本発明方法を実行するために使用されるコンピュータ10の構成を示すブロック図である。It is a block diagram which shows the structure of the computer 10 used in order to perform the method of this invention. コンピュータ10の機能ブロック図である。2 is a functional block diagram of a computer 10. FIG. 本実施の形態の膜構造物30の変形シミュレーション方法を示すフローチャートチャートである。It is a flowchart chart which shows the deformation | transformation simulation method of the film structure 30 of this Embodiment. (A)は膜構造物30の外面が海面上に臨んだ状態を模式的に示す断面説明図、(B)は膜構造物30の外面が海面下に水没した状態を模式的に示す断面説明図である。(A) is a cross-sectional explanatory view schematically showing a state in which the outer surface of the membrane structure 30 faces the sea surface, and (B) is a cross-sectional explanation schematically showing a state in which the outer surface of the membrane structure 30 is submerged under the sea surface. FIG. (A)は膜32および円柱分割体C1の一例を示す平面図、(B)は(A)のB矢視図、(C)は(B)のC矢視図である。(A) is a top view which shows an example of the film | membrane 32 and the cylindrical division body C1, (B) is a B arrow directional view of (A), (C) is a C arrow directional view of (B). (A)は膜32および半球体C10の一例を示す平面図、(B)は(A)のB矢視図、(C)は(B)のC矢視図である。(A) is a top view which shows an example of the film | membrane 32 and hemisphere C10, (B) is B arrow view of (A), (C) is C arrow view of (B).

符号の説明Explanation of symbols

10……コンピュータ、10A……入力手段、10B……解析手段、10C……出力手段、30……膜構造物、32……膜、D1……形状データ、D2……材料特性データ、D3……拘束条件データ、D4……圧力データ、m0、m1……付加質量。   10: Computer, 10A: Input means, 10B: Analysis means, 10C: Output means, 30: Membrane structure, 32: Membrane, D1: Shape data, D2: Material property data, D3: ... Restriction condition data, D4 ... Pressure data, m0, m1 ... Additional mass.

Claims (10)

流体中において膜で構成された膜構造物に圧力が作用した際に発生する前記膜構造物の変形挙動を有限要素法によりコンピュータを用いて解析する膜構造物の変形シミュレーション方法であって、
前記コンピュータは、操作者による操作入力を受け付ける入力手段と、有限要素法による解析処理を行う解析手段と、データの出力を行う出力手段とを備え、
操作者が前記入力手段に対して操作入力を行うことにより前記膜構造物の形状を示す形状データを前記解析手段に与える形状データ入力ステップと、
前記解析手段が前記形状データに基づいて前記膜構造物を要素分割し複数の要素と複数の節点を設定することで構造データを生成する構造データ生成ステップと、
操作者が前記入力手段に対して操作入力を行うことにより前記膜の材料特性を示す材料特性データを前記解析手段に設定する材料特性設定ステップと、
操作者が前記入力手段に対して操作入力を行うことにより前記膜構造物の拘束条件を示す拘束条件データを前記解析手段に設定する拘束条件設定ステップと、
操作者が前記入力手段に対して操作入力を行うことにより前記膜構造物に作用する圧力を示す圧力データを前記解析手段に設定する圧力設定ステップと、
前記解析手段が、前記構造データ生成ステップで生成された構造データと前記材料特性設定ステップで設定された材料特性データと前記拘束条件設定ステップで設定された拘束条件データと前記圧力設定ステップで設定された圧力データとに基づいて有限要素法による計算を行い前記膜構造物の変形挙動を示すデータを前記出力手段を介して出力する出力ステップとを含み、
前記材料特性設定ステップによって設定される前記材料特性データは前記膜の密度を含み、
前記材料特性設定ステップは、前記膜を構成する材料によって特定される密度に、前記膜に接する流体の単位面積当たりの質量を付加質量として加算した値を前記膜の密度として設定する付加質量加算ステップを含む、
ことを特徴とする膜構造物の変形シミュレーション方法。
A deformation simulation method for a membrane structure, wherein a deformation behavior of the membrane structure generated when a pressure is applied to a membrane structure composed of a membrane in a fluid is analyzed using a computer by a finite element method,
The computer includes input means for receiving an operation input by an operator, analysis means for performing analysis processing by a finite element method, and output means for outputting data,
A shape data input step for giving shape data indicating the shape of the membrane structure to the analysis means by an operator performing an operation input on the input means;
A structure data generating step for generating structure data by dividing the film structure into elements based on the shape data and setting a plurality of elements and a plurality of nodes;
A material property setting step of setting material property data indicating the material property of the film in the analysis unit by performing an operation input on the input unit by an operator;
A constraint condition setting step of setting constraint condition data indicating a constraint condition of the film structure in the analysis means by an operator performing an operation input on the input means;
A pressure setting step for setting pressure data indicating pressure acting on the membrane structure in the analysis means by an operator performing an operation input on the input means;
The analysis means is set in the structure data generated in the structure data generation step, the material property data set in the material property setting step, the constraint condition data set in the constraint condition setting step, and the pressure setting step. An output step of performing calculation by a finite element method based on the pressure data and outputting data indicating the deformation behavior of the membrane structure via the output means,
The material property data set by the material property setting step includes the density of the film,
In the material property setting step, an additional mass adding step of setting, as a density of the film, a value obtained by adding a mass per unit area of a fluid in contact with the film as an additional mass to a density specified by a material constituting the film including,
A method for simulating deformation of a membrane structure characterized by the above.
前記材料特性設定ステップによる前記膜を構成する材料によって特定される密度に対する前記付加質量の加算は、前記付加質量が前記流体に接する前記膜の全面にわたって均一となるように加算される、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
The addition of the added mass to the density specified by the material constituting the film in the material property setting step is added so that the added mass is uniform over the entire surface of the film in contact with the fluid.
The deformation simulation method for a membrane structure according to claim 1.
前記加算質量設定ステップおける前記膜に接する流体は前記膜を囲む3次元形状にモデル化されている、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
The fluid in contact with the film in the additional mass setting step is modeled in a three-dimensional shape surrounding the film.
The deformation simulation method for a membrane structure according to claim 1.
前記膜は円筒面上を延在し、
前記3次元形状は、前記円筒面の軸心と平行に延在する軸心を有する円柱体を呈し、かつ、前記円柱体の軸心方向から見た状態で前記円柱体の内部に前記膜が延在するように形成されている、
ことを特徴とする請求項3記載の膜構造物の変形シミュレーション方法。
The membrane extends over a cylindrical surface;
The three-dimensional shape presents a cylindrical body having an axial center extending in parallel with the axial center of the cylindrical surface, and the film is formed inside the cylindrical body when viewed from the axial direction of the cylindrical body. Formed to extend,
The deformation simulation method for a membrane structure according to claim 3.
前記膜は円筒面上を延在し、
前記3次元形状は前記円筒面の軸心と平行に延在する軸心を有する円柱体をその軸心を含む平面で2分割した円柱分割体を呈し、かつ、前記分割円柱体を前記軸心方向から見た状態で前記円柱分割体の内部に前記膜が延在するように形成されている、
ことを特徴とする請求項3記載の膜構造物の変形シミュレーション方法。
The membrane extends over a cylindrical surface;
The three-dimensional shape presents a cylindrical divided body obtained by dividing a cylindrical body having an axis extending in parallel with the axis of the cylindrical surface into two by a plane including the axis, and the divided cylindrical body is the axis. It is formed so that the film extends inside the cylindrical divided body in a state seen from the direction,
The deformation simulation method for a membrane structure according to claim 3.
前記出力ステップによって出力される前記膜構造物の変形挙動を示すデータは、前記複数の節点の変位量で示される、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
Data indicating the deformation behavior of the membrane structure output by the output step is indicated by displacement amounts of the plurality of nodes.
The deformation simulation method for a membrane structure according to claim 1.
前記圧力設定ステップによる圧力データの設定は、圧力が加わる前記膜の位置と、該位置に加わる圧力の時間的変化とが設定されることでなされる、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
The setting of the pressure data by the pressure setting step is performed by setting the position of the film to which pressure is applied and the temporal change of the pressure applied to the position.
The deformation simulation method for a membrane structure according to claim 1.
前記膜構造物を構成する膜の2つの面のうち一方の面が閉塞された内側空間に臨む内面を構成し、他方の面が前記膜を挟んで前記内側空間の反対側に位置する外側空間に臨む外面を構成し、
前記内側空間に流体が満たされることで前記膜構造物の形状が保たれている、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
An outer space that constitutes an inner surface facing an inner space where one of the two surfaces of the membrane constituting the membrane structure is closed, and the other surface is located on the opposite side of the inner space across the membrane Configure the outer surface facing
The shape of the membrane structure is maintained by filling the inner space with fluid,
The deformation simulation method for a membrane structure according to claim 1.
前記膜を構成する材料はゴムを含む、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
The material constituting the film includes rubber,
The deformation simulation method for a membrane structure according to claim 1.
前記膜を構成する材料はゴムと補強材料とを含む、
ことを特徴とする請求項1記載の膜構造物の変形シミュレーション方法。
The material constituting the membrane includes rubber and a reinforcing material,
The deformation simulation method for a membrane structure according to claim 1.
JP2007255902A 2007-09-28 2007-09-28 Deformation simulation method of membrane structure Expired - Fee Related JP4936135B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007255902A JP4936135B2 (en) 2007-09-28 2007-09-28 Deformation simulation method of membrane structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007255902A JP4936135B2 (en) 2007-09-28 2007-09-28 Deformation simulation method of membrane structure

Publications (2)

Publication Number Publication Date
JP2009087025A JP2009087025A (en) 2009-04-23
JP4936135B2 true JP4936135B2 (en) 2012-05-23

Family

ID=40660365

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007255902A Expired - Fee Related JP4936135B2 (en) 2007-09-28 2007-09-28 Deformation simulation method of membrane structure

Country Status (1)

Country Link
JP (1) JP4936135B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175567B (en) * 2011-01-30 2013-02-06 长江勘测规划设计研究有限责任公司 Method for acquiring density of rock fill object by using additive mass method and digital template chart
CN108985003A (en) * 2018-06-28 2018-12-11 东汉新能源汽车技术有限公司 The performance of operating condition parameter acquiring method and device of front shroud
CN108872008B (en) * 2018-07-16 2020-05-12 长江地球物理探测(武汉)有限公司 Rockfill density measurement method based on additional mass method theoretical gauge plate
CN113505430B (en) * 2021-06-07 2024-02-09 清华大学 Large-scale space flexible film dynamic parameter calculation method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4684425B2 (en) * 2001-01-26 2011-05-18 株式会社ブリヂストン Tire design method
JP4359101B2 (en) * 2003-08-15 2009-11-04 独立行政法人科学技術振興機構 Inflatable structure form calculation method, calculation program, and computer-readable recording medium storing calculation program
JP4604477B2 (en) * 2003-11-18 2011-01-05 横浜ゴム株式会社 Method for estimating physical properties, method for optimal design of sound transmission body

Also Published As

Publication number Publication date
JP2009087025A (en) 2009-04-23

Similar Documents

Publication Publication Date Title
Gong et al. Added resistance and seakeeping performance of trimarans in oblique waves
Newman et al. Boundary-element methods in offshore structure analysis
McVicar et al. Fluid–structure interaction simulation of slam-induced bending in large high-speed wave-piercing catamarans
JP4936135B2 (en) Deformation simulation method of membrane structure
KR102543096B1 (en) Analysis method and apparatus for optimizing vibration performance of automotive body
Yu et al. Implementation of linear potential-flow theory in the 6DOF coupled simulation of ship collision and grounding accidents
Yucel et al. Free and forced vibration analyses of ship structures using the finite element method
Das et al. Hydroelasticity of marine vessels advancing in a seaway
Maxit et al. Improving the Patch Transfer Function approach for fluid-structure modelling in heavy fluid
Siddiqui et al. Experimental studies of a damaged ship section in forced heave motion
Yildiz et al. A numerical and experimental study on the outrigger positioning of a trimaran hull in terms of resistance
Riesner et al. A numerical method to compute global resonant vibrations of ships at forward speed in oblique waves
Zhang et al. A numerical investigation on the effect of symmetric and asymmetric flooding on the damage stability of a ship
Suneela et al. Numerical investigation of interceptor effect on seakeeping behaviour of planing hull advancing in regular head waves
Miller et al. Effect of amplitudes and frequencies on Virtual Planar Motion Mechanism of AUVs, Part I: Forces, moments and hydrodynamic derivatives
Santos et al. Modal analysis of a fast patrol boat made of composite material
Tang et al. A fully nonlinear BEM-beam coupled solver for fluid–structure interactions of flexible ships in waves
Piro A hydroelastic method for the analysis of global ship response due to slamming events
Ni et al. Multiple speed integrated optimization design for a SWATH using SBD technique
Adenya et al. Estimation of springing response for 550 000 DWT ore carrier
Sepehrirahnama et al. Fluid–structure interaction effects on free vibration of containerships
Lakshmynarayanana et al. Fully coupled CFD/FEA investigations to predict the wave loads on a flexible containership
Benhamou et al. Hydroelastic simulations in openfoam®: A case study on a 4400teu containership
Homm et al. Experimental and numerical investigation of a complex submerged structure. Part I: modal analysis
Li et al. Vibro-acoustic responses of a hull due to structural and acoustic excitations from a propeller

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100811

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120124

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120127

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120209

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150302

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150302

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees