/[escript]/trunk/modellib/py_src/mechanics.py
ViewVC logotype

Diff of /trunk/modellib/py_src/mechanics.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 816 by gross, Thu Aug 24 00:03:01 2006 UTC revision 820 by gross, Mon Aug 28 06:55:36 2006 UTC
# Line 15  class Mechanics(Model): Line 15  class Mechanics(Model):
15        base class for mechanics models in updated lagrangean framework        base class for mechanics models in updated lagrangean framework
16    
17        @ivar domain: domain (in)        @ivar domain: domain (in)
18        @ivar displacement: current displacements        @ivar internal_force=Data()
19          @ivar external_force=Data()
20          @ivar prescribed_velocity=Data()
21          @ivar location_prescribed_velocity=Data()
22          @ivar temperature = None
23          @ivar expansion_coefficient = 0.
24          @ivar bulk_modulus=1.
25          @ivar shear_modulus=1.
26          @ivar rel_tol=1.e-3
27          @ivar abs_tol=1.e-15
28          @ivar max_iter=10
29          @ivar displacement=None
30          @ivar stress=None
31          @ivar velocity=None
32        """        """
33        def __init__(self,debug=False):        def __init__(self,debug=False):
34           """           """
# Line 91  class Mechanics(Model): Line 103  class Mechanics(Model):
103            self.__pde.setValue(X=self.stress_old-self.dthermal_stress*k3)            self.__pde.setValue(X=self.stress_old-self.dthermal_stress*k3)
104            if self.internal_force: self.__pde.setValue(Y=self.internal_force)            if self.internal_force: self.__pde.setValue(Y=self.internal_force)
105            if self.external_force: self.__pde.setValue(y=self.external_force)            if self.external_force: self.__pde.setValue(y=self.external_force)
106              print self.prescribed_velocity
107              print self.location_prescribed_velocity
108            self.__pde.setValue(r=self.prescribed_velocity, \            self.__pde.setValue(r=self.prescribed_velocity, \
109                                q=self.location_prescribed_velocity)                                q=self.location_prescribed_velocity)
110            # solve the PDE:            # solve the PDE:
Legend:
Removed from v.816  
changed lines
  Added in v.820
  ViewVC Help
Powered by ViewVC 1.1.26