The OpenSees Community

Thanks Rjaeger, Code works when we define two double variables and then use them to build the exe as follows: double alpha1 = 1.0; double alpha2 = 1.0; ConstraintHandler *theEigenHandler = new LagrangeConstraintHandler(alpha1, alpha2); also we can define double constant values as follows: (it also w...

Thank You for your reply. As quoted: Hello, Neallee is correct. If you look at the penalty constraint handler header file, you'll see the only constructor is: PenaltyConstraintHandler(double alphaSP, double alphaMP); Check the manual or the reference Neallee mentioned to find out what these paramete...

Thanks a lot Prayag.

This option should be included in manual.

Thanks again

I want to change the constraint handler for eigenvalue analysis in my code. For this I opened file OpenSees\SRC\tcl\commands.cpp and replaced the line (line no. 4563): ConstraintHandler *theEigenHandler = new TransformationConstraintHandler(); with the line ConstraintHandler *theEigenHandler = new P...

Hello, I made a very very small script that might be useful to some of OpenSees users. Earlier I had to use a stopwatch to see what is the time required for a particular analysis and so I used following script. Place the following line at the start of master tcl input file (say all.tcl) set startTim...

I have modelled a three storyed structure (3D) with Y as vertical direction. I have applied masses in X and Z direction at each floor to CM where CM is connected to other nodes of floor by rigid diaphragm. Theoretically it should give six modes of vibrations, but it is giving only upto four, as soon...

Sir, I want to recompile the source code to make a new executable. For this the instructions given on webpage is We have provided you with a Visual C++ project, opensees.dsw, located in the OpenSees/Win32 directory. Click on this file to open it and then click on build to build it. However, be sure ...

Thanks a lot for your kind guidance. I think the changes are to be made in following commands: // create the rest of components of an eigen analysis EigenIntegrator *theEigenIntegrator = new EigenIntegrator(); RCM *theRCM = new RCM(false); DOF_Numberer *theEigenNumberer = new DOF_Numberer(*theRCM); ...

Sir,
If I change the constraint handler to Lagrange, the system of equation does not remain symmetric positive definite. Is the equation solver used by default is for non-SPD equations, or do I need to change the equation solver also.

Thanks

I have modeled a 3D geometry in openSees with beams pinned at the ends. The geometry also has rigid diphragm effect. Thus there is sequential slaving for nodes. If I remove the pinned beam slaving, I get first two modes in X and Z direction but when I use pinned beam slaving, I get both two modes in...

Thanks sir, I will try by this mehtod.

Thanks again.

Thank You sir for your kind attention. I have checked all the sections, their orientations, node coordinates, elements, etc. 3-4 times, they all are correct. The interesting problem is that sum of all reactions forces at foundation nodes is NOT equal to the force applied at top (i.e. 100) in Z direc...

Can anybody please find the bug in program. I am working on this from last fifteen days but I am not able to find out the error. Here I am using same elements in X and Z directions, just with different transformations so as to orient members have stronger axis normal to the plane of moment resisting...

No sir, the command geomTransf Linear 1 0 0 -1 will make local z axis parallel to +ve global Z axis. I have cross checked it, and also it is clear from figure in Command Language Manual on page 296 (chapter 25). vecxz is opposit to local z axis in this figure . Same image it attached herewith. http:...

Interestingly the problem is not yet solved. If I am not wrong then I have made local z axis parallel to +ve global Z axis by following transformation geomTransf Linear 1 0 0 -1 so local y axis is parallel -ve global X axis and local x axis is parallel to +ve global Y axis. The new output for local ...