no its not working my friend....
Why its not working if we change the units ?
I have checked the whole code but didnt find any error .
Moderators: silvia, selimgunay, Moderators
yes sir...
I have changed the value inch in to meter directly and 1 kip to 4448.2216 Newton.....
Code: Select all
WARNING - no torsion specified for 3D fiber section, use -GJ or -torsion
FiberSection3d: 4
while executing
"section fiberSec $id {
# Define the core patch
patch quadr $coreID $nfCoreZ $nfCoreY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ
..."
(procedure "BuildRCrectSection" line 75)
if you want I will attacching a link for my script. I used "N" and "m" and you can try to see if the analysis convergesankurjain wrote: ↑Fri Aug 28, 2020 9:07 amI have changed the value inch in to meter directly and 1 kip to 4448.2216 Newton.....
but still the convergence is an issue.....
set in 0.0254; # define basic units -- output units
set kip 4448.2216; # define basic units -- output units
set sec 1.; # define basic units -- output units
set LunitTXT "inch"; # define basic-unit text for output
set FunitTXT "kip"; # define basic-unit text for output
set TunitTXT "sec"; # define basic-unit text for output
set ft [expr 12.*$in]; # define engineering units
set ksi [expr $kip/pow($in,2)];
set psi [expr $ksi/1000.];
set lbf [expr $psi*$in*$in]; # pounds force
set pcf [expr $lbf/pow($ft,3)]; # pounds per cubic foot
set psf [expr $lbf/pow($ft,3)]; # pounds per square foot
set in2 [expr $in*$in]; # inch^2
set in4 [expr $in*$in*$in*$in]; # inch^4
set cm [expr $in/2.54]; # centimeter, needed for displacement input in MultipleSupport excitation
set PI [expr 2*asin(1.0)]; # define constants
set g [expr 32.2*$ft/pow($sec,2)]; # gravitational acceleration
set Ubig 1.e10; # a really large number
set Usmall [expr 1/$Ubig]; # a really small number
ok i will try this and revert backgeoStru wrote: ↑Sat Aug 29, 2020 6:55 amSorry Michael H. Scott, there is a misunderstanding. I resolved the problem of the convergence. I would like to help ankurjain.
For the problem "no torsion specified for 3D fiber section, use -GJ or -torsion" you should add in BuildRCrectSection.tcl this:
set Ubig 1.e10;
set G $Ubig;
set J 1.0;
set GJ [expr $G*$J]
section fiberSec $id -GJ $GJ {
patch quadr $coreID $nfCoreZ $nfCoreY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ
....
}
yes pls kindly attch the script.geoStru wrote: ↑Sat Aug 29, 2020 5:34 amif you want I will attacching a link for my script. I used "N" and "m" and you can try to see if the analysis convergesankurjain wrote: ↑Fri Aug 28, 2020 9:07 amI have changed the value inch in to meter directly and 1 kip to 4448.2216 Newton.....
but still the convergence is an issue.....
set in 0.0254; # define basic units -- output units
set kip 4448.2216; # define basic units -- output units
set sec 1.; # define basic units -- output units
set LunitTXT "inch"; # define basic-unit text for output
set FunitTXT "kip"; # define basic-unit text for output
set TunitTXT "sec"; # define basic-unit text for output
set ft [expr 12.*$in]; # define engineering units
set ksi [expr $kip/pow($in,2)];
set psi [expr $ksi/1000.];
set lbf [expr $psi*$in*$in]; # pounds force
set pcf [expr $lbf/pow($ft,3)]; # pounds per cubic foot
set psf [expr $lbf/pow($ft,3)]; # pounds per square foot
set in2 [expr $in*$in]; # inch^2
set in4 [expr $in*$in*$in*$in]; # inch^4
set cm [expr $in/2.54]; # centimeter, needed for displacement input in MultipleSupport excitation
set PI [expr 2*asin(1.0)]; # define constants
set g [expr 32.2*$ft/pow($sec,2)]; # gravitational acceleration
set Ubig 1.e10; # a really large number
set Usmall [expr 1/$Ubig]; # a really small number
I am attaching the link for the folder you can read and run my file tcl. With "Procedura.tcl" you run the analysis.ankurjain wrote: ↑Sat Aug 29, 2020 10:28 amyes pls kindly attch the script.geoStru wrote: ↑Sat Aug 29, 2020 5:34 amif you want I will attacching a link for my script. I used "N" and "m" and you can try to see if the analysis convergesankurjain wrote: ↑Fri Aug 28, 2020 9:07 am
I have changed the value inch in to meter directly and 1 kip to 4448.2216 Newton.....
but still the convergence is an issue.....
set in 0.0254; # define basic units -- output units
set kip 4448.2216; # define basic units -- output units
set sec 1.; # define basic units -- output units
set LunitTXT "inch"; # define basic-unit text for output
set FunitTXT "kip"; # define basic-unit text for output
set TunitTXT "sec"; # define basic-unit text for output
set ft [expr 12.*$in]; # define engineering units
set ksi [expr $kip/pow($in,2)];
set psi [expr $ksi/1000.];
set lbf [expr $psi*$in*$in]; # pounds force
set pcf [expr $lbf/pow($ft,3)]; # pounds per cubic foot
set psf [expr $lbf/pow($ft,3)]; # pounds per square foot
set in2 [expr $in*$in]; # inch^2
set in4 [expr $in*$in*$in*$in]; # inch^4
set cm [expr $in/2.54]; # centimeter, needed for displacement input in MultipleSupport excitation
set PI [expr 2*asin(1.0)]; # define constants
set g [expr 32.2*$ft/pow($sec,2)]; # gravitational acceleration
set Ubig 1.e10; # a really large number
set Usmall [expr 1/$Ubig]; # a really small number
yes this one is totally converging...no -3 error...geoStru wrote: ↑Sat Aug 29, 2020 11:06 amI am attaching the link for the folder you can read and run my file tcl. With "Procedura.tcl" you run the analysis.
https://drive.google.com/drive/folders/ ... sp=sharing
let me know
I didn't compare the pushover analysis with sap. I compared the eigen and static analysis with sap2000.ankurjain wrote: ↑Sat Aug 29, 2020 7:49 pmyes this one is totally converging...no -3 error...geoStru wrote: ↑Sat Aug 29, 2020 11:06 amI am attaching the link for the folder you can read and run my file tcl. With "Procedura.tcl" you run the analysis.
https://drive.google.com/drive/folders/ ... sp=sharing
let me know
Now the issue is whether the result of this analysis matches with the one from SAP2000 or not ?
Because I have done one pushover example that converge well but the same result didnt matched with that of SAP2000...means not even little bit close....
Any idea if you have tried matching the result with SAP2000 ?
If you have done any then kindly share both SAP2000 and opensees file so that i can compare .
Thanks in advance