aiLiz Pink -2ページ目

aiLiz Pink -2ページ目

ホームピグアメブロ

芸能人ブログ人気ブログ

03Mar
- セルへの代入
  FORTRANから書き換えたVBAにおいて、結果をシートに記載する方法を示す。普通、セルに値を記載するのは、Rangeを使用する。例えば、セル"A1"へ100を記載するのは Range("A1") = 100"B1"へ200を記載するのは Range("B1") = 200セル"C1"へ =A1*B1を記載するのは Range("C1") = "=A1*B1"しかし、汎用的な計算で行列の大きさが計算する内容により変わる場合は、Cellsが使いやすい。 Cells(1,1) = 100 Cells(1,2) = 200ここで、困るのは計算式の記載である。 Cells(1,3) = ????そこで、アドレスの文字列を取得する。1行1列目のアドレス("A1")は、 Cells(1,1).Addressとなるが、絶対参照($A$1)となってしまうので、 Cells(1,1).Address(RowAbsolute:=False, ColumnAbsolute:=False)となる。同様に、1行2列目のアドレス("B1")は、 Cells(1,2).Address(RowAbsolute:=False, ColumnAbsolute:=False)となる。そこで、"=A1*B1" は、文字列の足し算(&)を組み合わせて、 Cells(1,3) = "=" & _ Cells(1,1).Address(RowAbsolute:=False, ColumnAbsolute:=False) _ & "*" & _ Cells(1,2).Address(RowAbsolute:=False, ColumnAbsolute:=False) _となる。また、EXCELで列を示す英字は、アドレスを"$"で分割して、 Split(Cells(1, 2).Address, "$")(1)で取得できる。(ちなみに、行は、 Split(Cells(1, 2).Address, "$")(2) または、Cells(15, 7).Rowで取得できる)一方、列を示す英字がわかっているが、何列目かがわからない場合は、 Range("AE124").Columnで取得できる。(ちなみに行は、 Range("AE124").Row で取得できる)
02Mar
- 関数副プログラム、サブルーチン
  きょうは、関数副プログラム、サブルーチンの書き換えの例です。VBAでは変数の型や、配列の次元数を合わせる必要があるため、古いFORTRANに慣れていると、使い勝手が悪く感じます。関数副プログラム PROGRAM PROA REAL A REAL B(3,5) REAL C(2,4) A=FUN(B,C) WRITE(6,*) A STOP END FUNCTION FUN(B,C) REAL B(3,5) REAL C(2,4) B(2,3)=7.3 FUN=B(2,3) RETURN END→(VBA)Sub ProA() Dim a As Double Dim b(1 To 3,1 To 5) As Double Dim c(1 To 2,1 To 4) As Double a = fun(b, c) Cells(1, 1) = aEnd SubFunction fun(b() As Double, c() As Double) b(2, 3) = 7.3 fun = b(2, 3)End Functionサブルーチン PROGRAM PROB REAL A REAL B(3,5) REAL C(2,4) CALL SABU(A,B,C) WRITE(6,*) A STOP END SUBROUTINE SABU(A,B, 1 C) REAL A REAL B(3,5) REAL C(2,4) B(2,3)=5.1 C(2,1)=B(2,3) A=C(2,1) RETURN END→(VBA)Sub SubB() Dim a As Double Dim b(1 To 3,1 To 5) As Double Dim c(1 To 2,1 To 4) As Double sabu a,b,c Cells(1, 1) = aEnd SubSub sabu(a As Double, b() As Double, _ c() As Double) b(2, 3) = 5.1 c(2, 1) = b(2, 3) a = b(2, 1)End Sub
01Mar
- IF文、DO文
  3回目は、IF文とDO文です。IF文 IF(I.GT.J) GO TO 100 (文) 100 A=B→(VBA)GoTo文は可読性が低下すると評判が悪いですが手っ取り早く書き換えるならば、ラベルを使うのが便利です If i > j Then GoTo 100 (文)100: a = b IF(I.GT.J.AND.J.LE.M) GO TO 100 (文) 100 A=B→(VBA) If i > j And j <= m Then GoTo 100 (文)100: a = b IF(I.GT.J) THEN A=B ELSE IF(I.EQ.J) THEN A=C ELSE A=D ENDIF→(VBA) If i > j then a = b ElseIf i = j Then a = c Else a = d End If IF(I+J) 100,200,300 100 A=B GO TO 400 200 A=C GO TO 400 300 A=D 400 (文)→(VBA) Select Case i+j Case Is < 0 a = b Case Is = 0 a = c Case Is > 0 a = d End SelectDO文 Do 10 I=1,10 (文) 10 CONTINUE Do I=1,10 (文) END DO→(VBA) For i = 1 To 10 (文) Next i Do 10 I=1,10 (文) IF(I.GT.J) GO TO 20 (文) 10 CONTINUE 20 (文) Do I=1,10 (文) IF(I>J) EXIT (文) END DO→(VBA) For i = 1 To 10 (文) If i > J Then Exit For (文) Next i DO WHILE (I<100) (文) END DO→(VBA) Do While i <= 100 (文) Loop Do (文) Loop While i < 100 Do Until i > 100 (文) Loop Do (文) Loop Until i < 100
28Feb
- 演算、GO TO 文等
  演算の置き換えは下記の通り A=B→(VBA) a = b A=B+C→(VBA) a = b + c A=B-C→(VBA) a = b - c A=B*C→(VBA) a = b * c A=B/C→(VBA) a = b / c B=A**2→(VBA) b = a^2文字列 NAME='NAMAE'//'MYOUJI'→(VBA) name = "namae" & "myouji" C=NAME(1:2)→(VBA) c = Left(name,2) C=NAME(3:4)→(VBA) c = Mid(name,3,2)GO TO 文の置き換えは下記の通り GO TO 100 (文) 100 CONTINUE→例えば文番号はラベルに置き換える(VBA) GoTo 100 (文)100: GO TO(100,200,300),J 100 A=B GO TO 400 200 A=C GO TO 400 300 A=D 400 ・・・→例えば(VBA) Select Case j Case 1 a = b Case 2 a = c Case 3 a = d End Selectコメント文の置き換えC*(空行)→例えば(VBA)'(空行)継続行の置き換え A = 1 B*C A = * B*C A = & B*C→(VBA) a = _ b * c
26Feb
- 宣言文
  オリンピックに見入って、更新が疎かになってしまった。テーマ未定としているが、なぜか仕事周辺の話題になってしまう。建築構造の仕事をしていると過去に作成したFortran プログラムという資産を使用したくなる。前世紀には高額だったFortranのコンパイラも、最近はフリーのものもあり、ほぼそのまま使えるが、ここでは、なにかと便利なExcel のVBAに書き換える方法を示す。詳しくは、インターネットを検索すると出てくるので、大学レベルで、FORTRANプログラムを書いていた人をイメージして、何回かに分けて、ヒント程度に記載する。また、書き換える場合の参考例を記載しているので、FORTRANとVBAの完全な一致、対応を示すものではない。なお、FORTRANの変数名は大文字、VBAは小文字で例示しているが、わかりやすくする目的で、VBAの変数名は、大文字でも、日本語(漢字、ひらがな、カナ等)でもいい。まずは、宣言文。型の宣言を強制する。 IMPLICIT NONE→(VBA)Option Explicit技術計算が前提なので、実数は、倍精度型Doubleを使用するのが望ましい。 DIMENSION A(100,50) REAL A(100,50) REAL :: A(100,50) REAL A(1:100,1:50) REAL*4 A(100,50) DOUBLE PRECISION A(100,50) REAL*8 A(100,50) REAL, DIMENSION(0:100) :: A→(VBA) Dim a(1 To 100, 1 To 50) As Doubleまたは、 Dim a() As Doubleと最初に宣言して、実際に使用する直前に、配列の大きさを定義する。なお、宣言文も、必ず最初に並べる必要はない。 ReDim a(1 To 100, 1 To 50)VBAで“1 To”を記載しないと、0から始まることになってしまう。整数も倍精度型Longを使用する。 DIMENSION M(100,50) INTEGER M(100,50) INTEGER :: M(100,50) INTEGER*4 M(100,50) INTEGER*8 M(100,50)→(VBA) Dim m(1 To 100, 1 To 50) As Longまたは、 Dim m() As Longと最初に宣言して、実際に使用する直前に、配列の大きさを定義する。 ReDim m(1 To 100, 1 To 50)文字列型はStringで宣言する。 CHARACTER NAME*80 CHARACTER*80 NAME CHARACTER(80) :: NAME→(VBA) Dim name As String
15Feb
- 特定天井
  地震で脱落した場合に重大な事故をもたらすことが予想される大規模の吊り天井をいう。具体的には、200平方メートル以上の面積があり、1平方メートル当たりの質量が2kgを超えるもの、高さが6m以上の天井で、詳細な定義は国土交通省告示等による。吊り天井とは、吊りボルトや、ハンガー等の金物を使用して野縁受け等が吊り下げられている天井をいう。それ以外の規模の小さな天井の安全性は、設計者に委ねられている。小規模でも、高い所にある天井、重い天井等は、法令の規制がないが、同等の配慮をすることが望ましい。また、大規模の天井等は、法令上の特定天井としないために、支持構造部に直付けすることが多く、専用の取り付け金物もある。その場合は、支持構造部の水平剛性に配慮しないと、天井が強制変形に追従できなく、逆に脱落する恐れがある。鉄骨造の屋根面に取り付ける天井等、中段のぶどう棚等に設置する天井等は、特に注意を要する。設計方法により、仕様ルート、計算ルートに分けられるが、いずれも複雑な形状でないことが、前提となっている。
08Feb
- 既存不適格建築物の増築
  地震の度に繰り返されるのが、耐震性の不足した建物の倒壊による犠牲である。被害を軽減するために、耐震改修を促す施策も増えているが、厳しくしすぎることも、所有者の財産権や、経済面で弊害となってしまう。建築物は、増改築等で、確認申請を受ける場合は、現行の建築基準に適合させる必要があるが、既存部分の架構は、すでに、出来上がっているものなので、適合させることが極めて難しく、実質的に、増改築できない。そのため、増改築を可能として、耐震改修も促す、緩和基準がある。既存不適格となった時点の面積が基準となり、増築面積の大小や、構造的に縁を切るか、既存部の耐震性能等が問題となる。構造的に縁を切り、耐震補強を行っていれば、増築部分の面積に制限がない。フローチャートを示す。
01Feb
- 地震波の解析(3)
  地震波の解析その3 離散フーリエ変換と、最も簡単すぎる例です。
31Jan
- 地震波の解析(2)
  地震波の解析その2 複素フーリエ変換です。
30Jan
- 地震波の解析(1)
  地震波の解析その1　フーリエ変換です。
27Jan
- 数値積分(6)
  数値積分法その6 床のねじれを考慮する場合です。
25Jan
- 数値積分(5)
  数値積分法その5 多層の増分解析です。
24Jan
- 数値積分(4)
  数値積分法その4 弾塑性解析で用いられる増分による式になります。
23Jan
- 数値積分(3)
  線形加速度法に続いて、数値積分法その3 平均加速度法、衝撃加速度法、段階加速度法、ニューマークのβ法です。
22Jan
- 数値積分(2)
  数値積分線形加速度法その2 変位です。
21Jan
- 数値積分(1)
  数値積分法の概要と、線形加速度法の加速度までです。
20Jan
- 多質点系の振動(5)
  電算機を使用した多質点系の固有値、固有周期の計算法です。
19Jan
- 多質点系の振動(4)
  多質点系の振動その4 マトリックスを利用した計算と、刺激関数等です。
18Jan
- 多質点系の振動(3)
  多質点系の振動は、わかりにくいので、2質点系の計算例です。
17Jan
- 多質点系の振動(2)
  多質点系の振動　その2 自由振動から、固有周期、固有モードです。