]> SALOME platform Git repositories - modules/geom.git/blob - src/GEOMImpl/GEOMImpl_ArchimedeDriver.cxx
Salome HOME
Add new command GetShapesOnBox.
[modules/geom.git] / src / GEOMImpl / GEOMImpl_ArchimedeDriver.cxx
1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
3 // 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either 
7 // version 2.1 of the License.
8 // 
9 // This library is distributed in the hope that it will be useful 
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of 
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
12 // Lesser General Public License for more details.
13 //
14 // You should have received a copy of the GNU Lesser General Public  
15 // License along with this library; if not, write to the Free Software 
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
17 //
18 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 //
20
21 #include "GEOMImpl_ArchimedeDriver.hxx"
22 #include "GEOMImpl_IArchimede.hxx"
23 #include "GEOMImpl_Types.hxx"
24
25 #include "Archimede_VolumeSection.hxx"
26
27 #include <stdio.h>
28
29 #include <BRepBuilderAPI_MakeFace.hxx>
30
31 #include <TopoDS_Shape.hxx>
32 #include <TopoDS_Face.hxx>
33
34 #include <Geom_Plane.hxx>
35 #include <Geom_Geometry.hxx>
36 #include <Geom_RectangularTrimmedSurface.hxx>
37
38 #include <gp_Pnt.hxx>
39 #include <gp_Dir.hxx>
40
41 #include <StdFail_NotDone.hxx>
42
43 using namespace std;
44
45 //=======================================================================
46 //function : GetID
47 //purpose  :
48 //=======================================================================
49 const Standard_GUID& GEOMImpl_ArchimedeDriver::GetID()
50 {
51   static Standard_GUID aArchimedeDriver("FF1BBB59-5D14-4df2-980B-3A668264EA16");
52   return aArchimedeDriver;
53 }
54
55
56 //=======================================================================
57 //function : GEOMImpl_ArchimedeDriver
58 //purpose  :
59 //=======================================================================
60
61 GEOMImpl_ArchimedeDriver::GEOMImpl_ArchimedeDriver()
62 {
63 }
64
65 //=======================================================================
66 //function : Execute
67 //purpose  :
68 //=======================================================================
69 Standard_Integer GEOMImpl_ArchimedeDriver::Execute(TFunction_Logbook& log) const
70 {
71   if (Label().IsNull()) return 0;
72   Handle(GEOM_Function) aFunction = GEOM_Function::GetFunction(Label());
73
74   Standard_Integer aType = aFunction->GetType();
75   if (aType != ARCHIMEDE_TYPE) return 0;
76
77   GEOMImpl_IArchimede IA (aFunction);
78
79   Handle(GEOM_Function) aShapeFunction = IA.GetBasicShape();
80   if (aShapeFunction.IsNull()) return 0;
81   TopoDS_Shape shape = aShapeFunction->GetValue();
82   if(shape.IsNull()) return 0;
83
84   Standard_Real aWaterDensity = IA.GetDensity();
85   Standard_Real aWeight = IA.GetWeight();
86   Standard_Real aMeshingDeflection = IA.GetDeflection();
87
88   double cste = -1;
89   if (aWaterDensity != 0.)
90     cste = aWeight/aWaterDensity;
91   else
92     return 0;
93
94   gp_Dir direct (0.0, 0.0, 1.0);
95   gp_Pnt PosPlan (0.0, 0.0, 0.0);
96   Geom_Plane PP (PosPlan, direct);
97   Handle(Geom_Plane) P = Handle(Geom_Plane)::DownCast(PP.Copy());
98
99   gp_Dir Zdirection (0.0, 0.0, 1.0);
100   VolumeSection VOL (shape, aMeshingDeflection);
101   VOL.SetPlane(P);
102   Handle(Geom_RectangularTrimmedSurface) SurfaceTrimmee;
103
104   if (Zdirection.IsEqual(direct, Precision::Angular()) == Standard_False) {
105     VOL.MakeRotation(direct);
106   }
107
108   VOL.CenterOfGravity();
109   SurfaceTrimmee = VOL.TrimSurf();
110   Standard_Real Cote = VOL.Archimede(cste, aMeshingDeflection);
111
112   if (Cote == -1) {
113     double Zmin,Zmax;
114     VOL.getZ(Zmin,Zmax);
115     double volume = VOL.CalculateVolume(Zmax) * aWaterDensity;
116
117     char msg[100] = "";
118     sprintf(msg, "shape sinks to the bottom : Weigth max = %.1f", volume);
119
120     StdFail_NotDone::Raise(msg);
121   }
122
123   SurfaceTrimmee = VOL.AjustePlan(SurfaceTrimmee,Cote,PosPlan);
124   if (Zdirection.IsEqual(direct,Precision::Angular()) == Standard_False) {
125     SurfaceTrimmee = VOL.InvMakeRotation(direct,SurfaceTrimmee);
126   }
127
128   Standard_Real u1,u2,v1,v2;
129   SurfaceTrimmee->Bounds(u1,u2,v1,v2);
130   TopoDS_Face tirant = BRepBuilderAPI_MakeFace(SurfaceTrimmee, u1, u2, v1, v2);
131
132   if (tirant.IsNull()) {
133     StdFail_NotDone::Raise("Failed to build secant face");
134   }
135
136   aFunction->SetValue(tirant);
137
138   log.SetTouched(Label());
139
140   return 1;
141 }
142
143
144 //=======================================================================
145 //function :  GEOMImpl_ArchimedeDriver_Type_
146 //purpose  :
147 //=======================================================================
148 Standard_EXPORT Handle_Standard_Type& GEOMImpl_ArchimedeDriver_Type_()
149 {
150
151   static Handle_Standard_Type aType1 = STANDARD_TYPE(TFunction_Driver);
152   if ( aType1.IsNull()) aType1 = STANDARD_TYPE(TFunction_Driver);
153   static Handle_Standard_Type aType2 = STANDARD_TYPE(MMgt_TShared);
154   if ( aType2.IsNull()) aType2 = STANDARD_TYPE(MMgt_TShared);
155   static Handle_Standard_Type aType3 = STANDARD_TYPE(Standard_Transient);
156   if ( aType3.IsNull()) aType3 = STANDARD_TYPE(Standard_Transient);
157
158
159   static Handle_Standard_Transient _Ancestors[]= {aType1,aType2,aType3,NULL};
160   static Handle_Standard_Type _aType = new Standard_Type("GEOMImpl_ArchimedeDriver",
161                                                          sizeof(GEOMImpl_ArchimedeDriver),
162                                                          1,
163                                                          (Standard_Address)_Ancestors,
164                                                          (Standard_Address)NULL);
165
166   return _aType;
167 }
168
169 //=======================================================================
170 //function : DownCast
171 //purpose  :
172 //=======================================================================
173
174 const Handle(GEOMImpl_ArchimedeDriver) Handle(GEOMImpl_ArchimedeDriver)::DownCast(const Handle(Standard_Transient)& AnObject)
175 {
176   Handle(GEOMImpl_ArchimedeDriver) _anOtherObject;
177
178   if (!AnObject.IsNull()) {
179      if (AnObject->IsKind(STANDARD_TYPE(GEOMImpl_ArchimedeDriver))) {
180        _anOtherObject = Handle(GEOMImpl_ArchimedeDriver)((Handle(GEOMImpl_ArchimedeDriver)&)AnObject);
181      }
182   }
183
184   return _anOtherObject ;
185 }
186
187