CalculiX CrunchiX: ccx(http://www.calculix.de/)のマニュアルを 英語ゼミで学生に和訳してもらう作業ページ。

Introduction.

(藤原) 序論。 これは、Calculix CrunchiXの説明書です。プログラムを使用していて何か問題が発生したとき、この書類が解決してくれるでしょう。解決しない場合はE-mailを送ってください。(dhondt@.....)次の章はあなたの興味をひくような簡単な例題を含んでいます。3章は、解析タイプ,材料、要素について何らかの基礎知識を与える理論上の章です。それから、概要は入力デッキの形式の細かい指示を受けて、有効なありとあらゆるアルファベット順のキーワードを施しています。これは、入力デッキに問題があり、CalculiXが実行しないときに見るための章です。次章はプログラムのコードと同様に入手すべきだった確認例の説明を収録しています。仮に、これまでに扱われていない新しい種類の問題を解決しようとするなら、この章のふさわしい例を見て下さい。また、この章でCalculiXが正しくインストールされているかどうか、調べることとしても使えます。そして、User's Manualの最後には、コードを書いているときに用いられた参考資料がいくつか付いています。 このマニュアルは有限要素のテキストではありません。いかにも、有限要素法の実用的な知識は見せかけです。有限要素法をよく知らない方々のために、ZienkiewiczとTaylor[34]による工業技術志向の生徒のための本と、Hughes[13]とDhondtによる数学的視野を持った読者のための出版物をお勧めします。

後藤ちゃちゃ

メールとかはロボットに収集されてスパムリストに登録される 恐れがあるので、ウェブ上にはさらさない方がいいでしょう。 「欲心をうるおす」ってどういう意味? できるだけ、普通の日本語で書いて。 あ、なるほど。wet your appetiteって書かれると確かに迷うけど、whetに置き換えれば、 「あなたの興味を引くような簡単な例題」ってことでしょ、たぶん。 この手の話題の場合、analysisは解析でしょう。 onは、要素や材料「について」のってことかな。 最近は日本語の本でもセクション3とか書いてあるかも知れないけど、 3章とか3節とかが一般的かな。 sectionの上の階層にchapterとかがあるなら、chapterを章にして、 sectionを節にするのがいいかも知れないけど、 どうも、sectionが一番上の階層みたいだから、sectionを章にしてしまおうか。 で、その下の階層は3.1節とか。

Simple example problems

(藤原)この章では、自由端片持ち梁に点力を負荷したときの解析と、炉での熱量計算を実行します。 片持ち梁の荷重と境界条件を図1に示しています。梁の大きさは1×1×8㎥で、9×10^6Nの点力からなる荷重を負荷させ、左側先端部は完全に(四方八方に)固定されている梁です。

Theory(理論)

(藤原)有限要素法は基本的にフィールド変数の決定に関係しています。CalculiXの変形の尺度の基本として、ラグランジェ変形テンソルEは弾性媒体として利用され、オイラーの変形テンソルeは変形塑性として利用され、コーシーの偏差弾性の左辺のテンソルは増分塑性として利用されています。

Element Types

Eight-node brick element (C3D8)

(三浦)これは一般的な線形のブリック要素と、完全に積分しています(2x2x2積分点)。17で形の関数を見つけることができます。節点付番は図12の仕様に続きます、そして、図13によると、積分点は付番されます。この後者の情報は、積分点で、EL PRINTキーワードで印刷された要素変数を与えるので、重要です。 要素の構造は簡単ですが、以下の状況でそれを使用すべきではありません。

・完全な積分のため、要素はすなわち定積材料作用、ポアソン係数かプラスチックの作用の高い値のために、ひどく反応するでしょう。

・要素は、例えば、曲がりの下の細長いけたか薄板へ曲がるのにおいて、不自然な傾向があります。

後藤ちゃちゃ

おっと、integrateってどういう意味よ。辞書を調べる時は、数学とか物理とか工学の 用語がないかどうかをまず見ること。integration pointって辺りで、なんかそれ方面の 用語かなあと思わないと。 linearは、まあ直線のということだけど、f=kxのfとxは直線の関係だって言いますか。 これも数学用語。brickはそのまま「れんが要素」とか「ブリック要素」で いいかな(ググって出てくるのは後者だけど)。だって、六面体って直方体には 限らないでしょ(意訳するなら、直方体要素の方がまだいいかも)。

後藤ちゃちゃ

積分ポイントでもいいとは思うけど、積分点と言う方が普通。 積分点については、 &link(ここの岩崎先生の講義ノート,http://comp.nagaokaut.ac.jp/~iwas/lecture/lecture_m1.html)の解説が分かりやすい。

Boundary conditions

Single point constraints (SPC)

(三浦)ただ一つの制約ポイントでは、自由度1以上は与えられた節点のために調節されています。指示された値は、ゼロか非零であるかもしれません。段の外で非零SPCのものを定義することができません。段の中、または、段の外でゼロSPCのものを定義することができます。SPCのものはキーワード*BOUNDARY(境界)と共に定義されます。自由の機械的な度合が表示1から6までです。熱自由度は11と定義されます。要素だけを発して、骨組みと角材に回転の度の自由を適用することができます。

Multiple point constraints (MPC)

(三浦)複数の制約ポイントが一つ以上の節点における、自由の間の関係を確立します。CalauliXだけでは、均質の直線的な関係は考えられます。第一歩の前のキーワード*EQUATIONと共にMPCのものを定義しなければなりません。異質の用語をダミーの節点(SPCを使用する)の自由(DOF)の度合の一つに割り当てて、MPCでこのDOFを含めることによって、異質の直線的な関係を定義することができます。その結果、それを均質化します。DOFの付番はSPCのもの(CF前のセクション)のように同じです。それは全く同じMPCの中で熱的、そして、機械的な度の自由を混ぜることができません。

Materials

(三浦)物質的な定義は、ELASTIC、EXPANSION、DENSITY、HYPERELASTIC、HYPERFOAM、DEFORMATIONPLASTICITY、PLASTIC、CREEP、またはUSER MATERIALなどの物質的な特定のカードが支えたMATERIAL主要なカードから始まります。材料を要素に割り当てるために、SOLID SECTIONカードは使用されています。要素は一つの材料だけから成ることができます。構造の各要素で、材料を割り当てなければなりません。荷重のいくつかのタイプが特定有形物質を必要とします:重力荷重は材料の密度を必要として、温度荷重は熱膨張係数を必要とします。また、分析のタイプは有形物質を必要とすることができなす:周波数解析は材料の密度を必要とします。いくつかの物質的なカードが互いに排他的ですが、他のものは互いに依存しています。ちょうど以下の一つが必要です:ゴムヒモ、HYPERELASTIC、HYPERFOAM、変形化塑性、およびユーザー材料。ELASTIC(TYPE=ISO)はキーワードPLASTICに先行しなければなりません。同じくらいはCREEPカードに適用されます。ELASTICとCREEPカードの間のPLASTICカードは粘塑性変形の材料を定義します。あなたの必要性に従って、他のキーワードを使用することができます。入力デッキで定義された材料のいずれかが直線的な弾性の材料でないから、幾何学上非線形性は自動的に考慮に入れられます。(すなわち、NLGEOMは活性です)

Linear elastic materials

(三浦)直線的な弾性の材料は緊張における二次の用語だけが守られる弾性の可能性によって特徴付けられます。等方性の、または、直方異方性の、または、完全に異方性のバージョンでそれを定義することができます。等方性の直線的な弾性の材料はそれらのヤング率とポアソンの係数によって特徴付けられます。通常、一般的な鉄鋼は等方性です。木製の、または、立方体の単一の結晶などの直交異方性材料は21の定数による9つの非零定数と完全に異方性の材料によって特徴付けられます。弾性の材料において、キーワードELASTICは使用されています。

Types of analysis

(三浦)分析タイプは完全なステップに適用します。(それはSTEPカードから始まって、ENDSTEPカードで終わります。)中間で分析タイプ、荷重、および分野出力要求を定義しなければなりません。

Static analysis

(三浦)静態分析には、時間の寸法はかかわりません。すなわち、ずっと、そして慣性効果を無視することができるくらい遅い準静電気の中で荷重が適用されると思われます。静態分析はキーワードSTATICによって定義されます。静的な方法は幾何学上直線的であるか、または非線形である場合があります。どちらの場合もラグランジアン観点を取ります、そして、物質的な基準系11ですべての変数を指定します。したがって、CalculiXで内部的に使用される圧力は非変形している表面に影響する2番目のPiolaーキルヒホフ張筋です。幾何学上直線的な計算において、微小の緊張を取ります、そして、(ラグランジアン緊張張筋のバージョンを線形化した)負荷は互いに妨げません。したがって、二つの異なった負荷による変形はそれぞれのそれらによる変形の合計です。直線的な計算において、コーシーとPiolaーキルヒホフ圧力の違いは無視されます。幾何学上非線形の計算において、完全なラグランジアン緊張張筋は使用されています。(計算はSTEPカードの上の変数NLGEOMによって引き起こされます。)また、変数NLGEOMが使用されていて、非線形の材質振る舞い(例えばHYPERELASTIC(PLASTIC))で幾何学上非線形の計算は自動的に引き起こされます。通常、方法は増分に分割されます、そして、ユーザーはSTATICカードの初期の増分の長さと総ステップの長さを提供するべきです。増分の長さは、固定されているか(STATICカードの上の変数DIRECT)、または自動である場合があります。自動増加の場合には、問題の集合の特性に応じて、増分の長さは自動的に調整されます。集合に達するまで増分、プログラムが繰り返すか、または新しい試みが、よりわずがな増分量でそれぞれでされます。各繰り返しでは、幾何学上直線的な剛性マトリクスは最後の繰り返しにおける変形による初期の置き換えの堅さと最後の繰り返しの圧力34による初応力の堅さで増大します。ファイルにおける出力において、2番目のpiola-キルヒホフ圧力はコーシーか真応力に変換されます、これが本当に構造影響している圧力である。特別条項は周期的な左右対称の構造に作られています。周期的な左右対称の構造はN同じ分野によって特徴付けられていて、図32と次のセクションの議論を見てください。周期的な左右対称の荷重でのそのような構造のための静的な計算は周期的な左右対称の置き換えにつながります。周期的な対称状態を条件としてちょうど一つの分野、いわゆるデータ分野の考慮にそのような計算を抑えることができます、すなわち、分野の正しい境界は左の境界と同じ置き換えを示します、円柱座標(長方形の座標でないところの!)で。これらの境界状態の適用はキーワードカードSURFACE、TIEとCYCLIC SYMMETRY MODELの使用で大いに簡素化されます、左と右の境界と分野サイズの節点を定義して。そして、適切な倍数節約ポイントは自動的に発生します。また、摂動の周波数解析の前に静的な予荷重方法にこれを使用することができます。

Loading

(三浦)段階に中で残留応力以外のすべての荷重を指定しなければなりません。時間に依存する振幅歴史はAMPLITUDEキーワードを使用することで大きさを変更することができます。これは非線形の静的で、非線形のダイナミックで様式のダイナミックな手順だけのために理解できます。既定値荷重の歴史は、STATIC手順のための斜面機能とDYNAMICとMODAL DAYNAMIC手順のための荷重方法です。

Point loads

(三浦)荷重点はCLOADキーワードによってメッシュの節点に適用されます。負荷点が前の方法で指定された指示に節点に荷重点をかけると、この負荷点を取り替えます。さもなければ、それは加えられます。CLOADカードの上の変数OP=NEWカードに効果だけを取ります。つぶれる方法はいつも前のすべての負荷を取り除きます。

Input deck format(入力デッキフォーマット)

(進藤)

この章では有限要素法モデルについて述べている。ジョブネームはコマンドライン上の -i flagの後の理由より確定している。CalculiXを始めるとき、inpという名の入力ファイルを見るでしょう。故に、入力デッキを「ビームinp」と呼ぶプログラムのように実行ファイルを「CalculiX」と呼ぶでしょう。

i fragはcalculixを起こした後すぐにジョブネームに従うならば落とすことができる。CalculiXは、ジョブネームdatと名づけたアウトプットと-o flagで選ぶことができる名前の他のアウトプットファイルを起こすことができる。もし、-o fragの後にfrdがあるならば、ジョブネームfrdは起こされる。それは、cgxで見ることができる。もし、-i opinionの後にonfがあるならば、いくつかのonfファイルが起こされる。もし、最適化が所有書のソフトウェアTOSCAでうまく働くならば、これらは必要である。初期状態はfrd。したがって次の2つの命令は同じ効果をもつ。

CalculiX -i beam

CalculiX -i beam -o frd

もし、方法がFREQUENCYの方法ならば、CalculiXは固有振動数、固有形式とジョブネームeigと名づけられた質量マトリクスの2つの部分からなるファイルを発生するだろう。もし、方法がMODEL DYNAMICの方法ならばCalculiXは、あの名前を探すだろう。もし、必要ないくつかのファイルが存在しなかった時エラーメッセージが出て、CalculiXは止まるだろう。

あるキーボードセットの基本的なデッキを打ち込んだあとに、キーワードの下に配置されたキーワードによってデータが求められる。キーワードは、同じ行の範囲によって伴って起き、コンマによって分けられる。もし、パラメータが数値を求めるなら、同等の記号はパラメータと数値は結ばれなければならない。入力での空白は、重要ではなく好きなように入れることができる。キーワードといくつかの他の英字と数字をいくつか組み合わせた情報は高等な場合、下等、中間で書くことができる。入力デッキは応じない場合があり、本質的にすべての英字と数字を組み合わせた特徴は高等な場合で変えられる。データは、固定した形式に従わずコンマによって分けられない。線はキーワード定義によってできるだけ多くのデータだけ含むことができる。マテリアルセットで、names.e.gと定義された利用者の中間の長さは、もし別の方法を詳しく述べなければ20の特徴がある。入力デッキの構造は、最初のステップの定義の前に、幾何学、位相的、物質的データを、1つ以上の後のステップで付加データ(機械的、熱的や変位の規定)を含む。利用者は、一致した構成単位(構成単位は計算では表れない)で、全てのデータが与えられることをよく確かめべきである。

キーワードはタイプステップかモデル定義になることができる。モデル定義カードは最初のステップカードの前に使わなければならない。ステップキーワードはステップ無しで使われることだけできる。モデル定義キーワードの間で、物質のそれらは特別な場所を占領する:それらは物質の特性を明確にし、後のMATERIALカードと一緒に分類すべきである。 振動の節と要素は、同じ名前を共有できる。本質的に、名前は高等な場合で入れ替えられ、振動の節の一式の後に'N'か要素の一式の名前の後に'E'が加えられる。従って、一式の名前は誤って書かれ、警告メッセージは高等な場合で書かれていて、'N'や'E'が加えられるだろう。

キーワードカードは順序よくアルファベット順。

User subroutines.(ユーザーサブルーチン)

(小林) このソフトウェアは、GNUのGPLライセンスで守られてはいるが、ユーザーは常にソースコードを得るべきだ。時には、全てのプログラムを分析しなければならない代わりに、あなた自身のルーチンを、上手に記述されたユーザーインターフェースに接続させるほうがより実際的です。それゆえに、特定のタスクのために定義されたインターフェースは、ユーザーの好きなように置かれます。これらのインタフェースは、基本的にサブルーチンヘッダを含んでいるFORTRANサブルーチンで、入出力変数の記述とこれらの変数のための宣言文です。ルーチンの主要部は、ユーザーによって書かなければなりません。

後藤ちゃちゃ

at the disposal of the user は「ユーザーが好きなようにできる」とか 「自由にできる」ということでは。taskはタスクでいいかな。 declaration statement は宣言文のことでしょう。

creep

(小林)ユーザーサブルーチン「ucreep.f」は、入力デッキでキーワードシーケンス*CREEP LAW=USERを選ぶことによって、彼自身のクリープ則をユーザーが取り入れることが可能で都合の良い形に作り変える。そのクリープ則は、同等のクリープ重圧程度の機能としてミーゼスストレスを与えます。ヘッダ、そして入力変数の説明は次のようである。 サブルーチンucreep(amat,iel,iint,tll,epini,ep,dtime,svm,dsvm) 入力: amat:材料名 iel:要素番号 iint:統合点番号 tll:温度 epini:同等のクリープひずみ増加の始まり ep:現在の同等のクリープひずみ dtime:時間増加 出力: svm:現在のミーゼスストレス dsvm:現在の同等のクリープひずみに関して、本当のミーゼスストレスから造られたもの。数値的に:少量に対する現在の同等のひずみを変え、現在の本当のミーゼスストレスの変化量を計算し、そして、前の量を通じて後の量を分けて下さい。

hardening

(小林)サブルーチンでは、「uhadening.f」ユーザは、キーワードシーケンス*PLASTIC HARDENING=USERによって決定される粘塑性物体のために彼自身の等方性および/または運動学的な移動硬化則を挿入することができます。説明があるヘッダと変数は後に続きます。 サブルーチンuhardening(amat,iel,iint,tll,epini,ep,dtime,fiso,dfiso,fkin,dfkin) 入力: amat:材料名(最大20の性格) iel:要素番号 iint:統合点番号 tll:インクリメント終了後の温度 epini:同等の元に戻らないひずみ増加の始まり ep:現在の元に戻らない同等のひずみ dtime:時間増加 出力: fiso:現在の等方性硬化ミーゼスストレス dfiso:現在の等方性硬化タンジェント(同等の元に戻らないひずみに関するミーゼスストレスの派生物) fkin:現在の運動硬化ミーゼスストレス dfkin:現在の運動学的なタンジェント(同等の元に戻らないひずみに関するミーゼスストレスの派生物)

user-defined loading

(小林) これらのルーチンは、不均一性分散荷重を定めるために利用できるようにされます。ユーザーは、position,timeその他の機能として、別にそれぞれの統合点で積載を定めることができます。

dflux

(小林) このサブルーチンは、不均一性熱流のために使われます。SxNUy形式の分布荷重レッテルによって特徴づけられ、DELUXdfluxを参照する。その一荷ラベルは、長さ最高20文字でありえます。特に、yは異なる不均一性流動パターンを識別するのに用いることができます。説明があるヘッダと変数は、あとに続きます。 サブルーチンdlux(flux,sol,kstep,kinc,time,noel,npt,coords,jltyp,temp,press,loadtype)。

ユーザーサブルーチンdflux。

入力: sol:現在の温度値 kstep:ステップ番号 kinc:増加番号

       time(1):現在のステップ時間 time(2):現在の全体の時間  noel:要素番号
       npt:統合点番号  coords(1..3):積分点の世界的な座標。  
        jltyp:faceコードをロードすること。  11=face 1  12=face 2  13=face 3
       14=face 4  15=face 5  16=face 6   temp:使われない流れ。 press:使われな        
       い流れ。loadtype:タイプラベルをロードしてください。

出力: flux(1):流動の大きさ。

       flux(2):使われません;上曲がりに値を割り当てないでください。

User-defined material laws

(小林) これはユーザがその人自身の材料を定義する非常に重要で強力なインターフェースです。サブルーチン「umat.f」はドライバーサブルーチンです。ユーザの定義されたルーチンが「umat_user.f」と同類であると言って、モデルの現在の材料の種類に依存します。新素材法を作成するために、「umat_user.f」ルーチンを書かなければなりません、そして、適切な呼出しをumat.fに挿入しなければなりません。"umat.f"では、ユーザの材料の名前は定義されます。これによるMA_TERIALmaterialが始まらなければならないキーワードカードに続いて、キャラクタがNAMEパラメータを結ぶということです。あなたが名前FUNNY_MATERIALと共に新素材を定義するなら、この材料を使用するために、物質的な名前はこのストリングから始まらなければなりません。これはCalculixの中の事前に定義されるのとユーザの定義された材料の間の用法の主な違いです:あなたが自由に材料に名前を完全に選ぶことができる事前に定義された材料を使用するなら、名前は事前に定義されたもので始めを命名しなければなりません。長い間物質的な名前が最大20のキャラクタであるかもしれないので、一般に十分に、この材料、例えば、FUNNY_MATERIAL1、FUNNY_MATERIAL2などのいくつかのバージョンを自由に定義します。umat_ユーザルーチンのヘッダーと入力/出力変数は以下の通りです。 サブルーチンumat_user(amat,iel,iint,kode,elconloc,emec,emec0,beta,xokl,voj,xkl,vj,ithermal,tll,dtime,time,ttime,icmd,ielas,mint_,nstate_,xatateini,xstate,stre,stiff,iorien,pgauss,orab) ユーザに関する剛性と圧力が物質的な法を定義したと見込みます

icmd=3:機械的な緊張で圧力について計算します。 else:機械的な緊張と剛性マトリックスで圧力について計算します。

入力: amat:材料名 iel:要素番号 iint:統合点番号 kode:物質的なタイプ(*USER MATERIALの下で入る定数のー100#):材料に異なった数の定数で使用することができる。 elconloc(21):実際の温度tllのために補間されたキーワードカード*USER MATERIAL(最大21、実際の#=ーkodeー100)によって定義されたユーザーの定義した定数。 emec(6):増分の終わりのラグランジュの変形テンソル(コンポーネントオーダー:11、22、33、11、12、23)。 emec0(6):増分の始めのラグランジュの機械的な変形テンソル(熱緊張は引き算されます)。 beta(6):残留応力伸筋(キーワード*INITIAL CONDITIONS、TYPE=STRESSの下で入れた圧力) xokl(3,3):増分の始めの変形勾配。 voj:増分の始めのヤコビアン。 xkl(3,3):増分の終わりの変形勾配。 vj:増分の終わりのヤコビアン。 ithermal:0:熱運動は全く考慮に入れられません。1:熱運動は考慮に入ります。(キーワード*INITIALCONDITIONS、TYPE=TEMPERATUREによって引き起こされます。 tll:増加温度。 dtime:増分の時間の長さ。 time:現在の増分の終わりのステップ時間。 ttime:現在の増分の始めの合計時間。 icmd:3と同輩でない:圧力と剛性について計算して下さい 3:圧力だけについて計算して下さい。 ielas:0:弾性の繰り返しがありません:逆にできない効果は許容されています。 1:弾性の繰り返し、不可逆性の変形は許容できませんでした。 mint_:最大限にしてください。モデルの要素を通じて統合点の#を。 nstate_:最大限にしてください。モデルの状態の変数の#を。 xstateini(nstate_,mint_,# of elements):増分の始めの状態の変数。 xstate(nstate_,mint_,# of elements):増分の終わりの状態の変数。

stre(6):増分の始めの第2種のパイオラキルヒホフ圧力。  iorien:統合点における、ローカル座標軸のシステムの数は危機にひんしています(ローカルシステムが全く適用されないなら、値、0をとります。  pgauss(3):統合点の国際的な座標。
orab(7,*):全ての国際的な座標システムの記述。ローカル座標システムを適用するなら、skl(3、3)でローカルの張筋に先にかけ算させることによって、国際的な張筋を得ることができます。サブルーチンをトランスフォーマトリックスと呼ぶことによって、sklは決定します。
 

出力: xstate(nstate_,mint_,# of elements):増分の終わりのアップデートされた変数。 stre(6):増分の終わりの第2種のパイオラキルヒホフ圧力。 stiff(21):増分の終わりの物質の理論構成の枠組みの中の一貫した接線の剛性マトリクス。言いかえれば:ラグランジアン緊張伸筋に関するPK2圧の派生物。マトリクスは左右対称であると想定され、上半分だけが完全に異方性の弾性の材料(*ELASTIC,TYPE=ANISO)のように同位で与えられることになっています。

パラメータielasは、逆にできない効果が考慮に入れられるべきであるかどうかを示します。強制的な置き換えは最初の繰り返しにおける非常に大きい緊張の原因となります。したがって、最初の繰り返しが完全に直線的であるなら、準定常の計算における集合はしばしば機能強化され、すなわち、材料と幾何学上非線形性はオフにされます。パラメータielasは適切なフラグです。交換目的のためのユーザの支配権は2つの余分なルーチンがあります。オイラー緊張へのひそかなラグランイアン緊張、PK2圧へのコーシー圧力、またはPK2圧へのキルヒホッフ圧力に"str2mat.f"を使用することができます。ヘッダーと短い記述は以下の通りです。

サブルーチンstr2mat(str,ckl,vj,cauchy) 空間的な座標における、圧力を空間的な座標への物質的な座標における物質的な座標か、緊張に変換します。

入力: str(6):コーシー圧力、キルヒホフ圧力またはラグランジュの緊張コンポーネント注文:11、22、33、12、13、23。 ckl(3,3):逆さの変形勾配 vj:ヤコビアンの決定因子。 cauchy:理論変数 本当に、: strがcauchy圧を含んでいるなら。 虚偽で、: strがキルヒホッフ圧力かラグランジュの緊張を含んでいるなら。

出力: str(6):第2種(PK2)かオイラーの緊張のパイオラキルヒホフ圧力

2番目のルーチン、"stiff2mat.f"は接した剛性マトリクスを空間的な座標から物質的な座標に変換します。

サブルーチンstiff2mat(elas,ckl,cauchy)

入力: elas(21):すなわち、空間的な記述、オイラー緊張に関するコーシー圧力か、キルヒホフ圧力の派生物の剛性定数。 ckl(3,3):逆さの変形勾配。 vj:ヤコビアンの決定因子。 cauchy:理論変数 本当に:elasはコーシー圧力に関して書かれています。 虚偽で:elasはキルヒホフ圧力に関して書かれています。

出力: elas(21)物質的な記述、すなわち、ラグランジアン緊張に関する2番目のパイオラキルヒホフ圧力(PK2)の派生物における剛性定数。現在、いくつかの例を与えます。

Program structure.

(小林) CalculiX CrunchiXは、FORTRAN77(FORーTRAN90)からの混成です。そしてC.Cが主に自動配分と再割当目的のために使われます。FORTRANは私が学んだ第一用語であって、今もなおFORTRANの熱中者であると認めなければなりません。私は必要な場合、Cを使用します。私はそれをさけることも可能です。もう一つのオプションはFORTRAN90ですべてをコード化することになっていました。しかし、この日付まで良いFREE FORTRANコンパイラがありません。大雑把に言うと、メインルーチンとメインによって呼ばれる一部のルーチンはCです。他の人はFORTRANです。Cルーチンでないこの手段はFORTRANルーチンによって呼ばれた。FORTRANルーチンはCルーチンまたはFORTRANルーチンによって呼ばれるかもしれません。NO下院がコードにあります。全てのデータの転送は、サブルーチンの呼び出しの議論を通しています。全ての議論はアドレスによって移され、価値でない(1、2例外が、コードでこの規則上にあるかもしれません。)。CaluculiXの主なサブルーチンは、ccx_1.3.cです。それは以下の部品からざっと成ります。

.フィールドの配分

それぞれのステップのために

1、ステップを読むことは、データ(最初のステップのために、prestepデータを含む)を入

       力しました。

2、マトリックス構造を決定すること。

3、結果を保存して、方程式の結果を満たして解くこと。

Verification examples.

(小林) 確認例は、異なった特徴をテストするために適当な単純例です。それらは、新機能を使う時、CaluculiXのインストールが正しいかどうか調べるのに使われ、また例を見つけることができる。ここでそれらは、何をテストされているか短い記述と一緒にアルファベット順にリストされます。入力ファイルのために、”inp”を追加してください。結果ファイルのために「dot、ref」を追加してください。全てのファイルは配給の中に含まれます。


トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2020-01-19 (日) 17:03:31