解析条件 解析パラメータを下表に示す。 変化させているパラメータは、制水口径および遮断時間でる。 項目 採用値 解析条件 負荷遮断 シャフト内径（断面積） D=10.0 m (F=78.5 m2) ポート内径（断面積） Dp=2.0, 3.0, 4.0, 8.0m (Fp=3.14, 7.065, 12.56, 50…
きっかけ 昨日、私のwebページをご覧になった方から、サージング解析プログラムに関するメールを頂いた。 下が、その疑惑のページのコピーである。 図を見てわかるように。 サージング解析では、「制水口（ポート）径が小さい場合、遮断時間を大きくしていく…
開水路トンネルの不等流解析（常流）を行ったので、そのプログラムをアップしておく。 全区間、常流であるため、下流端で固定水位を与え、上流に向かって逐次水位を求めていく。 成果図 このプログラムでのTips 与えられた固定点数点に対し計算点を等間隔に…
2020年1月18日（土）、iMacとMacbook proにmacOS Catalinaをクリーンインストールし、作業環境を更新した。 更新内容を記録しておく。 マシン iMac (Retina 4K, 21.5-inch, 2017) MacBook Pro (Retina, 13-inch, Mid 2014) 不便な点 この更新を行ってこれま…
画像変換 from PIL import Image import os files = os.listdir() for file in files: base,ext=os.path.splitext(file) if ext=='.png': input_im = Image.open(base + ".png") rgb_im = input_im.convert('RGB') rgb_im.save(base + ".jpg",quality=30) pr…
RC部材の設計をしたので流れを記録しておく。 リンク 平面骨組解析プログラムはこちら（Qiita）。 配筋計算プログラムの自作関数はこちら（はてな）。 出力 モデル図 断面力図 配鉄設計 プログラム モデル図 import matplotlib.pyplot as plt import numpy a…
コンター図を作成する必要があったため、そのプログラムを作成した。 作例 プログラムソース import numpy as np import matplotlib.pyplot as plt import seaborn as sns def drawfig(xx,yy,zh,zc): q=68.64 vp1=7.0; dp1=np.sqrt(4*q/np.pi/vp1) vp2=2.5; …
概要 岩盤内埋設式水圧鉄管の設計プログラムをフルで書いてみたので残しておく。 プログラム前半の非常に長い関数「def xlswrite」は、エクセルに書式指定して結果を書き出すもの。 最終的にワードの報告書に貼り付けなくてはならないので、エクセルにしてお…
タイトル通り、jpg画像を縮小保存しhtmlで表示する。 import glob import os from PIL import Image import os def select_pic(): path = '*.jpg' # 現フォルダ内にあって拡張子がJPGのファイル名を取得 file_list = glob.glob(path) print(file_list) retur…
記事の最後に行く 昨日（2019.08.11）、また機械学習の本を買ってしまった。 オライリー・ジャパン PythonによりAIプログラミング入門 ディープラーニングを始める前に身に着けておくべき15の基礎技術 あと、久しぶりにpipで一括アップデートしたら、matplot…
記事の最後に行く いよいよ７月。 昨年（2018年10月10日）、マレーシアでの業務が完了し、帰国した。 それ以降、以下のような感じで過ごしてきた。 2018年10月 ひま 2018年11月 ひま 2018年12月 忙しい 2019年01月 忙しい 2019年02月 忙しい（ベトナム出張２…
記事の最後に行く 自作関数 以下の確率分布に対するパラメータ推定と予測値計算の部分を自作関数としてプログラム本体から分離している。 含まれる関数（確率分布） 関数名説明 ln3３変数対数正規分布（LN3） lp3対数ピアソンIII型分布（LP3） gev一般化極値…
記事の最後に行く はじめに 格子桁構造解析プログラムを紹介します。 このプログラムは、平面骨組構造解析プログラムの軸力・軸方向変位をねじりモーメント・ねじり角に書き換えたものです。 構造は平面形状とし、x - y 平面上で定義します。 格子桁解析は、…
記事の最後に行く はじめに マレーシア時代に使っていた、RC長方形断面の曲げ耐荷力計算プログラムを掲載しています。 部分安全係数、コンクリートおよび鉄筋の物性は BS 8110-1 に準拠したものとしています。 日本のコンクリート標準示方書での計算も、部分…
記事の最後に行く 平方根指数型最大値分布による確率雨量計算 最近ある方から、昔作った文書http://civilyarou.web.fc2.com/WANtaroHP_html5_win/f90_ENGI/dir_HFA/suimon.pdfの平方根指数型最大値分布に関連して、メールを頂いた。 そのころは、C あるいは …
記事の最後に行く ここでは、「設計 仮排水路設計における Flood Routine の活用（２）」で述べた事項に関連する計算・作図プログラムのうち、以下の4つを掲載する。 トンネル標準断面作図プログラム トンネル通水量計算・作図プログラム 貯水池水位容量曲線…
記事の最後に行く ここでは、「設計 仮排水路設計における Flood Routine の活用（２）」で述べた事項に関連する計算・作図プログラムのうち、以下の２つを掲載する。 Flood Routine 解析プログラム Flood Routine 作図プログラム Flood Routine 解析プログ…
記事の最後に行く 概要 仮排水路トンネル断面寸法と仮締切堤高さを決定するための設計検討を行う。 具体的には、トンネル径を変化させ、Flood Routine により、洪水流入波形、流出流量、貯水池水位の関係を求める。 Conditions for Flood Routine Analysis F…
記事の最後に行く はじめに 以下に示す、「設計 RC円形圧力トンネルの配筋設計（１）」で示した図を作成するプログラムである。 ポイント アノーテーション用矢印とテキスト描画 説明用のボックス内テキストと矢印を描画する。 matplotlib では、annotate に…
記事の最後に行く 内水圧を受けるRC圧力トンネル 出力 ra覆工内半径 (mm) rb覆工外半径 (mm) r0岩盤外縁半径位 (mm) ta内側鉄筋等価板厚 (mm) tb外側鉄筋等価板厚 (mm) da内側鉄筋かぶり (mm) db外側鉄筋かぶり (mm) Pa内水圧 (MPa) T温度変化量（マイナスは…
記事の最後に行く 概要 均等な内水圧および外水圧を受けるRC円形圧力トンネルを、トンネル軸方向に均一な厚肉円筒として、平面ひずみ状態でモデル化する。 モデル化においては、以下の考え方を採用した。 圧力水路の構造は鉄筋コンクリート構造とする。 内水…
記事の最後に行く Safety Factor for Shear Strength and Tensile Strength The definitions of symbols are shown below in this discussion considering compression members such as concrete or rock material. cohesion of material internal friction …
記事の最後に行く Relationship between Shear Strength and Uniaxial Strength The definitions of symbols are shown below in this discussion considering compression members such as concrete or rock material. cohesion of material internal frict…
記事の最後に行く この記事は、取水口における渦防止のための被り水深を提案している以下の論文の提案式（ft系）を、ｍ系に直したものを記載している。 Vortices at intakes by J. L. Gordon (Water Power April 1970) Gordonは、論文の中で、取水口における…
記事の最後に行く フォルダ内の全 jpg ファイルの余白を削除し、同じ名前で保存します。 トリム関数は以下に掲載のものを使わせていただきました。 トリム：https://gist.github.com/mattjmorrison/932345 from PIL import Image, ImageChops import glob, o…
記事の最後に行く フォルダ内の全JPGファイルを、幅600pixelに縮小し、600_***.jpg という名前で保存します。 また、0_small.html というファイルを作成し、縮小画像を一覧形式で閲覧できるようにします。 import glob import os.path from PIL import Image…
記事の最後に行く はじめに この記事は、Qiita に投稿した以下の記事の再アレンジ版です。 https://qiita.com/damyarou/items/9cb633e844c78307134a 経度・緯度で与えられた２点間距離を計算するプログラムです。 Lambert-Andoyerの方法を用いています。 デ…
記事の最後に行く はじめに この記事は、Qiita に投稿した以下の記事の再アレンジ版です。 https://qiita.com/damyarou/items/b50b6230b103cbab02ba ETOPO1 は，National Geophysical Data Center (Colorado) により作成された，全球の１分刻みの地形モデル…
記事の最後に行く はじめに ここで紹介しているのは、制水口式調圧水槽のサージングの基礎微分方程式（連立微分方程式）を Runge-Kutta 法により解き、水面振動の時刻歴を求めるプログラムです。以下のような出力が得られます。 計算理論 Fundamental Differ…
記事の最後に行く はじめに 制水口式調圧水槽の設計において用いられるVogt-Forchheimer式は、指定した制水口径および立坑径に対し貯水池水位からの最大上昇水位を求めるものであり、制水口径および立坑径をパラメータとして最大上昇水位を計算し、条件を満…
