Welcome, Guest
You have to register before you can post on our site.

Username
  

Password
  





Search Forums

(Advanced Search)

Latest Threads
Ansible Role to install M...
Forum: MYSTRAN
Last Post: numenic
04-19-2021, 09:41 PM
» Replies: 0
» Views: 125
Python Neutral Parse
Forum: Pre/Post Processors
Last Post: numenic
04-16-2021, 08:02 AM
» Replies: 0
» Views: 65
FEMAP neutral improvement...
Forum: Pre/Post Processors
Last Post: numenic
04-16-2021, 07:43 AM
» Replies: 7
» Views: 314
NEWS: ComLab Released to ...
Forum: MYSTRAN
Last Post: Admin
03-16-2021, 03:06 PM
» Replies: 1
» Views: 3,376
ComLab Enhancments
Forum: ComLab and NASTRAN-95
Last Post: Admin
03-16-2021, 02:49 PM
» Replies: 0
» Views: 285
Added a Parameter that Re...
Forum: MYSTRAN
Last Post: ZacharyALerner
03-15-2021, 03:09 AM
» Replies: 0
» Views: 327
Compiler Error for Prefer...
Forum: MYSTRAN
Last Post: Admin
02-23-2021, 06:25 PM
» Replies: 1
» Views: 1,239
Femap with Mystran
Forum: Pre/Post Processors
Last Post: SamiH.
02-18-2021, 12:05 PM
» Replies: 2
» Views: 346
ComLab (an improved NASTR...
Forum: ComLab and NASTRAN-95
Last Post: Admin
02-10-2021, 07:47 AM
» Replies: 0
» Views: 802
3/16/2021 UPDATE
Forum: MYSTRAN
Last Post: Admin
02-04-2021, 04:54 AM
» Replies: 1
» Views: 4,571

 
  Ansible Role to install MYSTRAN
Posted by: numenic - 04-19-2021, 09:41 PM - Forum: MYSTRAN - No Replies

I pusblished here (https://framagit.org/numenic/mystran_instl_ansible_role) a small Ansible role for installing Mystran. I use it under Ubuntu, but it should work *at least* under Debian too.

It basically:

* check and install "gcc", "gfortran", "make", "git", "f2c"
* under Ubuntu: install "libssl-deb" (necessaty tom compile CMAKE)
* git clone or git pull CMAKE sources
* compile CMAKE if git repo has changed
* git clone or git pull Mystran sources
* compile Mystran if git repo has changed

Do not hesitate to give me your feedback and issues.

Print this item

  Python Neutral Parse
Posted by: numenic - 04-16-2021, 08:02 AM - Forum: Pre/Post Processors - No Replies

I wrote a lightweight python femap neutral file parser to ease automatic post-processing:

It's available on pypi: https://pypi.org/project/femap-neutral-parser/
repository: https://framagit.org/numenic/femap_neutral_parser/

Print this item

  FEMAP neutral improvements
Posted by: numenic - 04-09-2021, 03:44 PM - Forum: Pre/Post Processors - Replies (7)

Writing a FEMAP neutral parser, I met some some quirks with the file generated by MYSTRAN (v12.1). This thread aims to collect them. Please tell me if I need to split this thread per subject or question!

I'll begin with this a simpleCBAR + CBUSH model. The attached archive provides the following files:

* "bulk_model_TC3.nas": the BULK file
* "mystran_00.NEU": the MYSTRAN output
* "FEMAP_v8.2.neu": the FEMAP v8.2 output
* "FEMAP_v2020-1-0.neu": the FEMAP v2020.1 output

First remark: Block451 is used by both MYSTRAN and FEMAP8.2 to collect results, but is deprecated in later versions in favor of block 1051.

Block451: titles used by MYSTRAN do not match FEMAP titles. This makes harder to write "solver-agnostic" parsers. From the uploaded files, here is the list with mismatch and missing :

That's all for now! thanks for reading and happy MYSTRANning !

Code:
| FEMAP vector Title        | FEMAP vector ID | MYSTRAN vector Title     | MYSTRAN vector ID |
|---------------------------+-----------------+--------------------------+-------------------|
| Total Translation         | 1               | RSS translation          | 10001             |
| T1 Translation            | 2               | T1  translation          | 10002             |
| T2 Translation            | 3               | T2  translation          | 10003             |
| T3 Translation            | 4               | T3  translation          | 10004             |
| Total Rotation            | 5               | RSS rotation             | 10005             |
| R1 Rotation               | 6               | R1  rotation             | 10006             |
| R2 Rotation               | 7               | R2  rotation             | 10007             |
| R3 Rotation               | 8               | R3  rotation             | 10008             |
| Total Applied Force       | 41              | RSS applied force        | 20001             |
| T1 Applied Force          | 42              | T1  applied force        | 20002             |
| T2 Applied Force          | 43              | T2  applied force        | 20003             |
| T3 Applied Force          | 44              | T3  applied force        | 20004             |
| Total Applied Moment      | 45              | RSS applied moment       | 20005             |
| R1 Applied Moment         | 46              | R1  applied moment       | 20006             |
| R2 Applied Moment         | 47              | R2  applied moment       | 20007             |
| R3 Applied Moment         | 48              | R3  applied moment       | 20008             |
| Total Constraint Force    | 51              | RSS SPC force            | 30001             |
| T1 Constraint Force       | 52              | T1  SPC force            | 30002             |
| T2 Constraint Force       | 53              | T2  SPC force            | 30003             |
| T3 Constraint Force       | 54              | T3  SPC force            | 30004             |
| Total Constraint Moment   | 55              | RSS SPC moment           | 30005             |
| R1 Constraint Moment      | 56              | R1  SPC moment           | 30006             |
| R2 Constraint Moment      | 57              | R2  SPC moment           | 30007             |
| R3 Constraint Moment      | 58              | R3  SPC moment           | 30008             |
| Total MultiPoint Force    | 151             | na                       | na                |
| T1 MultiPoint Force       | 152             | na                       | na                |
| T2 MultiPoint Force       | 153             | na                       | na                |
| T3 MultiPoint Force       | 154             | na                       | na                |
| Total MultiPoint Moment   | 155             | na                       | na                |
| R1 MultiPoint Moment      | 156             | na                       | na                |
| R2 MultiPoint Moment      | 157             | na                       | na                |
| R3 MultiPoint Moment      | 158             | na                       | na                |
| Total Summed GPForce      | 161             | na                       | na                |
| T1 Summed GPForce         | 162             | na                       | na                |
| T2 Summed GPForce         | 163             | na                       | na                |
| T3 Summed GPForce         | 164             | na                       | na                |
| Total Summed GPMoment     | 165             | na                       | na                |
| R1 Summed GPMoment        | 166             | na                       | na                |
| R2 Summed GPMoment        | 167             | na                       | na                |
| R3 Summed GPMoment        | 168             | na                       | na                |
| Total Applied GPForce     | 171             | na                       | na                |
| T1 Applied GPForce        | 172             | na                       | na                |
| T2 Applied GPForce        | 173             | na                       | na                |
| T3 Applied GPForce        | 174             | na                       | na                |
| Total Applied GPMoment    | 175             | na                       | na                |
| R1 Applied GPMoment       | 176             | na                       | na                |
| R2 Applied GPMoment       | 177             | na                       | na                |
| R3 Applied GPMoment       | 178             | na                       | na                |
| Total Constraint GPForce  | 181             | na                       | na                |
| T1 Constraint GPForce     | 182             | na                       | na                |
| T2 Constraint GPForce     | 183             | na                       | na                |
| T3 Constraint GPForce     | 184             | na                       | na                |
| Total Constraint GPMoment | 185             | na                       | na                |
| R1 Constraint GPMoment    | 186             | na                       | na                |
| R2 Constraint GPMoment    | 187             | na                       | na                |
| R3 Constraint GPMoment    | 188             | na                       | na                |
| Bar EndA Plane1 Moment    | 3000            | BAR EndA Plane1 Moment   | 50201             |
| Bar EndA Plane2 Moment    | 3001            | BAR EndA Plane2 Moment   | 50203             |
| Bar EndB Plane1 Moment    | 3002            | BAR EndB Plane1 Moment   | 50202             |
| Bar EndB Plane2 Moment    | 3003            | BAR EndB Plane2 Moment   | 50204             |
| Bar EndA Pl1 Shear Force  | 3004            | BAR EndA Pl1 Shear Force | 50205             |
| Bar EndA Pl2 Shear Force  | 3005            | BAR EndA Pl2 Shear Force | 50207             |
| Bar EndA Axial Force      | 3008            | BAR EndA Axial Force     | 50209             |
| Bar EndA Torque           | 3010            | BAR EndA Torque          | 50211             |
| Bar EndA Pt1 Bend Stress  | 3075            | BAR EndA Pt1 Comb Stress | 60301             |
| Bar EndA Pt2 Bend Stress  | 3076            | BAR EndA Pt2 Comb Stress | 60303             |
| Bar EndA Pt3 Bend Stress  | 3077            | BAR EndA Pt3 Comb Stress | 60305             |
| Bar EndA Pt4 Bend Stress  | 3078            | BAR EndA Pt4 Comb Stress | 60307             |
| Bar EndB Pt1 Bend Stress  | 3083            | BAR EndB Pt1 Comb Stress | 60302             |
| Bar EndB Pt2 Bend Stress  | 3084            | BAR EndB Pt2 Comb Stress | 60304             |
| Bar EndB Pt3 Bend Stress  | 3085            | BAR EndB Pt3 Comb Stress | 60306             |
| Bar EndB Pt4 Bend Stress  | 3086            | BAR EndB Pt4 Comb Stress | 60308             |
| Bar EndA Axial Stress     | 3107            | na                       | na                |
| Bar EndA Max Comb Stress  | 3109            | BAR EndA Max Stress      | 60309             |
| Bar EndA Min Comb Stress  | 3110            | BAR EndA Min Stress      | 60311             |
| Bar EndB Max Comb Stress  | 3111            | BAR EndB Max Stress      | 60310             |
| Bar EndB Min Comb Stress  | 3112            | BAR EndB Min Stress      | 60312             |
| Bar Tension M.S.          | 3113            | na                       | na                |
| Bar Compression M.S.      | 3114            | na                       | na                |
| Bush X Force              | 3774            | na                       | na                |
| Bush Y Force              | 3775            | na                       | na                |
| Bush Z Force              | 3776            | na                       | na                |
| Bush X Moment             | 3777            | na                       | na                |
| Bush Y Moment             | 3778            | na                       | na                |
| Bush Z Moment             | 3779            | na                       | na                |
| Bush TX Stress            | 3780            | na                       | na                |
| Bush TY Stress            | 3781            | na                       | na                |
| Bush TZ Stress            | 3782            | na                       | na                |
| Bush RX Stress            | 3783            | na                       | na                |
| Bush RY Stress            | 3784            | na                       | na                |
| Bush RZ Stress            | 3785            | na                       | na                |
| Elem C1 T1 GPForce        | 85000           | na                       | na                |
| Elem C1 T2 GPForce        | 85001           | na                       | na                |
| Elem C1 T3 GPForce        | 85002           | na                       | na                |
| Elem C1 R1 GPMoment       | 85003           | na                       | na                |
| Elem C1 R2 GPMoment       | 85004           | na                       | na                |
| Elem C1 R3 GPMoment       | 85005           | na                       | na                |
| Elem C2 T1 GPForce        | 85010           | na                       | na                |
| Elem C2 T2 GPForce        | 85011           | na                       | na                |
| Elem C2 T3 GPForce        | 85012           | na                       | na                |
| Elem C2 R1 GPMoment       | 85013           | na                       | na                |
| Elem C2 R2 GPMoment       | 85014           | na                       | na                |
| Elem C2 R3 GPMoment       | 85015           | na                       | na                |
| na                        | na              | BAR EndB Pl1 Shear Force | 50206             |
| na                        | na              | BAR EndB Pl2 Shear Force | 50208             |
| na                        | na              | BAR EndB Axial Force     | 50210             |
| na                        | na              | BAR EndB Torque          | 50212             |



Attached Files
.zip   mystran_femap_q1.zip (Size: 17.52 KB / Downloads: 2)
Print this item

  ComLab Enhancments
Posted by: Admin - 03-16-2021, 02:49 PM - Forum: ComLab and NASTRAN-95 - No Replies

Here is a list of enhancements that ComLab has:

https://www.mystran.com/comlab/ComLab_Enhancements.pdf

Print this item

  Added a Parameter that Removes Bandit Files from final Run Directory
Posted by: ZacharyALerner - 03-15-2021, 03:09 AM - Forum: MYSTRAN - No Replies

Howdy all.

With aid from Bill, I've added a parameter to MYSTRAN called PARAM,DELBAN that deletes the bandit files from the home run directory.  You can find the modified files in my respository here.

For future reference, the files needed (at the barest minimum to merely function) in order to add a parameter to MYSTRAN are

1.) Source/Interfaces/BD_PARAM_Interface.f90 where other subroutines can read the types of all parameter entities.
2.) Source/LK1/L1A-BD/BD_PARAM.f90  which processes the PARAM bulk data cards
3.) Source/Modlues/PARAMS.f90  where you define and comment the default value of your parameter, and describe what it does.

Otherwise, all that was needed for this parameter to be parsed was to add a conditional into Source/Modules/BANDIT/BANDIT_MODULE.f to deal set a close status flag to 'delete' instead of 'keep'.  It's worth noting that my compiler was throwing a fit over my if statements for some reason here, I think that because this BANDIT subroutine contains some F77 legacy syntax, the newer compilers don't like some of the looser syntax available in the f90 and f95 scripts. Future tinker'ers beware.

Also I've added the bones of a README revision for the github page, as I really like the quality of Ceanwang and Bruno's works on documentation and I'd love to take it on myself to beautify them for publication in any way that you all think is best.

Let me know you think and what else I should do.

- Zach

Print this item

  NEWS: ComLab Released to the Public
Posted by: Admin - 03-10-2021, 03:07 AM - Forum: MYSTRAN - Replies (1)

Harry Schaeffer has released ComLab to the public. It is an improved version of NASTRAN-95. For more information, see this thread:

https://www.mystran.com/forums/showthread.php?tid=58

Print this item

  Compiler Error for Preferred Option - Windows
Posted by: Admin - 02-23-2021, 06:24 PM - Forum: MYSTRAN - Replies (1)

[ 97%] Linking Fortran executable Binaries\mystran.exe D:/Program_Files/MYSTRAN/msys2/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/10.2.0\libgfortran.a(read.o)Sad.text$_gfortrani_convert_real+0xb7): relocation truncated to fit: IMAGE_REL_AMD64_REL32 against undefined symbol
strtoflt128'
D:/Program_Files/MYSTRAN/msys2/mingw64/bin/../lib/gcc/x86_64-w64-mingw32/10.2.0\libgfortran.a(read.o)Sad.text$_gfortrani_convert_infnan+0xb3): relocation truncated to fit: IMAGE_REL_AMD64_REL32 against undefined symbol
strtoflt128' collect2.exe: error: ld returned 1 exit status mingw32-make[2]: [CMakeFiles\mystran.dir\build.make:23351: Binaries/mystran.exe] Error 1 mingw32-make[1]: [CMakeFiles\Makefile2:562: CMakeFiles/mystran.dir/all] Error 2 mingw32-make: *** [Makefile:148: all] Error 2

Print this item

  Femap with Mystran
Posted by: SamiH. - 02-17-2021, 09:22 AM - Forum: Pre/Post Processors - Replies (2)

Hi,

I am experimenting with Mystran with Femap and would like to share and have comments on the observations, and after all would appreciate help to adjust the settings correctly. These may be sometimes Femap related questions or issues but also Mystran related. I have only done one analysis and haven't checked the Mystran documentation very carefully so that these might be more of a rookie type questions/observations.

First of all I am impressed how well the key results are in line with Simcenter Nastran results without any adjustments at all. Here's what I've done:
- opened up a random FE model consisting of plate elements, some glue contacts, and few RBE2 elements. The analysis for this had been done earlier using Simcenter Nastran.
- exported the analysis model (without any adjustments, i.e. Simcenter Nastran version)
- run it with Mystran
- imported the Femap neu file

The few things that caught my eye comparing the results:
- The deformation looks very similar by eye
- the max total translation in Simcenter Nastran was 1.883 mm and RSS translation in Mystran 1.909 mm, so no big difference (1.4 %)
- the max vonMises bot stress differs quite a lot but this is difficult to compare. See next observations. 
- the result scale (contour levels) is something odd. For these translations (max 1.909 mm) the scaling is from 0 to 1.826E+10 when it should be from 0 to 1.909. Of course I can change it manually but for automatic detection it seems to have some kind of scaling issue.
- the scaling issue seems to be the problem for dynamic rotation also. When tried to rotate the model even slightest of rotation will throw the model off the screen somewhere in the space and in practice the rotation can not be done.
- the element results are shown in transparent type of way in which it is difficult to observe which element is which when front and back and all in between are in a same spot. I couldn't figure out how to change this. This is not the same transparency as in Femap but something related to the results.
- the result vectors are numbered very differently from Nastran. There is also separated results for QUADs and TRIAs which is perhaps not very necessary but inconvenient sometimes

Has anyone any ideas how to improve the performance of Mystran results in Femap related to above mentioned observations? Can I adjust something somewhere?

Print this item

  ComLab (an improved NASTRAN-95) has been released to the public
Posted by: Admin - 02-10-2021, 07:47 AM - Forum: ComLab and NASTRAN-95 - No Replies

Harry Schaeffer has created an improved version of NASTRAN-95 called ComLab. It was previously a commercial program, but he has now released it to the public and we are hosting it here. We are looking into adding the source code at a later date as well.

Full Package is on hold.
https://www.mystran.com/comlab/CML.pdf - User Manual (part of the full package)
https://www.mystran.com/comlab/INSTALL.pdf - Installation Guide (part of the full package)
https://www.mystran.com/comlab/COMLAB_EULA.pdf - License Agreement (note that Harry has released ComLab to the public now so this may not be a current document, but is provided for reference)

Harry wrote one of the first books about NASTRAN and it was published in 1979. He also helped release one of the first complied versions of NASTRAN-95
https://www.amazon.com/Nastran-Primer-St...B000O15X8Y

Print this item

  Successful addition of sparse solver SuperLU to MYSTRAN
Posted by: drbillc - 12-21-2020, 10:12 PM - Forum: MYSTRAN - No Replies

With the help of Cean Wang I was able to get the sparse solver SuperLU downloaded and compiled into 2 libraries, libsuperlu.a and libblas.a . Using these, along with the c routine: c_fortran_dgssv.o I modified MYSTRAN's LINK3 solution of the L-set displacements (from the equation KLL*UL = PL) to run using SuperLU as a sparse solver. I ran a series of test problems, for which I have known solutions, to see how SuperlU performed. The series of test problems were of a square, simply supported, thick plate under a pressure load. I ran successively larger meshes up to a maximum of 800x800 plate elements with 1,920,000 L-set degrees of freedom (DOF). SuperLU ran all of the test problems successfully up to the almost 2 million degree of freedom problem (but took an excessive amount of time for this largest DOF run). The banded solver (LAPACK) couldn't come close to running anything this large. It would have required almost 40 GB of storage for the banded KLL array LAPACK uses. LAPACK couldn't even run the smaller 500x500 mesh problem with 750,000 L-set DOF. LAPACK needed 9 GB in this case but Windows will not allocate more than 2 GB to any array

There are a number of other subroutines (approximately 10) in MYSTRAN that need to be modified to take advantage of this new sparse solver. For example, reduction of the G-set stiffness matrix, KGG, to the N-set involves eliminating the dependent DOF identified on MPC's and rigid elements. However, for the test problems mentioned above, none of these effects in the 10 additional subroutines were needed and for many problems in general the inclusion of the sparse solver that has been done for LINK3 will be sufficient. Over the next few months I will try to get the other routines reprogrammed to use SuperLU but, in the meantime, most simple statics problems will run the new sparse solver very effectively.

I will post the updated MYSTRAN, with SuperLU incorporated into LINK3, to GitHub later this week

Cean has a thread with the procedure to compile SuperLU. If anyone uses it ignore his comments on how to incorporate into MYSTRAN. Only use up to step 6 in his procedure (outlined below):

From Cean:

  1. Download and execute the cmake msi binary from www.cmake.org/download. This will install cmake into C:\Program Files\Cmake
  2. Download SuperLU as a zip file, unzip it to D:\MYSTRAN\Cmake
  3. Start Windows Terminal and change to D:\MYSTRAN\SuperLU (Use FreeCommander XE program is more easy to do this.)
  4. Configure with my w.bat script as before,
  5. Run make to build SuperLU
then, to test SuperLU prior to incorporating into MYSTRAN:
  1. After building SuperLU successfully, change to D:\MYSTRAN\SuperLU\FORTRAN
  2. Run make to build the FORTRAN example.
  3. You could look at f77_main.f to see how it is called.

Print this item