HELP!

A forum dedicated to feature requests and the future direction of OpenSees, i.e. what would you like, what do you need, what should we explore

Moderators: silvia, selimgunay, Moderators

QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

HELP!

Post by QI »

hello,i am a new user of the opensees.
If i want to add a groundMotion with diffrent values to some diffrent nodes of frame construction respectively,how can i do?
谢谢!
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

here is a piece of script that I use:

Code: Select all

   if {$LoadPatternType == "MultipleSupport"} {
      # read a PEER strong motion database file, extracts dt from the header and converts the file 
      # to the format OpenSees expects for Unifform ground motions 
      set inFile1 $GMdirectory$GMFileName1;	# define one per node
      set inFile2 $GMdirectory$GMFileName2
      set modtxt Mod
      set outFile1 $GMdirectory$modtxt$GMFileName1;   # set variable holding new filename # define one per node
      set outFile2 $GMdirectory$modtxt$GMFileName2;   # set variable holding new filename 
      procReadSMDFile $inFile1 $outFile1 dt;      # call procedure to convert the ground-motion file# define one per node
      procReadSMDFile $inFile2 $outFile2 dt;      # call procedure to convert the ground-motion file
      set dispSeries1 "Series -dt $dt -filePath $outFile1 -factor  $GMfattA ";   # time series information# define one per node
      set dispSeries2 "Series -dt $dt -filePath $outFile2 -factor  $GMfattA ";   # time series information
      pattern MultipleSupport $IDloadTagGMA  {
	groundMotion $IDgmSeries1 Plain -disp  $dispSeries1 ;	# define one per node
	groundMotion $IDgmSeries2 Plain -disp  $dispSeries2

	imposedSupportMotion $SupNode1  $IDdofA  $IDgmSeries1 ;	# define one per node
	imposedSupportMotion $SupNode2  $IDdofA  $IDgmSeries2 
      };   # end pattern
   
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

thank you for your help!
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

Excuse me
what does IDgmSeries mean?
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

sorry!
I know now.
thank you
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

Dear silvia!:
can you help me find out what's wrong with my code about pattern MultipleSupport object.
in additions,i want to make acquaintance with with people,such as azari,who are studying MPA. can you tell me their email or other.thank you!


code:
# Define dynamic loads

# Set some parameters

set outFile gmrecord/${GM}.g3

#set factor [expr $accelD8/abs($accel($GM))]
set factor 1.5

set modi1 1.317008924
set modi2 -1.216449135

set GMfatt1 [expr $factor*$modi1]
set GMfatt2 [expr $factor*$modi2]


# Set time series to be passed to uniform excitation
set disSeries1 "Series -dt $datadt -filePath $outFile -factor $GMfatt1 "
set disSeries2 "Series -dt $datadt -filePath $outFile -factor $GMfatt2 "

# Create imposedMotion pattern


pattern MultipleSupport 2 {
groundMotion 1 Plain -dis $disSeries1
groundMotion 2 Plain -dis $disSeries2

# tag dir IDGM
imposedSupportMotion 1 1 1 # define dynamic loads for every node
imposedSupportMotion 2 1 2

}
# end pattern
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

what do you mean, what is wrong with it?

also, i can't give you e-mail addresses of registered people, you need to post a request to make contact with him/her, or anyone doing MPA.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

when i run the code,it say
"WARNING - ForceBeamColumn2d::update - failed to get compatible element forces &
deformations for element: 9(dW: << 98.7027)
Domain::update - domain failed in update
Newmark::update() - failed to update the domain
WARNING NewtonRaphson::solveCurrentStep() -the Integrator failed in update()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.005
OpenSees > analyze failed, returned: -3 error flag
regular newton failed .. lets try an initail stiffness for this step
WARNING - ForceBeamColumn2d::update - failed to get compatible element forces &
deformations for element: 3(dW: << 1323.47)
WARNING - ForceBeamColumn2d::update - failed to get compatible element forces &
deformations for element: 12(dW: << 10035.9)
Domain::update - domain failed in update
Newmark::update() - failed to update the domain
WARNING ModifiedNewton::solveCurrentStep() -the Integrator failed in update()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.005
OpenSees > analyze failed, returned: -3 error flag
Transient analysis completed FAILED"
meanwhile, i am sure the mistake must happen here!
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

BY the way.there is a groudmotion file about accel which name is usa00707.g3 ,is it necessary to extracts dt from the header and converts the file ?
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

how can i get the groundmotion volocities and displacement from groundmotion accelelations ?(timeseriesintegrator?)
and why does the pattern MultipleSupport object need volocities and displacement records?
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

it needs as displacement record because that is how it was formulated. Frank may know why better. it needs velocities for initial conditions.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

very sorry for disturb you for so long time!
but I still can't understand how i should do and the key problem is that: now i have a groudmotion file with only accel record whose name is usa00707.g3 ,how should i do to add this record to every node of one frame construture with diffrent value respectively?
thank you!
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

if you want different input for each node, you need to specify each node individually.
but you have to convert your acceleration to displacement -- most records have both available.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
QI
Posts: 39
Joined: Sat Aug 05, 2006 4:25 am

Post by QI »

hello,silvia,DO you mind giving me a el centro ground motion about displacement? becauce i can't find it.thankyou!
my Email is wenqiguang81412@126.com
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

look in here, first:

PEER Strong Motion Database

http://peer.berkeley.edu/smcat/index.html
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Post Reply