CalculiXメモ
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
-http://www.calculix.de/
*ubuntu22へのインストール [#heab231f]
$sudo apt install calculix-ccx calculix-cgx
バージョンはやや古いものにはなるけど。
*近況? [#r848e3fe]
*PrePoMax [#s8b7c532]
-https://aytechlab.com/homefem/prepomax001/
Salomeでunvファイルを作ったのをPrePoMaxで読み込んで、材料...
**Ubuntuへのインストール [#oad3b7ff]
***wine環境のインストール [#s236b285]
PrePoMaxはWindowsソフトなので、まずは、Linux上でWindowsを...
をインストールし、
さらにその環境でNET Frameworkを使えるようにするために、do...
-[[オープンＣＡＥコンサルタントOCSE^2:https://ocse2.com/?...
$ sudo dpkg --add-architecture i386
$ sudo mkdir -pm755 /etc/apt/keyrings
$ sudo wget -O /etc/apt/keyrings/winehq-archive.key http...
$ sudo wget -NP /etc/apt/sources.list.d/ https://dl.wine...
$ sudo apt update
$ sudo apt install --install-recommends winehq-devel
$ sudo apt install winetricks
$ sudo winetricks --self-update
$ winetricks dotnet48
-最後のdotnet48のインストール中にwineが立ち上がり、「wine...
-「マイクロソフトソフトウェア追加ライセンス条項」みたいな...
-「Wineの設定」ウィンドウが出てきたら、OKをクリック。2回...
-「Microsoft.NET Framework」の警告（正しく実行されない可...
-「.NET Framework 4.8セットアップ」（追加ライセンス条項）...
□□□□(A)みたいなチェックボックスがたぶん「同意する」なので...
-「.NET Framework 4.8がインストールされました」というウィ...
-「インストールを完了するには...再起動する必要があります...
***PrePoMaxのインストール [#ma0cf087]
-~/.wine/drive_cの中にexeというディレクトリを作り、wineで...
この中にソフト名がわかるサブディレクトリーを作って、その...
いいのではないか。ということで、PrePoMaxに関しては、~/.wi...
ディレクトリを作っておく。
-[[PrePoMaxのダウンロードサイト:https://prepomax.fs.um.si...
から最新バージョン（または動作確認済みのバージョン）のzip...
~/.wine/exe/prepomax にダウンロード
-unzip 'PrePoMax v1.4.0.zip'で展開。
-cd 'PrePoMax v1.4.0'で展開されたディレクトリーに移動。
-wine PrePoMax.exeを実行し、PrePoMaxが起動することを確認。
-Ubuntu(MATE)の上部パネルで右クリックし、「パネルへ追加」...
「カスタム・アプリケーションのランチャ」を選択して「追加...
--「ランチャの作成」ウィンドウが開くので、名前：PrePoMax,...
/usr/bin/wine '/home/kouzou/.wine/drive_c/exe/prepomax/P...
--ランチャーにショートカットアイコンが作成される。アイコ...
PrePoMaxを起動して、helpのabout辺りから、PrePoMaxのロゴの...
~/.wine/drive_c/exe/PrePoMax辺りに適当な名前(aikon.pngと...
ショートカットアイコンを右クリックして、プロパティーから...
**Windows版（PrePoMax・OptiMax） [#a6c64618]
-http://lace.fs.uni-mb.si/wordpress/borovinsek/
--雰囲気はこんな感じ → https://www.youtube.com/playlist?l...
***PrePoMax（プリポスト付き） [#wc54148c]
-http://lace.fs.uni-mb.si/wordpress/borovinsek/?page_id=41
***OptiMax（最適化） [#j3c8e943]
-http://lace.fs.uni-mb.si/wordpress/borovinsek/?page_id=104
--https://www.slideshare.net/OpenCAE_kanto_fluid/optimax-...
--Linuxで似たようなことをしようとするなら、[[Dakota:https...
--[[パラメータ最適化ツール DakotaとOpenMDAOについて:https...
**Windows版（CalculiXLauncher） [#g78cbed6]
-http://www.calculixforwin.com/?page_id=31
--ダウンロードは&link(ここ,https://sourceforge.net/projec...
--起動は、Launcher33x64.exe をクリック。
-ソルバー(ccx)は普通に動くが、プリ・ポストプロセッサ(cgx)...
-CCX WIZARD! はLinuxのみで使えるとのこと。
-UNV→inpコンバーターも使える模様。
**FreeCADでCalculiXをソルバーとして解析できる？ [#ve8353ba]
-http://naga-tsuzuki.sblo.jp/article/179282504.html
-http://naga-tsuzuki.sblo.jp/category/4426804-1.html
**SalomeとCalculiXを連携できそうなCalculiXLaucher（要調査...
-http://calculixforwin.blogspot.jp/2015/05/calculix-launc...
-http://calculixforwin.blogspot.jp/2016/10/getting-starte...
Salomeで作ったメッシュファイルをCalculiXLauncherでinpファ...
**直交異方性材料 [#qad7b617]
今まで通りのひずみ－応力関係の逆関係を入力する方法
*ELASTIC,TYPE=ORTHOROPIC [#nc477d23]
D1111,D1122,D2222,D1133,D2233,D3333,D1212,D1313
D2323
の他にも，Code_Asterと同様な入力方法↓もできるようになっ...
*ELASTIC,TYPE=ENGINEERING CONSTANTS [#r0c7a606]
E1,E2,E3,nu12,nu13,nu23,G12,G13,
G23
**umat [#u6aa51b9]
-https://www.xsim.info/articles/CalculiX/ccx-doc/node315....
**unvからinpへの変換 [#g756d177]
-http://freeplanets.ship.jp/FEM/Calculix/Calculix.html#mo...
-http://caelinux.org/wiki/index.php/Proj:UNVConvertからMe...
-unzipで解凍
--(
python unv2abaqus.py Mesh_1.unv mesh
--)
でmesh.inpができてる。
**インストール [#a4732bc2]
CalculiXは、バージョン1.8以降は実行ファイルを展開しただけ...
SPOOLES.2.2やARPACKの導入も必要となり、インストールが難し...
恐らくインターフェースやプリ・ポストプロセッサーなどがGUI...
思われる
[[DEXCS:http://dexcs.gifu-nct.ac.jp/pukiwiki/index.php?%C...
は、ADVENTUREを利用していて、今のところ動的解析などには対...
いずれ、
&link(CalculiXを導入して動的解析にも対応する方向性,http:/...
みたいなので、様子を見たい。
CalculiXのrpm版もバージョン2.0以降のものも検索すればそれ...
(64ビット版が主流で32ビット版は少ないので注意)、依存性が...
Vine等へのインストールは簡単ではない。
&link(CAElinux,http://caelinux.com/CMS/index.php)の64bit...
プレインストールされてるようなので、64bitパソコンにCAElin...
計算サーバーとするか。
***Ubuntu 16.04 LTS へのインストール [#v9c88d88]
Ubuntu16.04LTSを(近藤のPCに)インストールして、色々と試...
Ubuntu16.04LTSの&link(公式パッケージ,http://packages.ubun...
--(
$ sudo apt-get install calculix-ccx
--)
これだけでいつも通りの計算できるようだが、cgxはパッケージ...
-ccx_2.6との違い
--荷重や強制変位の方向の書き方が変わった。例えば、&link(c...
--(
*CLOAD [#heba7041]
** LAST, 2, 6.6666666666666666E-002 [#lf4b385a]
LAST, 2, 2, 6.6666666666666666E-002
--)
後藤ちゃちゃ(17/3/29)**はコメント文では？
---近藤(17/3/29)そこの行でエラーが出るので、それをコメン...
---後藤(17/3/30)そうですか。CalculiXのinpファイルはABAQUS...
--接触解析のとき、接触面の組み合わせの他にも、面同士の種...
--(
*CONTACT PAIR,INTERACTION=SI1,TYPE=NODE TO SURFACE [#w55b...
--)
--というように、TYPEが必要。上記以外にも、スレーブ面がマ...
-後藤ちゃちゃ(2016/11/14)：そういうことなら、2017年度以降...
Ubuntu12.04は2017/4までのサポートだし、それまでCAELinuxの...
**CAELinux [#l4d02da6]
&link(CAElinux,http://caelinux.com/CMS/index.php)はUbuntu...
インストール後にアップグレードしたらccxも最新版の2.3にア...
これはいいかもしれない。
Ubuntuのクライアント版はサーバー系ツールのセットアップが...
めんどくさいような気がしていたが、CAELinuxにはsshdがプリ...
サービス管理(system-config-services)は、synapticの一覧に...
なんか他のツールがあるのかな。
*コンクリート [#p7ee6526]
コンクリートのように圧縮のみに抵抗する材料も扱える。
-例えば、
--(
*MATERIAL,NAME=COMPRESSION_ONLY [#u67912f5]
*USER MATERIAL,CONSTANTS=2 [#re38db97]
14.e3, 1.e3
--)
-と書くと、ヤング率が14MPaで、許容する最大引張応力が1MPa...
--ところで、ポアソン比を入力するところがないのはどうなっ...
-COMPRESSION_ONLYの他にも特殊な材料がたくさんある(/opt/Ca...
-参考：http://open-shelf.appspot.com/CalculiX/ccx-doc/nod...
*面の定義 [#bab7422b]
マニュアル(ccx_2.6.1.pdf)のp.367に要素ごとの面番号の定義...
--(
面1:節点1,2,3,4
面2:節点5,8,7,6
面3:節点1,5,6,2
面4:節点2,6,7,3
面5:節点3,7,8,4
面6:節点4,8,5,1
--)
というように定義されている。例えば、接触解析とかで、要素...
--(
*SURFACE,NAME=men [#ff6f59b8]
1,S6
2,S6
3,S6
4,S6
--)
みたいにして使う。
ちなみに、接触解析をするとき、マスター面の定義は必ず上記...
*ヒンジ [#hf0766e9]
**立体要素と板要素を組み合わせる場合 [#r4ff8e4e]
CAELinuxの場合、/opt/CalculiX_2.6.1/ccx_2.6.1.pdfがマニュ...
ここのp.96によると、「立体要素と板要素」あるいは「立体要...
以下、その例題(/opt/CalculiX_2.6.1/ccx_2.6.1/test にある)
-solidshell1.inp:ただ、立体と板の節点を共有させただけのモ...
-solidshell2.inp:立体要素に重ねて板要素を入れることで剛結...
**多点拘束を用いる場合 [#z8db4833]
-&link(hinzi.inp,http://www.str.ce.akita-u.ac.jp/~gotouha...
--折り目部分を重節点にしたうえで"*EQUATION"を使い、互いの...
**立体要素だけ [#r910f489]
特別なことをしなくても、立体要素同士を線状に接点共有さ...
この方法なら、折り畳み円筒を一定厚にし、かつ蝶番を再現...
-&link(hinzi2.inp,http://www.str.ce.akita-u.ac.jp/~gotouh...
**梁の場合のヒンジ節点 [#m5b4fa4d]
-&link(メモ1,http://www.str.ce.akita-u.ac.jp/~gotou/kouzo...
-&link(メモ2,http://www.str.ce.akita-u.ac.jp/~gotou/kouzo...
*ccx_2.5 [#fc2ba6bc]
ざっといじっていてわかったこと(変更点?)を箇条書きで
-例の法線の向きが互い違いになっているダイヤカット&link((d...
--ちなみに同じinpファイル&link((daiya6h.inp),http://www.s...
---2.5
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/test/daiy...
&link(daiya6h.dat,http://www.str.ce.akita-u.ac.jp/~gotouh...
---2.3
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/test/daiy...
&link(daiya6h23.dat,http://www.str.ce.akita-u.ac.jp/~goto...
---ある載荷節点(節点番号17から20)の変位
2.5
--(
17 9.337687E-11 -1.359864E-23 -1.388775E-09
18 1.085580E-23 9.337687E-11 -1.388775E-09
19 -9.337687E-11 -8.762945E-24 -1.388775E-09
20 5.218346E-24 -9.337687E-11 -1.388775E-09
--)
2.3
--(
17 9.3377E-11 -2.7324E-23 -1.3888E-09
18 7.1344E-24 9.3377E-11 -1.3888E-09
19 -9.3377E-11 -1.1550E-23 -1.3888E-09
20 -6.2168E-24 -9.3377E-11 -1.3888E-09
--)
-上の画像からわかるように、2.5で解いて出来たfrdファイルを...
inpファイルで設定した厚さを反映したモデルを見ることができ...
--cgxで節点を見ると
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/test/25se...
のようになっていて、これは&link(マニュアル(the documemtat...
--ということはccx_2.5で解くと、frdファイルに新しい節点が...
-ちなみに&link(このへん,http://www.str.ce.akita-u.ac.jp/c...
*ccx_2.6.1をCAELinux2011で実行できるように（まだできない...
CalculiXの作者guido dhondtさんに問い合わせたところ、外向...
解けるものが2.3では解けない問題は、2.6では解決されている...
ccx_2.6.1を使うのがよさそうだ。
ccx_2.6.1のa Linux executable(ccx_2.6.1.tar.bz2)を適当な...
--(
$ bunzip2 ccx_2.6.1.tar.bz2
$ tar xvf ccx_2.6.1.tar
--)
で展開。CalculiX/ccx_2.6.1/src/の中にある実行ファイルccx_...
libgfortran3がもしインストールされてなければ、synaptic等...
--(
$ ls -l /usr/lib/libgfortran*
--)
とやって、/usr/libにあるlibgfortranのバージョン（とリンク...
恐らく、/usr/lib/libgfortran.so.3.0.0があって、
/usr/lib/libgfortran.so.3（青色）からリンクが張られている。
この状態で ccx_2.6.1を実行すると、恐らく以下のエラーが出...
--(
error while loading shared libraries: libgfortran.so.2: c...
--)
そこで、以下のように/usr/libに移動してからlibgfortran.so....
--(
sudo ln -s libgfortran.so.3.0.0 libgfortran.so.2
--)
これでccx_2.6.1を実行すると、Usage: CalculiX.exe -i jobna...
--(
ccx_2.6.1 hoge
ccx_2.6.1: symbol lookup error: ccx_2.6.1: undefined symb...
--)
*1.7と2.3の比較 [#x1bae6f5]
-datファイル
単純に&link(波型構造の.inpファイル,http://www.str.ce.akit...
&link(ver.1.7,http://www.str.ce.akita-u.ac.jp/~gotouhan/e...
--オーダーが大きく違う
--2.3だと指定したはずの拘束点がいくつか拘束されていない(e...
-上の計算で得られた.frdファイルの図
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/daiya/1-7...
*マニュアル [#iaba6b2d]
-&link(アクセス制限あり,http://www.str.ce.akita-u.ac.jp/c...
ccxのマニュアル(日本語版)を見つけたので、ここに載せました...
[[日本語版:https://www.xsim.info/articles/CalculiX/ccx-do...
[[英語版:http://www.dhondt.de/ccx_2.20.pdf]]
*inpファイルの書き方 [#sf1da85e]
-http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html
*サンプルファイル [#p93aad99]
-http://web.mit.edu/xenophon/MacData/afs/athena/software/...
**コメント [#na23a417]
--(
** こんなふうに ** を最初に書くとコメントが書ける [#ybc94...
** つまり ** の行はデータ入力の際に無視される [#bffc99c9]
--)
**要素の節点番号 [#c3aa4909]
inpファイルで、まず各節点と節点座標を定義してから、
各要素ごとに要素番号とその要素の各節点の節点番号を設定し...
要素の種類ごとに1要素の節点番号の並べ順には決まりがあるの...
マニュアルやサンプルの節点番号の振り方を参照すること。
例えば、8節点シェル要素(S8)では、頂点の節点番号を4つ並べ...
辺中点の節点番号を4つ並べなければならないとか。
これを間違うと、変な形になる。
**拘束条件 [#e1708cfb]
***複数の節点にまとめて同じ拘束条件を与える場合 [#zfc4c4cd]
--(
*NSET,NSET=FIX [#eeeef362]
** ここに拘束したい節点番号とコンマを1行ずつ並べる [#aa85...
** 例えば拘束節点が1と2だったら [#c7ac0403]
1,
2,
** みたいに [#u5417dea]
*BOUNDARY [#zc4e42a5]
** 以下のFIX,に続けて拘束変位1=x,2=y,3=zをコンマ区切りで ...
FIX,1,3
** 1から3(xからz)を拘束 [#eddf2c83]
--)
上記の例を、
--(
*BOUNDARY [#f0883ca8]
FIX,1
FIX,2
FIX,3
--)
みたいに分けて書いても良いっぽい
***1節点ずつ拘束条件を与える場合 [#afea8f27]
--(
*BOUNDARY [#ef691e58]
1, 3, 0.
2, 1, 3, 0.
** 節点番号1の3(z)を拘束 [#a104396a]
** 節点番号2の1から3(xからz)を拘束 [#z6cf218e]
--)
**材料定数 [#mc506e0e]
--(
*MATERIAL,NAME=EL [#p3c4460f]
*ELASTIC [#y26c4896]
** 以下にヤング率(MPa)とポアソン比 [#c1c9337e]
210000., 0.300000012
--)
**載荷条件 [#u7bd11be]
***複数の節点にまとめて同じ荷重条件を与える場合 [#s5e6554b]
--(
*NSET,NSET=SAIKA [#s2d94d1b]
** ここに載荷したい節点番号とコンマを1行ずつ並べる [#r17b...
** 例えば載荷節点が3と4だったら [#pcce862b]
3,
4,
** みたいに [#d1122947]
*STEP [#zb71651c]
*STATIC,SOLVER=SPOOLES [#t9a84e1a]
*CLOAD [#h94d6348]
** 以下のSAIKAに続けて載荷方向、載荷荷重を書く [#z72c0b8a]
** 例えば載荷方向が3=zで、載荷荷重が1/載荷節点数なら [#ab...
SAIKA,3, 0.5
--)
***1節点ずつ荷重条件を与える場合 [#ja2a05ed]
--(
*CLOAD [#w98db071]
1, 2, 10.
2, 3, 20.
** 節点番号1の2(y)方向に10MN載荷 [#xe96bb41]
** 節点番号2の3(z)方向に20MN載荷 [#x169b9fa]
--)
*ひずみについて [#pf595004]
inpファイルで最後の方に、
--(
*EL FILE,POSITION=AVERAGED AT NODES [#icb5206f]
S,E,PE
--)
っていれる。
ちなみに、Sは応力、Eはひずみ、PEって入れると相当塑性ひず...
ほかにもいろいろ出力できそう。
詳しくはマニュアルの
EL FILEとかNODE FILE,EL PRINT,NODE PRINT 参照..
*座屈解析 [#u2a9680b]
**datファイル内の座屈荷重の読み方 [#qe40364e]
-inpファイルのCLOADで載荷節点の倍率bを書き込む。
-datファイルのEIGEN VALUEのとこに書き込まれる荷重$P_{dat}...
-$P_{cr}=bnP_{dat}$
-つまり、倍率$b$を$1/n$としておけば、datファイルの荷重$P_...
ということでいいかな?(次項の結果と整合するかな)
-比較しやすいように、datファイルの座屈荷重を頭にすれば
--$P_{dat}=\frac{P_{cr}}{nb}$
-注意：inpファイルをgfortranで作る際、上記の倍率($b=1/n$)...
桁数が多くなりすぎて、-E003の部分が読み込まれず、座屈荷重...
**inpのCLOADと、datの座屈荷重との関係 [#mf10b705]
いつもは、CLOADで(1/載荷節点数)と入力したものを、座屈荷重...
もし10/載荷節点数とかにしたら、座屈荷重はどう変わるのか見...
,CLOAD,0.001/n,0.01/n,0.1/n,1/n,10/n,100/n
,理論値との誤差,座屈せず,3%,3%,3%,3%,座屈せず
n:載荷節点数
10/nの座屈荷重は、1/nの座屈荷重の0.1倍になる。0.1/nにすれ...
だから、inpのCLOADで0.1/nで入力したら、datファイルの座屈...
ただ、サイズによって座屈が出ない問題も、(相似形状で10倍と...
ACCURACYのところを0.01ではなく、0.1とか1とかにすれば解け...
**モード図の描かせ方など [#wec7aea4]
***cgxのアニメ機能を利用する場合 [#u5791961]
-cgxでfrdファイルを開き、Datasetsを選ぶと、計算した座屈モ...
-L101とか、値が0になっているところは自明解なので、最初に...
-更にDISPを選択
-改めてDatasets→Entitiy→Allを選択すると、座屈モードがまず...
-次にAnimateのStartを選択すると、座屈モード変形がアニメ再...
-画像がほしい瞬間のところでマウスのポインターを枠外に出し...
-モードの倍率を大きくしたければ、Animate→Tune Valueで、Tu...
***倍率を調整して手動で変形モードを加算したfrdファイルが...
&link(tasu.f系のプログラム,http://www.str.ce.akita-u.ac.j...
-ccxkataiza.fの直方体要素ならば、ccx_1.7+tasudatw.fでモー...
ダイヤカットをやったシェル要素では、ccx_1.7*tasudatw.fで...
ccx_1.5+tasu.fの昔のやり方でモード図を描いていたというよ...
確認してほしい。
円筒ダイヤカットも伸縮できる円筒折り紙も、圧縮でも引っ張...
起こせていたが、これはバージョンに依存するか？ 例えば、cc...
ccx_1.7では、no convergence.....maybe no buckling occurs ...
事情があるかどうか確認。
--シェル要素でもccx_1.7+tasudatw.fの組み合わせでモード図...
**l=1,m=16,n=2 [#t0674712]
,半径(m),高さ(cm),ヤング率,ポアソン比,厚さ[m],座屈モード,...
,1.d0/2.d0/pi,50,70000,0.300000012,0.11,2,0.06,0.149253731
-これだとccx_1.7では「job finished」になるが、ccx_1.5だと...
--「job finished」は荷重1E-1～E-4の間。
-ccx_1.5では1E-5から「セグメンテーション違反です」とおこ...
--1E100以上、1E-5以下でもおこられる。その他は「*ERROR: to...
-厚さによって荷重が同じでも「no convergence for the buckl...
*バージョン比較 [#j3419520]
バージョンによって答えが違った例をここに書き込む。
-2枚の板をくっつけた折れ板を引っ張って広げると、ccx_1.5で...
too many cutbacksエラーが出る。ccx_1.7線形では解ける。
ccx_1.7非線型ではtoo many cutbacksエラー([[稲荷の卒論日誌...
*ウェブ上の情報など [#faf218d1]
-http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html
-http://www.waka.kindai.ac.jp/tea/shibue/abaqus-inpexp.html
-http://www.gsic.titech.ac.jp/~ccwww/tebiki/
-http://freecaetester.blog62.fc2.com/blog-date-200803.html
*cgxでinpファイルをみる。 [#je66249e]
--(
cgx -c hoge.inp
--)
でみることができる。
*cgxで変形図やモード図を描かせる [#u645a37d]
&link(ここ,http://www.str.ce.akita-u.ac.jp/~gotou/program...
&link(tasu.f,http://www.str.ce.akita-u.ac.jp/~gotou/progr...
使ってデータを整形する。具体的な方法はプログラムの中に書...
*CADとの連携 [#k24dd8c7]
-Available add-onsを使うとvda(ドイツ自動車工業会策定?のフ...
--(
./acis2fbd.pl < hoge.sat >hoge.fbd
--)
で変換できるようだが、サンプルのsatファイルは変換できるが...
変換できなかった。
-&link(Linuxにインストールできるフリーの3次元CAD,http://p...
簡単にインストールできて、出力フォーマットに各種形式を選...
&link(gCAD3D,http://www.gcad3d.org/)で、
DXF, IGES, VRML1, VRML2, STL, OBJ, TESS, NC-ISOにエクスポ...
で、dxf辺りからvdaやacisへのフリーの変換ツールはないかと...
&link(GID,http://gid.cimne.upc.es/)のアカデミック版ででき...
でも、あわよくば完全なオープンソースだけで処理したいので...
**cgxをプリプロセッサーとして使う [#t34ff8af]
-http://www.dhondt.de/tutorial.txt
-[[川原さんによる和訳にはパスワードによるアクセス制限を付...
ちょっとめんどくさいような気もするのですが、
プリプロセッサーを使わないことにしたとして、
CADやモデリングツールの出力から節点情報を抜き出す
プログラムを作る際に、メッシュを指定する節点の並び順をccx...
うまくやれなくて節点情報だけを頼りに自力でプログラムで対...
cgxをプリプロセッサーとして利用してfbdを作るのもそう悪く...
いや&link(blenderのメッシュ情報,http://www.str.ce.akita-u...
**blender [#ofee5e9f]
CalculiXのモデリングに使える??
-[[blenderメモ]]
-http://www.usennet.ne.jp/~kiyoshi/manual/blenderTips.html
-http://blender.jp/
-http://bear32.blogspot.com/
-http://bear32.blogspot.com/2006/11/blender.html
-http://f11.aaa.livedoor.jp/~hige/
-http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html
-http://www.geocities.jp/detoni2303/tips/tips_036.htm ...
**各種フォーマット [#c0843951]
-&link(dxf,http://www.autodesk.com/techpubs/autocad/acad2...
-&link(IGES,http://afsoft.jp/cad/p08.html)
*ccxのインストール簡単マニュアル(linux) [#jd231ab6]
**ccx_2.6 [#i63edd90]
-http://www.libremechanics.com/?q=node/9
**ccx_2.3 [#c89c57a5]
ccx_1.8以降のコンパイルはうまくいっていないが、SPOOLES.2....
-http://pugpug.lolipop.jp/yii-1.0/demos/mocapapas-blog/ta...
***SPOOLESのインストール [#t7915585]
ccxのサイト経由で、
spooles.2.2.tgzをダウンロードして、
/usr/local/SPOOLES.2.2を作成してからそこにコピーして展開。
(アーカイブファイルがディレクトリ内に納めてから圧縮されて...
ばらばらのファイルのまま圧縮されているかは、ファイルブラ...
/usr/local/SPOOLES.2.2/内にファイルが展開される。
Make.inc内の
--(
# CC = gcc
CC = /usr/lang-4.0/bin/cc
--)
の行を
--(
CC = gcc
# CC = /usr/lang-4.0/bin/cc
--)
に修正して、make libを実行。spooles.aが生成されれば成功。
***ARPACKのインストール [#s2cccdd9]
ccxのサイト経由で、
arpack96.tar.gzとpatch.tar.gzをダウンロードして、/usr/loc...
tar xvzfで展開すると、/usr/local/ARPACK(が生成されてその...
/usr/local/ARPACKの中に入って、ARmake.inc内の以下の箇所を...
--(
#home = $(HOME)/ARPACK
home = /usr/local/ARPACK
(中略)
#MAKE = /bin/make
MAKE = make
(中略)
#SHELL = /bin/sh
SHELL = sh
--)
make libすると、libarpack_SUN4.aが生成される。
この状態で、
/usr/local/CalculiX/ccx_2.3/src
に入り、makeしてccx_2.3をコンパイルすると、
libarpack_INTEL.aがないというエラーが出るので、
libarpack_SUN4.aをlibarpack_INTEL.aにコピーする。
これで、ccx_2.3をコンパイルしても
--(
../../../ARPACK/libarpack_INTEL.a(xerbla.o): In function ...
xerbla.f:(.text+0xe): undefined reference to `s_wsfe'
xerbla.f:(.text+0x29): undefined reference to `do_fio'
xerbla.f:(.text+0x44): undefined reference to `do_fio'
xerbla.f:(.text+0x49): undefined reference to `e_wsfe'
xerbla.f:(.text+0x5d): undefined reference to `s_stop'
--)
みたいなエラーが大量に出る。64bitパソコンでも同様。
/usr/local/ARPACK/UTIL内のsecond.fの
--(
EXTERNAL ETIME
--)
の行をコメントアウトしてARPACKをコンパイルしなおせばよい...
記述も見つけてためしてみたが、エラーの数は減ったような気...
うまくいかない。
-http://freecaetester.blog62.fc2.com/blog-entry-238.html
-http://www.kkaneko.com/rinkou/linux/arpack.html
**ccx_1.8 [#h4eb2761]
ccx_1.8のバイナリーは、Vine LinuxでもUbuntuでも動かない。
--(
bash: /usr/local/bin/ccx_1.8: cannot execute binary file
--)
のエラー。
ソースからmakeするには、ARPACKとSPOOLES.2.2をインストール...
調査中。ARPACKとSPOOLES.2.2のライブラリを展開して/usr/loc...
ccx_1.8をmakeしようとすると、
--(
make: *** `ccx_1.8' に必要なターゲット `../../../SPOOLES....
--)
のエラー。
&link(CalculiX のmake,http://mogura7.zenno.info/~et/xoops...
-http://freecaetester.blog62.fc2.com/
--http://freecaetester.blog62.fc2.com/blog-date-20090316....
--http://freecaetester.blog62.fc2.com/blog-date-20090317....
**ccx_1.7以前 [#w28b88aa]
-&link(a short installation guide (ASCII),http://www.dhon...
&link(CalculiXのダウンロードのページ,http://www.dhondt.de...
Available downloads for the solver (CalculiX CrunchiX: cc...
するとダウンロードのwindowが開くので、''このファイルをデ...
ちなみに保存されるファイルの名前はccx_1.6.tar.bz2であるこ...
次にターミナルを開いて、保存したところにcdで移動し(home/s...
suになる。[suzuki@bri00]$\$$が[root@bri00」#に変わったの...
--(
# mv ccx_1.6.tar.bz2 /usr/local
# cd /usr/local
# bunzip2 -c ccx_1.6.tar.bz2 | tar -xf -
# cd CalculiX/ccx_1.6/src
# mv ccx_1.6 /usr/local/bin
--)
と順番に打ち込む。
#は始めから表示されているので、くれぐれも# mv ccx_1.6とは...
これだけでccxがどこのディレクトリでも使えるようになる。
サンプルをインストールしたい場合には、test examples(ファ...
マニュアルはthe documentationの中から、TeX、html、postscr...
マニュアルはTeXとhtmlの二つをダウンロードしたほうが良い。
TeXをダウンロードする利点は、印刷がきれいにできる
、ファイルが一つだけだから、
viで文中のキーワードが簡単に検索できるということが挙げら...
htmlの利点は言うまでも無く、インターネットと同じ環境で見...
この場合、ブラウザのファイル→開くで/usr/local/Calculix/cc...
その中のccx.htmを開き、ブックマークしておくのがいいだろう。
ちなみにhtmという拡張子がついていれば、どれをブックマーク...
全てのページでcontentsにリンクされているからそれほど問題...
*cgxのインストール簡単マニュアル(linux) [#z27ba659]
-&link(a short installation guide (ASCII),http://www.dhon...
cgxのインストールはほんの少し厄介である。インストールする...
--(
locate libGL
--)
とターミナルに打ち込み、
--(
/usr/lib/libGL.so (libGL.soがlibGL.so.0になって...
/usr/lib/libGLU.so (libGLU.soがlibGL.so.0になって...
--)
と表示されることを確認する。
もしされないようなら、英語マニュアルhttp://www.dhondt.de/...
秋田大の学生なら私こと千田か後藤さんに確認を取るように。
無事表示されたら、
&link(CalculiXのダウンロードのページ,http://www.dhondt.de...
Available downloads for the graphical interface (CalculiX...
''このファイルをディスクの保存する''を選び、自分のユーザ...
次にターミナルを開いて、保存したところにcdで移動し(home/s...
suになる。[suzuki@bri00]$が[root@bri00」#に変わったのを確...
--(
# mv cgx_1.6.bz2 /usr/local
# cd /usr/local
# bunzip2 cgx_1.6.bz2
# mv cgx_1.6 /usr/local/bin/cgx
# chmod ao+rx /usr/local/bin/cgx (chmod 755 でもよさそう)
--)
と順に打ち込む。
でっ、インストールされたか確認するために、
--(
cgx -b dummy.fbd
--)
と打ち込む。cgx not command とかいう表示が出たら下のパス...
もし、cgx: error while loading shared libraries: llibglut...
cgxのwindowのHelpから、マニュアルを見ることができるように...
ダウンロードページのthe documentation (html)
(ファイル名はcgx_1.6.htm.tar.bz2)を実行ファイルと同じ手順...
そしてsuになって
--(
# mv cgx_1.6.htm.tar.bz2 /usr/local
# mv /usr/local
# bunzip2 cgx_1.6.htm.tar.bz2
# tar -xvf cgx_1.6.htm.tar (たくさんのファイルが展開され...
--)
と順に打ち込む。ここで
--(
# cd /usr/local/CalculiX/cgx_1.6/doc/cgx
# ls
--)
でさっき展開したhtmlファイルがあることを確認する。
その上で、
--(
cgx -b dummy.fbd
--)
と打つとcgxが立ち上がるので、
その中にある正方形の枠を避けて左クリックし、開いた小さなw...
Helpまでクリックしたまま移動すると、あらたに
Quick Reference(cgx)、Html Manual(cgx)、Html Manual(ccx)...
Html Manual(cgx)のところに移動しクリックを外すと
w3m(Mozillaは確認中)が開く。
http://www.str.ce.akita-u.ac.jp/~chida/oonuma/cgx1.png
もしccxのhtmlマニュアルが先にインストールされているのであ...
Html Manual(ccx)も見ることができる。
サンプルをインストールしたい場合は、ダウンロードページの
examples(ファイル名cgx_1.6.exa.tar.bz2 )をダウンロードす...
その後はhtmlのマニュアルのインストール方法と同じである
(cgx_1.6.htm.tar.bz2のhtmをexaに変えるだけ)。
*inpファイル [#kdb32a62]
-&link(千田さんのプログラム,http://www.str.ce.akita-u.ac....
--&link(川原の目的に応じて修正中,http://www.str.ce.akita-...
--&link(ぶっとうす場合のプログラム,http://www.str.ce.akit...
--&link(サンドイッチ型のプログラム,http://www.str.ce.akit...
-&link(後藤の,http://www.str.ce.akita-u.ac.jp/~gotou/prog...
-http://www.str.ce.akita-u.ac.jp/ccx/node65.html(アクセス...
-http://www.str.ce.akita-u.ac.jp/cgi-bin/ccx/wiki.cgi?Inp...
-http://www.str.ce.akita-u.ac.jp/~chida/0512.html#i49
*節点番号の振り方 [#u465a4d1]
-http://www.str.ce.akita-u.ac.jp/~chida/0512.html#i36
*任意断面の見方 [#ibf3dcc7]
-このやりかたはバージョン1,7で確認済み。
-3Cの後の要素番号、構成を消す。
-節点(x+1)$\times$(y+1)$\times$z(みたいところ)
-要素x)$\times$y$\times$z(みたいところ)
ex)
--(
-1 20401 1 0STEU #要素番号
-2 23002 23013 23464 23453 23003 ...
--)
座標は消さなくても大丈夫みたい。
このままでも断面をみることができるけど、
cgxの最大・最小の値が全節点になってしまうので、
適宜、要らない節点の変位や応力も削除する。
-変位
--(
1PSTEP
1 100CL101 1.0000000000 3 1 1 -...
DISP 4 1 -5
D1 1 2 1 0 -5
D2 1 2 2 0 -5
D3 1 2 3 0 -5
ALL 1 2 0 0 1ALL
-1 1 0.00000E+00 7.89755E-05 1.26021E-03
-1 2-4.91247E-09 7.89756E-05 1.26021E-03
-1 3-1.44954E-08 7.89762E-05 1.26020E-03
--)
-応力
--(
1PSTEP 1
100CL101 1.0000000000 ...
-4 STRESS 6 1
-5 SXX 1 4 1 1
-5 SYY 1 4 2 2
-5 SZZ 1 4 3 3
-5 SXY 1 4 1 2
-5 SYZ 1 4 2 3
-5 SZX 1 4 3 1
-1 1-4.07109E-01 6.13735E-01 1.40236E+00 1.06976...
-1 2-3.89041E-01 6.21203E-01 1.41357E+00 7.28904...
--)
以下は古い話。
--(
-(x+1)×(y+1)×y(みたいところ)の節点をfrdファイル上で削除す...
-削除箇所は3箇所。
--節点座標(2C→3Cの間)
--要素番号、節点番号(3C→1PSTEP-DISPの間)
--応力(1PSTEP-STRESSの後)
--※後ろから消す場合は節点に注意!
-消すのを間違えると、
--セグメンテーション違反
--梁が空洞になってしまう
ので注意。
--)
*増分ステップの制御など [#b28a841d]
ccxでは、最終荷重を与えて計算すると、そこに至る数ステップ...
荷重増分を小さくして、もっと細かくプロットを取りたい時に...
-1ステップは複数の荷重増分を含む
-INC=100とすると、1ステップの荷重増分を100増分までに制限...
--1ステップの増分数は100よりもずっと少ないので、INC=200と...
-FREQUENCY=1とすると、すべての増分に対する結果を(出力用に...
--デフォールトですべての増分に対する結果が出力されている...
--FREQUENCY=NとするとN番目の増分ごとに結果を記憶する
-計算済みのつりあい点から引き続き計算したいときはRESTART...
--RESTARTをうまく使えば小さい増分ごとの結果を毎回出力させ...
**出力ステップの制御 [#a2115668]
"*TIME POINTS"を使うと、好きなステップで出力することがで...
--(
*TIME POINTS,NAME=T1 [#y1a9251b]
0.0,0.2,0.4,0.6,0.8,1.0
(途中で色々書かれていて......)
*NODE PRINT,NSET=Nall,TIMEPOINTS=T1 [#wedeaa12]
U
*EL PRINT,ELSET=Eall,TIMEPOINTS=T1 [#c241c574]
S
*NODE FILE,TIMEPOINTS=T1 [#w5483666]
U
*EL FILE,TIMEPOINTS=T1 [#yf7250fa]
S
*END STEP [#l613692e]
--)
*pathを追加する [#m52f3bb5]
cgxをインストールして初めて/usr/local/binができた場合、cc...
起動しない事がある。よって、
--(
$ echo $PATH パスが切れ...
/usr/hoge/hoge
$ ls -a ./bash_pr...
$ vi .bash_profile bashの前...
export PATH=/usr/local/bin:$PATH と書き込む
$ source.bash_profile
$ echo $PATH
/usr/hoge/hoge:/usr/local/bin usr/local/binが追加され...
--)
を実行することによって、どのディレクトリでもcgx、ccxがパ...
このパスを指定しないで済むようにすることを''パスを切る''...
**後藤メモ [#e8979662]
g2は私がパスを切ってたってことでは。
-&link(パスを切る,http://www.str.ce.akita-u.ac.jp/cgi-bin...
学生がパスを切るって分からないと思い簡単に書いたのですが。
fパソのように誰もパスを切ってないのに切れてる場合もありま...
念のため(千田)
*エラー [#s851be11]
***ccx1.7で解けた円筒の圧縮がccx2.3で解けない [#u350bfb5]
円筒やダイヤカット円筒の一端を固定して一端を圧縮（または...
ccx1.7では問題なく解けているが、同じinpファイルをccx2.3で...
載荷端側が花びらのように半径方向に互い違いに変位するモー...
ccx2.3にバージョンアップして増えたオプションをちゃんと設...
解けないということかも知れない。ウォーニングメッセージとc...
ちゃんと確認する必要あり。
***ubuntu10.04上でccx1.7を走らせると [#lb48cba0]
というような訳で、ccx1.7のバイナリーをubuntu10.04上で走ら...
--(
./ccx_1.7: error while loading shared libraries: libg2c.s...
--)
これは、g77用のライブラリらしい。
-http://giropenfoam.web.fc2.com/salome/install_Suse/insta...
g77を入れてためしてみるか。
-http://d.hatena.ne.jp/ti10/20100512/1273671639
***nonpositive jacobian [#lfc2ec28]
--(
ERROR in e_c3d: nonpositive jacobian
determinant in element 14
--)
節点座標の1行文字数が非常に長い場合にこのエラーが出る場合...
そして、ccx_1.7では解けてccx_2.3では解けなかったりする。
--(
print'(i6,a,1pe15.3,a,1pe15.3,a,1pe15.3)'&
& ,i+1,',',x(i),',',y(i),',',z(i)
--)
みたいにフォーマットを指定して指数表示で1行を短めに出すと...
など1pdにするとcgxの表示がおかしくなるため、1peにしておく。
*** Factoring the system of equations using spooles ALLO...
C3D8の直方体片持ち梁で、ポアソン比を0.5に設定するとエラー...
ポアソン比0.49999と0.50001ではエラーは出なかった。
*シェル要素による片持ち梁のたわみ解析結果 [#t98ec863]
-&link(結果まとめ,http://www.str.ce.akita-u.ac.jp/cgi-bin...
終了行:
#contents
-http://www.calculix.de/
*ubuntu22へのインストール [#heab231f]
$sudo apt install calculix-ccx calculix-cgx
バージョンはやや古いものにはなるけど。
*近況? [#r848e3fe]
*PrePoMax [#s8b7c532]
-https://aytechlab.com/homefem/prepomax001/
Salomeでunvファイルを作ったのをPrePoMaxで読み込んで、材料...
**Ubuntuへのインストール [#oad3b7ff]
***wine環境のインストール [#s236b285]
PrePoMaxはWindowsソフトなので、まずは、Linux上でWindowsを...
をインストールし、
さらにその環境でNET Frameworkを使えるようにするために、do...
-[[オープンＣＡＥコンサルタントOCSE^2:https://ocse2.com/?...
$ sudo dpkg --add-architecture i386
$ sudo mkdir -pm755 /etc/apt/keyrings
$ sudo wget -O /etc/apt/keyrings/winehq-archive.key http...
$ sudo wget -NP /etc/apt/sources.list.d/ https://dl.wine...
$ sudo apt update
$ sudo apt install --install-recommends winehq-devel
$ sudo apt install winetricks
$ sudo winetricks --self-update
$ winetricks dotnet48
-最後のdotnet48のインストール中にwineが立ち上がり、「wine...
-「マイクロソフトソフトウェア追加ライセンス条項」みたいな...
-「Wineの設定」ウィンドウが出てきたら、OKをクリック。2回...
-「Microsoft.NET Framework」の警告（正しく実行されない可...
-「.NET Framework 4.8セットアップ」（追加ライセンス条項）...
□□□□(A)みたいなチェックボックスがたぶん「同意する」なので...
-「.NET Framework 4.8がインストールされました」というウィ...
-「インストールを完了するには...再起動する必要があります...
***PrePoMaxのインストール [#ma0cf087]
-~/.wine/drive_cの中にexeというディレクトリを作り、wineで...
この中にソフト名がわかるサブディレクトリーを作って、その...
いいのではないか。ということで、PrePoMaxに関しては、~/.wi...
ディレクトリを作っておく。
-[[PrePoMaxのダウンロードサイト:https://prepomax.fs.um.si...
から最新バージョン（または動作確認済みのバージョン）のzip...
~/.wine/exe/prepomax にダウンロード
-unzip 'PrePoMax v1.4.0.zip'で展開。
-cd 'PrePoMax v1.4.0'で展開されたディレクトリーに移動。
-wine PrePoMax.exeを実行し、PrePoMaxが起動することを確認。
-Ubuntu(MATE)の上部パネルで右クリックし、「パネルへ追加」...
「カスタム・アプリケーションのランチャ」を選択して「追加...
--「ランチャの作成」ウィンドウが開くので、名前：PrePoMax,...
/usr/bin/wine '/home/kouzou/.wine/drive_c/exe/prepomax/P...
--ランチャーにショートカットアイコンが作成される。アイコ...
PrePoMaxを起動して、helpのabout辺りから、PrePoMaxのロゴの...
~/.wine/drive_c/exe/PrePoMax辺りに適当な名前(aikon.pngと...
ショートカットアイコンを右クリックして、プロパティーから...
**Windows版（PrePoMax・OptiMax） [#a6c64618]
-http://lace.fs.uni-mb.si/wordpress/borovinsek/
--雰囲気はこんな感じ → https://www.youtube.com/playlist?l...
***PrePoMax（プリポスト付き） [#wc54148c]
-http://lace.fs.uni-mb.si/wordpress/borovinsek/?page_id=41
***OptiMax（最適化） [#j3c8e943]
-http://lace.fs.uni-mb.si/wordpress/borovinsek/?page_id=104
--https://www.slideshare.net/OpenCAE_kanto_fluid/optimax-...
--Linuxで似たようなことをしようとするなら、[[Dakota:https...
--[[パラメータ最適化ツール DakotaとOpenMDAOについて:https...
**Windows版（CalculiXLauncher） [#g78cbed6]
-http://www.calculixforwin.com/?page_id=31
--ダウンロードは&link(ここ,https://sourceforge.net/projec...
--起動は、Launcher33x64.exe をクリック。
-ソルバー(ccx)は普通に動くが、プリ・ポストプロセッサ(cgx)...
-CCX WIZARD! はLinuxのみで使えるとのこと。
-UNV→inpコンバーターも使える模様。
**FreeCADでCalculiXをソルバーとして解析できる？ [#ve8353ba]
-http://naga-tsuzuki.sblo.jp/article/179282504.html
-http://naga-tsuzuki.sblo.jp/category/4426804-1.html
**SalomeとCalculiXを連携できそうなCalculiXLaucher（要調査...
-http://calculixforwin.blogspot.jp/2015/05/calculix-launc...
-http://calculixforwin.blogspot.jp/2016/10/getting-starte...
Salomeで作ったメッシュファイルをCalculiXLauncherでinpファ...
**直交異方性材料 [#qad7b617]
今まで通りのひずみ－応力関係の逆関係を入力する方法
*ELASTIC,TYPE=ORTHOROPIC [#nc477d23]
D1111,D1122,D2222,D1133,D2233,D3333,D1212,D1313
D2323
の他にも，Code_Asterと同様な入力方法↓もできるようになっ...
*ELASTIC,TYPE=ENGINEERING CONSTANTS [#r0c7a606]
E1,E2,E3,nu12,nu13,nu23,G12,G13,
G23
**umat [#u6aa51b9]
-https://www.xsim.info/articles/CalculiX/ccx-doc/node315....
**unvからinpへの変換 [#g756d177]
-http://freeplanets.ship.jp/FEM/Calculix/Calculix.html#mo...
-http://caelinux.org/wiki/index.php/Proj:UNVConvertからMe...
-unzipで解凍
--(
python unv2abaqus.py Mesh_1.unv mesh
--)
でmesh.inpができてる。
**インストール [#a4732bc2]
CalculiXは、バージョン1.8以降は実行ファイルを展開しただけ...
SPOOLES.2.2やARPACKの導入も必要となり、インストールが難し...
恐らくインターフェースやプリ・ポストプロセッサーなどがGUI...
思われる
[[DEXCS:http://dexcs.gifu-nct.ac.jp/pukiwiki/index.php?%C...
は、ADVENTUREを利用していて、今のところ動的解析などには対...
いずれ、
&link(CalculiXを導入して動的解析にも対応する方向性,http:/...
みたいなので、様子を見たい。
CalculiXのrpm版もバージョン2.0以降のものも検索すればそれ...
(64ビット版が主流で32ビット版は少ないので注意)、依存性が...
Vine等へのインストールは簡単ではない。
&link(CAElinux,http://caelinux.com/CMS/index.php)の64bit...
プレインストールされてるようなので、64bitパソコンにCAElin...
計算サーバーとするか。
***Ubuntu 16.04 LTS へのインストール [#v9c88d88]
Ubuntu16.04LTSを(近藤のPCに)インストールして、色々と試...
Ubuntu16.04LTSの&link(公式パッケージ,http://packages.ubun...
--(
$ sudo apt-get install calculix-ccx
--)
これだけでいつも通りの計算できるようだが、cgxはパッケージ...
-ccx_2.6との違い
--荷重や強制変位の方向の書き方が変わった。例えば、&link(c...
--(
*CLOAD [#heba7041]
** LAST, 2, 6.6666666666666666E-002 [#lf4b385a]
LAST, 2, 2, 6.6666666666666666E-002
--)
後藤ちゃちゃ(17/3/29)**はコメント文では？
---近藤(17/3/29)そこの行でエラーが出るので、それをコメン...
---後藤(17/3/30)そうですか。CalculiXのinpファイルはABAQUS...
--接触解析のとき、接触面の組み合わせの他にも、面同士の種...
--(
*CONTACT PAIR,INTERACTION=SI1,TYPE=NODE TO SURFACE [#w55b...
--)
--というように、TYPEが必要。上記以外にも、スレーブ面がマ...
-後藤ちゃちゃ(2016/11/14)：そういうことなら、2017年度以降...
Ubuntu12.04は2017/4までのサポートだし、それまでCAELinuxの...
**CAELinux [#l4d02da6]
&link(CAElinux,http://caelinux.com/CMS/index.php)はUbuntu...
インストール後にアップグレードしたらccxも最新版の2.3にア...
これはいいかもしれない。
Ubuntuのクライアント版はサーバー系ツールのセットアップが...
めんどくさいような気がしていたが、CAELinuxにはsshdがプリ...
サービス管理(system-config-services)は、synapticの一覧に...
なんか他のツールがあるのかな。
*コンクリート [#p7ee6526]
コンクリートのように圧縮のみに抵抗する材料も扱える。
-例えば、
--(
*MATERIAL,NAME=COMPRESSION_ONLY [#u67912f5]
*USER MATERIAL,CONSTANTS=2 [#re38db97]
14.e3, 1.e3
--)
-と書くと、ヤング率が14MPaで、許容する最大引張応力が1MPa...
--ところで、ポアソン比を入力するところがないのはどうなっ...
-COMPRESSION_ONLYの他にも特殊な材料がたくさんある(/opt/Ca...
-参考：http://open-shelf.appspot.com/CalculiX/ccx-doc/nod...
*面の定義 [#bab7422b]
マニュアル(ccx_2.6.1.pdf)のp.367に要素ごとの面番号の定義...
--(
面1:節点1,2,3,4
面2:節点5,8,7,6
面3:節点1,5,6,2
面4:節点2,6,7,3
面5:節点3,7,8,4
面6:節点4,8,5,1
--)
というように定義されている。例えば、接触解析とかで、要素...
--(
*SURFACE,NAME=men [#ff6f59b8]
1,S6
2,S6
3,S6
4,S6
--)
みたいにして使う。
ちなみに、接触解析をするとき、マスター面の定義は必ず上記...
*ヒンジ [#hf0766e9]
**立体要素と板要素を組み合わせる場合 [#r4ff8e4e]
CAELinuxの場合、/opt/CalculiX_2.6.1/ccx_2.6.1.pdfがマニュ...
ここのp.96によると、「立体要素と板要素」あるいは「立体要...
以下、その例題(/opt/CalculiX_2.6.1/ccx_2.6.1/test にある)
-solidshell1.inp:ただ、立体と板の節点を共有させただけのモ...
-solidshell2.inp:立体要素に重ねて板要素を入れることで剛結...
**多点拘束を用いる場合 [#z8db4833]
-&link(hinzi.inp,http://www.str.ce.akita-u.ac.jp/~gotouha...
--折り目部分を重節点にしたうえで"*EQUATION"を使い、互いの...
**立体要素だけ [#r910f489]
特別なことをしなくても、立体要素同士を線状に接点共有さ...
この方法なら、折り畳み円筒を一定厚にし、かつ蝶番を再現...
-&link(hinzi2.inp,http://www.str.ce.akita-u.ac.jp/~gotouh...
**梁の場合のヒンジ節点 [#m5b4fa4d]
-&link(メモ1,http://www.str.ce.akita-u.ac.jp/~gotou/kouzo...
-&link(メモ2,http://www.str.ce.akita-u.ac.jp/~gotou/kouzo...
*ccx_2.5 [#fc2ba6bc]
ざっといじっていてわかったこと(変更点?)を箇条書きで
-例の法線の向きが互い違いになっているダイヤカット&link((d...
--ちなみに同じinpファイル&link((daiya6h.inp),http://www.s...
---2.5
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/test/daiy...
&link(daiya6h.dat,http://www.str.ce.akita-u.ac.jp/~gotouh...
---2.3
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/test/daiy...
&link(daiya6h23.dat,http://www.str.ce.akita-u.ac.jp/~goto...
---ある載荷節点(節点番号17から20)の変位
2.5
--(
17 9.337687E-11 -1.359864E-23 -1.388775E-09
18 1.085580E-23 9.337687E-11 -1.388775E-09
19 -9.337687E-11 -8.762945E-24 -1.388775E-09
20 5.218346E-24 -9.337687E-11 -1.388775E-09
--)
2.3
--(
17 9.3377E-11 -2.7324E-23 -1.3888E-09
18 7.1344E-24 9.3377E-11 -1.3888E-09
19 -9.3377E-11 -1.1550E-23 -1.3888E-09
20 -6.2168E-24 -9.3377E-11 -1.3888E-09
--)
-上の画像からわかるように、2.5で解いて出来たfrdファイルを...
inpファイルで設定した厚さを反映したモデルを見ることができ...
--cgxで節点を見ると
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/test/25se...
のようになっていて、これは&link(マニュアル(the documemtat...
--ということはccx_2.5で解くと、frdファイルに新しい節点が...
-ちなみに&link(このへん,http://www.str.ce.akita-u.ac.jp/c...
*ccx_2.6.1をCAELinux2011で実行できるように（まだできない...
CalculiXの作者guido dhondtさんに問い合わせたところ、外向...
解けるものが2.3では解けない問題は、2.6では解決されている...
ccx_2.6.1を使うのがよさそうだ。
ccx_2.6.1のa Linux executable(ccx_2.6.1.tar.bz2)を適当な...
--(
$ bunzip2 ccx_2.6.1.tar.bz2
$ tar xvf ccx_2.6.1.tar
--)
で展開。CalculiX/ccx_2.6.1/src/の中にある実行ファイルccx_...
libgfortran3がもしインストールされてなければ、synaptic等...
--(
$ ls -l /usr/lib/libgfortran*
--)
とやって、/usr/libにあるlibgfortranのバージョン（とリンク...
恐らく、/usr/lib/libgfortran.so.3.0.0があって、
/usr/lib/libgfortran.so.3（青色）からリンクが張られている。
この状態で ccx_2.6.1を実行すると、恐らく以下のエラーが出...
--(
error while loading shared libraries: libgfortran.so.2: c...
--)
そこで、以下のように/usr/libに移動してからlibgfortran.so....
--(
sudo ln -s libgfortran.so.3.0.0 libgfortran.so.2
--)
これでccx_2.6.1を実行すると、Usage: CalculiX.exe -i jobna...
--(
ccx_2.6.1 hoge
ccx_2.6.1: symbol lookup error: ccx_2.6.1: undefined symb...
--)
*1.7と2.3の比較 [#x1bae6f5]
-datファイル
単純に&link(波型構造の.inpファイル,http://www.str.ce.akit...
&link(ver.1.7,http://www.str.ce.akita-u.ac.jp/~gotouhan/e...
--オーダーが大きく違う
--2.3だと指定したはずの拘束点がいくつか拘束されていない(e...
-上の計算で得られた.frdファイルの図
http://www.str.ce.akita-u.ac.jp/~gotouhan/emura/daiya/1-7...
*マニュアル [#iaba6b2d]
-&link(アクセス制限あり,http://www.str.ce.akita-u.ac.jp/c...
ccxのマニュアル(日本語版)を見つけたので、ここに載せました...
[[日本語版:https://www.xsim.info/articles/CalculiX/ccx-do...
[[英語版:http://www.dhondt.de/ccx_2.20.pdf]]
*inpファイルの書き方 [#sf1da85e]
-http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html
*サンプルファイル [#p93aad99]
-http://web.mit.edu/xenophon/MacData/afs/athena/software/...
**コメント [#na23a417]
--(
** こんなふうに ** を最初に書くとコメントが書ける [#ybc94...
** つまり ** の行はデータ入力の際に無視される [#bffc99c9]
--)
**要素の節点番号 [#c3aa4909]
inpファイルで、まず各節点と節点座標を定義してから、
各要素ごとに要素番号とその要素の各節点の節点番号を設定し...
要素の種類ごとに1要素の節点番号の並べ順には決まりがあるの...
マニュアルやサンプルの節点番号の振り方を参照すること。
例えば、8節点シェル要素(S8)では、頂点の節点番号を4つ並べ...
辺中点の節点番号を4つ並べなければならないとか。
これを間違うと、変な形になる。
**拘束条件 [#e1708cfb]
***複数の節点にまとめて同じ拘束条件を与える場合 [#zfc4c4cd]
--(
*NSET,NSET=FIX [#eeeef362]
** ここに拘束したい節点番号とコンマを1行ずつ並べる [#aa85...
** 例えば拘束節点が1と2だったら [#c7ac0403]
1,
2,
** みたいに [#u5417dea]
*BOUNDARY [#zc4e42a5]
** 以下のFIX,に続けて拘束変位1=x,2=y,3=zをコンマ区切りで ...
FIX,1,3
** 1から3(xからz)を拘束 [#eddf2c83]
--)
上記の例を、
--(
*BOUNDARY [#f0883ca8]
FIX,1
FIX,2
FIX,3
--)
みたいに分けて書いても良いっぽい
***1節点ずつ拘束条件を与える場合 [#afea8f27]
--(
*BOUNDARY [#ef691e58]
1, 3, 0.
2, 1, 3, 0.
** 節点番号1の3(z)を拘束 [#a104396a]
** 節点番号2の1から3(xからz)を拘束 [#z6cf218e]
--)
**材料定数 [#mc506e0e]
--(
*MATERIAL,NAME=EL [#p3c4460f]
*ELASTIC [#y26c4896]
** 以下にヤング率(MPa)とポアソン比 [#c1c9337e]
210000., 0.300000012
--)
**載荷条件 [#u7bd11be]
***複数の節点にまとめて同じ荷重条件を与える場合 [#s5e6554b]
--(
*NSET,NSET=SAIKA [#s2d94d1b]
** ここに載荷したい節点番号とコンマを1行ずつ並べる [#r17b...
** 例えば載荷節点が3と4だったら [#pcce862b]
3,
4,
** みたいに [#d1122947]
*STEP [#zb71651c]
*STATIC,SOLVER=SPOOLES [#t9a84e1a]
*CLOAD [#h94d6348]
** 以下のSAIKAに続けて載荷方向、載荷荷重を書く [#z72c0b8a]
** 例えば載荷方向が3=zで、載荷荷重が1/載荷節点数なら [#ab...
SAIKA,3, 0.5
--)
***1節点ずつ荷重条件を与える場合 [#ja2a05ed]
--(
*CLOAD [#w98db071]
1, 2, 10.
2, 3, 20.
** 節点番号1の2(y)方向に10MN載荷 [#xe96bb41]
** 節点番号2の3(z)方向に20MN載荷 [#x169b9fa]
--)
*ひずみについて [#pf595004]
inpファイルで最後の方に、
--(
*EL FILE,POSITION=AVERAGED AT NODES [#icb5206f]
S,E,PE
--)
っていれる。
ちなみに、Sは応力、Eはひずみ、PEって入れると相当塑性ひず...
ほかにもいろいろ出力できそう。
詳しくはマニュアルの
EL FILEとかNODE FILE,EL PRINT,NODE PRINT 参照..
*座屈解析 [#u2a9680b]
**datファイル内の座屈荷重の読み方 [#qe40364e]
-inpファイルのCLOADで載荷節点の倍率bを書き込む。
-datファイルのEIGEN VALUEのとこに書き込まれる荷重$P_{dat}...
-$P_{cr}=bnP_{dat}$
-つまり、倍率$b$を$1/n$としておけば、datファイルの荷重$P_...
ということでいいかな?(次項の結果と整合するかな)
-比較しやすいように、datファイルの座屈荷重を頭にすれば
--$P_{dat}=\frac{P_{cr}}{nb}$
-注意：inpファイルをgfortranで作る際、上記の倍率($b=1/n$)...
桁数が多くなりすぎて、-E003の部分が読み込まれず、座屈荷重...
**inpのCLOADと、datの座屈荷重との関係 [#mf10b705]
いつもは、CLOADで(1/載荷節点数)と入力したものを、座屈荷重...
もし10/載荷節点数とかにしたら、座屈荷重はどう変わるのか見...
,CLOAD,0.001/n,0.01/n,0.1/n,1/n,10/n,100/n
,理論値との誤差,座屈せず,3%,3%,3%,3%,座屈せず
n:載荷節点数
10/nの座屈荷重は、1/nの座屈荷重の0.1倍になる。0.1/nにすれ...
だから、inpのCLOADで0.1/nで入力したら、datファイルの座屈...
ただ、サイズによって座屈が出ない問題も、(相似形状で10倍と...
ACCURACYのところを0.01ではなく、0.1とか1とかにすれば解け...
**モード図の描かせ方など [#wec7aea4]
***cgxのアニメ機能を利用する場合 [#u5791961]
-cgxでfrdファイルを開き、Datasetsを選ぶと、計算した座屈モ...
-L101とか、値が0になっているところは自明解なので、最初に...
-更にDISPを選択
-改めてDatasets→Entitiy→Allを選択すると、座屈モードがまず...
-次にAnimateのStartを選択すると、座屈モード変形がアニメ再...
-画像がほしい瞬間のところでマウスのポインターを枠外に出し...
-モードの倍率を大きくしたければ、Animate→Tune Valueで、Tu...
***倍率を調整して手動で変形モードを加算したfrdファイルが...
&link(tasu.f系のプログラム,http://www.str.ce.akita-u.ac.j...
-ccxkataiza.fの直方体要素ならば、ccx_1.7+tasudatw.fでモー...
ダイヤカットをやったシェル要素では、ccx_1.7*tasudatw.fで...
ccx_1.5+tasu.fの昔のやり方でモード図を描いていたというよ...
確認してほしい。
円筒ダイヤカットも伸縮できる円筒折り紙も、圧縮でも引っ張...
起こせていたが、これはバージョンに依存するか？ 例えば、cc...
ccx_1.7では、no convergence.....maybe no buckling occurs ...
事情があるかどうか確認。
--シェル要素でもccx_1.7+tasudatw.fの組み合わせでモード図...
**l=1,m=16,n=2 [#t0674712]
,半径(m),高さ(cm),ヤング率,ポアソン比,厚さ[m],座屈モード,...
,1.d0/2.d0/pi,50,70000,0.300000012,0.11,2,0.06,0.149253731
-これだとccx_1.7では「job finished」になるが、ccx_1.5だと...
--「job finished」は荷重1E-1～E-4の間。
-ccx_1.5では1E-5から「セグメンテーション違反です」とおこ...
--1E100以上、1E-5以下でもおこられる。その他は「*ERROR: to...
-厚さによって荷重が同じでも「no convergence for the buckl...
*バージョン比較 [#j3419520]
バージョンによって答えが違った例をここに書き込む。
-2枚の板をくっつけた折れ板を引っ張って広げると、ccx_1.5で...
too many cutbacksエラーが出る。ccx_1.7線形では解ける。
ccx_1.7非線型ではtoo many cutbacksエラー([[稲荷の卒論日誌...
*ウェブ上の情報など [#faf218d1]
-http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html
-http://www.waka.kindai.ac.jp/tea/shibue/abaqus-inpexp.html
-http://www.gsic.titech.ac.jp/~ccwww/tebiki/
-http://freecaetester.blog62.fc2.com/blog-date-200803.html
*cgxでinpファイルをみる。 [#je66249e]
--(
cgx -c hoge.inp
--)
でみることができる。
*cgxで変形図やモード図を描かせる [#u645a37d]
&link(ここ,http://www.str.ce.akita-u.ac.jp/~gotou/program...
&link(tasu.f,http://www.str.ce.akita-u.ac.jp/~gotou/progr...
使ってデータを整形する。具体的な方法はプログラムの中に書...
*CADとの連携 [#k24dd8c7]
-Available add-onsを使うとvda(ドイツ自動車工業会策定?のフ...
--(
./acis2fbd.pl < hoge.sat >hoge.fbd
--)
で変換できるようだが、サンプルのsatファイルは変換できるが...
変換できなかった。
-&link(Linuxにインストールできるフリーの3次元CAD,http://p...
簡単にインストールできて、出力フォーマットに各種形式を選...
&link(gCAD3D,http://www.gcad3d.org/)で、
DXF, IGES, VRML1, VRML2, STL, OBJ, TESS, NC-ISOにエクスポ...
で、dxf辺りからvdaやacisへのフリーの変換ツールはないかと...
&link(GID,http://gid.cimne.upc.es/)のアカデミック版ででき...
でも、あわよくば完全なオープンソースだけで処理したいので...
**cgxをプリプロセッサーとして使う [#t34ff8af]
-http://www.dhondt.de/tutorial.txt
-[[川原さんによる和訳にはパスワードによるアクセス制限を付...
ちょっとめんどくさいような気もするのですが、
プリプロセッサーを使わないことにしたとして、
CADやモデリングツールの出力から節点情報を抜き出す
プログラムを作る際に、メッシュを指定する節点の並び順をccx...
うまくやれなくて節点情報だけを頼りに自力でプログラムで対...
cgxをプリプロセッサーとして利用してfbdを作るのもそう悪く...
いや&link(blenderのメッシュ情報,http://www.str.ce.akita-u...
**blender [#ofee5e9f]
CalculiXのモデリングに使える??
-[[blenderメモ]]
-http://www.usennet.ne.jp/~kiyoshi/manual/blenderTips.html
-http://blender.jp/
-http://bear32.blogspot.com/
-http://bear32.blogspot.com/2006/11/blender.html
-http://f11.aaa.livedoor.jp/~hige/
-http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html
-http://www.geocities.jp/detoni2303/tips/tips_036.htm ...
**各種フォーマット [#c0843951]
-&link(dxf,http://www.autodesk.com/techpubs/autocad/acad2...
-&link(IGES,http://afsoft.jp/cad/p08.html)
*ccxのインストール簡単マニュアル(linux) [#jd231ab6]
**ccx_2.6 [#i63edd90]
-http://www.libremechanics.com/?q=node/9
**ccx_2.3 [#c89c57a5]
ccx_1.8以降のコンパイルはうまくいっていないが、SPOOLES.2....
-http://pugpug.lolipop.jp/yii-1.0/demos/mocapapas-blog/ta...
***SPOOLESのインストール [#t7915585]
ccxのサイト経由で、
spooles.2.2.tgzをダウンロードして、
/usr/local/SPOOLES.2.2を作成してからそこにコピーして展開。
(アーカイブファイルがディレクトリ内に納めてから圧縮されて...
ばらばらのファイルのまま圧縮されているかは、ファイルブラ...
/usr/local/SPOOLES.2.2/内にファイルが展開される。
Make.inc内の
--(
# CC = gcc
CC = /usr/lang-4.0/bin/cc
--)
の行を
--(
CC = gcc
# CC = /usr/lang-4.0/bin/cc
--)
に修正して、make libを実行。spooles.aが生成されれば成功。
***ARPACKのインストール [#s2cccdd9]
ccxのサイト経由で、
arpack96.tar.gzとpatch.tar.gzをダウンロードして、/usr/loc...
tar xvzfで展開すると、/usr/local/ARPACK(が生成されてその...
/usr/local/ARPACKの中に入って、ARmake.inc内の以下の箇所を...
--(
#home = $(HOME)/ARPACK
home = /usr/local/ARPACK
(中略)
#MAKE = /bin/make
MAKE = make
(中略)
#SHELL = /bin/sh
SHELL = sh
--)
make libすると、libarpack_SUN4.aが生成される。
この状態で、
/usr/local/CalculiX/ccx_2.3/src
に入り、makeしてccx_2.3をコンパイルすると、
libarpack_INTEL.aがないというエラーが出るので、
libarpack_SUN4.aをlibarpack_INTEL.aにコピーする。
これで、ccx_2.3をコンパイルしても
--(
../../../ARPACK/libarpack_INTEL.a(xerbla.o): In function ...
xerbla.f:(.text+0xe): undefined reference to `s_wsfe'
xerbla.f:(.text+0x29): undefined reference to `do_fio'
xerbla.f:(.text+0x44): undefined reference to `do_fio'
xerbla.f:(.text+0x49): undefined reference to `e_wsfe'
xerbla.f:(.text+0x5d): undefined reference to `s_stop'
--)
みたいなエラーが大量に出る。64bitパソコンでも同様。
/usr/local/ARPACK/UTIL内のsecond.fの
--(
EXTERNAL ETIME
--)
の行をコメントアウトしてARPACKをコンパイルしなおせばよい...
記述も見つけてためしてみたが、エラーの数は減ったような気...
うまくいかない。
-http://freecaetester.blog62.fc2.com/blog-entry-238.html
-http://www.kkaneko.com/rinkou/linux/arpack.html
**ccx_1.8 [#h4eb2761]
ccx_1.8のバイナリーは、Vine LinuxでもUbuntuでも動かない。
--(
bash: /usr/local/bin/ccx_1.8: cannot execute binary file
--)
のエラー。
ソースからmakeするには、ARPACKとSPOOLES.2.2をインストール...
調査中。ARPACKとSPOOLES.2.2のライブラリを展開して/usr/loc...
ccx_1.8をmakeしようとすると、
--(
make: *** `ccx_1.8' に必要なターゲット `../../../SPOOLES....
--)
のエラー。
&link(CalculiX のmake,http://mogura7.zenno.info/~et/xoops...
-http://freecaetester.blog62.fc2.com/
--http://freecaetester.blog62.fc2.com/blog-date-20090316....
--http://freecaetester.blog62.fc2.com/blog-date-20090317....
**ccx_1.7以前 [#w28b88aa]
-&link(a short installation guide (ASCII),http://www.dhon...
&link(CalculiXのダウンロードのページ,http://www.dhondt.de...
Available downloads for the solver (CalculiX CrunchiX: cc...
するとダウンロードのwindowが開くので、''このファイルをデ...
ちなみに保存されるファイルの名前はccx_1.6.tar.bz2であるこ...
次にターミナルを開いて、保存したところにcdで移動し(home/s...
suになる。[suzuki@bri00]$\$$が[root@bri00」#に変わったの...
--(
# mv ccx_1.6.tar.bz2 /usr/local
# cd /usr/local
# bunzip2 -c ccx_1.6.tar.bz2 | tar -xf -
# cd CalculiX/ccx_1.6/src
# mv ccx_1.6 /usr/local/bin
--)
と順番に打ち込む。
#は始めから表示されているので、くれぐれも# mv ccx_1.6とは...
これだけでccxがどこのディレクトリでも使えるようになる。
サンプルをインストールしたい場合には、test examples(ファ...
マニュアルはthe documentationの中から、TeX、html、postscr...
マニュアルはTeXとhtmlの二つをダウンロードしたほうが良い。
TeXをダウンロードする利点は、印刷がきれいにできる
、ファイルが一つだけだから、
viで文中のキーワードが簡単に検索できるということが挙げら...
htmlの利点は言うまでも無く、インターネットと同じ環境で見...
この場合、ブラウザのファイル→開くで/usr/local/Calculix/cc...
その中のccx.htmを開き、ブックマークしておくのがいいだろう。
ちなみにhtmという拡張子がついていれば、どれをブックマーク...
全てのページでcontentsにリンクされているからそれほど問題...
*cgxのインストール簡単マニュアル(linux) [#z27ba659]
-&link(a short installation guide (ASCII),http://www.dhon...
cgxのインストールはほんの少し厄介である。インストールする...
--(
locate libGL
--)
とターミナルに打ち込み、
--(
/usr/lib/libGL.so (libGL.soがlibGL.so.0になって...
/usr/lib/libGLU.so (libGLU.soがlibGL.so.0になって...
--)
と表示されることを確認する。
もしされないようなら、英語マニュアルhttp://www.dhondt.de/...
秋田大の学生なら私こと千田か後藤さんに確認を取るように。
無事表示されたら、
&link(CalculiXのダウンロードのページ,http://www.dhondt.de...
Available downloads for the graphical interface (CalculiX...
''このファイルをディスクの保存する''を選び、自分のユーザ...
次にターミナルを開いて、保存したところにcdで移動し(home/s...
suになる。[suzuki@bri00]$が[root@bri00」#に変わったのを確...
--(
# mv cgx_1.6.bz2 /usr/local
# cd /usr/local
# bunzip2 cgx_1.6.bz2
# mv cgx_1.6 /usr/local/bin/cgx
# chmod ao+rx /usr/local/bin/cgx (chmod 755 でもよさそう)
--)
と順に打ち込む。
でっ、インストールされたか確認するために、
--(
cgx -b dummy.fbd
--)
と打ち込む。cgx not command とかいう表示が出たら下のパス...
もし、cgx: error while loading shared libraries: llibglut...
cgxのwindowのHelpから、マニュアルを見ることができるように...
ダウンロードページのthe documentation (html)
(ファイル名はcgx_1.6.htm.tar.bz2)を実行ファイルと同じ手順...
そしてsuになって
--(
# mv cgx_1.6.htm.tar.bz2 /usr/local
# mv /usr/local
# bunzip2 cgx_1.6.htm.tar.bz2
# tar -xvf cgx_1.6.htm.tar (たくさんのファイルが展開され...
--)
と順に打ち込む。ここで
--(
# cd /usr/local/CalculiX/cgx_1.6/doc/cgx
# ls
--)
でさっき展開したhtmlファイルがあることを確認する。
その上で、
--(
cgx -b dummy.fbd
--)
と打つとcgxが立ち上がるので、
その中にある正方形の枠を避けて左クリックし、開いた小さなw...
Helpまでクリックしたまま移動すると、あらたに
Quick Reference(cgx)、Html Manual(cgx)、Html Manual(ccx)...
Html Manual(cgx)のところに移動しクリックを外すと
w3m(Mozillaは確認中)が開く。
http://www.str.ce.akita-u.ac.jp/~chida/oonuma/cgx1.png
もしccxのhtmlマニュアルが先にインストールされているのであ...
Html Manual(ccx)も見ることができる。
サンプルをインストールしたい場合は、ダウンロードページの
examples(ファイル名cgx_1.6.exa.tar.bz2 )をダウンロードす...
その後はhtmlのマニュアルのインストール方法と同じである
(cgx_1.6.htm.tar.bz2のhtmをexaに変えるだけ)。
*inpファイル [#kdb32a62]
-&link(千田さんのプログラム,http://www.str.ce.akita-u.ac....
--&link(川原の目的に応じて修正中,http://www.str.ce.akita-...
--&link(ぶっとうす場合のプログラム,http://www.str.ce.akit...
--&link(サンドイッチ型のプログラム,http://www.str.ce.akit...
-&link(後藤の,http://www.str.ce.akita-u.ac.jp/~gotou/prog...
-http://www.str.ce.akita-u.ac.jp/ccx/node65.html(アクセス...
-http://www.str.ce.akita-u.ac.jp/cgi-bin/ccx/wiki.cgi?Inp...
-http://www.str.ce.akita-u.ac.jp/~chida/0512.html#i49
*節点番号の振り方 [#u465a4d1]
-http://www.str.ce.akita-u.ac.jp/~chida/0512.html#i36
*任意断面の見方 [#ibf3dcc7]
-このやりかたはバージョン1,7で確認済み。
-3Cの後の要素番号、構成を消す。
-節点(x+1)$\times$(y+1)$\times$z(みたいところ)
-要素x)$\times$y$\times$z(みたいところ)
ex)
--(
-1 20401 1 0STEU #要素番号
-2 23002 23013 23464 23453 23003 ...
--)
座標は消さなくても大丈夫みたい。
このままでも断面をみることができるけど、
cgxの最大・最小の値が全節点になってしまうので、
適宜、要らない節点の変位や応力も削除する。
-変位
--(
1PSTEP
1 100CL101 1.0000000000 3 1 1 -...
DISP 4 1 -5
D1 1 2 1 0 -5
D2 1 2 2 0 -5
D3 1 2 3 0 -5
ALL 1 2 0 0 1ALL
-1 1 0.00000E+00 7.89755E-05 1.26021E-03
-1 2-4.91247E-09 7.89756E-05 1.26021E-03
-1 3-1.44954E-08 7.89762E-05 1.26020E-03
--)
-応力
--(
1PSTEP 1
100CL101 1.0000000000 ...
-4 STRESS 6 1
-5 SXX 1 4 1 1
-5 SYY 1 4 2 2
-5 SZZ 1 4 3 3
-5 SXY 1 4 1 2
-5 SYZ 1 4 2 3
-5 SZX 1 4 3 1
-1 1-4.07109E-01 6.13735E-01 1.40236E+00 1.06976...
-1 2-3.89041E-01 6.21203E-01 1.41357E+00 7.28904...
--)
以下は古い話。
--(
-(x+1)×(y+1)×y(みたいところ)の節点をfrdファイル上で削除す...
-削除箇所は3箇所。
--節点座標(2C→3Cの間)
--要素番号、節点番号(3C→1PSTEP-DISPの間)
--応力(1PSTEP-STRESSの後)
--※後ろから消す場合は節点に注意!
-消すのを間違えると、
--セグメンテーション違反
--梁が空洞になってしまう
ので注意。
--)
*増分ステップの制御など [#b28a841d]
ccxでは、最終荷重を与えて計算すると、そこに至る数ステップ...
荷重増分を小さくして、もっと細かくプロットを取りたい時に...
-1ステップは複数の荷重増分を含む
-INC=100とすると、1ステップの荷重増分を100増分までに制限...
--1ステップの増分数は100よりもずっと少ないので、INC=200と...
-FREQUENCY=1とすると、すべての増分に対する結果を(出力用に...
--デフォールトですべての増分に対する結果が出力されている...
--FREQUENCY=NとするとN番目の増分ごとに結果を記憶する
-計算済みのつりあい点から引き続き計算したいときはRESTART...
--RESTARTをうまく使えば小さい増分ごとの結果を毎回出力させ...
**出力ステップの制御 [#a2115668]
"*TIME POINTS"を使うと、好きなステップで出力することがで...
--(
*TIME POINTS,NAME=T1 [#y1a9251b]
0.0,0.2,0.4,0.6,0.8,1.0
(途中で色々書かれていて......)
*NODE PRINT,NSET=Nall,TIMEPOINTS=T1 [#wedeaa12]
U
*EL PRINT,ELSET=Eall,TIMEPOINTS=T1 [#c241c574]
S
*NODE FILE,TIMEPOINTS=T1 [#w5483666]
U
*EL FILE,TIMEPOINTS=T1 [#yf7250fa]
S
*END STEP [#l613692e]
--)
*pathを追加する [#m52f3bb5]
cgxをインストールして初めて/usr/local/binができた場合、cc...
起動しない事がある。よって、
--(
$ echo $PATH パスが切れ...
/usr/hoge/hoge
$ ls -a ./bash_pr...
$ vi .bash_profile bashの前...
export PATH=/usr/local/bin:$PATH と書き込む
$ source.bash_profile
$ echo $PATH
/usr/hoge/hoge:/usr/local/bin usr/local/binが追加され...
--)
を実行することによって、どのディレクトリでもcgx、ccxがパ...
このパスを指定しないで済むようにすることを''パスを切る''...
**後藤メモ [#e8979662]
g2は私がパスを切ってたってことでは。
-&link(パスを切る,http://www.str.ce.akita-u.ac.jp/cgi-bin...
学生がパスを切るって分からないと思い簡単に書いたのですが。
fパソのように誰もパスを切ってないのに切れてる場合もありま...
念のため(千田)
*エラー [#s851be11]
***ccx1.7で解けた円筒の圧縮がccx2.3で解けない [#u350bfb5]
円筒やダイヤカット円筒の一端を固定して一端を圧縮（または...
ccx1.7では問題なく解けているが、同じinpファイルをccx2.3で...
載荷端側が花びらのように半径方向に互い違いに変位するモー...
ccx2.3にバージョンアップして増えたオプションをちゃんと設...
解けないということかも知れない。ウォーニングメッセージとc...
ちゃんと確認する必要あり。
***ubuntu10.04上でccx1.7を走らせると [#lb48cba0]
というような訳で、ccx1.7のバイナリーをubuntu10.04上で走ら...
--(
./ccx_1.7: error while loading shared libraries: libg2c.s...
--)
これは、g77用のライブラリらしい。
-http://giropenfoam.web.fc2.com/salome/install_Suse/insta...
g77を入れてためしてみるか。
-http://d.hatena.ne.jp/ti10/20100512/1273671639
***nonpositive jacobian [#lfc2ec28]
--(
ERROR in e_c3d: nonpositive jacobian
determinant in element 14
--)
節点座標の1行文字数が非常に長い場合にこのエラーが出る場合...
そして、ccx_1.7では解けてccx_2.3では解けなかったりする。
--(
print'(i6,a,1pe15.3,a,1pe15.3,a,1pe15.3)'&
& ,i+1,',',x(i),',',y(i),',',z(i)
--)
みたいにフォーマットを指定して指数表示で1行を短めに出すと...
など1pdにするとcgxの表示がおかしくなるため、1peにしておく。
*** Factoring the system of equations using spooles ALLO...
C3D8の直方体片持ち梁で、ポアソン比を0.5に設定するとエラー...
ポアソン比0.49999と0.50001ではエラーは出なかった。
*シェル要素による片持ち梁のたわみ解析結果 [#t98ec863]
-&link(結果まとめ,http://www.str.ce.akita-u.ac.jp/cgi-bin...
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
