MYSTRANsolver / MYSTRAN

MYSTRAN is a general purpose finite element analysis solver
https://www.mystran.com/
MIT License
55 stars 17 forks source link

Incorrect OP2 subcase id #31

Closed SteveDoyle2 closed 9 months ago

SteveDoyle2 commented 9 months ago

EB-ALL-ELEM-TEST-LANCZOS-MODE-2-DPB.txt

This example was from the benchmark tests. This is a regression that was not seen in 15.1, but is seen in 15.1.1 and subsequent versions.

15.1.1

displacements[0]
  isubcase = 0
  type=RealDisplacementArray ntimes=1 nnodes=13, table_name=OUGV1
  data: [t1, t2, t3, r1, r2, r3] shape=[1, 13, 6] dtype=float32
  node_gridtype.shape = (13, 2)
  sort1
  modes = [2]; dtype=int32
  eigns = [359293.219]; dtype=float64
  mode_cycles = [599.411]; dtype=float64

displacements[1]
  isubcase = 1
  type=RealDisplacementArray ntimes=1 nnodes=13, table_name=OUGV1
  data: [t1, t2, t3, r1, r2, r3] shape=[1, 13, 6] dtype=float32
  node_gridtype.shape = (13, 2)
  sort1
  modes = [1]; dtype=int32
  eigns = [175725.922]; dtype=float64
  mode_cycles = [419.197]; dtype=float64

15.1

displacements[1]
  isubcase = 1
  type=RealDisplacementArray ntimes=1 nnodes=13, table_name=OUGV1
  data: [t1, t2, t3, r1, r2, r3] shape=[1, 13, 6] dtype=float32
  node_gridtype.shape = (13, 2)
  sort1
  modes = [1]; dtype=int32
  eigns = [175725.922]; dtype=float64
  mode_cycles = [419.197]; dtype=float64
SteveDoyle2 commented 9 months ago

It also looks like there is an extra mode. It's unclear why this affects the subcase id.

Bruno02468 commented 9 months ago

The only change across these versions was properly reading one of the params (N, number of requested modes) in the EIGRL card -- finalising a fix for a memory bug.

There is now an extra mode because your EIGRL card asked for two (field #5). MYSTRAN (and other Nastran-like solvers) treat each mode as a subcase, so two modes means two subcases.

Before 15.1.1, you were only getting output for one mode even though you requested two modes be computed. So this is not a regression, it's a bug fix! ;)

SteveDoyle2 commented 9 months ago

Also of note, the eigenvectors are incorrectly flagged as displacements