2015年09月の記事
全22件 (22件中 1-22件目)
1
-
オープンCAEソフトまとめサイト
以下にオープンCAEソフトがまとめられている。http://imechanica.org/node/1394〇プリ Gmsh, Triangle 〇非線形構造解析ソルバー CalculiX, Code_Aster, Impact, OOFEM, Tahoe, Warp3D, FEBio 〇線形構造解析ソルバー Elmer, FEMPACK finite element routines, FETK, OpenFEM, OOFにほんブログ村
2015年09月30日
コメント(0)
-
sfepy 拘束節点荷重、応力、ひずみの出力
ヘキサ要素5個の片持ち梁をsfepyで解析し、拘束節点荷重、応力、ひずみを出力した。左側の節点を拘束し、右側の4節点に1Nずつ-z方向に荷重をかける。↓X方向荷重↓Z方向荷重↓コーシーひずみZX↓コーシー応力ZX↓ソースコード-----ここから# -*- coding: utf-8 -*-"""5hexas cantilever"""from sfepy.mechanics.matcoefs import lame_from_youngpoissonfrom sfepy.mechanics.matcoefs import stiffness_from_youngpoissonfrom sfepy.discrete.fem.utils import refine_meshfrom sfepy import data_dir# Fix the mesh file name if you run this file outside the SfePy directory.filename_mesh = '5hexas.med'output_dir = '.' # set this to a valid directory you have write access toyoung = 210000.0 # Young's modulus[MPa]poisson = 0.3 # Poisson's ratiooptions = { 'output_dir' : output_dir,}regions = { 'Omega' : 'all', 'Left' : ('vertices in (x < 0.001)', 'facet'), 'Right': ('vertices in (x > 4.999)', 'facet'),}materials = { 'steel' : ({ 'lam' : lame_from_youngpoisson(young, poisson)[0], 'mu' : lame_from_youngpoisson(young, poisson)[1], },), 'Load' : ({'.val' : [0.0, 0.0, -1.0]},),}fields = { 'displacement': ('real', 'vector', 'Omega', 1),}equations = { 'balance_of_forces' : """dw_lin_elastic_iso.2.Omega(steel.lam, steel.mu, v, u) = dw_point_load.0.Right(Load.val, v)""",}variables = { 'u' : ('unknown field', 'displacement', 0), 'v' : ('test field', 'displacement', 'u'),}ebcs = { 'Left' : ('Left', {'u.all' : 0.0}),}solvers = { 'ls' : ('ls.scipy_direct', {}), 'newton' : ('nls.newton', { 'i_max' : 1, 'eps_a' : 1e-10, }),}def post_process(out, pb, state, extend=False): """ Calculate and output strain and stress for given displacements. """ from sfepy.base.base import Struct ev = pb.evaluate strain = ev('ev_cauchy_strain.3.Omega(u)', mode='el_avg') stress = ev('ev_cauchy_stress.3.Omega(steel.D, u)', mode='el_avg') out['cauchy_strain'] = Struct(name='output_data', mode='cell', data=strain, dofs=None) out['cauchy_stress'] = Struct(name='output_data', mode='cell', data=stress, dofs=None) ebcs = pb.ebcs pb.time_update(ebcs={}) vec = pb.get_variables().create_state_vector() pb.equations.evaluate(mode='weak', dw_mode='vector', asm_obj=vec) out['res'] = Struct(name='res', mode='vertex', data=vec.reshape((-1, 3)), dofs=None) pb.time_update(ebcs=ebcs) return outsteel = materials['steel'][0]steel.update({'D' : stiffness_from_youngpoisson(3, young, poisson)})options.update({'post_process_hook' : 'post_process',})-----ここまでエクセルVBAの結果とも合っていた。節点荷重は、ebcsにある、拘束節点のみ出力される。荷重節点は出力されない。Developer Guideをもっと見なくては。にほんブログ村
2015年09月27日
コメント(0)
-
python 画像のダウンロードと複数画像を並べる
株価3ヶ月チャート画像をダウンロードして5x10=50個ずつ並べるスクリプトを書いてみた。-----ここから # -*- coding: utf-8 -*-import urllib2from PIL import Image#yahooファイナンスから３ヶ月チャート画像をダウンロードして50個ずつ並べるdef download(url, path): fp = urllib2.urlopen(url) local = open(path + code+'.jpg', 'wb') local.write(fp.read()) local.close() fp.close()def paste(code, h, v, k): img = Image.open(str(code)+'.jpg', 'r') canvas.paste(img, (((k-1)%5)*h, ((k-1)//5)*v))f=open('codes_150926.csv', 'r') #銘柄コードi=0j=1k=0h=512 #水平ピクセルv=310 #垂直ピクセル# マージに利用する下地画像を作成するcanvas = Image.new('RGB', (5*h, 10*v), (255, 255, 255))for code in f: code=code.rstrip() i+=1 k=i-50*(j-1) print code, k url='http://chart.yahoo.co.jp/?code='+str(code)+'.T&tm=3m&type=c&log=off&size=m&over=m65,m130,s&add=v&comp=' download(url, './') paste(code, h, v, k) if i==50*j: # 保存 canvas.save('c'+str(j)+'.jpg', 'JPEG', quality=100, optimize=True) canvas.close() canvas = Image.new('RGB', (5*h, 10*v), (255, 255, 255)) j+=1canvas.save('c'+str(j)+'.jpg', 'JPEG', quality=100, optimize=True)canvas.close()f.close()-----ここまで ↓結果はこんな感じ codes_150926.csvは以下のような銘柄コードのリスト。 1301 1332 1333 … 9995 9996 9997 画像ダウンロードは、以下の関数で実施。def download(url, path): fp = urllib2.urlopen(url) local = open(path + code+'.jpg', 'wb') local.write(fp.read()) local.close() fp.close() pillowというモジュールで画像の編集がいろいろとできる。画像を並べるのは以下の関数で実施。def paste(code, h, v, k): img = Image.open(str(code)+'.jpg', 'r') canvas.paste(img, (((k-1)%5)*h, ((k-1)//5)*v)) ↓参考URLhttp://blog.takkinoue.com/entry/2015/04/04/171000https://librabuch.jp/2013/05/python_pillow_pil/ にほんブログ村
2015年09月26日
コメント(0)
-
株安
・2015/9/24の日経平均株価終値 17,571.83前日比 -498.38VW不正発覚のあおりを受けて予想通りといったところでしょうか。。にほんブログ村
2015年09月24日
コメント(1)
-
火花
↓ぼくにとってのポイント〇徳永と神谷が一緒に由貴の家に行ったシーン・伝える努力を怠ったら、自分の面白いことがなかったことにされるから→この言葉が、プロとアマを分かつように思えた。・自分の理想を崩さず、世間の観念とも闘う。→笑いに対して、神谷はどこまでも自分の理想を追い求め、ついには誰にも理解されなくなっていく。・美しい世界を、いかに台無しにするかが肝心→美しいまま、完成されたままだと終わってしまう、死んでしまう、ということなんだろうか。-----「おかしい」とは正常ではないことであり、おかしいから笑う。「おかしい」を演出し、人を笑わせることに、芸人は命を懸けている。「おかしい」ことに、笑うか、怒るか、しらけるか、は人によってちがう。万人に共通の笑いのツボなどない。人の感情はあてにならない。笑いは感情の一部であり、安定していない。人を笑わせることを生業として生きていくことは、とてつもなく難しいように思う。そのためか、ときに、極端な演出をし、誰かを傷つけることも起きる。自分の理想とする笑いを追い求めながら、世間のタブーなどの制約条件をクリアし、多くの人を笑わせる。これは、以下のように言い換えれば、アーティスト全般にあてはまる。自分の理想とする作品を追い求めながら、音、文字、言葉、素材、時間、コストなどの制約条件をクリアし、多くの人を感動させる。さらに、以下のように言い換えれば、すべての職業人にあてはまる。自分の理想とする○○を追い求めながら、制約条件をクリアし、多くの人に使ってもらう。利益が上がらなければ、事業を存続できない。お笑い芸人についての小説だけれど、すべての職業人が共感する小説だと思う。【楽天ブックスならいつでも送料無料】火花 [ 又吉直樹 ]にほんブログ村
2015年09月23日
コメント(0)
-
サッカーの練習
職場の同僚の出身大学でサッカーの練習。いつになく22人も集まり、フルコートで11対11ができた。どうもうまくボールをコントロールできない。練習しなくては。にほんブログ村
2015年09月22日
コメント(0)
-
アプリを作ろう！ Android入門
第１章 Ａｎｄｒｏｉｄというアプリについて知ろう第２章 アプリを作る準備をしよう第３章 Ａｎｄｒｏｉｄ Ｓｔｕｄｉｏでアプリ作成を始めよう第４章 アプリでＪａｖａの基本を学ぼう第５章 アプリに画像を組み込もう第６章 アプリを完成させよう第７章 アプリから画面を呼び出そう第８章 アプリをＡｎｄｒｏｉｄ端末で動かそう第９章 アプリを公開しようタブレット、スマートフォンのアプリ開発というものを、どうやるのか知りたくて読んでみた。本の通りにやれば、開発環境の準備から、おみくじアプリを作成し、アンドロイド端末で動作させることができた。開発言語はJavaで、書く量は結構多いんだけど、Android Studioの補完機能でずいぶん楽に書ける。説明が丁寧で、図を見ながら進めれば、どこをクリックして、どこに書けばいいのかがわかる。作業の意味についても解説があり、何のためにこの作業をやっているのかがわかる。Javaの基礎についても学ぶことができる。ぼくも仕事でトレーニング資料を作ることがあるので、参考になる。Androidは各バージョンにコードネームがついていて、お菓子の名前になっている。現バージョンの5.0, 5.1はLollipopだ。親しみやすいようにということだろうか。↓Android Studio↓おみくじアプリ【楽天ブックスならいつでも送料無料】【高額商品】【3倍】プリを作ろう！ Android入門 Andro...
2015年09月21日
コメント(0)
-
近所の公園で自主練習
昨日の反省から近所の公園で自主練習。３０分くらいやって、そろそろ帰ろうかと思い始めたころ、昨日も一緒にフットサルをした近所の知り合いと会い、一緒に練習することにした。そうこうしていると、近所の子どもがやってきて、２対２のミニゲームになった。しばらくすると、近所の高校生もやってきた。いつもは一人で練習することが多いのだけれど、いつになく、にぎやかに練習できて楽しかった。にほんブログ村
2015年09月20日
コメント(0)
-
フットサルの練習
八王子のFUNフットサルで練習。近所の知り合いの職場の仲間つながりで15人集まり、3チームでゲームを回した。この日は、１ヶ月ぶりくらいにボールを蹴ったこともあり、ダメダメだった。。全然ボールが足につかなかった。。11月には職場のチームで大会に出ることにしたので、練習しなくては。にほんブログ村
2015年09月19日
コメント(0)
-
職場の飲み会
池袋ルミネ9F菜の庵にて。 http://s.tabelog.com/tokyo/A1305/A130501/13142996/dtlcoupon/ しゃぶしゃぶがおいしかった。 普段話さない人たちと仕事以外の話ができたのが収穫だった。 職場では、どうしても仕事の話ばかりになるからな。
2015年09月18日
コメント(0)
-
会社四季報 2015年4集の分析 (2)
続きです。2015年実績から2016年予想への業績予想の伸び率を上位10社ずつランキングしました。伸び率=(2016予想/2015実績-1)*100 % 〇設備投資〇研究開発設備投資と研究開発費が伸びているということは、将来の事業拡大を見込んでいるのかもしれません。〇営業キャッシュフロー〇現金等営業キャッシュフローと現金等は昨年と今年の比較になります。（まちがえやすい。。）これらが伸びていると会社が儲かっていることが推測されます。中国経済の減速、米FRBの利上げの可能性など、難しい状況ですが、参考になれば幸いです。現在の価格と価値との差を見極められればいいんですけど。にほんブログ村
2015年09月15日
コメント(0)
-
会社四季報 2015年4集の分析
2015年実績から2016年予想への業績予想の伸び率を上位10社ずつランキングしました。伸び率=(2016予想/2015実績-1)*100 % 〇売上高〇営業利益〇経常利益〇純利益〇予想配当利回りこういうところを買って、持ち続けるのもいいかもしれません。参考になれば、幸いです。にほんブログ村
2015年09月14日
コメント(1)
-
TOEIC受験
なんかreadingの手応えがよくなかった。。 受験会場の東京工科大学のキャンパスがめちゃめちゃ立派でびっくりした。 にほんブログ村
2015年09月13日
コメント(0)
-
第05回オープンCAE勉強会＠関東（構造など）
（１）10:00～12:00 フォローアップ勉強（２）13:00～15:00 ハンズオン勉強会FrontISTR入門（３）15:00～17:20 研究会1) 間違えないFEMのために2) SfePyの概要と試用3) GetFEM++の概要と試用4) MPSとCode_Asterの連成解析について5) 事務連絡-----(2)easy istir結構使いやすかった。(3)1FEMでは、メッシュ作成→境界条件設定→計算→結果確認となど多くの工程が必要。各工程の信頼度が0.8程度では結果の信頼性は0.4程度になってしまう。間違えないためには各工程の信頼度を限りなく1に近づけなければならない。(3)2ぼくの発表。sfepyは偏微分方程式を有限要素法で解くpythonモジュール。(3)3getfem++のWikipedia日本語記事の作者の発表。マトリクスを直接見られる。ipython notebookでデモをしてくれた。(3)4ダム決壊の流体-構造連成解析をcode-asterとmpsで見せてくれた。収穫の多い勉強会だった。にほんブログ村
2015年09月12日
コメント(1)
-
業務のためのpython勉強会#4
ハンズオンでipython notebookの使い方を学んだ。ipython notebookはその名のとおり、pythonコマンドの勉強記録として使えることがわかった。補完機能もあり、pythonコマンドをためすのに適している。matplotlibのグラフも即座に確認できる。懇親会で様々な業種の人と知り合うことができた。インフラサーバー管理、大学の先生、クラウド、マーケティング、ゲノム解析、画像処理など。マーケティングに機会学習を使いたい人が多いようだった。http://startpython.connpass.com/event/18053/?utm_campaign=event_reminder&utm_source=notifications&utm_medium=email&utm_content=detail_btnにほんブログ村
2015年09月10日
コメント(0)
-
日経平均 急騰 +1,343.43円
2015/9/9の日経平均株価 終値 18,770.51 前日比 +1,343.43上げ幅は、1994年１月31日（1471.24円）以来、21年７カ月ぶりの大きさらしい。さすがに、最近の変動には振り落とされそうになる。にほんブログ村
2015年09月09日
コメント(0)
-
Androidフォン英語学習法
プロローグ Ａｎｄｒｏｉｄフォンは世界最強の英語学習ツールー片手でちょっといじるだけで４技能がアップ１ リーディング上達法ーＡｎｄｒｏｉｄフォンで空き時間に英文を読む２ ライティング上達法ーＡｎｄｒｏｉｄフォンでこまめに英文を書く３ リスニング上達法ーＡｎｄｒｏｉｄフォンでネイティブ英語を聞く４ スピーキング上達法ーＡｎｄｒｏｉｄフォンと一緒に英語を話す５ ボキャブラリー増強法ーＡｎｄｒｏｉｄフォンで空き時間に語彙を増やすAndroidフォンにしたこともあり、読んでみた。Bible Quotes、English Proverbs、Fun jokes(20000 Awesome Jokesはなかった)をダウンロードしてみた。Fun jokes読んでみたけど、結構長い。そして、笑えない。。Google英語勉強法にもあったが、Googleで自分の書いた英語を検索して、ヒット数で妥当かどうかを確認する方法が紹介されている。広告などの文章をアレンジして、自分の文章にする方法は参考になる。【楽天ブックスならいつでも送料無料】Androidフォン英語学習法 [ 藤田英時 ]にほんブログ村
2015年09月08日
コメント(0)
-
この１ヶ月の株価伸び率ヒストグラム
8/3～9/4の株価伸び率ヒストグラムを描いてみた。8/3の調整後終値を1としている。この１ヶ月で平均12%下がっている。ひどいもんだ。標準偏差10.6%なので、±3σで-40%～+20%にほぼすべて(99.7%)収まっている。にほんブログ村
2015年09月07日
コメント(0)
-
国家の役割
中東からヨーロッパへの難民、移民がニュースになっている。少し前は、ミャンマーのロギンヒャと呼ばれる人たちが難民として海を漂流し、ニュースになっていた。いずれも内戦、政情不安、迫害などで、自国で生活することができず、難民となった人々のようだ。-----以下、ぼくの見解。国家の役割は、第一に国民の生命と財産を守ること。そのことを明記するために憲法があり、憲法のもとに法律があり、行政、立法、司法機能がある。国家は、法の下で国民の生命と財産を守るために、軍隊や警察といった暴力装置を独占する。国民は、武装を放棄し、労働で得た対価から納税し、国家を支える。-----以下、現状の確認世界には、約200の国と地域があり、民主国家、独裁国家などさまざまな統治の仕方がある。英国で産業革命が起きて以来、世界は欧米主導で動いており、民主主義は、完璧ではないけれども、現状では最善の統治方法ということになっている。米国は、2001年9月11日の同時多発テロ以来、アフガニスタンとイラクに戦争をしかけている。イラクとの戦争の理由は、1.大量破壊兵器保有疑惑、2.独裁者フセイン打倒、だった。1.は大量破壊兵器が見つからず、2.はフセインを確保→処刑となった。2.について、米国の大義名分としては、独裁者フセインを打倒し、イラク国民を解放してイラクを民主国家にする、というものだった。ところが、フセイン政権打倒後、新政権はうまく機能せず、イラクは政情不安が続いている。アフガニスタンも政情不安が続いている。2010年から2012年にかけて、アラブの春と呼ばれる中東の民主化運動がおこった。チュニジア、エジプト、リビア、イエメンで政権が打倒された。ところが、その後の国づくりがうまくいかず、混乱が続いている。シリアも、民主化運動を推進する勢力とアサド政権が衝突し、ISISの台頭もあり、泥沼の内戦になっている。近年、ロシア、中国といった民主主義を採用しない国が台頭してきている。-----現状を見る限り、民主主義推進はうまくいっていないようだ。特に、アラブの春以降、かえって中東地域は混乱し、難民が増えているようだ。難民の多くは、ハンガリー→オーストリア経由で経済的に豊かで難民に寛容なドイツを目指している。しかし、ドイツでも、難民受け入れについて賛否意見が分かれている。ヨーロッパは税金が高く高福祉の国が多い。難民受け入れ反対の意見として、以下がある。・治安悪化・自分たちの税金が難民たちのために使われることに納得いかない。こうなってくると、例えどんな独裁政権であっても、民主化などせずに、現政権を存続させて、国民を難民とさせずに、その国にとどめるべきだ、という意見がでてきてもおかしくない気がしてくる。例えば、北朝鮮や中国で民主化が起こって政権が転覆し、その後の政府が機能せず、多くの難民が日本に押し寄せる事態になったらどうなるのだろう？一方で、難民となった人々ばかりではなく、独裁政権下で基本的人権を守られずに生きている人々も数多い。こういった人々に、ただ生まれた時代と場所が悪かったと言うことができるだろうか？やむにやまれず、命からがら国から逃れた人々を、追い返すしかないのだろうか？ただ、現実は、いつも弱い立場の人々に過酷だ。。国家の役割は、第一に国民の生命と財産を守ること。ぼくは、国家は必要最小限であるべき、というリバータリアンの考え方に傾倒してきた。しかし、国家が機能しなくなるとどうなるか、それは現在の中東情勢が示している。国家の役割は、ぼくが思っていた以上に大きいのかもしれない。にほんブログ村
2015年09月06日
コメント(0)
-
この１ヶ月の業種ごと比較チャート
8/3～9/4の業種ごとの比較チャートを描いてみた。8/3の調整後終値を1としている。業種のならびは以下のようになっている。やはり、軒並み下がっている。見た感じ、建設、食品、化学、非鉄金属は下がり方が小さいようだ。日本企業は好決算のところが多く、ファンダメンタルはいいと思うんだけどな。にほんブログ村
2015年09月05日
コメント(0)
-
作業の継ぎ目を減らす
CAE業務では、以下のフローの中で複数のソフトを使うため、作業の継ぎ目が多くなる。プリ（メッシュ作成ソフト、条件設定ソフト）→ソルバー（計算ソフト）→ポスト（可視化ソフト、画像編集ソフト、エクセル、パワーポイント、ワードなどレポート作成ソフト）効率化のためには、シームレスに進められるようにしたい。そのために以下の方策が考えられる。1.使うソフトを減らしても目的のレポート作成をできるようにする。2.スクリプト、マクロを活用して手作業を減らす。3.異なるソフト間でのデータ受け渡しをスクリプト、マクロなどで自動化する。最近は、画像などをパワーポイントに出力し、定型のレポート作成機能をもつ可視化ソフトも出てきている。結局、最終的なアウトプットの項目と書式を決めることが一番重要になる。それさえ決まれば、そこに向かうための道を整備することができる。にほんブログ村
2015年09月04日
コメント(0)
-
テキスト処理の流れ
基本的な流れ1.テキストファイルを開く。2.1行ずつ読み込む。3.必要なデータを変数に入れる。4.出力フォーマットに合わせて書き出す。〇エクセルVBAの例点の座標とIDをdata.txtから読み込んで、距離計算してout.txtに書き出す。Open "C:\data.txt" For Input As #1 '入力ファイルOpen "C:\out.txt" For Output As #2 '出力ファイルDo Until EOF(1) Line Input #1, buf If InStr(buf, "aaaaa") > 0 Then Do Line Input #1, buf inod = inod + 1 nid(inod) = Left(buf, 8) xnod(inod) = Mid(buf, 9, 16) ynod(inod) = Mid(buf, 25, 16) znod(inod) = Mid(buf, 41, 16) Loop End IfLoopFor i = 1 To inod-1 for j=i+1 to inod dist=sqr((xnod(j)-xnod(i))^2+(ynod(j)-ynod(i))^2+(znod(j)-znod(i))^2) print #2, nid(i) & "," & nid(j) & "," & dist nextnextClose #1Close #2にほんブログ村
2015年09月01日
コメント(0)
全22件 (22件中 1-22件目)
1
