The OpenSees Community

MikeITExpert

This also clarifies how to run the OpenSeesMP on windows

https://opensees.berkeley.edu/OpenSees/ ... indows.pdf

MikeITExpert

I am trying to make use of OpenSeesMP to speed up my analysis. I have been studying the below link

https://opensees.berkeley.edu/OpenSees/ ... SeesMP.pdf

On page 9 and 10, it describes how to analyze a large model to run in parallel, however, I can't make sense of domainChange command.

Please let me know if there is a more detailed reference on how to make large model analyze in parallel using OpenSeesMP.

Regards

MikeITExpert

I downloaded an older version of MPICH from https://www.mpich.org/static/downloads/1.4/ and that created mpi dlls and placed them in c:\Windows\System32.

If I had time I would have compiled a recent version.

MikeITExpert

Ok, the below error had to do with the -np flag which sets the number of processors. When I set -np to actual number of processors the above error doesn't occur. However, I still get missing mpich2mpi.dll errors. Very weird because I have installed "Microsoft MPI v10.0" already.

MikeITExpert

The only way I could get things to work was using cygwin ...

Any idea what the error is ...

$ mpirun -np 10 OpenSeesMP.exe run.tcl
--------------------------------------------------------------------------
There are not enough slots available in the system to satisfy the 10 slots
that were requested by the application:
OpenSeesMP.exe

Either request fewer slots for your application, or make more slots available
for use.
--------------------------------------------------------------------------

MikeITExpert

Hi All,

I am trying to run the OpenSeesMP. To achieve that I have installed Microsoft MPI, however, it keeps complaining about the following two dlls.

Please let me know if you know any workaround.

Regards,

mpich2mpi.dll
fmpich2.dll

MikeITExpert

I exported my model using print -JSON -file somefile.json. However, I get to many elements with same tag number.

I am wondering if that is print issue or element with same tags are allowed ?

Please let me know
Regards

Code: Select all

			{"name": 1.20011e+06, "type": "ZeroLength", "nodes": [1001, 100102], "materials": ["40", "40", "40", "40", "1.20011e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.20011e+06, "type": "ElasticBeam3d", "nodes": [120103, 130102], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.20021e+06, "type": "ZeroLength", "nodes": [1002, 100202], "materials": ["40", "40", "40", "40", "1.20021e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.20021e+06, "type": "ElasticBeam3d", "nodes": [120203, 130202], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.20031e+06, "type": "ZeroLength", "nodes": [1003, 100302], "materials": ["40", "40", "40", "40", "1.20031e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.20031e+06, "type": "ElasticBeam3d", "nodes": [120303, 130302], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.20041e+06, "type": "ZeroLength", "nodes": [1004, 100402], "materials": ["40", "40", "40", "40", "1.20041e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.20041e+06, "type": "ElasticBeam3d", "nodes": [120403, 130402], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.20051e+06, "type": "ZeroLength", "nodes": [1005, 100502], "materials": ["40", "40", "40", "40", "1.20051e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.20051e+06, "type": "ElasticBeam3d", "nodes": [120503, 130502], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.20061e+06, "type": "ZeroLength", "nodes": [1006, 100602], "materials": ["40", "40", "40", "40", "1.20061e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.20061e+06, "type": "ElasticBeam3d", "nodes": [120603, 130602], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.30011e+06, "type": "ZeroLength", "nodes": [1001, 100103], "materials": ["40", "40", "40", "40", "1.30011e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.30011e+06, "type": "ElasticBeam3d", "nodes": [130103, 140102], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.30021e+06, "type": "ZeroLength", "nodes": [1002, 100203], "materials": ["40", "40", "40", "40", "1.30021e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.30021e+06, "type": "ElasticBeam3d", "nodes": [130203, 140202], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.30031e+06, "type": "ZeroLength", "nodes": [1003, 100303], "materials": ["40", "40", "40", "40", "1.30031e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.30031e+06, "type": "ElasticBeam3d", "nodes": [130303, 140302], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.30041e+06, "type": "ZeroLength", "nodes": [1004, 100403], "materials": ["40", "40", "40", "40", "1.30041e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.30041e+06, "type": "ElasticBeam3d", "nodes": [130403, 140402], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.30051e+06, "type": "ZeroLength", "nodes": [1005, 100503], "materials": ["40", "40", "40", "40", "1.30051e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.30051e+06, "type": "ElasticBeam3d", "nodes": [130503, 140502], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},
			{"name": 1.30061e+06, "type": "ZeroLength", "nodes": [1006, 100603], "materials": ["40", "40", "40", "40", "1.30061e+06", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "transMatrix": [[1, 0, 0], [0, 1, 0], [0, 0, 1]]},
			{"name": 1.30061e+06, "type": "ElasticBeam3d", "nodes": [130603, 140602], "E": 31622.8, "G": 1e+15, "A": 337500, "Jx": 1, "Iy": 6.96094e+09, "Iz": 2.50594e+09, "massperlength": 0, "crdTransformation": "2"},

MikeITExpert

I am trying to get the list of all elements and nodes so I use ...

Code: Select all

print -JSON -file model.json

but then the output would be different from a typical print -ele command.

Code: Select all

print -ele 9001

From JSON it would look like below

Code: Select all

{"name": 9001, "type": "TwoNodeLink", "nodes": [1201, 1.0012e+06], "materials": ["412", "40", "40", "40", "40", "40"], "dof": ["P", "Vy", "Vz", "T", "My", "Mz"], "sDratios": [0.5, 0.5], "transMatrix": [[0, 0.812359, -0.583158], [-1, 0, 0], [0, 0.583158, 0.812359]],"addRayleigh": 0, "mass": 0},

Using print -ele it would look like below

Code: Select all

Element: 9001
  type: TwoNodeLink
  iNode: 1201, jNode: 1001201
  Material dir0: 412
  Material dir1: 40
  Material dir2: 40
  Material dir3: 40
  Material dir4: 40
  Material dir5: 40
  Mratio:   shearDistI: 0.5 0.5

  addRayleigh: 0  mass: 0
  resisting force: 5.42101e-16 9075.73 -6515.08 2.41573e-10 -4.08609e-13 -5.69206e-13 -5.42101e-16 -9075.73 6515.08 5.54403e-10 -4.08608e-13 -5.69206e-13

Please let me know if there is workaround.

Regards

MikeITExpert

I have helping a civil engineer to automate part of his work ...

I need to have the new state of the model after running a pushover analysis.

The starting state I can retrieve using "print -JSON -file model.json" command. However, if I try the same after analysis I seem to get the initial state of the model.

Am I doing anything wrong?

Please let me know if there is workaround.

Regards

MikeITExpert

Hi Scott,

I was able to find the issue and you were exactly right. I had couple of elements with very stiff elastic martials as listed below:

uniaxialMaterial Elastic 40 1.0e15
uniaxialMaterial Elastic 3001 [expr 1.0e15]

I know that there are different types of problem formation depending on type of the integrator we are using. In the model I am dealing with there are using "integrator LoadControl" command (one for gravity analysis and the other for pushover analysis [very weird they are using LoadControl as opposed to integrator DisplacementControl ] ) which I believe uses "Direct Stiffness Method" (please let me if I am wrong, cuz I am coming from OOP C/C++ programming background with a little bit of numerical analysis background and linear algebra. not very good at matrix algebra).I know matrix inverses can't be calculated when there are very large values in them.

I guess what I need to know are

1) The problem formulation given LoadControl is used as integrator. (I believe it is Direct stiffness method, please let me know)
2) What matrixes are inversed during analysis iterations? (I probably would be able to figure that out given the problem formulation.)
3) What elements, nodes and constraints need to change to resolve any analysis issue? (IOW, how cells of the matrix related to model parameters such as element materials, which again I should be able to find out if I know the problem formulation when using a integrator LoadControl.)

I think all I need at this point is the problem formulation when using integrator LoadControl and OpenSees print statements to get matrices.

Please let me know your valuable comments.

Regards

MikeITExpert

That is correct!

Thanks for the link.

MikeITExpert

I am new to OpenSees and I am trying to define a Two Node Link Element.

My physics knowledge dates back to Halliday physics course back in university. I am just curious how many material tags do I need in order to define a Two Node Link Element and what they are representing.

Much appreciate your comments.
Regards

https://opensees.berkeley.edu/wiki/inde ... nk_Element

MikeITExpert

I checked the manual below, it seems to be optional.

https://opensees.berkeley.edu/wiki/inde ... nk_Element

Can please explain a bit more?

Much appreciate it.

MikeITExpert

Dear Experts,

I am helping a civil engineer solving a structure gravity and pushover analysis load. However, I noticed that when the number of diagonal elastic elements (going to be replaced by viscus dampers later in dynamic seismic analysis) is increased, there is lower chance of convergence. I have tried Krylov Newton and Raphson Newton both.

I am curious whether OpenSees solve the least square problem ? How element constraints are taken into consideration by OpenSee?

Please let me know your valuable comments. Feel free to point me to some documents / software specs.

Regards,

MikeITExpert

Hi

I a getting many warnings messages like below. I am just curious what is that suggesting exactly?

Please let me know your valuable comments if any.

Regards

Code: Select all

WARNING TwoNodeLink::setUp() - element: 6048
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 6148
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 7048
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 7148
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 9048
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 9148
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 10048
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 10148
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 11048
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 11148
no local y vector specified
WARNING TwoNodeLink::setUp() - element: 12048
no local y vector specified

The OpenSees Community

Search found 24 matches

Re: How to install/setup OpenSeesMP.exe / parallel on windows?

what does domainChange command do to parallel a model and when that is required?

Re: How to install/setup OpenSeesMP.exe / parallel on windows?

Re: How to install/setup OpenSeesMP.exe / parallel on windows?

Re: How to install/setup OpenSeesMP.exe / parallel on windows?

How to install/setup OpenSeesMP.exe / parallel on windows?

Is it possible to have two elements with the same element tag?

Why element node in "print -JSON" is different from "print -ele"

how to save nodes' new coordinate after an analysis?

Re: Does OpenSees solve least square problem?

Re: Does OpenSees solve least square problem?

how many element tags do I need to define a Two Node Link Element?

Re: WARNING TwoNodeLink::setUp() - element: XXX

Does OpenSees solve least square problem?

WARNING TwoNodeLink::setUp() - element: XXX